Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header

 !goto proc

!endif

slib_title=Geometric mean of data

slib_parms=1\

 ,[data] or [data],[weight]

slib_author=Adeline Grelot

slib_out=geometric mean,total weight

slib_comment=The numbers must be strictly positive. \

 Data and weight can be recognized under many formats, \

 in particular as a matrix (with data before weight).

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