From: eap <eap@opencascade.com>
Date: Fri, 3 Oct 2008 10:38:21 +0000 (+0000)
Subject: 19978PAL20462 4.x Computation of the mesh via NETGEN is failed.
X-Git-Tag: V5_1_0_for_TC1~3
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=ce350d72ee1023518a3d78d0e0db7f8ae1f99729;p=plugins%2Fnetgenplugin.git

19978PAL20462 4.x Computation of the mesh via NETGEN is failed.
 fix PrepareOCCgeometry() to add subshape with right orientation
---

diff --git a/src/NETGENPlugin/NETGENPlugin_Mesher.cxx b/src/NETGENPlugin/NETGENPlugin_Mesher.cxx
index 6389d3c..f9de58a 100644
--- a/src/NETGENPlugin/NETGENPlugin_Mesher.cxx
+++ b/src/NETGENPlugin/NETGENPlugin_Mesher.cxx
@@ -212,15 +212,8 @@ void NETGENPlugin_Mesher::PrepareOCCgeometry(netgen::OCCGeometry&     occgeo,
 
   occgeo.shape = shape;
   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
 
   // 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:;
         }
       }