Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header
  !goto proc
!endif

slib_author=
slib_parms=6\
,[zone de dessin x min,max et y min,max]\
,coordonnees des 2 points D=(xb,yb) et F=(xe,ye)\
,(optionnel) sens de la parabole up pour dessus du segment, down en dessous

slib_title=

slib_out=3ème point x,y de la courbe par la méthode de Bezier, \
  puis les points (éventuellement) modifiés de départ et d'arrivée de la courbe.\
  doit pouvoir être amélioré !

slib_comment=[zone de dessin xmin et xmax,zone de dessin ymin et ymax],\
  coordonnées xb et yb du point D,coordonnées xe et ye du point F,\
  (optionnel) sens de la parabole up dessus et down dessous

slib_example=[-10,10,-10,10],8,8,-7,0,down\
[-10,10,-10,10],8,8,-7,0,up\
[-10,10,-10,10],8,0,-7,8,down\
[-10,10,-10,10],8,0,-7,8,up\
[-10,10,-10,10],8,0,-7,0,down\
[-10,10,-10,10],8,0,-7,0,up\
[-10,10,-10,10],-7,0,8,8,down\
[-10,10,-10,10],-7,0,8,8,up\
[-10,10,-10,10],-7,8,8,0,down\
[-10,10,-10,10],-7,8,8,0,up\
[-10,10,-10,10],-7,0,8,0,down\
[-10,10,-10,10],-7,0,8,0,up\
[-10,10,-10,10],8,0,8,8,down\
[-10,10,-10,10],8,0,8,8,up\
[-10,10,-10,10],8,8,8,0,down\
[-10,10,-10,10],8,8,8,0,up\

!exit

:proc

!!--- init
!reset slib_out slib_med
!distribute items $wims_read_parm into slib_zone,slib_xb,slib_yb,slib_xe,slib_ye,slib_pos
!set slib_zone=!declosing $slib_zone
!set slib_xmin=$($slib_zone[1])
!set slib_xmax=$($slib_zone[2])
!set slib_ymin=$($slib_zone[3])
!set slib_ymax=$($slib_zone[4])
slib_xy=$slib_xb,$slib_yb,$slib_xe,$slib_ye
!!slib_xyinitial=$slib_xy
!!distribute items $slib_xy into slib_xb,slib_yb,slib_xe,slib_ye
!!--- calculs pour trouver down ou up
!if $slib_pos=
  !set slib_L=$[sqrt((($slib_xe-$slib_xb)^2)+(($slib_ye-$slib_yb)^2))]
  !set slib_xm=$[($slib_xe+$slib_xb)/2]
  !set slib_ym=$[($slib_ye+$slib_yb)/2]
  !set slib_m=$[max($[$slib_xmax-$slib_xmin],$[$slib_ymax-$slib_ymin])]
  !set slib_b=$[min($[$slib_xmax-$slib_xm],$[$slib_ymax-$slib_ym])]
  !set slib_d=$[min($[$slib_xm-$slib_xmin],$[$slib_ym-$slib_ymin])]
  !set slib_s=$[max($[($slib_m-$slib_L)/4],$[$slib_m/10])]
  !if $slib_s < $slib_b
    !set slib_pos=up
  !else
    !set slib_pos=down
  !endif
!endif
!!on change un peu les extrémités
!distribute item 0,0,0,0, into slib_signxb,slib_signyb,slib_signxe,slib_signye
!set slib_ss=1
!if $slib_pos=up
  !set slib_ss=1
!endif
!if $slib_pos=down
  !set slib_ss=-1
!endif
!if $slib_xb < $slib_xe and $slib_yb < $slib_ye
    !distribute item 0,$slib_ss,0,($slib_ss) into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb < $slib_xe and $slib_yb > $slib_ye
  !distribute item 0,($slib_ss),0,($slib_ss) into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb < $slib_xe and $slib_yb = $slib_ye
  !distribute item 0,$slib_ss,0,$slib_ss into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb > $slib_xe and $slib_yb < $slib_ye
  !distribute item 0,$slib_ss,0,$slib_ss into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb > $slib_xe and $slib_yb > $slib_ye
  !distribute item 0,$slib_ss,0,$slib_ss into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb > $slib_xe and $slib_yb = $slib_ye
  !distribute item 0,$slib_ss,0,$slib_ss into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb = $slib_xe and $slib_yb > $slib_ye
  !distribute item $slib_ss,0,$slib_ss,0 into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!if $slib_xb = $slib_xe and $slib_yb < $slib_ye
  !distribute item $slib_ss,0,$slib_ss,0 into slib_signxb,slib_signyb,slib_signxe,slib_signye
!endif
!set slib_med=!getopt med in $wims_read_parm
!if $(slib_med[1])=0
  !distribute items 0,0 into slib_signxb,slib_signyb
!endif
!if $(slib_med[2])=0
  !distribute items 0,0 into slib_signxe,slib_signye
!endif
!set slib_tmp=0.15
slib_tmp=($slib_signxb)*$slib_tmp,($slib_signyb)*$slib_tmp,($slib_signxe)*$slib_tmp,($slib_signye)*$slib_tmp
slib_xy=!exec pari [$slib_xy]+3*[$slib_tmp]

!distribute items $slib_xy into slib_xb,slib_yb,slib_xe,slib_ye
!!--- calculs preliminaires
!set slib_L=$[sqrt((($slib_xe-$slib_xb)^2)+(($slib_ye-$slib_yb)^2))]

!set slib_xm=$[($slib_xe+$slib_xb)/2]
!set slib_ym=$[($slib_ye+$slib_yb)/2]

!set slib_m=$[max($[$slib_xmax-$slib_xmin],$[$slib_ymax-$slib_ymin])]
!set slib_b=$[min($[$slib_xmax-$slib_xm],$[$slib_ymax-$slib_ym])]
!set slib_d=$[min($[$slib_xm-$slib_xmin],$[$slib_ym-$slib_ymin])]

!set slib_s=$[max($[($slib_m-$slib_L)/4],$[$slib_m/10])]

!if $slib_pos issametext up
  !set slib_h=$[min($slib_s,$slib_b)]
!else
  !if $slib_pos issametext down
    !set slib_h=$[-min($slib_s,$slib_d)]
  !endif
!endif

!!--- determination XC,YC

!if $[abs(($slib_xb-$slib_xe))] < 1.E-6
  !set slib_xc=$[$slib_xb+$slib_h]
  !set slib_yc=$slib_ym
!else
  !set slib_r=$[($slib_ye-($slib_yb))/($slib_xe-($slib_xb))]
  !set slib_xc=$[-($slib_h*$slib_r/(1+($slib_r)^2)^0.5) + (($slib_xb+$slib_xe)/2)]
  !set slib_yc=$[($slib_h/((1+($slib_r)^2)^0.5)) + (($slib_yb+$slib_ye)/2)]
!endif

slib_out=$slib_xc,$slib_yc,$slib_xy

!!for debug
!!!slib_out=!exec canvasdraw\
!!!size 400,400\
!!!xrange $(slib_zone[1,2])\
!!!yrange $(slib_zone[3,4])\
!!!grid 1,1,grey\
!!!linewidth 4\
!!!points red, $slib_xyinitial\
!!!linewidth 1\
!!!curvedarrows blue,$(slib_out[3,4,1,2,5,6])