X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FDriverDAT%2FDriverDAT_W_SMDS_Mesh.cxx;h=403dad4f5634748cf6a0d38af896efa95deb158b;hp=400daa548c30b9d97c3a350ed6303aac1dd5ae9a;hb=d8bc11015235c68631ac85d7e8cf278edc9809d1;hpb=b0a908c0d20341651771d0249fb10882f54b2aad diff --git a/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx b/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx index 400daa548..403dad4f5 100644 --- a/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx +++ b/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx @@ -1,4 +1,4 @@ -// Copyright (C) 2007-2014 CEA/DEN, EDF R&D, OPEN CASCADE +// Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE // // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS @@ -38,11 +38,11 @@ Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() Status aResult = DRS_OK; int nbNodes, nbCells; - //int i; - + char *file2Read = (char *)myFile.c_str(); FILE* aFileId = fopen(file2Read, "w+"); - if (aFileId < 0) { + if ( !aFileId ) + { fprintf(stderr, ">> ERREUR : ouverture du fichier %s \n", file2Read); return DRS_FAIL; } @@ -50,12 +50,12 @@ Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() /**************************************************************************** * NOMBRES D'OBJETS * ****************************************************************************/ - + /* Combien de noeuds ? */ nbNodes = myMesh->NbNodes(); - + /* Combien de mailles, faces ou aretes ? */ - int /*nb_of_nodes,*/ nb_of_edges, nb_of_faces, nb_of_volumes; + int nb_of_edges, nb_of_faces, nb_of_volumes; nb_of_edges = myMesh->NbEdges(); nb_of_faces = myMesh->NbFaces(); nb_of_volumes = myMesh->NbVolumes(); @@ -63,29 +63,32 @@ Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() SCRUTE(nb_of_edges); SCRUTE(nb_of_faces); SCRUTE(nb_of_volumes); - - fprintf(stdout, "%d %d\n", nbNodes, nbCells); + + //fprintf(stdout, "%d %d\n", nbNodes, nbCells); fprintf(aFileId, "%d %d\n", nbNodes, nbCells); - + /**************************************************************************** * ECRITURE DES NOEUDS * ****************************************************************************/ - + SMDS_NodeIteratorPtr itNodes=myMesh->nodesIterator(); - while(itNodes->more()){ + while(itNodes->more()) + { const SMDS_MeshNode * node = itNodes->next(); - fprintf(aFileId, "%d %e %e %e\n", node->GetID(), node->X(), node->Y(), node->Z()); + fprintf(aFileId, "%d %.14e %.14e %.14e\n", node->GetID(), node->X(), node->Y(), node->Z()); } - + /**************************************************************************** * ECRITURE DES ELEMENTS * ****************************************************************************/ /* Ecriture des connectivites, noms, numeros des mailles */ - + SMDS_EdgeIteratorPtr itEdges=myMesh->edgesIterator(); - while(itEdges->more()){ + while(itEdges->more()) + { const SMDS_MeshElement * elem = itEdges->next(); - switch (elem->NbNodes()) { + switch (elem->NbNodes()) + { case 2: fprintf(aFileId, "%d %d ", elem->GetID(), 102); break; @@ -94,38 +97,56 @@ Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() break; } SMDS_ElemIteratorPtr it=elem->nodesIterator(); - while(it->more()) + while(it->more()) fprintf(aFileId, "%d ", it->next()->GetID()); fprintf(aFileId, "\n"); } - + SMDS_FaceIteratorPtr itFaces=myMesh->facesIterator(); - while(itFaces->more()){ + while(itFaces->more()) + { const SMDS_MeshElement * elem = itFaces->next(); if ( elem->IsPoly() ) fprintf(aFileId, "%d %d ", elem->GetID(), 400+elem->NbNodes()); else fprintf(aFileId, "%d %d ", elem->GetID(), 200+elem->NbNodes()); SMDS_ElemIteratorPtr it=elem->nodesIterator(); - while(it->more()) + while(it->more()) fprintf(aFileId, "%d ", it->next()->GetID()); fprintf(aFileId, "\n"); } SMDS_VolumeIteratorPtr itVolumes=myMesh->volumesIterator(); - while(itVolumes->more()){ + const SMDS_MeshVolume* v; + while(itVolumes->more()) + { const SMDS_MeshElement * elem = itVolumes->next(); if ( elem->IsPoly() ) + { fprintf(aFileId, "%d %d ", elem->GetID(), 500+elem->NbNodes()); + + if (( v = myMesh->DownCast< SMDS_MeshVolume >( elem ))) + { + std::vector quant = v->GetQuantities(); + if ( !quant.empty() ) + { + fprintf(aFileId, "%d %d ", (int)quant.size(), quant[0]); + for ( size_t i = 1; i < quant.size(); ++i ) + fprintf(aFileId, "%d ", quant[i]); + } + } + } else + { fprintf(aFileId, "%d %d ", elem->GetID(), 300+elem->NbNodes()); + } SMDS_ElemIteratorPtr it=elem->nodesIterator(); - while(it->more()) + while(it->more()) fprintf(aFileId, "%d ", it->next()->GetID()); fprintf(aFileId, "\n"); } - + fclose(aFileId); return aResult;