Subversion Repositories wimsdev

Rev

Rev 4330 | Blame | Compare with Previous | Last modification | View Log | RSS feed

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

slib_author=Adeline, Grelot; Bernadette, PERRIN-RIOU

slib_example=[2,-7,6]\
[2,1;-7,1;6,4]\
[2,-7,6],[1,1,4]\
[2,-7,6;1,1,4]
slib_require=pari
!exit

:proc

!readproc slib/stat/dataproc $wims_read_parm
!if $slib_cnt=0
  slib_out=0
!else
  !if $slib_weight =$empty
    slib_out=!sum x^2 for x in $slib_data
    slib_out=$[sqrt($slib_out/$slib_cnt)]
  !else
    slib_cnt=!itemcnt $slib_data
    slib_out=!exec pari print(sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x])^2)/$slib_tw))

  !endif weight
!endif

slib_out=!trim $slib_out