1 // Copyright (C) 2007-2008 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.
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
25 #include "UNV2411_Structure.hxx"
26 #include "UNV_Utilities.hxx"
30 using namespace UNV2411;
33 static int MYDEBUG = 1;
35 static int MYDEBUG = 0;
38 static string _label_dataset = "2411";
40 UNV2411::TRecord::TRecord():
42 disp_coord_sys_num(0),
46 void UNV2411::Read(std::ifstream& in_stream, TDataSet& theDataSet)
49 EXCEPTION(runtime_error,"ERROR: Input file not good.");
52 * adjust the \p istream to our
55 if(!beginning_of_dataset(in_stream,_label_dataset))
56 EXCEPTION(runtime_error,"ERROR: Could not find "<<_label_dataset<<" dataset!");
59 * always 3 coordinates in the UNV file, no matter
60 * which dimensionality libMesh is in
64 for(; !in_stream.eof();){
67 // end of dataset is reached
72 in_stream>>aRec.exp_coord_sys_num;
73 in_stream>>aRec.disp_coord_sys_num;
74 in_stream>>aRec.color;
80 for(int d = 0; d < 3; d++){
82 aRec.coord[d] = D_to_e(num_buf);
85 theDataSet.insert(TDataSet::value_type(aLabel,aRec));
90 void UNV2411::Write(std::ofstream& out_stream, const TDataSet& theDataSet)
92 if(!out_stream.good())
93 EXCEPTION(runtime_error,"ERROR: Output file not good.");
96 * Write beginning of dataset
99 out_stream<<" "<<_label_dataset<<"\n";
101 TDataSet::const_iterator anIter = theDataSet.begin();
102 for(; anIter != theDataSet.end(); anIter++){
103 const TNodeLab& aLabel = anIter->first;
104 const TRecord& aRec = anIter->second;
106 sprintf(buf, "%10d%10d%10d%10d\n",
108 aRec.exp_coord_sys_num,
109 aRec.disp_coord_sys_num,
114 sprintf(buf, "%25.16E%25.16E%25.16E\n",
123 * Write end of dataset