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_parms=1\
 5 
slib_parms=1\
6 
 ,[data] or [data],[weight]
 6 
 ,[data] or [data],[weight]
7 
slib_author=Adeline, Grelot
 7 
slib_author=Adeline, Grelot
Line 14... Line 14...
14 
 
 14 
 
15 
:proc
 15 
:proc
16 
!readproc slib/stat/dataproc $wims_read_parm
 16 
!readproc slib/stat/dataproc $wims_read_parm
17 
 
 17 
 
18 
!if $slib_cnt=0
 18 
!if $slib_cnt=0
19 
slib_out=1
 19 
  slib_out=1
20 
!else
 20 
!else
21 
!if $slib_weight =$empty
 21 
  !if $slib_weight =$empty
22 
slib_data=!product x for x in $slib_data
 22 
    slib_data=!product x for x in $slib_data
23 
slib_out=$[($slib_data)^(1/$slib_cnt)],$slib_cnt
 23 
    slib_out=$[($slib_data)^(1/$slib_cnt)],$slib_cnt
24 
!else
 24 
    !else
25 
slib_out=!exec pari print(prod(x=1,$slib_cnt,[$slib_data][x]^[$slib_weight][x])^(1/$slib_tw))
 25 
    slib_out=!exec pari print(prod(x=1,$slib_cnt,[$slib_data][x]^[$slib_weight][x])^(1/$slib_tw))
26 
slib_out=$slib_out,$slib_tw
 26 
    slib_out=$slib_out,$slib_tw
27 
 
 - 
 
28 
!endif weight
 27 
  !endif weight
29 
!endif
 28 
!endif
30 
slib_out =!trim $slib_out
 29 
slib_out =!trim $slib_out