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 |