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