fluid-lin-solve   xtlang

(bind-func fluid-lin-solve
  (lambda (b:i64 x:double* x0:double* a c iter:i64 Ny:i64 N:i64)
    (let ((cRecip (/ 1.0 c))
    (k 0)
    (m 0)
    (j 0)
    (i 0))
      (dotimes (k iter)
  (dotimes (j (- Ny 2))
    (dotimes (i (- N 2))
      (pset! x (+ (+ i 1) (* (+ j 1) N))
       (* cRecip
          (+ (pref x0 (+ (+ i 1) (* (+ j 1) N)))
       (* a (+ (pref x (+ (+ i 2) (* (+ j 1) N)))
         (pref x (+ i (* (+ j 1) N)))
         (pref x (+ (+ i 1) (* (+ j 2) N)))
         (pref x (+ (+ i 1) (* j N)))))))))))
      (fluid-set-boundary b x Ny N))
    1))