puts "============" puts "BUC60901" puts "============" puts "" ################################################################ ## We should get the edges forming a closed intersection line but we see breaks. ################################################################ set n 3 set dx 1 set sx 20 set sz 10 set xc [dval (2*$dx+2*$sx)*$n/2] set yc 0 set zc [dval $sz*10] set nxc 0.5 set nyc -0.3 set nzc 1 set r [dval $xc*0.9] set x 0 set y [dval -$r*1.2] set z 0 set pol "polyline w $x $y $z" for {set i 1} {$i <= $n} {incr i} { set x [dval $x+$dx] set pol "$pol $x $y $z" set x [dval $x+$sx] set z [dval $z+$sz] set pol "$pol $x $y $z" set x [dval $x+$dx] set pol "$pol $x $y $z" set x [dval $x+$sx] set z [dval $z-$sz] set pol "$pol $x $y $z" } eval $pol circle c $xc $yc $zc $nxc $nyc $nzc $r mkedge c c set lcyl [dval -$zc*2] set lcrg [dval -$y*2] prism cyl c 0 0 $lcyl prism crg w 0 $lcrg 0 bsection result crg cyl set nb_info [nbshapes result] regexp {VERTEX +: +([-0-9.+eE]+)} $nb_info full nb_vertices regexp {EDGE +: +([-0-9.+eE]+)} $nb_info full nb_edges if { $nb_vertices != $nb_edges } { puts "Section is opened: Faulty shape" } else { puts "Section is closed: OK" } set length 350.404 set 2dviewer 0