0021059: EDF 1666 SMESH: 515rc1 Regression in Netgen

   Make NETGEN 2d work with pre-computed sub-meshes similarly to work w/o them

diff --git a/src/NETGENPlugin/NETGENPlugin_Mesher.cxx b/src/NETGENPlugin/NETGENPlugin_Mesher.cxx
index 76093a6c1263840f52dd87f58a1f7749e7a0ccaa..89e9ca829a7e454f4310d7d7b5aa7c011c3082d1 100644 (file)
--- a/src/NETGENPlugin/NETGENPlugin_Mesher.cxx
+++ b/src/NETGENPlugin/NETGENPlugin_Mesher.cxx
@@ -493,6 +493,10 @@ bool NETGENPlugin_Mesher::fillNgMesh(const netgen::OCCGeometry&     occgeom,
           seg.si = faceID;                   // = geom.fmap.FindIndex (face);
           seg.edgenr = ngMesh.GetNSeg() + 1; // segment id
           ngMesh.AddSegment (seg);
+
+          netgen::Point3d ngP1(p1.node->X(), p1.node->Y(), p1.node->Z());
+          netgen::Point3d ngP2(p2.node->X(), p2.node->Y(), p2.node->Z());
+          ngMesh.RestrictLocalH( netgen::Center( ngP1,ngP2), Dist(ngP1,ngP2));
 #ifdef DUMP_SEGMENTS
           cout << "Segment: " << seg.edgenr << " on SMESH face " << helper.GetMeshDS()->ShapeToIndex( face ) << endl
                << "\tface index: " << seg.si << endl