OpenCascade/occt
1
0
Fork 0
mirror of https://git.dev.opencascade.org/repos/occt.git synced 2025-04-03 17:56:21 +03:00
Code Packages Releases Activity

0028076: Test bugs moddata_2 bug469 fails when FPE are enabled

Browse Source 
The reason of FPE exception has been eliminated.
This commit is contained in:
nbv 2016-11-09 11:22:12 +03:00 committed by apn
parent bf7b2cebbe
commit 326b3e283a
1 changed files with 20 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
src/ApproxInt/ApproxInt_KnotTools.cxx
View File

@ -57,17 +57,33 @@ static Standard_Real EvalCurv(const Standard_Real dim,
mp += p*p;
}
}
//mp *= 2.; //P(j,i) = -P(i,j);
mp = Sqrt(mp);
//
Standard_Real q = 0.;
for(i = 0; i < dim; ++i)
{
q += V1[i]*V1[i];
}
q = Sqrt(q);
if (q < 1 / Precision::Infinite())
{
// Indeed, if q is small then we can
// obtain equivocation of "0/0" type.
// In this case, local curvature can be
// not equal to 0 or Infinity.
// However, it is good solution to insert
// knot in the place with such singularity.
// Therefore, we need imitation of curvature
// jumping. Return of Precision::Infinite() is
// enough for it.
return Precision::Infinite();
}
q = Min(q, Precision::Infinite());
q *= q*q;
//
Standard_Real curv = mp / (q*q*q);
Standard_Real curv = Sqrt(mp / q);
return curv;
}