(bind-func Cmag (lambda (num:gsl_complex*) (let ((real (aref (tref num 0) 0)) (imag (aref (tref num 0) 1))) (sqrt (+ (* real real) (* imag imag))))))