Rev 4158 | Rev 4351 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 4158 | Rev 4349 | ||
|---|---|---|---|
| Line 20... | Line 20... | ||
| 20 | !readproc slib/matrix/givenrank $slib_n,$slib_size,$slib_n,$slib_range |
20 | !readproc slib/matrix/givenrank $slib_n,$slib_size,$slib_n,$slib_range |
| 21 | !!slib_matrix=!translate internal ; to $\ |
21 | !!slib_matrix=!translate internal ; to $\ |
| 22 | $ in $slib_out |
22 | $ in $slib_out |
| 23 | 23 | ||
| 24 | slib_matrix=$slib_out |
24 | slib_matrix=$slib_out |
| 25 | slib_out=!exec pari {slib_A=matid($slib_size);slib_n=$slib_n; |
25 | slib_out=!exec pari {slib_A=matid($slib_size);slib_n=$slib_n; \ |
| 26 | slib_M= Mat([$slib_matrix]);\\ |
26 | slib_M= Mat([$slib_matrix]);\\ |
| 27 | for(slib_j=1 , slib_n,slib_vv=slib_M[slib_j,]; slib_A=slib_A*(matid($slib_size)-2/norml2(slib_vv)*slib_vv~*slib_vv)); |
27 | for(slib_j=1 , slib_n,slib_vv=slib_M[slib_j,]; slib_A=slib_A*(matid($slib_size)-2/norml2(slib_vv)*slib_vv~*slib_vv));\ |
| 28 | print(slib_A)} |
28 | print(slib_A)} |
| 29 | 29 | ||
| 30 | slib_out= !trim $slib_out |
30 | slib_out= !trim $slib_out |
| 31 | slib_out=!line -1 of $slib_out |
31 | slib_out=!line -1 of $slib_out |
| 32 | 32 | ||