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

0025223: Wrong result done by projection algorithm

Browse Source 
class ProjLib_ComputeApprox
method
 ProjLib_ComputeApprox::ProjLib_ComputeApprox
   (const Handle(Adaptor3d_HCurve)   & C,
    const Handle(Adaptor3d_HSurface) & S,
    const Standard_Real              Tol )

Using Precision::PConfusion() as a 2d tolerance when shifting the projected curve.

Test case for issue CR25223
This commit is contained in:
emv 2014-09-18 14:26:14 +04:00 committed by bugmaster
parent 53c888d174
commit 6aa39ac3ad
2 changed files with 96 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/ProjLib/ProjLib_ComputeApprox.cxx
View File

@ -1072,9 +1072,14 @@ ProjLib_ComputeApprox::ProjLib_ComputeApprox
}
if (F.UCouture || (F.VCouture && SType == GeomAbs_Sphere))
{
gp_Pnt2d P2d = F.Value( UFirst );
number = (Standard_Integer) (Floor((P2d.X()-u)/M_PI + Epsilon(M_PI)));
du = -number*M_PI;
Standard_Real aNbPer;
gp_Pnt2d P2d = F.Value(UFirst);
du = u - P2d.X();
du = (du < 0) ? (du - Precision::PConfusion()) :
(du + Precision::PConfusion());
modf(du/M_PI, &aNbPer);
number = (Standard_Integer)aNbPer;
du = number*M_PI;
}
if (!myBSpline.IsNull())

88
tests/bugs/moddata_3/bug25223 Normal file
View File

@ -0,0 +1,88 @@
puts "================"
puts "OCC25223"
puts "================"
puts ""
#######################################################################
# Wrong result done by projection algorithm
#######################################################################
restore [locate_data_file bug25223_surface] s
restore [locate_data_file bug25223_curve] c
project c2d c s
# I.
#Parameters : 2.50050886855066 2.71802559491834
regexp {Parameters : ([-0-9.+eE]+) ([-0-9.+eE]+)} [dump c] full a1 a2
set good_a1 2.50050886855066
set good_a2 2.71802559491834
set c_tol 1.e-13
if { [expr abs(${a1} - ${good_a1}) ] > ${c_tol} } {
puts "Error: invalid a1 parameter"
}
if { [expr abs(${a2} - ${good_a2}) ] > ${c_tol} } {
puts "Error: invalid a2 parameter"
}
# II.
set u1 2.6
cvalue c $u1 x y z
vertex v1 x y z
point p1 x y z
mkface f s
set proj_fp [projponf f p1]
#proj dist = 9.8455689542838467e-014
regexp {proj dist = ([-0-9.+eE]+) uvproj = \(([-0-9.+eE]+) ([-0-9.+eE]+)\); pproj = \(([-0-9.+eE]+) ([-0-9.+eE]+) ([-0-9.+eE]+)\)} ${proj_fp} full dist uproj vproj proj1 proj2
puts "dist=${dist}"
puts "uproj=${uproj}"
puts "vproj=${vproj}"
puts "proj1=${proj1}"
puts "proj2=${proj2}"
set tolmax_f [tolmax f]
regexp {max tol = ([-0-9.+eE]+)} ${tolmax_f} full CMP_TOL
puts "CMP_TOL=${CMP_TOL}"
set good_dist 9.8455689542838467e-14
set good_uproj 0.66864242394263962
set good_vproj 0.53857983853386682
set good_proj1 336.79267755053786
set good_proj2 -256.45869438720723
if { [expr abs(${dist} - ${good_dist}) ] > ${CMP_TOL} } {
puts "Error: invalid dist"
}
if { [expr abs(${uproj} - ${good_uproj}) ] > ${CMP_TOL} } {
puts "Error: invalid uproj"
}
if { [expr abs(${vproj} - ${good_vproj}) ] > ${CMP_TOL} } {
puts "Error: invalid vproj"
}
if { [expr abs(${proj1} - ${good_proj1}) ] > ${CMP_TOL} } {
puts "Error: invalid proj1"
}
if { [expr abs(${proj2} - ${good_proj2}) ] > ${CMP_TOL} } {
puts "Error: invalid proj2"
}
# III.
2dcvalue c2d $u1 u v
svalue s u v x y z
vertex v2 x y z
distmini d v1 v2
#the distance value is : 3.99072821948773e-08
regexp {([-0-9.+eE]+)$} [dump d_val] full dist
regexp { +Vertex +: +Min +[-0-9.+eE]+ +Max +([-0-9.+eE]+)} [ maxtolerance d ] full toler
set good_dist 3.99072821948773e-08
if { [expr abs( ${dist} - ${good_dist} )] > ${toler} } {
puts "Error : the distanse is ${dist}. It is bad value"
}
smallview
fit
set only_screen_axo 1