- fprintf(stdout,"\n(**************************)\n");
- fprintf(stdout,"(* ELEMENTS DU MAILLAGE : *)\n");
- fprintf(stdout,"(**************************)");
-
- fprintf(stdout,"%d elements\n",nbCells);
-
- for (i=0; i<nbCells; i++) {
- fscanf(myFileId,"%d %d",&intNumMaille,&ValElement);
- Degre=abs(ValElement/100);
- nbNoeuds=ValElement-(Degre*100);
-
- // Recuperation des noeuds de la maille
- for (j=0; j<nbNoeuds; j++) {
- fscanf(myFileId,"%d",&NoeudMaille);
- NoeudsMaille[j]=NoeudMaille;
- }
-
- // Analyse des cas de cellules
- switch (ValElement) {
- case 102 : ;
- case 103 : ;
- {
- ValDegre=3;
- nbNoeuds=2;
- ok = myMesh->AddEdgeWithID(NoeudsMaille[0],NoeudsMaille[1],intNumMaille);
- break;
- }
- case 204 : ;
- case 208 : ;
- {
- ValDegre=9;
- nbNoeuds=4;
- ok = myMesh->AddFaceWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],intNumMaille);
- break;
- }
- case 203 : ;
- case 206 : ;
- {
- ValDegre=5;
- nbNoeuds=3;
- ok = myMesh->AddFaceWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],intNumMaille);
- break;
- }
- case 308 : ;
- case 320 : ;
- {
- ValDegre=12;
- nbNoeuds=8;
- if (ValElement==320) {
- //A voir, correspondance VTK
- NoeudsMaille[4]=NoeudsMaille[8];
- NoeudsMaille[5]=NoeudsMaille[9];
- NoeudsMaille[6]=NoeudsMaille[10];
- NoeudsMaille[7]=NoeudsMaille[11];
- }
- ok = myMesh->AddVolumeWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],NoeudsMaille[4],NoeudsMaille[5],NoeudsMaille[6],NoeudsMaille[7],intNumMaille);
- break;
- }
- case 304 : ;
- case 310 : ;
- {
- ValDegre=10;
- nbNoeuds=4;
- if (ValElement==310)
- NoeudsMaille[3]=NoeudsMaille[6];
- ok = myMesh->AddVolumeWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],intNumMaille);
- break;
- }
- case 306 : ;
- case 315 : ;
- {
- ValDegre=12;
- nbNoeuds=8;
- if (ValElement==315) {
- NoeudsMaille[3]=NoeudsMaille[6];
- NoeudsMaille[4]=NoeudsMaille[7];
- NoeudsMaille[5]=NoeudsMaille[8];
- }
- NoeudsMaille[7]=NoeudsMaille[5];
- NoeudsMaille[6]=NoeudsMaille[5];
- NoeudsMaille[5]=NoeudsMaille[4];
- NoeudsMaille[4]=NoeudsMaille[3];
- NoeudsMaille[3]=NoeudsMaille[2];
- ok = myMesh->AddVolumeWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],NoeudsMaille[4],NoeudsMaille[5],intNumMaille);
- break;
- }
-
- }
- }
+ fprintf(stdout, "\n(**************************)\n");
+ fprintf(stdout, "(* ELEMENTS DU MAILLAGE : *)\n");
+ fprintf(stdout, "(**************************)");
+
+ fprintf(stdout, "%d elements\n", nbCells);
+
+ for (i = 0; i < nbCells; i++)
+ {
+ fscanf(myFileId, "%d %d", &intNumMaille, &ValElement);
+ Degre = abs(ValElement / 100);
+ nbNoeuds = ValElement - (Degre * 100);
+
+ // Recuperation des noeuds de la maille
+ for (j = 0; j < nbNoeuds; j++)
+ {
+ fscanf(myFileId, "%d", &NoeudMaille);
+ NoeudsMaille[j] = NoeudMaille;
+ }
+
+ // Analyse des cas de cellules
+ switch (ValElement)
+ {
+ case 102:;
+ case 103:;
+ {
+ ValDegre = 3;
+ nbNoeuds = 2;
+ ok = myMesh->AddEdgeWithID(NoeudsMaille[0], NoeudsMaille[1],
+ intNumMaille);
+ break;
+ }
+ case 204:;
+ case 208:;
+ {
+ ValDegre = 9;
+ nbNoeuds = 4;
+ ok = myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
+ NoeudsMaille[2], NoeudsMaille[3], intNumMaille);
+ break;
+ }
+ case 203:;
+ case 206:;
+ {
+ ValDegre = 5;
+ nbNoeuds = 3;
+ ok = myMesh->AddFaceWithID(NoeudsMaille[0], NoeudsMaille[1],
+ NoeudsMaille[2], intNumMaille);
+ break;
+ }
+ case 308:;
+ case 320:;
+ {
+ ValDegre = 12;
+ nbNoeuds = 8;
+ if (ValElement == 320)
+ {
+ //A voir, correspondance VTK
+ NoeudsMaille[4] = NoeudsMaille[8];
+ NoeudsMaille[5] = NoeudsMaille[9];
+ NoeudsMaille[6] = NoeudsMaille[10];
+ NoeudsMaille[7] = NoeudsMaille[11];
+ }
+ ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
+ NoeudsMaille[2], NoeudsMaille[3], NoeudsMaille[4],
+ NoeudsMaille[5], NoeudsMaille[6], NoeudsMaille[7],
+ intNumMaille);
+ break;
+ }
+ case 304:;
+ case 310:;
+ {
+ ValDegre = 10;
+ nbNoeuds = 4;
+ if (ValElement == 310)
+ NoeudsMaille[3] = NoeudsMaille[6];
+ ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
+ NoeudsMaille[2], NoeudsMaille[3], intNumMaille);
+ break;
+ }
+ case 306:;
+ case 315:;
+ {
+ ValDegre = 12;
+ nbNoeuds = 8;
+ if (ValElement == 315)
+ {
+ NoeudsMaille[3] = NoeudsMaille[6];
+ NoeudsMaille[4] = NoeudsMaille[7];
+ NoeudsMaille[5] = NoeudsMaille[8];
+ }
+ NoeudsMaille[7] = NoeudsMaille[5];
+ NoeudsMaille[6] = NoeudsMaille[5];
+ NoeudsMaille[5] = NoeudsMaille[4];
+ NoeudsMaille[4] = NoeudsMaille[3];
+ NoeudsMaille[3] = NoeudsMaille[2];
+ ok = myMesh->AddVolumeWithID(NoeudsMaille[0], NoeudsMaille[1],
+ NoeudsMaille[2], NoeudsMaille[3], NoeudsMaille[4],
+ NoeudsMaille[5], intNumMaille);
+ break;
+ }
+
+ }
+ }