Subversion Repositories wimsdev

Rev

Rev 4330 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4330 Rev 7692
Line 4... Line 4...
4 
 
 4 
 
5 
slib_parms=2\
 5 
slib_parms=2\
6 
  2,m\
 6 
  2,m\
7 
  2*m,N\
 7 
  2*m,N\
8 
  R, K
 8 
  R, K
9
9 
 
10 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
 10 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
11 
 
 11 
 
12 
slib_example= 5,10\
 12 
slib_example= 5,10\
13 
7,10,Q
 13 
7,10,Q
14 
slib_require=pari
 14 
slib_require=pari
Line 17... Line 17...
17 
!distribute item $wims_read_parm into slib_mm,slib_NN,slib_K
 17 
!distribute item $wims_read_parm into slib_mm,slib_NN,slib_K
18 
 
 18 
 
19 
!default slib_mm=2
 19 
!default slib_mm=2
20 
!default slib_NN=$[2*$slib_mm]
 20 
!default slib_NN=$[2*$slib_mm]
21 
!default slib_K= R
 21 
!default slib_K= R
22 
 
 22 
 
23 
!if $slib_mm = 1
 23 
!if $slib_mm = 1
24 
 slib_out = 1
 24 
 slib_out = 1
25 
!else
 25 
!else
26 
 slib_cc = !values 1/$slib_mm for x = 1 to $[$slib_mm-1]
 26 
 slib_cc = !values 1/$slib_mm for x = 1 to $[$slib_mm-1]
27 
 
 27 
 
28 
 !readproc slib/stat/multinomial 1, $[$slib_NN], $slib_cc
 28 
 !readproc slib/stat/multinomial 1, $[$slib_NN], $slib_cc
29 
 slib_vv= $slib_out
 29 
 slib_vv= $slib_out
30 
 
 30 
 
31 
 !if $slib_K issametext Q
 31 
 !if $slib_K issametext Q
32 
  slib_out=!exec pari slib_a=Vec([$slib_vv]);vector($slib_mm, x,slib_a[x]/$slib_NN)
32 
  slib_out=!exec pari slib_a=Vec([$slib_vv]);vector($slib_mm, x,slib_a[x]/$slib_NN)
33 
 !else
 33 
 !else
34 
  slib_out=!values x/$slib_NN for x in $slib_vv
 34 
  slib_out=!values x/$slib_NN for x in $slib_vv
35 
 !endif
35 
 !endif
36 
!endif
 36 
!endif
37 
 
 - 
 
38 
 
 37 
 
39 
slib_out =!trim $slib_out
 38 
slib_out =!trim $slib_out