src/DriverDAT/DriverDAT_R_SMDS_Mesh.cxx

   1 //  SMESH DriverDAT : driver to read and write 'dat' files
   2 //
   3 //  Copyright (C) 2003  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.
  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.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
  21 //
  22 //
  23 //
  24 //  File   : DriverDAT_R_SMDS_Mesh.cxx
  25 //  Module : SMESH
  26
  27 using namespace std;
  28 #include "DriverDAT_R_SMDS_Mesh.h"
  29
  30 #include "utilities.h"
  31
  32 DriverDAT_R_SMDS_Mesh::DriverDAT_R_SMDS_Mesh() {
  33 ;
  34 }
  35
  36 DriverDAT_R_SMDS_Mesh::~DriverDAT_R_SMDS_Mesh() {
  37 ;
  38 }
  39
  40 void DriverDAT_R_SMDS_Mesh::SetMesh(Handle(SMDS_Mesh)& aMesh) {
  41   myMesh = aMesh;
  42 }
  43
  44 void DriverDAT_R_SMDS_Mesh::SetFile(string aFile) {
  45   myFile = aFile;
  46 }
  47
  48 void DriverDAT_R_SMDS_Mesh::SetFileId(FILE* aFileId) {
  49   myFileId = aFileId;
  50 }
  51
  52 void DriverDAT_R_SMDS_Mesh::SetMeshId(int aMeshId) {
  53   myMeshId = aMeshId;
  54 }
  55
  56 void DriverDAT_R_SMDS_Mesh::Add() {
  57   ;
  58 }
  59
  60 void DriverDAT_R_SMDS_Mesh::Read() {
  61
  62   int i,j;
  63   int nbNodes,nbCells;
  64   int intNumPoint;
  65   float coordX, coordY, coordZ;
  66   int nbNoeuds;
  67
  68   int intNumMaille,Degre;
  69   int ValElement;
  70   int ValDegre;
  71   int NoeudsMaille[20];
  72   int NoeudMaille;
  73
  74   bool ok;
  75
  76   MESSAGE("in DriverDAT_R_SMDS_Mesh::Read()");
  77   /****************************************************************************
  78   *                      OUVERTURE DU FICHIER EN LECTURE                      *
  79   ****************************************************************************/
  80   char* file2Read = (char*)myFile.c_str();
  81   myFileId = fopen(file2Read,"r");
  82   if (myFileId < 0)
  83     {
  84       fprintf(stderr,">> ERREUR : ouverture du fichier %s \n",file2Read);
  85       exit(EXIT_FAILURE);
  86     }
  87
  88   fscanf(myFileId,"%d %d\n",&nbNodes,&nbCells);
  89
  90   /****************************************************************************
  91   *                       LECTURE DES NOEUDS                                  *
  92   ****************************************************************************/
  93   fprintf(stdout,"\n(************************)\n");
  94   fprintf(stdout,"(* NOEUDS DU MAILLAGE : *)\n");
  95   fprintf(stdout,"(************************)\n");
  96
  97   for (i=0;i<nbNodes;i++) {
  98     fscanf(myFileId,"%d %e %e %e\n",&intNumPoint,&coordX,&coordY,&coordZ);
  99     ok = myMesh->AddNodeWithID(coordX,coordY,coordZ,intNumPoint);
 100   }
 101
 102   fprintf(stdout,"%d noeuds\n",myMesh->NbNodes());
 103   /****************************************************************************
 104   *                       LECTURE DES ELEMENTS                                *
 105   ****************************************************************************/
 106   fprintf(stdout,"\n(**************************)\n");
 107   fprintf(stdout,"(* ELEMENTS DU MAILLAGE : *)\n");
 108   fprintf(stdout,"(**************************)");
 109
 110   fprintf(stdout,"%d elements\n",nbCells);
 111
 112   for (i=0; i<nbCells; i++) {
 113     fscanf(myFileId,"%d %d",&intNumMaille,&ValElement);
 114     Degre=abs(ValElement/100);
 115     nbNoeuds=ValElement-(Degre*100);
 116
 117     // Recuperation des noeuds de la maille
 118     for (j=0; j<nbNoeuds; j++) {
 119       fscanf(myFileId,"%d",&NoeudMaille);
 120       NoeudsMaille[j]=NoeudMaille;
 121     }
 122
 123     // Analyse des cas de cellules
 124     switch (ValElement) {
 125     case 102 : ;
 126     case 103 : ;
 127       {
 128       ValDegre=3;
 129       nbNoeuds=2;
 130       ok = myMesh->AddEdgeWithID(NoeudsMaille[0],NoeudsMaille[1],intNumMaille);
 131       break;
 132     }
 133     case 204 : ;
 134     case 208 : ;
 135       {
 136       ValDegre=9;
 137       nbNoeuds=4;
 138       ok = myMesh->AddFaceWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],intNumMaille);
 139       break;
 140     }
 141     case 203 : ;
 142     case 206 : ;
 143       {
 144       ValDegre=5;
 145       nbNoeuds=3;
 146       ok = myMesh->AddFaceWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],intNumMaille);
 147       break;
 148     }
 149     case 308 : ;
 150     case 320 : ;
 151       {
 152       ValDegre=12;
 153       nbNoeuds=8;
 154       if (ValElement==320) {
 155         //A voir, correspondance VTK
 156         NoeudsMaille[4]=NoeudsMaille[8];
 157         NoeudsMaille[5]=NoeudsMaille[9];
 158         NoeudsMaille[6]=NoeudsMaille[10];
 159         NoeudsMaille[7]=NoeudsMaille[11];
 160       }
 161       ok = myMesh->AddVolumeWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],NoeudsMaille[4],NoeudsMaille[5],NoeudsMaille[6],NoeudsMaille[7],intNumMaille);
 162      break;
 163     }
 164     case 304 : ;
 165     case 310 : ;
 166       {
 167       ValDegre=10;
 168       nbNoeuds=4;
 169       if (ValElement==310)
 170         NoeudsMaille[3]=NoeudsMaille[6];
 171       ok = myMesh->AddVolumeWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],intNumMaille);
 172      break;
 173     }
 174     case 306 : ;
 175     case 315 : ;
 176       {
 177       ValDegre=12;
 178       nbNoeuds=8;
 179       if (ValElement==315) {
 180         NoeudsMaille[3]=NoeudsMaille[6];
 181         NoeudsMaille[4]=NoeudsMaille[7];
 182         NoeudsMaille[5]=NoeudsMaille[8];
 183       }
 184       NoeudsMaille[7]=NoeudsMaille[5];
 185       NoeudsMaille[6]=NoeudsMaille[5];
 186       NoeudsMaille[5]=NoeudsMaille[4];
 187       NoeudsMaille[4]=NoeudsMaille[3];
 188       NoeudsMaille[3]=NoeudsMaille[2];
 189       ok = myMesh->AddVolumeWithID(NoeudsMaille[0],NoeudsMaille[1],NoeudsMaille[2],NoeudsMaille[3],NoeudsMaille[4],NoeudsMaille[5],intNumMaille);
 190       break;
 191     }
 192
 193     }
 194   }
 195
 196   /****************************************************************************
 197   *                      FERMETURE DU FICHIER                      *
 198   ****************************************************************************/
 199   fclose(myFileId);
 200
 201 }