(bind-func fluid-diffuse (lambda (b:i64 x:double* x0:double* diff:double dt:double iter Ny N) (let ((a:double (* dt diff (i64tod (- N 2))))) (fluid-lin-solve b x x0 a (+ 1.0 (* 6.0 a)) iter Ny N))))