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