OpenCascade/occt
1
0
Fork 0
mirror of https://git.dev.opencascade.org/repos/occt.git synced 2025-08-19 13:40:49 +03:00
Code Activity

0030062: Modeling Algorithms - The shape healing is instable in test "parasolid/doc_3/E3"

Browse Source 
A solution to fix the instability was created.
This commit is contained in:
abk
2018-08-13 17:43:39 +03:00
parent 47ba172e98
commit c7c6b408cc
1 changed files with 23 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
src/ShapeAnalysis/ShapeAnalysis_Wire.cxx
View File

@@ -937,13 +937,29 @@ Standard_Boolean ShapeAnalysis_Wire::CheckDegenerated (const Standard_Integer nu
return Standard_False; return Standard_False;
} }
//added by rln 18/12/97 CSR# CTS18544 entity 13638 {
//the situation when degenerated edge already exists but flag is not set //added by rln 18/12/97 CSR# CTS18544 entity 13638
//(i.e. the parametric space is closed) //the situation when degenerated edge already exists but flag is not set
GeomAdaptor_Surface& Ads = mySurf->Adaptor3d()->ChangeSurface(); //(i.e. the parametric space is closed)
Standard_Real max = Max ( Ads.UResolution(myPrecision), GeomAdaptor_Surface& Ads = mySurf->Adaptor3d()->ChangeSurface();
Ads.VResolution(myPrecision) ); Standard_Real max = Max ( Ads.UResolution(myPrecision),
if ( p2d1.Distance (p2d2) /*Abs (par1 - par2)*/ <= max + gp::Resolution() ) return Standard_False; Ads.VResolution(myPrecision) );
// The constant is used to avoid numerical instability when the
// computational error for the coordinates of "p2d1" and "p2d2" is more than
// the distance between these points and even more than "max".
// The value of the constant should be not less than the relative
// computational error for the coordinates.
// The value of the constant is selected based on the case described in OCCT
// Mantiss issue 30062.
const Standard_Real aRelativePrecisionP2d1P2d2 = 2e-6;
const Standard_Real aCoordinateMax =
Max(Max(Abs(p2d1.X()), Abs(p2d1.Y())), Max(Abs(p2d2.X()), Abs(p2d2.Y())));
if (p2d1.Distance (p2d2) /*Abs (par1 - par2)*/ <= (max + gp::Resolution()) +
aCoordinateMax * aRelativePrecisionP2d1P2d2) return Standard_False;
}
//#84 rln p2d1 = aP2d.XY() + par1 * theDir2d.XY(); //#84 rln p2d1 = aP2d.XY() + par1 * theDir2d.XY();
//#84 rln p2d2 = aP2d.XY() + par2 * theDir2d.XY(); //#84 rln p2d2 = aP2d.XY() + par2 * theDir2d.XY();