(bind-func vmedian (lambda (buf:float* len:i64) (let ((h (vsort buf len)) (idx (dtoi64 (ceil (/ (i64tod len) 2.0))))) (if (evenp len) (/ (+ (pref h idx) (pref h (- idx 1))) 2.0) (pref h idx)))))