Rev 4161 | Rev 4352 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
||
3 | !endif |
||
4161 | bpr | 4 | |
4330 | bpr | 5 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
4161 | bpr | 6 | |
20 | reyssat | 7 | slib_example=[2,3,5],5 |
3265 | bpr | 8 | slib_require=pari |
20 | reyssat | 9 | !exit |
10 | |||
11 | :proc |
||
12 | slib_out= |
||
13 | |||
14 | !distribute item $wims_read_parm into slib_data,slib_degree |
||
15 | !default slib_degree= 1 |
||
16 | slib_data= !declosing $slib_data |
||
17 | slib_n= !itemcnt $slib_data |
||
18 | slib_u= !exec pari gamma($slib_degree*0.5)^(-1) |
||
19 | !for slib_i in $slib_data |
||
20 | !if $slib_i > 0 |
||
21 | slib_val=$[0.5*exp(-0.5*$slib_i)*(0.5*$slib_i)^( $slib_degree/2-1)*$slib_u] |
||
22 | !else |
||
23 | slib_val=0 |
||
24 | !endif |
||
25 | slib_out=!append item $slib_val to $slib_out |
||
26 | !next i |
||
27 | slib_out=!nonempty items $slib_out |
||
28 | slib_out= !trim $slib_out |