ForkMaster/smesh
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/smesh.git synced 2024-11-15 10:08:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

avoid FPE exception

Browse Source
This commit is contained in:
ptv 2009-09-15 07:51:41 +00:00
parent 7a2dff55c5
commit 4c604aa712
1 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/StdMeshers/StdMeshers_MEFISTO_2D.cxx
View File

@ -53,6 +53,7 @@
#include <Geom_Curve.hxx> #include <Geom_Curve.hxx>
#include <Geom2d_Curve.hxx> #include <Geom2d_Curve.hxx>
#include <Geom_Surface.hxx> #include <Geom_Surface.hxx>
#include <Precision.hxx>
#include <TopExp.hxx> #include <TopExp.hxx>
#include <TopExp_Explorer.hxx> #include <TopExp_Explorer.hxx>
#include <TopTools_ListIteratorOfListOfShape.hxx> #include <TopTools_ListIteratorOfListOfShape.hxx>
@ -303,7 +304,7 @@ bool StdMeshers_MEFISTO_2D::Evaluate(SMESH_Mesh & aMesh,
TopoDS_Face F = TopoDS::Face(aShape.Oriented(TopAbs_FORWARD)); TopoDS_Face F = TopoDS::Face(aShape.Oriented(TopAbs_FORWARD));
double aLen = 0.0; double aLen = 0.0;
double NbSeg = 0; int NbSeg = 0;
bool IsQuadratic = false; bool IsQuadratic = false;
bool IsFirst = true; bool IsFirst = true;
TopExp_Explorer exp(F,TopAbs_EDGE); TopExp_Explorer exp(F,TopAbs_EDGE);
@ -331,16 +332,20 @@ bool StdMeshers_MEFISTO_2D::Evaluate(SMESH_Mesh & aMesh,
P1 = P2; P1 = P2;
} }
} }
if ( NbSeg > 0 )
aLen = aLen/NbSeg; // middle length aLen = aLen/NbSeg; // middle length
else
aLen = Precision::Infinite();
_edgeLength = DBL_MAX; _edgeLength = Precision::Infinite();
double tmpLength = Min( _edgeLength, aLen ); double tmpLength = Min( _edgeLength, aLen );
GProp_GProps G; GProp_GProps G;
BRepGProp::SurfaceProperties(aShape,G); BRepGProp::SurfaceProperties(aShape,G);
double anArea = G.Mass(); double anArea = G.Mass();
int nbFaces = (int) ( anArea/(tmpLength*tmpLength*sqrt(3.)/4) ); int nbFaces = Precision::IsInfinite( tmpLength ) ? 0 :
(int)( anArea/(tmpLength*tmpLength*sqrt(3.)/4) );
int nbNodes = (int) ( nbFaces*3 - (NbSeg-1)*2 ) / 6; int nbNodes = (int) ( nbFaces*3 - (NbSeg-1)*2 ) / 6;
std::vector<int> aVec(SMDSEntity_Last); std::vector<int> aVec(SMDSEntity_Last);