Salome HOME
fix biquadratic quadratic pentahedron(18 nodes prism)
[modules/smesh.git] / src / MEDWrapper / Base / MED_Utilities.cxx
1 // Copyright (C) 2007-2016  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, or (at your option) any later version.
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 #include "MED_Utilities.hxx"
23 #include "MED_Common.hxx"
24
25 using namespace std;
26
27 #ifdef _DEBUG_
28 static int MYDEBUG = 0;
29 #else
30 // static int MYDEBUG = 0;
31 #endif
32
33
34 int MED::PrefixPrinter::myCounter = 0;
35
36 MED::PrefixPrinter::PrefixPrinter(bool theIsActive):
37   myIsActive(theIsActive)
38 {
39   if(myIsActive)
40     myCounter++;
41   MSG(MYDEBUG,"MED::PrefixPrinter::PrefixPrinter(...)- "<<myCounter);
42 }
43
44 MED::PrefixPrinter::~PrefixPrinter()
45 {
46   if(myIsActive){
47     myCounter--;
48     if(myCounter < 0)
49       EXCEPTION(runtime_error,"PrefixPrinter::~PrefixPrinter() - myCounter("<<myCounter<<") < 0");
50   }
51 }
52
53 string MED::PrefixPrinter::GetPrefix()
54 {
55   if(myCounter){
56     if(myCounter < 0)
57       EXCEPTION(runtime_error,"PrefixPrinter::~PrefixPrinter() - myCounter("<<myCounter<<") < 0");
58     return string(myCounter*2,' ');
59   }
60   return "";
61 }
62
63
64 static MED::TEntity2GeomSet Entity2GeomSet;
65
66 bool InitEntity2GeomSet()
67 {
68   using namespace MED;
69
70   TGeomSet& aGeomARETESet = Entity2GeomSet[eARETE];
71   aGeomARETESet.insert(eSEG2);
72   aGeomARETESet.insert(eSEG3);
73
74   TGeomSet& aGeomFACESet = Entity2GeomSet[eFACE];
75   aGeomFACESet.insert(eTRIA3);
76   aGeomFACESet.insert(eQUAD4);
77   aGeomFACESet.insert(eTRIA6);
78   aGeomFACESet.insert(eTRIA7);
79   aGeomFACESet.insert(eQUAD8);
80   aGeomFACESet.insert(eQUAD9);
81   aGeomFACESet.insert(ePOLYGONE);
82   aGeomFACESet.insert(ePOLYGON2);
83
84   TGeomSet& aGeomMAILLESet = Entity2GeomSet[eMAILLE];
85   aGeomMAILLESet.insert(ePOINT1);
86   aGeomMAILLESet.insert(aGeomARETESet.begin(),aGeomARETESet.end());
87   aGeomMAILLESet.insert(aGeomFACESet.begin(),aGeomFACESet.end());
88   aGeomMAILLESet.insert(eTETRA4);
89   aGeomMAILLESet.insert(ePYRA5);
90   aGeomMAILLESet.insert(ePENTA6);
91   aGeomMAILLESet.insert(eHEXA8);
92   aGeomMAILLESet.insert(eOCTA12);
93   aGeomMAILLESet.insert(eTETRA10);
94   aGeomMAILLESet.insert(ePYRA13);
95   aGeomMAILLESet.insert(ePENTA15);
96   aGeomMAILLESet.insert(ePENTA18);
97   aGeomMAILLESet.insert(eHEXA20);
98   aGeomMAILLESet.insert(eHEXA27);
99   aGeomMAILLESet.insert(ePOLYEDRE);
100
101   /* This combination allows reading nb of models of structure elements */
102   Entity2GeomSet[eSTRUCT_ELEMENT].insert(eAllGeoType); 
103
104   return true;
105 }
106
107 static bool anIsInited = InitEntity2GeomSet();
108
109 const MED::TEntity2GeomSet& MED::GetEntity2GeomSet()
110 {
111   return Entity2GeomSet;
112 }
113
114