From 58bcf172e925fb697d9e2b0dcc1e632c68421b3a Mon Sep 17 00:00:00 2001 From: eap Date: Fri, 3 Oct 2008 10:36:29 +0000 Subject: [PATCH] PAL20462 4.x Computation of the mesh via NETGEN is failed. fix PrepareOCCgeometry() to add subshape with right orientation --- src/NETGENPlugin/NETGENPlugin_Mesher.cxx | 25 ++++++++++++------------ 1 file changed, 12 insertions(+), 13 deletions(-) diff --git a/src/NETGENPlugin/NETGENPlugin_Mesher.cxx b/src/NETGENPlugin/NETGENPlugin_Mesher.cxx index 6dfafeb..f9de58a 100644 --- a/src/NETGENPlugin/NETGENPlugin_Mesher.cxx +++ b/src/NETGENPlugin/NETGENPlugin_Mesher.cxx @@ -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:; } } -- 2.39.2