]> SALOME platform Git repositories - modules/smesh.git/blob - src/DriverUNV/UNV2411_Structure.cxx
Salome HOME
Cosmetic changes
[modules/smesh.git] / src / DriverUNV / UNV2411_Structure.cxx
1 // Copyright (C) 2007-2020  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
23 #include <fstream>
24 #include <cstdio>
25
26 #include "UNV2411_Structure.hxx"
27 #include "UNV_Utilities.hxx"
28
29 using namespace std;
30 using namespace UNV;
31 using namespace UNV2411;
32
33 static string _label_dataset = "2411";
34
35 UNV2411::TRecord::TRecord():
36   exp_coord_sys_num(1),
37   disp_coord_sys_num(1),
38   color(11)//(0) -  0019936: EDF 794 SMESH : Export UNV : Node color and group id
39 {
40   coord[1] = coord[2] = 0.0; // prepare to e.g. 2D mesh
41 }
42
43 void UNV2411::Read(std::ifstream& in_stream, TDataSet& theDataSet)
44 {
45   if(!in_stream.good())
46     EXCEPTION(runtime_error,"ERROR: Input file not good.");
47
48   /*
49    * adjust the \p istream to our
50    * position
51    */
52   if(!beginning_of_dataset(in_stream,_label_dataset))
53     EXCEPTION(runtime_error,"ERROR: Could not find "<<_label_dataset<<" dataset!");
54
55   /**
56    * always 3 coordinates in the UNV file, no matter
57    * which dimensionality libMesh is in
58    */
59   int dim = 3;
60   std::string num_buf;
61
62   // Issue 22638. Find out space dimension to read a 2D mesh from a file
63   // generated by SIMAIL from Simulog
64   if ( !in_stream.eof() )
65   {
66     int where  = in_stream.tellg();
67
68     TRecord aRec;
69     in_stream >> aRec.label ;
70     if ( aRec.label == -1 ) return; // dataset end 
71
72     dim = 0;
73     num_buf = read_line( in_stream );
74     for ( size_t i = 0; i < num_buf.size(); )
75     {
76       // skip spaces
77       while ( i < num_buf.size() && num_buf[i] == ' ' )
78         ++i;
79
80       dim += ( i < num_buf.size() );
81
82       // skip non-spaces
83       while ( i < num_buf.size() && num_buf[i] != ' ' )
84         ++i;
85     }
86     if ( dim == 0 )
87       return;
88
89     in_stream.seekg( where, in_stream.beg );
90   }
91
92   // read the rest records
93   while ( !in_stream.eof() )
94   {
95     TRecord aRec;
96     in_stream >> aRec.label ;
97     if ( aRec.label == -1 ) {
98       // end of dataset is reached
99       break;
100     }
101
102     in_stream>>aRec.exp_coord_sys_num;
103     in_stream>>aRec.disp_coord_sys_num;
104     in_stream>>aRec.color;
105
106     /*
107      * take care of the
108      * floating-point data
109      */
110     for(int d = 0; d < dim; d++){
111       in_stream>>num_buf;
112       aRec.coord[d] = D_to_e(num_buf);
113     }
114
115     theDataSet.push_back(aRec);
116   }
117 }
118
119
120 void UNV2411::Write(std::ofstream& out_stream, const TDataSet& theDataSet)
121 {
122   if(!out_stream.good())
123     EXCEPTION(runtime_error,"ERROR: Output file not good.");
124   
125   /*
126    * Write beginning of dataset
127    */
128   out_stream<<"    -1\n";
129   out_stream<<"  "<<_label_dataset<<"\n";
130
131   TDataSet::const_iterator anIter = theDataSet.begin();
132   for(; anIter != theDataSet.end(); anIter++)
133   {
134     const TRecord& aRec = *anIter;
135     char buf[78];
136     sprintf(buf, "%10d%10d%10d%10d\n", 
137             aRec.label,
138             aRec.exp_coord_sys_num,
139             aRec.disp_coord_sys_num,
140             aRec.color);
141     out_stream<<buf;
142
143     // the coordinates
144     sprintf(buf, "%25.16E%25.16E%25.16E\n", 
145             aRec.coord[0],
146             aRec.coord[1],
147             aRec.coord[2]);
148     out_stream<<buf;
149   }
150   
151   
152   /*
153    * Write end of dataset
154    */
155   out_stream<<"    -1\n";
156 }