Subversion Repositories wimsdev

Rev

Rev 12676 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1
!if $wims_read_parm!=slib_header
13587 bpr 2
  !goto proc
20 reyssat 3
!endif
4158 bpr 4
 
4351 bpr 5
slib_author=Bernadette, Perrin-Riou
4158 bpr 6
 
20 reyssat 7
slib_example=[0,1,0,0;1,0,0,0;0,0,0,1;0,0,1,0]
3265 bpr 8
slib_require=pari
20 reyssat 9
!exit
10
 
11
:proc
12
 
13
slib_G= $wims_read_parm
14
slib_G=!declosing $slib_G
8560 bpr 15
slib_G=!translate internal $\
16
$ to ; in $slib_G
20 reyssat 17
 
18
slib_n=!exec pari slib_N=Mat([$slib_G]);matsize(slib_N)[1]
19
 
8560 bpr 20
slib_G= !exec pari slib_N=Mat([$slib_G]); matrix($slib_n,$slib_n,i,j, max(slib_N[i,j],slib_N[j,i]))
20 reyssat 21
 
22
slib_cnt=0
23
slib_liste=
24
slib_matrix=
25
slib_out=
26
!for slib_i=1 to $slib_n
12676 bpr 27
  !if $slib_i notitemof $slib_liste
28
    !readproc slib/graph/connexcomponent [$slib_G],$slib_i
29
    slib_liste=!append item $slib_out to $slib_liste
30
    !advance slib_cnt
31
    slib_matrix=!append line $slib_out to $slib_matrix
32
    slib_out=
33
  !endif
7692 bpr 34
!next slib_i
20 reyssat 35
 
36
slib_matrix=!translate internal $\
8560 bpr 37
$ to ; in $slib_matrix
20 reyssat 38
 
39
slib_out= $slib_cnt,[$slib_matrix]