ForkMaster/smesh
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/smesh.git synced 2024-11-15 10:08:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

0021096: EDF 1729 SMESH: Create a Projection 1D-2D algorithm

Browse Source 
fix merging nodes in case 1D-2D projection of the colsed face (care
   of nodes on the seam edge)
This commit is contained in:
eap 2012-04-10 10:06:23 +00:00
parent a17fd6ce6b
commit ef93a614a6
1 changed files with 51 additions and 28 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

79
src/StdMeshers/StdMeshers_Projection_2D.cxx
View File

@ -608,9 +608,16 @@ namespace {
{ {
list< TopoDS_Edge > tgtEdges; list< TopoDS_Edge > tgtEdges;
StdMeshers_FaceSidePtr srcWire = srcWires[iW]; StdMeshers_FaceSidePtr srcWire = srcWires[iW];
TopTools_IndexedMapOfShape edgeMap; // to detect seam edges
for ( int iE = 0; iE < srcWire->NbEdges(); ++iE ) for ( int iE = 0; iE < srcWire->NbEdges(); ++iE )
{
tgtEdges.push_back( TopoDS::Edge( shape2ShapeMap( srcWire->Edge( iE )))); tgtEdges.push_back( TopoDS::Edge( shape2ShapeMap( srcWire->Edge( iE ))));
// reverse a seam edge encountered for the second time
const int oldExtent = edgeMap.Extent();
edgeMap.Add( tgtEdges.back() );
if ( oldExtent == edgeMap.Extent() )
tgtEdges.back().Reverse();
}
tgtWires[ iW ].reset( new StdMeshers_FaceSide( tgtFace, tgtEdges, tgtMesh, tgtWires[ iW ].reset( new StdMeshers_FaceSide( tgtFace, tgtEdges, tgtMesh,
/*theIsForward = */ true, /*theIsForward = */ true,
/*theIgnoreMediumNodes = */false)); /*theIgnoreMediumNodes = */false));
@ -941,8 +948,8 @@ bool StdMeshers_Projection_2D::Compute(SMESH_Mesh& theMesh, const TopoDS_Shape&
{ {
SMESH_subMesh* sm = smIt->next(); SMESH_subMesh* sm = smIt->next();
SMESHDS_SubMesh* smDS = sm->GetSubMeshDS(); SMESHDS_SubMesh* smDS = sm->GetSubMeshDS();
if ( !sm->IsMeshComputed() ) if ( smDS->NbNodes() == 0 )
break; continue;
//if ( !is1DComputed && sm->GetSubShape().ShapeType() == TopAbs_EDGE ) //if ( !is1DComputed && sm->GetSubShape().ShapeType() == TopAbs_EDGE )
//break; //break;
@ -966,7 +973,7 @@ bool StdMeshers_Projection_2D::Compute(SMESH_Mesh& theMesh, const TopoDS_Shape&
bool isOld = isOldNode( node ); bool isOld = isOldNode( node );
if ( !isOld && isSeam ) { // new node on a seam edge if ( !isOld && isSeam ) { // new node on a seam edge
if ( seamNodes.find( node ) != seamNodes.end()) if ( seamNodes.count( node ) )
continue; // node is already in the map continue; // node is already in the map
} }
@ -992,40 +999,56 @@ bool StdMeshers_Projection_2D::Compute(SMESH_Mesh& theMesh, const TopoDS_Shape&
node->GetPosition()->GetTypeOfPosition()); node->GetPosition()->GetTypeOfPosition());
} }
} }
if ( u2nodesMaps[ NEW_NODES ].size() != u2nodesMaps[ OLD_NODES ].size() ) const bool mergeNewToOld =
( u2nodesMaps[ NEW_NODES ].size() == u2nodesMaps[ OLD_NODES ].size() );
const bool mergeSeamToNew =
( u2nodesMaps[ NEW_NODES ].size() == u2nodesOnSeam.size() );
if ( !mergeNewToOld )
{ {
if ( u2nodesMaps[ NEW_NODES ].size() == 0 && // if ( u2nodesMaps[ NEW_NODES ].size() == 0 &&
sm->GetSubShape().ShapeType() == TopAbs_EDGE && // sm->GetSubShape().ShapeType() == TopAbs_EDGE &&
helper.IsDegenShape( sm->GetId() ) ) // helper.IsDegenShape( sm->GetId() ) )
// NPAL15894 (tt88bis.py) - project mesh built by NETGEN_1d_2D that // // NPAL15894 (tt88bis.py) - project mesh built by NETGEN_1d_2D that
// does not make segments/nodes on degenerated edges // // does not make segments/nodes on degenerated edges
continue; // continue;
if ( u2nodesMaps[ OLD_NODES ].size() == 0 && // if ( u2nodesMaps[ OLD_NODES ].size() == 0 &&
sm->GetSubShape().ShapeType() == TopAbs_VERTEX ) // sm->GetSubShape().ShapeType() == TopAbs_VERTEX )
// old nodes are optional on vertices in the case of 1D-2D projection // // old nodes are optional on vertices in the case of 1D-2D projection
continue; // continue;
RETURN_BAD_RESULT("Different nb of old and new nodes on shape #"<< sm->GetId() <<" "<< //RETURN_BAD_RESULT
u2nodesMaps[ OLD_NODES ].size() << " != " << MESSAGE("Different nb of old and new nodes on shape #"<< sm->GetId() <<" "<<
u2nodesMaps[ NEW_NODES ].size()); u2nodesMaps[ OLD_NODES ].size() << " != " <<
u2nodesMaps[ NEW_NODES ].size());
} }
if ( isSeam && u2nodesMaps[ OLD_NODES ].size() != u2nodesOnSeam.size() ) { if ( isSeam && !mergeSeamToNew ) {
RETURN_BAD_RESULT("Different nb of old and seam nodes " << //RETURN_BAD_RESULT
u2nodesMaps[ OLD_NODES ].size() << " != " << u2nodesOnSeam.size()); MESSAGE("Different nb of old and seam nodes " <<
u2nodesMaps[ OLD_NODES ].size() << " != " << u2nodesOnSeam.size());
} }
// Make groups of nodes to merge // Make groups of nodes to merge
u_oldNode = u2nodesMaps[ OLD_NODES ].begin(); u_oldNode = u2nodesMaps[ OLD_NODES ].begin();
u_newNode = u2nodesMaps[ NEW_NODES ].begin(); u_newNode = u2nodesMaps[ NEW_NODES ].begin();
newEnd = u2nodesMaps[ NEW_NODES ].end(); newEnd = u2nodesMaps[ NEW_NODES ].end();
u_newOnSeam = u2nodesOnSeam.begin(); u_newOnSeam = u2nodesOnSeam.begin();
for ( ; u_newNode != newEnd; ++u_newNode, ++u_oldNode ) { if ( mergeNewToOld )
groupsOfNodes.push_back( list< const SMDS_MeshNode* >() ); for ( ; u_newNode != newEnd; ++u_newNode, ++u_oldNode )
groupsOfNodes.back().push_back( u_oldNode->second ); {
groupsOfNodes.back().push_back( u_newNode->second ); groupsOfNodes.push_back( list< const SMDS_MeshNode* >() );
if ( isSeam ) groupsOfNodes.back().push_back( u_oldNode->second );
groupsOfNodes.back().push_back( (u_newOnSeam++)->second ); groupsOfNodes.back().push_back( u_newNode->second );
} if ( mergeSeamToNew )
groupsOfNodes.back().push_back( (u_newOnSeam++)->second );
}
else if ( mergeSeamToNew )
for ( ; u_newNode != newEnd; ++u_newNode, ++u_newOnSeam )
{
groupsOfNodes.push_back( list< const SMDS_MeshNode* >() );
groupsOfNodes.back().push_back( u_newNode->second );
groupsOfNodes.back().push_back( u_newOnSeam->second );
}
} }
// Merge // Merge