1 // SMESH DriverDAT : driver to read and write 'dat' files
3 // Copyright (C) 2003 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.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
24 // File : DriverDAT_R_SMDS_Mesh.cxx
27 #include "DriverDAT_R_SMDS_Mesh.h"
28 #include "utilities.h"
34 * Factory function which will be called by SMESHDriver
36 void * SMESH_createDATMeshReader()
38 return new DriverDAT_R_SMDS_Mesh();
42 DriverDAT_R_SMDS_Mesh::DriverDAT_R_SMDS_Mesh()
47 DriverDAT_R_SMDS_Mesh::~DriverDAT_R_SMDS_Mesh()
52 void DriverDAT_R_SMDS_Mesh::SetMesh(SMDS_Mesh * aMesh)
57 void DriverDAT_R_SMDS_Mesh::SetFile(string aFile)
62 void DriverDAT_R_SMDS_Mesh::SetFileId(FILE * aFileId)
67 void DriverDAT_R_SMDS_Mesh::SetMeshId(int aMeshId)
72 void DriverDAT_R_SMDS_Mesh::Add()
77 void DriverDAT_R_SMDS_Mesh::Read()
83 float coordX, coordY, coordZ;
86 int intNumMaille, Degre;
94 MESSAGE("in DriverDAT_R_SMDS_Mesh::Read()");
95 /****************************************************************************
96 * OUVERTURE DU FICHIER EN LECTURE *
97 ****************************************************************************/
98 char *file2Read = (char *)myFile.c_str();
99 myFileId = fopen(file2Read, "r");
102 fprintf(stderr, ">> ERREUR : ouverture du fichier %s \n", file2Read);
106 fscanf(myFileId, "%d %d\n", &nbNodes, &nbCells);
108 /****************************************************************************
109 * LECTURE DES NOEUDS *
110 ****************************************************************************/
111 fprintf(stdout, "\n(************************)\n");
112 fprintf(stdout, "(* NOEUDS DU MAILLAGE : *)\n");
113 fprintf(stdout, "(************************)\n");
115 for (i = 0; i < nbNodes; i++)
117 fscanf(myFileId, "%d %e %e %e\n", &intNumPoint, &coordX, &coordY,
119 ok = myMesh->AddNodeWithID(coordX, coordY, coordZ, intNumPoint);
122 fprintf(stdout, "%d noeuds\n", myMesh->NbNodes());
123 /****************************************************************************
124 * LECTURE DES ELEMENTS *
125 ****************************************************************************/
126 fprintf(stdout, "\n(**************************)\n");
127 fprintf(stdout, "(* ELEMENTS DU MAILLAGE : *)\n");
128 fprintf(stdout, "(**************************)");
130 fprintf(stdout, "%d elements\n", nbCells);
132 for (i = 0; i < nbCells; i++)
134 fscanf(myFileId, "%d %d", &intNumMaille, &ValElement);
135 Degre = abs(ValElement / 100);
136 nbNoeuds = ValElement - (Degre * 100);
138 // Recuperation des noeuds de la maille
139 for (j = 0; j < nbNoeuds; j++)
141 fscanf(myFileId, "%d", &NoeudMaille);
142 NoeudsMaille[j] = NoeudMaille;
145 // Analyse des cas de cellules
153 ok = myMesh->AddEdgeWithID(NoeudsMaille[0], NoeudsMaille[1],
162 ok = myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
163 NoeudsMaille[2], NoeudsMaille[3], intNumMaille);
171 ok = myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
172 NoeudsMaille[2], intNumMaille);
180 if (ValElement == 320)
182 //A voir, correspondance VTK
183 NoeudsMaille[4] = NoeudsMaille[8];
184 NoeudsMaille[5] = NoeudsMaille[9];
185 NoeudsMaille[6] = NoeudsMaille[10];
186 NoeudsMaille[7] = NoeudsMaille[11];
188 ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
189 NoeudsMaille[2], NoeudsMaille[3], NoeudsMaille[4],
190 NoeudsMaille[5], NoeudsMaille[6], NoeudsMaille[7],
199 if (ValElement == 310)
200 NoeudsMaille[3] = NoeudsMaille[6];
201 ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
202 NoeudsMaille[2], NoeudsMaille[3], intNumMaille);
210 if (ValElement == 315)
212 NoeudsMaille[3] = NoeudsMaille[6];
213 NoeudsMaille[4] = NoeudsMaille[7];
214 NoeudsMaille[5] = NoeudsMaille[8];
216 NoeudsMaille[7] = NoeudsMaille[5];
217 NoeudsMaille[6] = NoeudsMaille[5];
218 NoeudsMaille[5] = NoeudsMaille[4];
219 NoeudsMaille[4] = NoeudsMaille[3];
220 NoeudsMaille[3] = NoeudsMaille[2];
221 ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
222 NoeudsMaille[2], NoeudsMaille[3], NoeudsMaille[4],
223 NoeudsMaille[5], intNumMaille);
230 /****************************************************************************
231 * FERMETURE DU FICHIER *
232 ****************************************************************************/