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 }