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=119f8f9a0efc023b73656c18763c0b6d759efa79;hp=a86d1b2954c189c842503a854a0ff485098ca8d2;hb=499f29d24922cec66e41b41a0039a954993bc6df;hpb=c98d9fcd7f02c1f1f5c24dd3e709ed75228d66c4 diff --git a/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx b/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx index a86d1b295..119f8f9a0 100644 --- a/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx +++ b/src/DriverDAT/DriverDAT_W_SMDS_Mesh.cxx @@ -1,4 +1,4 @@ -// Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE +// Copyright (C) 2007-2022 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 @@ -32,16 +32,30 @@ using namespace std; +//================================================================================ +/*! + * \brief Write mesh data to a file + * \param [in] renumber - if true, renumber nodes and cell starting from 1 + * \return Driver_Mesh::Status - Ok or not + */ +//================================================================================ + Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() { Kernel_Utils::Localizer loc; Status aResult = DRS_OK; - int nbNodes, nbCells; + smIdType nbNodes, nbCells; +#if defined(WIN32) && defined(UNICODE) + std::wstring file2Read = Kernel_Utils::utf8_decode_s(myFile); + FILE* aFileId = _wfopen(file2Read.c_str(), L"w+"); +#else char *file2Read = (char *)myFile.c_str(); FILE* aFileId = fopen(file2Read, "w+"); - if ( !aFileId ) { +#endif + if ( !aFileId ) + { fprintf(stderr, ">> ERREUR : ouverture du fichier %s \n", file2Read); return DRS_FAIL; } @@ -54,7 +68,7 @@ Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() nbNodes = myMesh->NbNodes(); /* Combien de mailles, faces ou aretes ? */ - int nb_of_edges, nb_of_faces, nb_of_volumes; + smIdType nb_of_edges, nb_of_faces, nb_of_volumes; nb_of_edges = myMesh->NbEdges(); nb_of_faces = myMesh->NbFaces(); nb_of_volumes = myMesh->NbVolumes(); @@ -64,67 +78,111 @@ Driver_Mesh::Status DriverDAT_W_SMDS_Mesh::Perform() SCRUTE(nb_of_volumes); //fprintf(stdout, "%d %d\n", nbNodes, nbCells); - fprintf(aFileId, "%d %d\n", nbNodes, nbCells); + fprintf(aFileId, "%ld %ld\n", static_cast< long >( nbNodes ), static_cast< long >( nbCells )); /**************************************************************************** * ECRITURE DES NOEUDS * ****************************************************************************/ + std::vector< size_t > nodeNumByID; + if ( myRenumber && myMesh->HasNumerationHoles() ) + nodeNumByID.resize( myMesh->MaxNodeID() + 1 ); + + smIdType num; SMDS_NodeIteratorPtr itNodes=myMesh->nodesIterator(); - while(itNodes->more()){ + for ( num = 1; itNodes->more(); ++num ) + { const SMDS_MeshNode * node = itNodes->next(); - fprintf(aFileId, "%d %.14e %.14e %.14e\n", node->GetID(), node->X(), node->Y(), node->Z()); + if ( !myRenumber ) + num = node->GetID(); + + fprintf(aFileId, "%d %.14e %.14e %.14e\n", num, node->X(), node->Y(), node->Z()); + + if ( !nodeNumByID.empty() ) + nodeNumByID[ node->GetID() ] = num; } - + /**************************************************************************** * ECRITURE DES ELEMENTS * ****************************************************************************/ /* Ecriture des connectivites, noms, numeros des mailles */ - - SMDS_EdgeIteratorPtr itEdges=myMesh->edgesIterator(); - while(itEdges->more()){ + + num = 1; + for ( SMDS_EdgeIteratorPtr itEdges = myMesh->edgesIterator(); itEdges->more(); ++num ) + { const SMDS_MeshElement * elem = itEdges->next(); - switch (elem->NbNodes()) { - case 2: - fprintf(aFileId, "%d %d ", elem->GetID(), 102); - break; - case 3: - fprintf(aFileId, "%d %d ", elem->GetID(), 103); - break; + if ( !myRenumber ) + num = elem->GetID(); + + fprintf(aFileId, "%d %d ", num, 100 + elem->NbNodes()); + + for ( SMDS_ElemIteratorPtr it = elem->nodesIterator(); it->more(); ) + { + smIdType nodeID = it->next()->GetID(); + if ( !nodeNumByID.empty() ) + nodeID = nodeNumByID[ nodeID ]; + fprintf(aFileId, "%ld ", static_cast< long >( nodeID )); } - SMDS_ElemIteratorPtr it=elem->nodesIterator(); - while(it->more()) - fprintf(aFileId, "%d ", it->next()->GetID()); fprintf(aFileId, "\n"); } - - SMDS_FaceIteratorPtr itFaces=myMesh->facesIterator(); - while(itFaces->more()){ + + for ( SMDS_FaceIteratorPtr itFaces = myMesh->facesIterator(); itFaces->more(); ++num ) + { 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()) - fprintf(aFileId, "%d ", it->next()->GetID()); + if ( !myRenumber ) + num = elem->GetID(); + + fprintf(aFileId, "%d %d ", num, (elem->IsPoly() ? 400 : 200 ) + elem->NbNodes() ); + + for( SMDS_ElemIteratorPtr it = elem->nodesIterator(); it->more(); ) + { + smIdType nodeID = it->next()->GetID(); + if ( !nodeNumByID.empty() ) + nodeID = nodeNumByID[ nodeID ]; + fprintf(aFileId, "%ld ", static_cast< long >( nodeID )); + } fprintf(aFileId, "\n"); } - SMDS_VolumeIteratorPtr itVolumes=myMesh->volumesIterator(); - while(itVolumes->more()){ + + const SMDS_MeshVolume* v; + for ( SMDS_VolumeIteratorPtr itVolumes=myMesh->volumesIterator(); itVolumes->more(); ++num ) + { const SMDS_MeshElement * elem = itVolumes->next(); + if ( !myRenumber ) + num = elem->GetID(); + if ( elem->IsPoly() ) - fprintf(aFileId, "%d %d ", elem->GetID(), 500+elem->NbNodes()); + { + fprintf(aFileId, "%d %d ", num, 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()) - fprintf(aFileId, "%d ", it->next()->GetID()); + { + fprintf(aFileId, "%d %d ", num, 300 + elem->NbNodes()); + } + + for( SMDS_ElemIteratorPtr it = elem->nodesIterator(); it->more(); ) + { + smIdType nodeID = it->next()->GetID(); + if ( !nodeNumByID.empty() ) + nodeID = nodeNumByID[ nodeID ]; + fprintf(aFileId, "%ld ", static_cast< long >( nodeID )); + } fprintf(aFileId, "\n"); } - + fclose(aFileId); return aResult;