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

PAL20462 4.x Computation of the mesh via NETGEN is failed.

Browse Source 
fix PrepareOCCgeometry() to add subshape with right orientation
This commit is contained in:
eap 2008-10-03 10:36:29 +00:00
parent c9722be8d0
commit 58bcf172e9
1 changed files with 12 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
src/NETGENPlugin/NETGENPlugin_Mesher.cxx
View File

@ -214,13 +214,6 @@ void NETGENPlugin_Mesher::PrepareOCCgeometry(netgen::OCCGeometry& occgeo,
occgeo.changed = 1;
//occgeo.BuildFMap();
//rnv to fix 19978 issue
bool includeSelf = true;
if ( shape.ShapeType() == TopAbs_FACE) {
occgeo.fmap.Add( shape );
includeSelf = false;
}
// fill maps of shapes of occgeo with not yet meshed subshapes
// get root submeshes
@ -237,16 +230,22 @@ void NETGENPlugin_Mesher::PrepareOCCgeometry(netgen::OCCGeometry& occgeo,
list< SMESH_subMesh* >::iterator rootIt = rootSM.begin(), rootEnd = rootSM.end();
for ( ; rootIt != rootEnd; ++rootIt ) {
SMESH_subMesh * root = *rootIt;
SMESH_subMeshIteratorPtr smIt = root->getDependsOnIterator(/*includeSelf=*/includeSelf,
SMESH_subMeshIteratorPtr smIt = root->getDependsOnIterator(/*includeSelf=*/true,
/*complexShapeFirst=*/true);
// to find a right orientation of subshapes (PAL20462)
TopTools_IndexedMapOfShape subShapes;
TopExp::MapShapes(root->GetSubShape(), subShapes);
while ( smIt->more() ) {
SMESH_subMesh* sm = smIt->next();
if ( sm->IsEmpty() ) {
switch ( sm->GetSubShape().ShapeType() ) {
case TopAbs_FACE : occgeo.fmap.Add( sm->GetSubShape() ); break;
case TopAbs_EDGE : occgeo.emap.Add( sm->GetSubShape() ); break;
case TopAbs_VERTEX: occgeo.vmap.Add( sm->GetSubShape() ); break;
case TopAbs_SOLID :occgeo.somap.Add( sm->GetSubShape() ); break;
TopoDS_Shape shape = sm->GetSubShape();
if ( shape.ShapeType() != TopAbs_VERTEX )
shape = subShapes( subShapes.FindIndex( shape ));// - shape->index->oriented shape
switch ( shape.ShapeType() ) {
case TopAbs_FACE : occgeo.fmap.Add( shape ); break;
case TopAbs_EDGE : occgeo.emap.Add( shape ); break;
case TopAbs_VERTEX: occgeo.vmap.Add( shape ); break;
case TopAbs_SOLID :occgeo.somap.Add( shape ); break;
default:;
}
}