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

Patch for issue #27842

Browse Source 
Task GEOM-04-009: Common operation fails on Set-56.
This commit is contained in:
nbv
2016-09-08 09:59:33 +03:00
parent 630d1a7f28
commit e93246f00a
1 changed files with 16 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

28
src/IntWalk/IntWalk_PWalking_1.gxx
View File

@@ -1813,6 +1813,7 @@ Standard_Boolean IntWalk_PWalking::
{
const Standard_Integer aNbIterMAX = 60;
const Standard_Real aTol = 1.0e-14;
const Standard_Real aTolNul = 1.0 / Precision::Infinite();
Handle(Geom_Surface) aS1, aS2;
switch(theASurf1->GetType())
@@ -1839,6 +1840,14 @@ Standard_Boolean IntWalk_PWalking::
return Standard_True;
}
// I.e. if theU1 = 0.0 then Epsilon(theU1) = DBL_MIN (~1.0e-308).
// Work with this number is impossible: there is a dangerous to
// obtain Floating-point-overflow. Therefore, we limit this value.
const Standard_Real aMinAddValU1 = Max(Epsilon(theU1), aTolNul);
const Standard_Real aMinAddValV1 = Max(Epsilon(theV1), aTolNul);
const Standard_Real aMinAddValU2 = Max(Epsilon(theU2), aTolNul);
const Standard_Real aMinAddValV2 = Max(Epsilon(theV2), aTolNul);
Standard_Boolean aStatus = Standard_False;
gp_Pnt aP1, aP2;
@@ -1864,21 +1873,16 @@ Standard_Boolean IntWalk_PWalking::
while(flRepeat)
{
Standard_Real anAdd = aGradFu*aSTEPuv;
Standard_Real aPARu = (anAdd >= 0.0)?
(theU1 - Max(anAdd, Epsilon(theU1))) :
(theU1 + Max(-anAdd, Epsilon(theU1)));
const Standard_Real aPARu = theU1 - Sign(Max(Abs(anAdd), aMinAddValU1), anAdd);
anAdd = aGradFv*aSTEPuv;
Standard_Real aPARv = (anAdd >= 0.0)?
(theV1 - Max(anAdd, Epsilon(theV1))) :
(theV1 + Max(-anAdd, Epsilon(theV1)));
const Standard_Real aPARv = theV1 - Sign(Max(Abs(anAdd), aMinAddValV1), anAdd);
anAdd = aGradFU*aStepUV;
Standard_Real aParU = (anAdd >= 0.0)?
(theU2 - Max(anAdd, Epsilon(theU2))) :
(theU2 + Max(-anAdd, Epsilon(theU2)));
const Standard_Real aParU = theU2 - Sign(Max(Abs(anAdd), aMinAddValU2), anAdd);
anAdd = aGradFV*aStepUV;
Standard_Real aParV = (anAdd >= 0.0)?
(theV2 - Max(anAdd, Epsilon(theV2))) :
(theV2 + Max(-anAdd, Epsilon(theV2)));
const Standard_Real aParV = theV2 - Sign(Max(Abs(anAdd), aMinAddValV2), anAdd);
gp_Pnt aPt1, aPt2;