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

0021131: EDF 1167 NETGENPLUGIN: SIGSEGV when trying to mesh a shape with degenerated edges

Browse Source 
prevent access to NULL curve of the degenerated edge
This commit is contained in:
eap 2011-01-21 17:45:09 +00:00
parent 7aacc6b391
commit c4fe195950
1 changed files with 19 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
src/NETGENPlugin/NETGENPlugin_Mesher.cxx
View File

@ -1416,17 +1416,28 @@ namespace
const int nb = 1000; const int nb = 1000;
Standard_Real u1, u2; Standard_Real u1, u2;
Handle(Geom_Curve) curve = BRep_Tool::Curve(edge, u1, u2); Handle(Geom_Curve) curve = BRep_Tool::Curve(edge, u1, u2);
Standard_Real delta = (u2-u1)/nb; if ( curve.IsNull() )
for(int i=0; i<nb; i++)
{ {
Standard_Real u = u1 + delta*i; TopoDS_Iterator vIt( edge );
gp_Pnt p = curve->Value(u); if ( !vIt.More() ) return;
gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( vIt.Value() ));
netgen::Point3d pi(p.X(), p.Y(), p.Z()); netgen::Point3d pi(p.X(), p.Y(), p.Z());
mesh.RestrictLocalH(pi, size); mesh.RestrictLocalH(pi, size);
double resultSize = mesh.GetH(pi); }
if ( resultSize - size > 0.1*size ) else
// netgen does restriction iff oldH/newH > 1.2 (localh.cpp:136) {
mesh.RestrictLocalH(pi, resultSize/1.201); Standard_Real delta = (u2-u1)/nb;
for(int i=0; i<nb; i++)
{
Standard_Real u = u1 + delta*i;
gp_Pnt p = curve->Value(u);
netgen::Point3d pi(p.X(), p.Y(), p.Z());
mesh.RestrictLocalH(pi, size);
double resultSize = mesh.GetH(pi);
if ( resultSize - size > 0.1*size )
// netgen does restriction iff oldH/newH > 1.2 (localh.cpp:136)
mesh.RestrictLocalH(pi, resultSize/1.201);
}
} }
} }