]> SALOME platform Git repositories - modules/smesh.git/blob - src/MEDWrapper/Base/MED_CoordUtils.cxx
Salome HOME
CMake:
[modules/smesh.git] / src / MEDWrapper / Base / MED_CoordUtils.cxx
1 // Copyright (C) 2007-2013  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22 //
23 #include "MED_CoordUtils.hxx"
24 #include "MED_Utilities.hxx"
25  
26 namespace MED
27 {
28
29   enum ECoordName{eX, eY, eZ, eNone};
30
31   template<ECoordName TCoordId>
32   TFloat 
33   GetCoord(const TCCoordSlice& theCoordSlice)
34   {
35     return theCoordSlice[TCoordId];
36   }
37
38   template<>
39   TFloat 
40   GetCoord<eNone>(const TCCoordSlice& theCoordSlice)
41   {
42     return 0.0;
43   }
44   
45   TGetCoord
46   aXYZGetCoord[3] = {
47     &GetCoord<eX>, 
48     &GetCoord<eY>, 
49     &GetCoord<eZ>
50   };
51   
52   TGetCoord
53   aXYGetCoord[3] = {
54     &GetCoord<eX>, 
55     &GetCoord<eY>, 
56     &GetCoord<eNone>
57   };
58   
59   TGetCoord
60   aYZGetCoord[3] = {
61     &GetCoord<eNone>,
62     &GetCoord<eX>, 
63     &GetCoord<eY>
64   };
65   
66   TGetCoord 
67   aXZGetCoord[3] = {
68     &GetCoord<eX>, 
69     &GetCoord<eNone>,
70     &GetCoord<eY>
71   };
72   
73   
74   TGetCoord 
75   aXGetCoord[3] = {
76     &GetCoord<eX>, 
77     &GetCoord<eNone>,
78     &GetCoord<eNone>
79   };
80   
81   TGetCoord
82   aYGetCoord[3] = {
83     &GetCoord<eNone>,
84     &GetCoord<eX>, 
85     &GetCoord<eNone>
86   };
87
88   TGetCoord
89   aZGetCoord[3] = {
90     &GetCoord<eNone>,
91     &GetCoord<eNone>,
92     &GetCoord<eX>
93   };
94
95   
96   //---------------------------------------------------------------
97   TCoordHelper
98   ::TCoordHelper(TGetCoord* theGetCoord):
99     myGetCoord(theGetCoord)
100   {}
101
102   TFloat 
103   TCoordHelper
104   ::GetCoord(TCCoordSlice& theCoordSlice, 
105              TInt theCoordId)
106   {
107     return (*myGetCoord[theCoordId])(theCoordSlice);
108   }
109
110
111   //---------------------------------------------------------------
112   PCoordHelper
113   GetCoordHelper(PNodeInfo theNodeInfo)
114   {
115     PCoordHelper aCoordHelper;
116     {
117       PMeshInfo aMeshInfo = theNodeInfo->GetMeshInfo();
118       TInt aMeshDimension = aMeshInfo->GetDim();
119       bool anIsDimPresent[3] = {false, false, false};
120       for(int iDim = 0; iDim < aMeshDimension; iDim++){
121         // PAL16857(SMESH not conform to the MED convention) ->
122         // 1D - always along X
123         // 2D - always in XOY plane
124         anIsDimPresent[iDim] = iDim < aMeshDimension;
125 //      std::string aName = theNodeInfo->GetCoordName(iDim);
126 //         if ( aName.size() > 1 ) // PAL12148, aName has size 8 or 16
127 //           aName = aName.substr(0,1);
128 //      if(aName == "x" || aName == "X")
129 //        anIsDimPresent[eX] = true;
130 //      else if(aName == "y" || aName == "Y")
131 //        anIsDimPresent[eY] = true;
132 //      else if(aName == "z" || aName == "Z")
133 //        anIsDimPresent[eZ] = true;
134       }
135
136       switch(aMeshDimension){
137       case 3:
138         aCoordHelper.reset(new TCoordHelper(aXYZGetCoord));
139         break;
140       case 2:
141         if(anIsDimPresent[eY] && anIsDimPresent[eZ])
142           aCoordHelper.reset(new TCoordHelper(aYZGetCoord));
143         else if(anIsDimPresent[eX] && anIsDimPresent[eZ])
144           aCoordHelper.reset(new TCoordHelper(aXZGetCoord));
145         else
146           aCoordHelper.reset(new TCoordHelper(aXYGetCoord));
147         break;
148       case 1:
149         if(anIsDimPresent[eY])
150           aCoordHelper.reset(new TCoordHelper(aYGetCoord));
151         else if(anIsDimPresent[eZ])
152           aCoordHelper.reset(new TCoordHelper(aZGetCoord));
153         else
154           aCoordHelper.reset(new TCoordHelper(aXGetCoord));
155         break;
156       }
157     }
158     return aCoordHelper;
159   }
160 }