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

0028692: Projection failed (projponf)

Browse Source 
Tolerance criterion is improved:
P.SquareDistance(p1) < Tol -> P.SquareDistance(p1) < Tol*Tol
This commit is contained in:
ifv 2017-05-16 17:37:13 +03:00 committed by bugmaster
parent b493028765
commit 8c745be5e0
2 changed files with 21 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/Extrema/Extrema_ExtPElS.cxx
View File

@ -313,9 +313,10 @@ void Extrema_ExtPElS::Perform(const gp_Pnt& P,
const gp_Torus& S,
const Standard_Real Tol)
{
const Standard_Real tol2 = Tol*Tol;
myDone = Standard_False;
myNbExt = 0;
// Projection of P in plane XOY ...
gp_Ax3 Pos = S.Position();
gp_Pnt O = Pos.Location();
@ -325,7 +326,7 @@ void Extrema_ExtPElS::Perform(const gp_Pnt& P,
// Calculation of extrema ...
gp_Vec OPp (O,Pp);
Standard_Real R2 = OPp.SquareMagnitude();
if (R2 < Tol * Tol) { return; }
if (R2 < tol2) { return; }
gp_Vec myZ = Pos.XDirection()^Pos.YDirection();
Standard_Real U1 = gp_Vec(Pos.XDirection()).AngleWithRef(OPp,myZ);
@ -338,8 +339,8 @@ void Extrema_ExtPElS::Perform(const gp_Pnt& P,
gp_Pnt O1 = O.Translated(OO1);
gp_Pnt O2 = O.Translated(OO2);
if(O1.SquareDistance(P) < Tol) { return; }
if(O2.SquareDistance(P) < Tol) { return; }
if(O1.SquareDistance(P) < tol2) { return; }
if(O2.SquareDistance(P) < tol2) { return; }
Standard_Real V1 = OPp.AngleWithRef(gp_Vec(O1,P),OPp.Crossed(OZ));
if (V1 > -ExtPElS_MyEps && V1 < ExtPElS_MyEps) { V1 = 0.; }

16
tests/bugs/modalg_6/bug28692 Normal file
View File

@ -0,0 +1,16 @@
puts "========"
puts "OCC28692"
puts "========"
puts ""
####################################################################################
## Projection failed (projponf)
####################################################################################
restore [locate_data_file bug28692_shapes.brep] bi3
point p -0.0153 0 0.017
set dist -1
regexp {proj dist = ([-0-9.+eE]+)} [projponf bi3 p] full dist
if { $dist < 0 || $dist > 1.e-7 } {
puts "Error : Projection failed "
}