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