Subversion Repositories wimsdev

Rev

Rev 7692 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1
!if $wims_read_parm!=slib_header
13585 bpr 2
  !goto proc
20 reyssat 3
!endif
4349 bpr 4
slib_author=Gang, XIAO
4158 bpr 5
 
20 reyssat 6
slib_example= 3,6
3265 bpr 7
slib_require=pari
20 reyssat 8
!exit
9
 
10
:proc
11
!distribute items $wims_read_parm into slib_size,slib_range,slib_ring
12
!default slib_ring=Z
13
!bound slib_size between integer 2 and 20 default 2
14
!bound slib_range between integer 1 and 1000000 default 2
15
 
16
!if $slib_ring=Z
13585 bpr 17
  slib_sh1=!shuffle $slib_size
18
  slib_sh2=!shuffle $slib_size
19
  slib_ran1=(random(2*$slib_range)*(2*random(2)-1))
20
  slib_ran2=(random($slib_range)*(4*random(2)-2))
21
  slib_ran3=(random($slib_range)*(4*random(2)-2)+2*random(2)-1)
22
  slib_coef=if(x==y,$slib_ran3,if(x>y,$slib_ran2,$slib_ran1))
23
  slib_M1=vecextract(matrix($slib_size,$slib_size,x,y,$slib_coef),[$slib_sh1],[$slib_sh2])
24
  slib_M2=matrix($slib_size,$slib_size,x,y,$slib_ran1)
25
  slib_out=!exec pari M2=$slib_M2;if(matdet(M2)<>0,print(M2),print($slib_M1))
26
  !exit
20 reyssat 27
!endif