puts "========="
puts "OCC23948"
puts "========="
puts ""
###########################################################
# Wrong intersection between a surface of revolution and a plane.
###########################################################

restore [locate_data_file bug22766_f1] f1
restore [locate_data_file bug22766_f2] f2

bopcurves f1 f2

mksurface s1 f1
mksurface s2 f2

dlog reset
dlog on
xdistcs c_1 s2 0 1 10 
set Log1 [dlog get]

set List1 [split ${Log1} {TD= \t\n}]

set L1 [llength ${List1}]
set L2 10
set L3 5
set N [expr (${L1} - ${L2})/${L3} + 1]
set Tolerance 1.0e-3
set D_good 0.

for {set i 1} {${i} <= ${N}} {incr i} {
   set j1 [expr ${L2} + (${i}-1)*${L3}]
   set j2 [expr ${j1} + 2]
   set T [lindex ${List1} ${j1}]
   set D [lindex ${List1} ${j2}]
   #puts "i=${i} j1=${j1} j2=${j2} T=${T} D=${D}"
   if { [expr abs(${D} - ${D_good})] > ${Tolerance} } {
      puts "Error: i=${i} T=${T} D=${D}"
   }
}

smallview
fit
set only_screen_axo 1 
