]> SALOME platform Git repositories - plugins/netgenplugin.git/commitdiff
Salome HOME
Fix of 20124 issue (EDF SMESH 920 : "Max element volume" hypothesis is not used)...
authorakl <akl@opencascade.com>
Mon, 19 Jan 2009 16:09:31 +0000 (16:09 +0000)
committerakl <akl@opencascade.com>
Mon, 19 Jan 2009 16:09:31 +0000 (16:09 +0000)
src/NETGENPlugin/NETGENPlugin_NETGEN_3D.cxx

index 0a22ea3d64e2dea767b7e1a960fc9f4c5bcb04d2..00d447d6f663a87cff4b5bd5472f57e63a6e6221 100644 (file)
@@ -276,7 +276,7 @@ bool NETGENPlugin_NETGEN_3D::Compute(SMESH_Mesh&         aMesh,
   int Netgen_NbOfNodes = 0;
   int Netgen_param2ndOrder = 0;
   double Netgen_paramFine = 1.;
-  double Netgen_paramSize = _maxElementVolume;
+  double Netgen_paramSize = pow( 72, 1/6. ) * pow( _maxElementVolume, 1/3. );
 
   double Netgen_point[3];
   int Netgen_triangle[3];
@@ -500,7 +500,7 @@ bool NETGENPlugin_NETGEN_3D::Compute(SMESH_Mesh& aMesh,
   int Netgen_NbOfNodes = 0;
   int Netgen_param2ndOrder = 0;
   double Netgen_paramFine = 1.;
-  double Netgen_paramSize = _maxElementVolume;
+  double Netgen_paramSize = pow( 72, 1/6. ) * pow( _maxElementVolume, 1/3. );
   
   double Netgen_point[3];
   int Netgen_triangle[3];
@@ -601,8 +601,8 @@ bool NETGENPlugin_NETGEN_3D::Compute(SMESH_Mesh& aMesh,
     {
       Ng_GetPoint( Netgen_mesh, nodeIndex, Netgen_point );
       SMDS_MeshNode * node = aHelper->AddNode(Netgen_point[0],
-                                             Netgen_point[1],
-                                             Netgen_point[2]);
+                                             Netgen_point[1],
+                                             Netgen_point[2]);
       nodeVec.at(nodeIndex) = node;
     }
 
@@ -611,9 +611,9 @@ bool NETGENPlugin_NETGEN_3D::Compute(SMESH_Mesh& aMesh,
     {
       Ng_GetVolumeElement(Netgen_mesh, elemIndex, Netgen_tetrahedron);
       aHelper->AddVolume (nodeVec.at( Netgen_tetrahedron[0] ),
-                                                 nodeVec.at( Netgen_tetrahedron[1] ),
-                                                 nodeVec.at( Netgen_tetrahedron[2] ),
-                                                 nodeVec.at( Netgen_tetrahedron[3] ));
+                         nodeVec.at( Netgen_tetrahedron[1] ),
+                         nodeVec.at( Netgen_tetrahedron[2] ),
+                         nodeVec.at( Netgen_tetrahedron[3] ));
     }
   }