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

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