puts "============"
puts "OCCGH648"
puts "============"
puts ""
#######################################################################
#    spiral helix: comphelix2, D1 = 100, D4 = 20, 1 part, pitche = 1; 10; 1, number of turns = 45; 1; 45
#######################################################################

set BugNumber OCCGH648

set np 3
set D1 100
set D4 20

set P1 1
set P2 10
set P3 1

set N1 45
set N2 1
set N3 45

set H1 [expr $P1*$N1]
set H2 [expr $P2*$N2]
set H3 [expr $P3*$N3]

set D2 [expr ($D1 - ($D1-$D4)*($H1)/($H1+$H2+$H3))]
set D3 [expr ($D1 - ($D1-$D4)*($H1+$H2)/($H1+$H2+$H3))]

set mistake 0
set result ""
if [catch { set info_result [comphelix2 result ${np} ${D1} ${D2} ${D3} ${D4} ${P1} ${P2} ${P3} ${N1} ${N2} ${N3}]} ] {
  puts "${info_result}"
  puts "Faulty ${BugNumber} : approximation algorithm is broken"
  set mistake 1
} elseif { [regexp {ErrorStatus} ${info_result}] } {
  puts "${info_result}"
  puts "Faulty ${BugNumber} : helix is broken"
  set mistake 1
}

if { ${mistake} == 0 } {
  set d 0.1
  set x [expr $D1/2]
  set y 0
  set z 0
  set dx 0
  set dy 1
  set dz 0

  circle c $x $y $z $dx $dy $dz $d
  mkedge e c
  wire s2 e

  renamevar result s1


set square 10778.1
set volume 538.896

set nb_v_good 92
set nb_e_good 183
set nb_w_good 93
set nb_f_good 93
set nb_sh_good 1
set nb_sol_good 1
set nb_compsol_good 0
set nb_compound_good 0
set nb_shape_good 463

CheckSweep
}
set 2dviewer 0
