Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header

 !goto proc

!endif

slib_title=Generates a random invertible matrix

slib_parms=3\

2,matrix size\

2,coefficient range\

Z,coefficient ring: Z for the time being,Q,R,C,G (Gauss ring), n (n is an integer for Z/nZ)

slib_author=Gang Xiao

slib_out=the generated matrix

slib_example= 3,6

slib_require=pari

!exit

:proc

!distribute items $wims_read_parm into slib_size,slib_range,slib_ring

!default slib_ring=Z

!bound slib_size between integer 2 and 20 default 2

!bound slib_range between integer 1 and 1000000 default 2

!if $slib_ring=Z

 slib_sh1=!shuffle $slib_size

 slib_sh2=!shuffle $slib_size

 slib_ran1=(random(2*$slib_range)*(2*random(2)-1))

 slib_ran2=(random($slib_range)*(4*random(2)-2))

 slib_ran3=(random($slib_range)*(4*random(2)-2)+2*random(2)-1)

 slib_coef=if(x==y,$slib_ran3,if(x>y,$slib_ran2,$slib_ran1))

 slib_M1=vecextract(matrix($slib_size,$slib_size,x,y,$slib_coef),[$slib_sh1],[$slib_sh2])

 slib_M2=matrix($slib_size,$slib_size,x,y,$slib_ran1)

 slib_out=!exec pari M2=$slib_M2;if(matdet(M2)<>0,print(M2),print($slib_M1))

 !exit

!endif