Rev 20 | Rev 4158 | 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 | slib_title=concatenation |
||
| 5 | slib_parms=1\ |
||
| 6 | ,list of square matrices under matrix format\ |
||
| 7 | |||
| 8 | slib_author=Bernadette PERRIN-RIOU |
||
| 9 | slib_out=block diagonal matrix with input matrices on the diagonal |
||
| 10 | slib_comment= |
||
| 11 | slib_example=[1,2;3,4], [1,2;3,4],[1] |
||
| 3265 | bpr | 12 | slib_require=pari |
| 20 | reyssat | 13 | !exit |
| 14 | |||
| 15 | :proc |
||
| 16 | slib_cnt=!itemcnt $wims_read_parm |
||
| 17 | !for slib_i=1 to $slib_cnt |
||
| 18 | slib_mat$slib_i = !item $slib_i of $wims_read_parm |
||
| 19 | !next |
||
| 20 | |||
| 21 | slib_out= |
||
| 22 | |||
| 23 | slib_wm= |
||
| 24 | slib_wn= |
||
| 25 | !!inutile ? |
||
| 26 | slib_data= |
||
| 27 | |||
| 28 | !for slib_i=1 to $slib_cnt |
||
| 29 | slib_N=!declosing $(slib_mat$slib_i) |
||
| 30 | slib_N = !translate internal ; to $\ |
||
| 31 | $ in $slib_N |
||
| 32 | slib_N= !translate internal $\ |
||
| 33 | $ to ; in $slib_N |
||
| 34 | slib_data=!append item [$slib_N] to $slib_data |
||
| 35 | !next slib_i |
||
| 36 | |||
| 37 | slib_out=!exec pari {slib_D=[$slib_data]; slib_R=Mat(slib_D[1]);\ |
||
| 38 | slib_a=matsize(slib_R)[1];\ |
||
| 39 | for(slib_i=2,$slib_cnt, slib_C= Mat(slib_D[slib_i]); slib_b=matsize(slib_C)[1];\ |
||
| 40 | slib_R=mattranspose(concat(slib_R, matrix(slib_a,slib_b)));\ |
||
| 41 | slib_S= mattranspose(concat(matrix(slib_b,slib_a), slib_C));\ |
||
| 42 | slib_R=concat(slib_R,slib_S)~ ; slib_a=slib_a+slib_b);\ |
||
| 43 | print(slib_R)} |
||
| 44 | |||
| 45 | slib_out=!trim $slib_out |