Rev 4349 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 4349 | Rev 13585 | ||
---|---|---|---|
Line 1... | Line 1... | ||
1 | !if $wims_read_parm!=slib_header |
1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
2 | !goto proc |
3 | !endif |
3 | !endif |
4 | 4 | ||
5 | slib_author=Gang, XIAO |
5 | slib_author=Gang, XIAO |
6 | 6 | ||
7 | slib_example=3, 5 |
7 | slib_example=3, 5 |
Line 13... | Line 13... | ||
13 | !default slib_ring=Z |
13 | !default slib_ring=Z |
14 | !bound slib_size between integer 2 and 20 default 2 |
14 | !bound slib_size between integer 2 and 20 default 2 |
15 | !bound slib_range between integer 1 and 1000000 default 2 |
15 | !bound slib_range between integer 1 and 1000000 default 2 |
16 | 16 | ||
17 | !if $slib_ring=Z |
17 | !if $slib_ring=Z |
18 | slib_sh1=!shuffle $slib_size |
18 | slib_sh1=!shuffle $slib_size |
19 | slib_sh2=!shuffle $slib_size |
19 | slib_sh2=!shuffle $slib_size |
20 | slib_ran1=(random($slib_range)*(2*random(2)-1)) |
20 | slib_ran1=(random($slib_range)*(2*random(2)-1)) |
21 | slib_ran2=(2*random(2)-1) |
21 | slib_ran2=(2*random(2)-1) |
22 | slib_coefu=if(x==y,$slib_ran2,if(x>y,0,$slib_ran1)) |
22 | slib_coefu=if(x==y,$slib_ran2,if(x>y,0,$slib_ran1)) |
23 | slib_coefd=if(x==y,$slib_ran2,if(x<y,0,$slib_ran1)) |
23 | slib_coefd=if(x==y,$slib_ran2,if(x<y,0,$slib_ran1)) |
24 | slib_triu=matrix($slib_size,$slib_size,x,y,$slib_coefu) |
24 | slib_triu=matrix($slib_size,$slib_size,x,y,$slib_coefu) |
25 | slib_trid=vecextract(matrix($slib_size,$slib_size,x,y,$slib_coefd),[$slib_sh1],[$slib_sh2]) |
25 | slib_trid=vecextract(matrix($slib_size,$slib_size,x,y,$slib_coefd),[$slib_sh1],[$slib_sh2]) |
26 | slib_out=!exec pari print($slib_triu*$slib_trid) |
26 | slib_out=!exec pari print($slib_triu*$slib_trid) |
27 | !exit |
27 | !exit |
28 | !endif |
28 | !endif |
29 | 29 |