mirror of
https://git.dev.opencascade.org/repos/occt.git
synced 2025-04-03 17:56:21 +03:00
3306fdd954
The algorithm has been improved for the cases when the intersection line goes through the cone apex. <!break> 1. All special points are put to the ALine forcefully (if they are true intersection point). Currently this step has not been implemented yet. 2. Now the tolerance of IntPatch_Point (put into ALine) is computed in order to cover the distance between it and the correspond ALine. 3. Test cases have been created. 4. Procedure of trimming IntAna_Curve has been improved. 5. Criterion when the discriminant of IntAna_Curve can be considered to be equal to 0 has been improved. 6. Methods IntAna_Curve::FindParameter(...) (and IntPatch_ALine::FindParameter(...)) currently returns list of all parameters corresponding the given point (IntAna_Curve can be self-interfered curve). Before the fix, this method always returned only one (randomly chosen) parameter. 7. Interfaces of the following methods have been changed: IntAna_Curve::FindParameter(...), IntPatch_ALine::FindParameter(...), IntPatch_ALine::ChangeVertex(...), IntPatch_SpecialPoints::AddPointOnUorVIso(...), IntPatch_SpecialPoints::AddSingularPole(...), IntPatch_WLineTool::ExtendTwoWLines(). 8. Following methods have been added: IntAna_Quadric::SpecialPoints(...), IntPatch_ALineToWLine::GetSectionRadius(...), IntPatch_SpecialPoints::ProcessSphere(...), IntPatch_SpecialPoints::ProcessCone(...), IntPatch_SpecialPoints::GetTangentToIntLineForCone(...). ------------------ 1) tests/boolean/volumemaker/C5 tests/boolean/volumemaker/C6 tests/boolean/volumemaker/E7 They are real IMPROVEMENTS. In the FIX (in compare with MASTER), section result between pairs of faces f2&f6 (C5), f3&f7 (C6) and f1&f5 (E7) is closed. Separated test cases have been created in order to focus on the problem with section. Bug #28503 has been fixed. Correction in test cases.
55 lines
1.1 KiB
Plaintext
55 lines
1.1 KiB
Plaintext
puts "========"
|
|
puts "0029824: Intersection of cylinder and sphere is incorrect"
|
|
puts "========"
|
|
puts ""
|
|
|
|
restore [locate_data_file OCC20964_revsolid.brep] b1
|
|
restore [locate_data_file OCC20964_sphere.brep] b2
|
|
explode b1 f
|
|
explode b2 f
|
|
donly b1_7 b2_1
|
|
|
|
bop b1_7 b2_1
|
|
bopsection result
|
|
|
|
checknbshapes result -edge 3
|
|
|
|
foreach a [explode result e] {
|
|
mkcurve cc $a
|
|
bounds cc U1 U2
|
|
|
|
if {[dval U2-U1] < 1.0e-9} {
|
|
puts "Error: Wrong curve's range!"
|
|
}
|
|
|
|
set aStep [expr [dval U2-U1]/100.0 ]
|
|
set isFirst 1
|
|
set aCosPrev 1.0
|
|
dset dx1 0.0
|
|
dset dy1 0.0
|
|
dset dz1 0.0
|
|
for {set aU [dval U1]} {$aU <= [dval U2]} {set aU [expr $aU + $aStep]} {
|
|
cvalue cc $aU xx yy zz dx dy dz
|
|
|
|
if {!$isFirst} {
|
|
set m1 [module dx1 dy1 dz1]
|
|
set m2 [module dx dy dz]
|
|
set aCos [dval dx1*dx+dy1*dy+dz1*dz]
|
|
set aCos [expr $aCos/($m1*$m2)]
|
|
|
|
if {abs([expr $aCos - $aCosPrev]) > 0.05} {
|
|
puts "Error: It seems that the curve $a change it direction at the point $aU. Please recheck."
|
|
break
|
|
}
|
|
|
|
set aCosPrev $aCos
|
|
}
|
|
|
|
set isFirst 0
|
|
|
|
dset dx1 dx
|
|
dset dy1 dy
|
|
dset dz1 dz
|
|
}
|
|
}
|