]> SALOME platform Git repositories - plugins/netgenplugin.git/commitdiff
Salome HOME
Fix regression of test grid 2d_mesh_NETGEN/D2
authoreap <eap@opencascade.com>
Wed, 2 Nov 2011 08:56:59 +0000 (08:56 +0000)
committereap <eap@opencascade.com>
Wed, 2 Nov 2011 08:56:59 +0000 (08:56 +0000)
  fix restricting local H, which was broken by a defaut netgen::mparam.minh

src/NETGENPlugin/NETGENPlugin_NETGEN_2D_ONLY.cxx

index 06a5bc42383839eaac857aef255c058ece2bdcf9..40acc112e120c67f7b734947cc85e70c607cc507 100644 (file)
@@ -301,15 +301,14 @@ static TError addSegmentsToMesh(netgen::Mesh&                    ngMesh,
       ngMesh.AddSegment (seg);
       {
         // restrict size of elements near the segment
-        netgen::Point3d ngP1(n->X(), n->Y(), n->Z());
-        n = uvPtVec[ i+1 ].node;
-        netgen::Point3d ngP2(n->X(), n->Y(), n->Z());
+        SMESH_TNodeXYZ np1( n ), np2( uvPtVec[ i+1 ].node );
         // get an average size of adjacent segments to avoid sharp change of
         // element size (regression on issue 0020452, note 0010898)
         int iPrev = SMESH_MesherHelper::WrapIndex( i-1, nbSegments );
         int iNext = SMESH_MesherHelper::WrapIndex( i+1, nbSegments );
         double avgH = ( segLen[ iPrev ] + segLen[ i ] + segLen[ iNext ]) / 3;
-        ngMesh.RestrictLocalH( netgen::Center( ngP1,ngP2), avgH );
+
+        NETGENPlugin_Mesher::RestrictLocalSize( ngMesh, 0.5*(np1+np2), avgH );
       }
 #ifdef DUMP_SEGMENTS
         cout << "Segment: " << seg.edgenr << endl