Rev 4330 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 4330 | Rev 13586 | ||
---|---|---|---|
Line 1... | Line 1... | ||
1 | !if $wims_read_parm!=slib_header |
1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
2 | !goto proc |
3 | !endif |
3 | !endif |
4 | 4 | ||
5 | slib_author=Adeline, Grelot; Bernadette, PERRIN-RIOU |
5 | slib_author=Adeline, Grelot; Bernadette, PERRIN-RIOU |
6 | 6 | ||
7 | slib_example=[2,-7,6]\ |
7 | slib_example=[2,-7,6]\ |
Line 13... | Line 13... | ||
13 | 13 | ||
14 | :proc |
14 | :proc |
15 | 15 | ||
16 | !readproc slib/stat/dataproc $wims_read_parm |
16 | !readproc slib/stat/dataproc $wims_read_parm |
17 | !if $slib_cnt=0 |
17 | !if $slib_cnt=0 |
18 | slib_out=0 |
18 | slib_out=0 |
19 | !else |
19 | !else |
20 | !if $slib_weight =$empty |
20 | !if $slib_weight =$empty |
21 | slib_out=!sum x^2 for x in $slib_data |
21 | slib_out=!sum x^2 for x in $slib_data |
22 | slib_out=$[sqrt($slib_out/$slib_cnt)] |
22 | slib_out=$[sqrt($slib_out/$slib_cnt)] |
23 | !else |
23 | !else |
24 | slib_cnt=!itemcnt $slib_data |
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)) |
25 | slib_out=!exec pari print(sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x])^2)/$slib_tw)) |
26 | 26 | ||
27 | !endif weight |
27 | !endif weight |
28 | !endif |
28 | !endif |
29 | 29 | ||
30 | slib_out=!trim $slib_out |
30 | slib_out=!trim $slib_out |