ForkMaster/netgenplugin
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

0019296: Too small Max element Area

Browse Source
This commit is contained in:
ptv 2009-10-01 06:28:39 +00:00
parent fb14d85a93
commit 347223da4b
1 changed files with 14 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
src/NETGENPlugin/NETGENPlugin_NETGEN_2D_ONLY.cxx
View File

@ -39,6 +39,7 @@
#include "StdMeshers_LengthFromEdges.hxx" #include "StdMeshers_LengthFromEdges.hxx"
#include "StdMeshers_QuadranglePreference.hxx" #include "StdMeshers_QuadranglePreference.hxx"
#include <Precision.hxx>
#include <Standard_ErrorHandler.hxx> #include <Standard_ErrorHandler.hxx>
#include <Standard_Failure.hxx> #include <Standard_Failure.hxx>
@ -492,17 +493,28 @@ bool NETGENPlugin_NETGEN_2D_ONLY::Evaluate(SMESH_Mesh& aMesh,
// compute edge length // compute edge length
double ELen = 0; double ELen = 0;
if (_hypLengthFromEdges || !_hypLengthFromEdges && !_hypMaxElementArea) { if (_hypLengthFromEdges || !_hypLengthFromEdges && !_hypMaxElementArea) {
ELen = fullLen / nb1d; if ( nb1d > 0 )
ELen = fullLen / nb1d;
} }
if ( _hypMaxElementArea ) { if ( _hypMaxElementArea ) {
double maxArea = _hypMaxElementArea->GetMaxArea(); double maxArea = _hypMaxElementArea->GetMaxArea();
ELen = sqrt(2. * maxArea/sqrt(3.0)); ELen = sqrt(2. * maxArea/sqrt(3.0));
} }
if ( ELen < Precision::Confusion() ) {
SMESH_subMesh *sm = aMesh.GetSubMesh(F);
if ( sm ) {
SMESH_ComputeErrorPtr& smError = sm->GetComputeError();
smError.reset( new SMESH_ComputeError(COMPERR_ALGO_FAILED,"Submesh can not be evaluated.\nToo small element length",this));
}
return false;
}
GProp_GProps G; GProp_GProps G;
BRepGProp::SurfaceProperties(F,G); BRepGProp::SurfaceProperties(F,G);
double anArea = G.Mass(); double anArea = G.Mass();
int nbFaces = (int) ( anArea / ( ELen*ELen*sqrt(3.) / 4 ) ); int nbFaces = 0;
if ( ELen > Precision::Confusion() )
nbFaces = (int) ( anArea / ( ELen*ELen*sqrt(3.) / 4 ) );
int nbNodes = (int) ( ( nbFaces*3 - (nb1d-1)*2 ) / 6 + 1 ); int nbNodes = (int) ( ( nbFaces*3 - (nb1d-1)*2 ) / 6 + 1 );
std::vector<int> aVec(SMDSEntity_Last); std::vector<int> aVec(SMDSEntity_Last);
for(int i=SMDSEntity_Node; i<SMDSEntity_Last; i++) aVec[i]=0; for(int i=SMDSEntity_Node; i<SMDSEntity_Last; i++) aVec[i]=0;