ForkMaster/netgenplugin
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

0021405: EDF 1965 SMESH:

Browse Source 
In fillNgMesh(), add a node on VERTEX only if a VERTEX is shared by a not meshed EDGE
This commit is contained in:
eap 2011-11-01 10:25:17 +00:00
parent 9c7a8090fd
commit 20b845d011
1 changed files with 17 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
src/NETGENPlugin/NETGENPlugin_Mesher.cxx
View File

@ -860,9 +860,23 @@ bool NETGENPlugin_Mesher::fillNgMesh(const netgen::OCCGeometry& occgeom,
case TopAbs_VERTEX: { // VERTEX
// --------------------------
SMDS_NodeIteratorPtr nodeIt = smDS->GetNodes();
if ( nodeIt->more() )
ngNodeId( nodeIt->next(), ngMesh, nodeNgIdMap );
// issue 21405. Add node only if a VERTEX is shared by a not meshed EDGE,
// else netgen removes a free node and nodeVector becomes invalid
PShapeIteratorPtr ansIt = helper.GetAncestors( sm->GetSubShape(),
*sm->GetFather(),
TopAbs_EDGE );
bool toAdd = false;
while ( const TopoDS_Shape* e = ansIt->next() )
{
SMESH_subMesh* eSub = helper.GetMesh()->GetSubMesh( *e );
if (( toAdd = eSub->IsEmpty() )) break;
}
if ( toAdd )
{
SMDS_NodeIteratorPtr nodeIt = smDS->GetNodes();
if ( nodeIt->more() )
ngNodeId( nodeIt->next(), ngMesh, nodeNgIdMap );
}
break;
}
default:;