Subversion Repositories wimsdev

Rev

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