Rev 4328 | 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=Adeline, Grelot; Bernadette, PERRIN-RIOU |
| 4161 | bpr | 6 | |
| 20 | reyssat | 7 | slib_example=[2,-7,6]\ |
| 8 | [2,1;-7,1;6,4]\ |
||
| 9 | [2,-7,6],[1,1,4]\ |
||
| 10 | [2,-7,6;1,1,4] |
||
| 3265 | bpr | 11 | slib_require=pari |
| 20 | reyssat | 12 | !exit |
| 13 | |||
| 14 | :proc |
||
| 15 | |||
| 16 | !readproc slib/stat/dataproc $wims_read_parm |
||
| 17 | !if $slib_cnt=0 |
||
| 18 | slib_out=0 |
||
| 19 | !else |
||
| 20 | !if $slib_weight =$empty |
||
| 21 | slib_out=!sum x^2 for x in $slib_data |
||
| 22 | slib_out=$[sqrt($slib_out/$slib_cnt)] |
||
| 23 | !else |
||
| 24 | slib_cnt=!itemcnt $slib_data |
||
| 25 | slib_out=!exec pari print(sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x])^2)/$slib_tw)) |
||
| 26 | |||
| 27 | !endif weight |
||
| 28 | !endif |
||
| 29 | |||
| 30 | slib_out=!trim $slib_out |