]> SALOME platform Git repositories - modules/smesh.git/blob - src/DriverDAT/DriverDAT_R_SMDS_Mesh.cxx
Salome HOME
PAL8238 - improvement for hypothesis "number of segments"
[modules/smesh.git] / src / DriverDAT / DriverDAT_R_SMDS_Mesh.cxx
1 //  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS 
3 // 
4 //  This library is free software; you can redistribute it and/or 
5 //  modify it under the terms of the GNU Lesser General Public 
6 //  License as published by the Free Software Foundation; either 
7 //  version 2.1 of the License. 
8 // 
9 //  This library is distributed in the hope that it will be useful, 
10 //  but WITHOUT ANY WARRANTY; without even the implied warranty of 
11 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
12 //  Lesser General Public License for more details. 
13 // 
14 //  You should have received a copy of the GNU Lesser General Public 
15 //  License along with this library; if not, write to the Free Software 
16 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA 
17 // 
18 //  See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org 
19
20 #include <stdio.h>
21
22 #include "DriverDAT_R_SMDS_Mesh.h"
23 #include "SMDS_Mesh.hxx"
24
25 #include "utilities.h"
26
27 using namespace std;
28
29 Driver_Mesh::Status DriverDAT_R_SMDS_Mesh::Perform()
30 {
31   Status aResult = DRS_OK;
32
33   int i, j;
34   int nbNodes, nbCells;
35   int intNumPoint;
36   float coordX, coordY, coordZ;
37   int nbNoeuds;
38   
39   int intNumMaille, Degre;
40   int ValElement;
41   int ValDegre;
42   int NoeudsMaille[20];
43   int NoeudMaille;
44   
45   bool ok;
46   
47   MESSAGE("in DriverDAT_R_SMDS_Mesh::Read()");
48   /****************************************************************************
49    *                      OUVERTURE DU FICHIER EN LECTURE                      *
50    ****************************************************************************/
51   char *file2Read = (char *)myFile.c_str();
52   FILE* aFileId = fopen(file2Read, "r");
53   if (aFileId < 0) {
54     fprintf(stderr, ">> ERREUR : ouverture du fichier %s \n", file2Read);
55     return DRS_FAIL;
56   }
57   
58   fscanf(aFileId, "%d %d\n", &nbNodes, &nbCells);
59   
60   /****************************************************************************
61    *                       LECTURE DES NOEUDS                                  *
62    ****************************************************************************/
63   fprintf(stdout, "\n(************************)\n");
64   fprintf(stdout, "(* NOEUDS DU MAILLAGE : *)\n");
65   fprintf(stdout, "(************************)\n");
66   
67   for (i = 0; i < nbNodes; i++){
68     fscanf(aFileId, "%d %e %e %e\n", &intNumPoint, &coordX, &coordY, &coordZ);
69     ok = myMesh->AddNodeWithID(coordX, coordY, coordZ, intNumPoint);
70   }
71   
72   fprintf(stdout, "%d noeuds\n", myMesh->NbNodes());
73   /****************************************************************************
74    *                       LECTURE DES ELEMENTS                                *
75    ****************************************************************************/
76   fprintf(stdout, "\n(**************************)\n");
77   fprintf(stdout, "(* ELEMENTS DU MAILLAGE : *)\n");
78   fprintf(stdout, "(**************************)");
79   
80   fprintf(stdout, "%d elements\n", nbCells);
81   
82   for (i = 0; i < nbCells; i++) {
83     fscanf(aFileId, "%d %d", &intNumMaille, &ValElement);
84     Degre = abs(ValElement / 100);
85     nbNoeuds = ValElement - (Degre * 100);
86     
87     // Recuperation des noeuds de la maille
88     for (j = 0; j < nbNoeuds; j++) {
89       fscanf(aFileId, "%d", &NoeudMaille);
90       NoeudsMaille[j] = NoeudMaille;
91     }
92     
93     // Analyse des cas de cellules
94     switch (ValElement) {
95     case 102:
96     case 103:
97       ValDegre = 3;
98       nbNoeuds = 2;
99       ok = myMesh->AddEdgeWithID(NoeudsMaille[0], NoeudsMaille[1], 
100                                  intNumMaille);
101       break;
102     case 204:
103     case 208:
104       ValDegre = 9;
105       nbNoeuds = 4;
106       ok = myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
107                                  NoeudsMaille[2], NoeudsMaille[3], 
108                                  intNumMaille);
109       break;
110     case 203:
111     case 206:
112       ValDegre = 5;
113       nbNoeuds = 3;
114       ok = myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
115                                  NoeudsMaille[2], intNumMaille);
116       break;
117     case 308:
118     case 320:
119       ValDegre = 12;
120       nbNoeuds = 8;
121       if (ValElement == 320){
122         //A voir, correspondance VTK
123         NoeudsMaille[4] = NoeudsMaille[8];
124         NoeudsMaille[5] = NoeudsMaille[9];
125         NoeudsMaille[6] = NoeudsMaille[10];
126         NoeudsMaille[7] = NoeudsMaille[11];
127       }
128       ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
129                                    NoeudsMaille[2], NoeudsMaille[3], 
130                                    NoeudsMaille[4], NoeudsMaille[5], 
131                                    NoeudsMaille[6], NoeudsMaille[7],
132                                    intNumMaille);
133       break;
134     case 304:
135     case 310:
136       ValDegre = 10;
137       nbNoeuds = 4;
138       if (ValElement == 310)
139         NoeudsMaille[3] = NoeudsMaille[6];
140       ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
141                                    NoeudsMaille[2], NoeudsMaille[3], 
142                                    intNumMaille);
143       break;
144     case 306:
145     case 315:
146       ValDegre = 12;
147       nbNoeuds = 8;
148       if (ValElement == 315) {
149         NoeudsMaille[3] = NoeudsMaille[6];
150         NoeudsMaille[4] = NoeudsMaille[7];
151         NoeudsMaille[5] = NoeudsMaille[8];
152       }
153       NoeudsMaille[7] = NoeudsMaille[5];
154       NoeudsMaille[6] = NoeudsMaille[5];
155       NoeudsMaille[5] = NoeudsMaille[4];
156       NoeudsMaille[4] = NoeudsMaille[3];
157       NoeudsMaille[3] = NoeudsMaille[2];
158       ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
159                                    NoeudsMaille[2], NoeudsMaille[3], 
160                                    NoeudsMaille[4], NoeudsMaille[5], 
161                                    intNumMaille);
162                                 break;
163     }
164   }
165   /****************************************************************************
166    *                      FERMETURE DU FICHIER                      *
167    ****************************************************************************/
168   fclose(aFileId);
169   return aResult;
170 }