occt/tests/bugs/moddata_2/bug22758

puts "========================"
puts "OCC22758"
puts "========================"
puts ""
#######################################################################
# Problem in BRepExtrema_DistShapeShape
#
# This script tries to reproduce the bug in BRepExtrema_DistanceSS class reported
# by Rob Bacnrach on OCCT Forum: http://www.opencascade.org/org/forum/thread_23040
# (when first shape is edge and second is vertex, the points are returned
# for the solutions are swapped)
#
#######################################################################

set BugNumber OCC22758

bsplinecurve c1 2   2 0 3 1 3       0 2 0 1  1 1 0 1  2 2 0 1
mkedge e c1
vertex v 0 0 1

set res [distmini r e v]

set status 0
# result should contain one edge, starting on e and ending on v
if { [llength $res] != 2 } {
    puts "Error: result has structure different from expected: $res"
    set status 1
}

# get start and end vertices -- these should be solutions on e1 and e2
set sol [explode [lindex $res 1] v]

# check distances in correct order
set d1 [dval [lindex [distmini d1 [lindex $sol 0] e] 0]]
set d2 [dval [lindex [distmini d2 [lindex $sol 1] v] 0]]
puts "Distances from solutions to relevant objecte: $d1, $d2"
if { $d1 > 1e-7 || $d2 > 1e-7 } {
    puts "Error: distances are non-zero!"
    set status 1
}

# check distances in swapped order
set d1 [dval [lindex [distmini d1 [lindex $sol 0] v] 0]]
set d2 [dval [lindex [distmini d2 [lindex $sol 1] e] 0]]
puts "Distances from solutions to swapped objecte: $d1, $d2"
if { $d1 < 0.1 || $d2 < 0.1 } {
    puts "Error: swapped distances are zero, i.e. solutions are swapped!"
    set status 1
}

# Resume
puts ""
if { ${status} == 1 } {
   puts "Faulty ${BugNumber}"
} else {
   puts "OK ${BugNumber}"
}