]> SALOME platform Git repositories - modules/smesh.git/blob - src/DriverDAT/DriverDAT_R_SMDS_Mesh.cxx
Salome HOME
Copyright update 2021
[modules/smesh.git] / src / DriverDAT / DriverDAT_R_SMDS_Mesh.cxx
1 // Copyright (C) 2007-2021  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 <stdio.h>
24
25 #include "DriverDAT_R_SMDS_Mesh.h"
26 #include "SMDS_Mesh.hxx"
27
28 #include "utilities.h"
29
30 #include <Basics_Utils.hxx>
31
32 using namespace std;
33
34 Driver_Mesh::Status DriverDAT_R_SMDS_Mesh::Perform()
35 {
36   Kernel_Utils::Localizer loc;
37   Status aResult = DRS_OK;
38
39   int i, j;
40   int nbNodes, nbCells;
41   int intNumPoint;
42   float coordX, coordY, coordZ;
43   int nbNoeuds;
44   
45   int intNumMaille, Degre;
46   int ValElement;
47   int NoeudsMaille[20];
48   int NoeudMaille;
49   
50   /****************************************************************************
51    *                      OUVERTURE DU FICHIER EN LECTURE                      *
52    ****************************************************************************/
53 #if defined(WIN32) && defined(UNICODE)
54   std::wstring file2Read = Kernel_Utils::utf8_decode_s(myFile);
55   FILE* aFileId = _wfopen(file2Read.c_str(), L"r");
56 #else
57   char *file2Read = (char *)myFile.c_str();
58   FILE* aFileId = fopen(file2Read, "r");
59 #endif
60   if ( !aFileId ) {
61     fprintf(stderr, ">> ERREUR : ouverture du fichier %s \n", file2Read);
62     return DRS_FAIL;
63   }
64   
65   fscanf(aFileId, "%d %d\n", &nbNodes, &nbCells);
66   
67   /****************************************************************************
68    *                       LECTURE DES NOEUDS                                  *
69    ****************************************************************************/
70   fprintf(stdout, "\n(************************)\n");
71   fprintf(stdout, "(* NOEUDS DU MAILLAGE : *)\n");
72   fprintf(stdout, "(************************)\n");
73   
74   for (i = 0; i < nbNodes; i++){
75     fscanf(aFileId, "%d %e %e %e\n", &intNumPoint, &coordX, &coordY, &coordZ);
76     myMesh->AddNodeWithID(coordX, coordY, coordZ, intNumPoint);
77   }
78   
79   fprintf(stdout, "%d noeuds\n", myMesh->NbNodes());
80   /****************************************************************************
81    *                       LECTURE DES ELEMENTS                                *
82    ****************************************************************************/
83   fprintf(stdout, "\n(**************************)\n");
84   fprintf(stdout, "(* ELEMENTS DU MAILLAGE : *)\n");
85   fprintf(stdout, "(**************************)");
86   
87   fprintf(stdout, "%d elements\n", nbCells);
88   
89   for (i = 0; i < nbCells; i++) {
90     fscanf(aFileId, "%d %d", &intNumMaille, &ValElement);
91     Degre = abs(ValElement / 100);
92     nbNoeuds = ValElement - (Degre * 100);
93     
94     // Recuperation des noeuds de la maille
95     for (j = 0; j < nbNoeuds; j++) {
96       fscanf(aFileId, "%d", &NoeudMaille);
97       NoeudsMaille[j] = NoeudMaille;
98     }
99     
100     // Analyse des cas de cellules
101     switch (ValElement) {
102     case 102:
103     case 103:
104       nbNoeuds = 2;
105       myMesh->AddEdgeWithID(NoeudsMaille[0], NoeudsMaille[1], 
106                                  intNumMaille);
107       break;
108     case 204:
109     case 208:
110       nbNoeuds = 4;
111       myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
112                                  NoeudsMaille[2], NoeudsMaille[3], 
113                                  intNumMaille);
114       break;
115     case 203:
116     case 206:
117       nbNoeuds = 3;
118       myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
119                                  NoeudsMaille[2], intNumMaille);
120       break;
121     case 308:
122     case 320:
123       nbNoeuds = 8;
124       if (ValElement == 320){
125         //A voir, correspondance VTK
126         NoeudsMaille[4] = NoeudsMaille[8];
127         NoeudsMaille[5] = NoeudsMaille[9];
128         NoeudsMaille[6] = NoeudsMaille[10];
129         NoeudsMaille[7] = NoeudsMaille[11];
130       }
131       myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
132                                    NoeudsMaille[2], NoeudsMaille[3], 
133                                    NoeudsMaille[4], NoeudsMaille[5], 
134                                    NoeudsMaille[6], NoeudsMaille[7],
135                                    intNumMaille);
136       break;
137     case 304:
138     case 310:
139       nbNoeuds = 4;
140       if (ValElement == 310)
141         NoeudsMaille[3] = NoeudsMaille[6];
142       myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
143                                    NoeudsMaille[2], NoeudsMaille[3], 
144                                    intNumMaille);
145       break;
146     case 306:
147     case 315:
148       nbNoeuds = 8;
149       if (ValElement == 315) {
150         NoeudsMaille[3] = NoeudsMaille[6];
151         NoeudsMaille[4] = NoeudsMaille[7];
152         NoeudsMaille[5] = NoeudsMaille[8];
153       }
154       NoeudsMaille[7] = NoeudsMaille[5];
155       NoeudsMaille[6] = NoeudsMaille[5];
156       NoeudsMaille[5] = NoeudsMaille[4];
157       NoeudsMaille[4] = NoeudsMaille[3];
158       NoeudsMaille[3] = NoeudsMaille[2];
159       myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
160                                    NoeudsMaille[2], NoeudsMaille[3], 
161                                    NoeudsMaille[4], NoeudsMaille[5], 
162                                    intNumMaille);
163                                 break;
164     }
165   }
166   /****************************************************************************
167    *                      FERMETURE DU FICHIER                      *
168    ****************************************************************************/
169   fclose(aFileId);
170   return aResult;
171 }