1 // Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
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.
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.
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
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
26 #include "UNV2411_Structure.hxx"
27 #include "UNV_Utilities.hxx"
31 using namespace UNV2411;
33 static string _label_dataset = "2411";
35 UNV2411::TRecord::TRecord():
37 disp_coord_sys_num(1),
38 color(11)//(0) - 0019936: EDF 794 SMESH : Export UNV : Node color and group id
40 coord[1] = coord[2] = 0.0; // prepare to e.g. 2D mesh
43 void UNV2411::Read(std::ifstream& in_stream, TDataSet& theDataSet)
46 EXCEPTION(runtime_error,"ERROR: Input file not good.");
49 * adjust the \p istream to our
52 if(!beginning_of_dataset(in_stream,_label_dataset))
53 EXCEPTION(runtime_error,"ERROR: Could not find "<<_label_dataset<<" dataset!");
56 * always 3 coordinates in the UNV file, no matter
57 * which dimensionality libMesh is in
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() )
66 int where = in_stream.tellg();
69 in_stream >> aRec.label ;
70 if ( aRec.label == -1 ) return; // dataset end
73 num_buf = read_line( in_stream );
74 for ( size_t i = 0; i < num_buf.size(); )
77 while ( i < num_buf.size() && num_buf[i] == ' ' )
80 dim += ( i < num_buf.size() );
83 while ( i < num_buf.size() && num_buf[i] != ' ' )
89 in_stream.seekg( where, in_stream.beg );
92 // read the rest records
93 while ( !in_stream.eof() )
96 in_stream >> aRec.label ;
97 if ( aRec.label == -1 ) {
98 // end of dataset is reached
102 in_stream>>aRec.exp_coord_sys_num;
103 in_stream>>aRec.disp_coord_sys_num;
104 in_stream>>aRec.color;
108 * floating-point data
110 for(int d = 0; d < dim; d++){
112 aRec.coord[d] = D_to_e(num_buf);
115 theDataSet.push_back(aRec);
120 void UNV2411::Write(std::ofstream& out_stream, const TDataSet& theDataSet)
122 if(!out_stream.good())
123 EXCEPTION(runtime_error,"ERROR: Output file not good.");
126 * Write beginning of dataset
129 out_stream<<" "<<_label_dataset<<"\n";
131 TDataSet::const_iterator anIter = theDataSet.begin();
132 for(; anIter != theDataSet.end(); anIter++)
134 const TRecord& aRec = *anIter;
136 sprintf(buf, "%10d%10d%10d%10d\n",
138 aRec.exp_coord_sys_num,
139 aRec.disp_coord_sys_num,
144 sprintf(buf, "%25.16E%25.16E%25.16E\n",
153 * Write end of dataset