Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header
 !goto proc
!endif

slib_parms=1\
 ,[data] or [data],[weight]
slib_author=Adeline Grelot

slib_example=[3,5,6]\
[3,5,6],[4,1,1]\
[3,5,6;4,1,1]
slib_require=pari
!exit

:proc
!readproc slib/stat/dataproc $wims_read_parm

!if $slib_cnt=0
slib_out=1
!else
!if $slib_weight =$empty
slib_data=!product x for x in $slib_data
slib_out=$[($slib_data)^(1/$slib_cnt)],$slib_cnt
!else
slib_out=!exec pari print(prod(x=1,$slib_cnt,[$slib_data][x]^[$slib_weight][x])^(1/$slib_tw))
slib_out=$slib_out,$slib_tw

!endif weight
!endif
slib_out =!trim $slib_out