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