Subversion Repositories wimsdev

Rev

Rev 4351 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1
!if $wims_read_parm!=slib_header
2
!goto proc
3
!endif
4
 
4351 bpr 5
slib_author=Bernadette, Perrin-Riou
4158 bpr 6
 
20 reyssat 7
slib_example=[1,2;3,4], [1,2;3,4],[1]
3265 bpr 8
slib_require=pari
20 reyssat 9
!exit
10
 
11
:proc
12
slib_cnt=!itemcnt $wims_read_parm
13
!for slib_i=1 to $slib_cnt
14
  slib_mat$slib_i = !item $slib_i of $wims_read_parm
7692 bpr 15
!next
20 reyssat 16
 
17
slib_out=
18
 
19
slib_wm=
20
slib_wn=
7692 bpr 21
!!inutile ?
20 reyssat 22
slib_data=
23
 
24
!for slib_i=1 to $slib_cnt
25
slib_N=!declosing $(slib_mat$slib_i)
26
slib_N = !translate internal ; to $\
27
$ in $slib_N
28
slib_N= !translate internal $\
29
$ to ; in $slib_N
30
slib_data=!append item [$slib_N] to $slib_data
31
!next slib_i
32
 
33
slib_out=!exec pari {slib_D=[$slib_data]; slib_R=Mat(slib_D[1]);\
34
slib_a=matsize(slib_R)[1];\
35
for(slib_i=2,$slib_cnt, slib_C= Mat(slib_D[slib_i]); slib_b=matsize(slib_C)[1];\
36
slib_R=mattranspose(concat(slib_R, matrix(slib_a,slib_b)));\
37
slib_S= mattranspose(concat(matrix(slib_b,slib_a), slib_C));\
38
slib_R=concat(slib_R,slib_S)~ ; slib_a=slib_a+slib_b);\
39
print(slib_R)}
40
 
41
slib_out=!trim $slib_out