Subversion Repositories wimsdev

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

slib_author=Bernadette, Perrin-Riou

slib_example=5,2,[a,b,c,d],blue,red\
5,2,[a,b,c,d],blue,red,html

slib_require=pari
!exit

:proc

!distribute item $wims_read_parm into slib_n,slib_s,slib_text,slib_color1,slib_color2,slib_option

!default slib_s=2
slib_color1=!declosing $slib_color1
slib_color2=!declosing $slib_color2
!default slib_color1=blue
!default slib_color2=green
!!default slib_text= !values x for x=1 to $slib_n

slib_text=!declosing $slib_text
!if [ isin $slib_s
  slib_G=!declosing $slib_s
!else
  slib_liste=
  !for slib_i = 1 to $slib_n
    slib_liste1=!values 1 for x=1 to $slib_s
    slib_liste0=!values 0 for x=1 to $[$slib_n-$slib_i-$slib_s]
    slib_liste2=!shuffle $slib_liste1,$slib_liste0
    slib_liste2=1,$slib_liste2
    slib_liste2=!item 1 to $[$slib_n-$slib_i+1] of $slib_liste2
    slib_liste=!append item $slib_liste2 to $slib_liste
  !next slib_i
  slib_liste=!nonempty items $slib_liste
  slib_G= !exec pari slib_n =$slib_n;slib_liste =[$slib_liste];slib_M=matrix(slib_n,slib_n,slib_i,slib_j, if(slib_i < slib_j,slib_liste[slib_n*(slib_n-1)/2-(slib_n-slib_i+1)*(slib_n-slib_i)/2 +slib_j-slib_i]));for(slib_i = 2,slib_n-1,if(slib_M[,slib_i]==0,slib_M[1,slib_i]=1);if(slib_M[slib_i,]==0,slib_M[slib_i,slib_n]=1));slib_M[1,slib_n]=0; slib_M

!endif
slib_GG=!translate internal ; to $\
$ in $slib_G

##########################
slib_niveau=!exec pari slib_n=$slib_n;slib_M=Mat([$slib_G]);slib_niveau=matrix(slib_n,slib_n);slib_niveau[1,1]=1;slib_degp=vector(slib_n,i,sum(x=1,slib_n,slib_M[x,i]));for(slib_k=1,slib_n-1,for(slib_i=1,slib_n,if(slib_niveau[slib_k,slib_i]!=0,for(slib_j=1,slib_n,if(slib_M[slib_i,slib_j]!=0,slib_degp[slib_j]-=1;if(slib_degp[slib_j]==0,slib_niveau[slib_k+1,slib_j]=1));););););slib_niveau
slib_height=!exec pari slib_n=$slib_n;slib_niveau=Mat([$slib_niveau]);slib_k=slib_n;for(slib_i=1,slib_n,if(slib_niveau[slib_i,]==0,slib_k=slib_i-1;break));print(slib_k);

slib_width =!exec pari slib_niveau=Mat([$slib_niveau]); vecmax(vector($slib_height,i, sum(j=1,$slib_n, slib_niveau[i,j])))

#hauteur des paraboles sur lesquelles sont mis les points
slib_listeb= !values ($slib_width/2+1-x)*0.98111111/$slib_width for x= 1 to $[ceil($slib_width/2)]
slib_listeb =!makelist x,-x for x in $slib_listeb
slib_niveau =!exec pari slib_niveau=Mat([$slib_niveau]); matrix($slib_height,$slib_n,i,j,slib_niveau[i,j])
slib_position= !exec pari matrix($slib_n,2)
slib_position=!translate internal ; to $\
$ in $slib_position
slib_niveau1=!translate internal ; to $\
$ in $slib_niveau
slib_sign=1
!for slib_i=1 to $slib_height
  slib_eps1=!item $slib_i of $slib_eps
  slib_ligne= !line $slib_i of $slib_niveau1
  slib_test=0
  slib_total =!replace internal , by + in $slib_ligne
  slib_total=$[$slib_total]
    !if $slib_total=1
      slib_sign=$[-$slib_sign]
    !endif

  slib_dist= $[round($slib_width/$slib_total)]
  !for slib_j = 1 to $slib_n
    slib_pt=!item $slib_j of $slib_ligne
    !if $slib_pt <>0
      !advance slib_test
      slib_bb=!item $[($slib_test)*$slib_dist] of $slib_listeb
      slib_position=!replace internal line number $slib_j by $[$slib_i-1], $[$slib_sign*$slib_bb*$slib_i*($slib_height -$slib_i)] in $slib_position
    !endif
  !next slib_j
!next slib_i
#dessin
####################
slib_yrange1= $[-$slib_height^2/8-0.5]
slib_yrange2=$[$slib_height^2/8+0.5]
slib_dessin = xrange -0.5, $[$slib_height]\
yrange $slib_yrange1, $slib_yrange2\
parallel -0.5,$slib_yrange1,-0.5,$slib_yrange2,1,0, $[$slib_height+2],yellow
####################

!for slib_i = 0 to $[floor($slib_height/2+2)]
  slib_dessin=!append line pointfill 2*$slib_i,0,3,3,yellow to $slib_dessin
!next slib_i
# laisser en premier les lignes correspondant aux sommets
!for slib_i = 1 to $slib_n
  slib_point=!line $slib_i of $slib_position
  slib_dessin =!append line disk $slib_point, 6, $slib_color2 to $slib_dessin
  !for slib_j = $slib_i to $slib_n
    slib_ligne=!line $slib_i of $slib_GG
    slib_test=!item $slib_j of $slib_ligne
    !if $slib_test <> 0
      slib_point1=!line $slib_j of $slib_position
      slib_interm= !exec pari ([$slib_point]+2*[$slib_point1])/3
      slib_dessin = !append line segment $slib_point, $slib_point1, $slib_color1 to $slib_dessin
      slib_dessin = !append line arrow $slib_point, $slib_interm, 10, $slib_color1 to $slib_dessin
    !endif
  !next slib_j
  !if $slib_text<>$empty
    slib_texti=!item $slib_i of $slib_text
    slib_dessin=!append line text black,$slib_point, medium, $slib_texti to $slib_dessin
  !endif
!next slib_i

slib_out = $slib_dessin
!if url iswordof $slib_option or html iswordof $slib_option
  insdraw_size=450,250
  !insdraw $slib_out
  slib_out=$ins_url
  !if html iswordof $slib_option
    slib_out=<img src="$ins_url" alt=""/>
  !else
    slib_out=$slib_out,[$slib_G],[$slib_niveau], [$slib_position],$insdraw_size
  !endif
!else
  slib_out=$slib_dessin,[$slib_G],[$slib_niveau], [$slib_position]
!endif