Subversion Repositories wimsdev

Rev

Rev 12676 | Rev 13103 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 12676 Rev 12678
Line 9... Line 9...
9 
slib_require=pari
 9 
slib_require=pari
10 
!exit
 10 
!exit
11 
 
 11 
 
12 
:proc
 12 
:proc
13 
!reset slib_m, slib_n
 13 
!reset slib_m, slib_n
14 
!distribute item $wims_read_parm into slib_n, slib_m
 14 
!distribute item $wims_read_parm into slib_n, slib_m, slib_option
15 
!default slib_m =!randint $[floor(1.5*$slib_n)], $[2*$slib_n]
 15 
!default slib_m =!randint $[floor(1.5*$slib_n)], $[2*$slib_n]
16 
slib_sommet=!values x for x=1 to $slib_n
 16 
slib_sommet=!values x for x=1 to $slib_n
17 
slib_u=!exec pari matrix($slib_n,$slib_n,slib_i,slib_j,if(slib_j ==slib_i,slib_i))
 17 
slib_u=!exec pari matrix($slib_n,$slib_n,slib_i,slib_j,if(slib_j==slib_i,slib_i))
18 
 
 18 
 
19 
slib_u=!translate internal ; to $\
 19 
slib_u=!translate internal ; to $\
20 
$  in $slib_u
 20 
$  in $slib_u
21 
slib_S=1
 21 
slib_S=1
22 
slib_euler=$slib_S
 22 
slib_euler=$slib_S
Line 61... Line 61...
61 
!for slib_i in $slib_euler
 61 
!for slib_i in $slib_euler
62 
  slib_pos=!positionof item $slib_i in $slib_p
 62 
  slib_pos=!positionof item $slib_i in $slib_p
63 
  slib_euler1=!append item $slib_pos to $slib_euler1
 63 
  slib_euler1=!append item $slib_pos to $slib_euler1
64 
!next slib_i
 64 
!next slib_i
65 
slib_euler=$slib_euler1
 65 
slib_euler=$slib_euler1
66 
slib_G=!exec pari slib_p = [$slib_p]; slib_M=Mat([$slib_u]); slib_M=matrix($slib_n,$slib_n,slib_i,slib_j,if(slib_i !=slib_j && slib_M[slib_i,slib_j] != 0,1)) ;\
 66 
slib_G=!exec pari slib_p = [$slib_p]; slib_M=Mat([$slib_u]); slib_M=matrix($slib_n,$slib_n,slib_i,slib_j,if(slib_i !=slib_j && slib_M[slib_i,slib_j] != 0,1)) ;matrix($slib_n,$slib_n,slib_i,slib_j,slib_M[slib_p[slib_i],slib_p[slib_j]])
- 
 
 67 
 
67 
matrix($slib_n,$slib_n,slib_i,slib_j,slib_M[slib_p[slib_i],slib_p[slib_j]])
 68 
!if $slib_option=1 and $(slib_euler[1])!=$(slib_euler[-1])
- 
 
 69 
  slib_G=!exec pari G={[$slib_G]};\
- 
 
 70 
    slib_e=[$slib_euler];\
- 
 
 71 
    G[slib_e[1],slib_e[$slib_m]]=1;\
- 
 
 72 
    G[slib_e[$slib_m],slib_e[1]]=1; G
- 
 
 73 
  slib_euler=$slib_euler, $(slib_euler[1])
- 
 
 74 
!endif
- 
 
 75 
 
68 
 
 76 
 
69 
slib_out=[$slib_G],[$slib_euler]
 77 
slib_out=[$slib_G],[$slib_euler]