]> SALOME platform Git repositories - modules/smesh.git/blob - src/DriverUNV/UNV2417_Structure.cxx
Salome HOME
PAL15429 Computation of the mesh, based on "014.brep" via Tetrahedron(NETGEN), is...
[modules/smesh.git] / src / DriverUNV / UNV2417_Structure.cxx
1 //  Copyright (C) 2007-2008  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 #include "UNV2417_Structure.hxx"
23 #include "UNV_Utilities.hxx"
24
25 #include <fstream>      
26 #include <iomanip>
27
28 using namespace std;
29 using namespace UNV;
30 using namespace UNV2417;
31
32 #ifdef _DEBUG_
33 static int MYDEBUG = 0;
34 #else
35 static int MYDEBUG = 0;
36 #endif
37
38
39 static string _group_labels[] = {"2417", "2429", "2430", "2432",
40                                  "2435", "2452", "2467", "2477"};
41 #define NBGROUP 8
42
43 static string _label_dataset = "2467";
44
45 void UNV2417::Read(std::ifstream& in_stream, TDataSet& theDataSet)
46 {
47   if(!in_stream.good())
48     EXCEPTION(runtime_error,"ERROR: Input file not good.");
49
50   std::string olds, news;
51   
52   while(true){
53     in_stream >> olds >> news;
54     /*
55      * a "-1" followed by a number means the beginning of a dataset
56      * stop combing at the end of the file
57      */
58     while( ((olds != "-1") || (news == "-1") ) && !in_stream.eof() ){     
59       olds = news;
60       in_stream >> news;
61     }
62     if(in_stream.eof())
63       return;
64     for (int i = 0; i < NBGROUP; i++) {
65       if (news == _group_labels[i]) {
66         ReadGroup(news, in_stream, theDataSet);
67       }
68     }
69   }
70 }
71
72
73
74 void UNV2417::ReadGroup(const std::string& myGroupLabel, std::ifstream& in_stream, TDataSet& theDataSet)
75 {
76   TGroupId aId;
77   for(; !in_stream.eof();){
78     in_stream >> aId ;
79     if(aId == -1){
80       // end of dataset is reached
81       break;
82     }
83
84     int n_nodes;
85     TRecord aRec;
86     int aTmp;
87     in_stream>>aTmp; // miss not necessary values
88     in_stream>>aTmp;
89     in_stream>>aTmp;
90     in_stream>>aTmp;
91     in_stream>>aTmp;
92     in_stream>>aTmp;
93     in_stream>>n_nodes;
94
95     std::getline(in_stream, aRec.GroupName, '\n'); // Finalise previous reading
96     std::getline(in_stream, aRec.GroupName, '\n');
97
98     int aElType;
99     int aElId;
100     int aNum;
101     for(int j=0; j < n_nodes; j++){
102       in_stream>>aElType;
103       in_stream>>aElId;
104       if ((myGroupLabel.compare("2435") == 0) ||
105           (myGroupLabel.compare("2452") == 0) ||
106           (myGroupLabel.compare("2467") == 0) ||
107           (myGroupLabel.compare("2477") == 0)) {
108         in_stream>>aTmp;
109         in_stream>>aTmp;
110       }
111       switch (aElType) {
112       case 7: // Nodes
113         aNum = aRec.NodeList.size();
114         aRec.NodeList.resize(aNum + 1);
115         aRec.NodeList[aNum] = aElId;
116         break;
117       case 8: // Elements
118         aNum = aRec.ElementList.size();
119         aRec.ElementList.resize(aNum + 1);
120         aRec.ElementList[aNum] = aElId;
121         break;
122       }
123     }
124     theDataSet.insert(TDataSet::value_type(aId,aRec));
125   }
126
127 }
128
129
130 void UNV2417::Write(std::ofstream& out_stream, const TDataSet& theDataSet)
131 {
132   if(!out_stream.good())
133     EXCEPTION(runtime_error,"ERROR: Output file not good.");
134   
135   /*
136    * Write beginning of dataset
137    */
138   out_stream<<"    -1\n";
139   out_stream<<"  "<<_label_dataset<<"\n";
140
141   TDataSet::const_iterator anIter = theDataSet.begin();
142   for(; anIter != theDataSet.end(); anIter++){
143     const TGroupId& aLabel = anIter->first;
144     const TRecord& aRec = anIter->second;
145     int aNbNodes = aRec.NodeList.size();
146     int aNbElements = aRec.ElementList.size();
147     int aNbRecords = aNbNodes + aNbElements;
148
149     out_stream<<std::setw(10)<<aLabel;  /* group ID */
150     out_stream<<std::setw(10)<<0;  
151     out_stream<<std::setw(10)<<0;
152     out_stream<<std::setw(10)<<0;
153     out_stream<<std::setw(10)<<0;
154     out_stream<<std::setw(10)<<0;
155     out_stream<<std::setw(10)<<0;
156     out_stream<<std::setw(10)<<aNbRecords<<std::endl; 
157
158     out_stream<<aRec.GroupName<<std::endl;
159     int aRow = 0;
160     int i;
161     for (i = 0; i < aNbNodes; i++) {
162       if (aRow == 2) {
163         out_stream<<std::endl; 
164         aRow = 0;
165       }
166       out_stream<<std::setw(10)<<7;
167       out_stream<<std::setw(10)<<aRec.NodeList[i];
168       out_stream<<std::setw(10)<<0;
169       out_stream<<std::setw(10)<<0;
170       aRow++;
171     }
172     for (i = 0; i < aNbElements; i++) {
173       if (aRow == 2) {
174         out_stream<<std::endl; 
175         aRow = 0;
176       }
177       out_stream<<std::setw(10)<<8;
178       out_stream<<std::setw(10)<<aRec.ElementList[i];
179       out_stream<<std::setw(10)<<0;
180       out_stream<<std::setw(10)<<0;
181       aRow++;
182     }
183     out_stream<<std::endl; 
184   }
185
186   /*
187    * Write end of dataset
188    */
189   out_stream<<"    -1\n";
190 }