(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))))