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