fluid_advect   xtlang

(bind-func fluid_advect
  (lambda (b:i64 d:double* d0:double* velocx:double* velocy:double* dt:double Ny:i64 N:i64)
    (let ((n_1 (i64tod (- N 1)))
          (dtx (* dt n_1))
          (dty (* dt n_1))
          (jfloat 0.0) (ifloat 0.0)
          (s0 0.0) (s1 0.0) (t0 0.0) (t1 0.0)
          (i0 0.0) (i0i 0) (i1 0.0) (i1i 0)
          (j0 0.0) (j0i 0) (j1 0.0) (j1i 0)
          (j:i64 0) (i:i64 0)
          (tmp1 0.0) (tmp2 0.0) (tmp3 0.0)
          (x 0.0) (y 0.0) (z 0.0)
          (Nfloat (i64tod N)))
      (dotimes (j (- Ny 1))
        (set! jfloat (+ jfloat 1.0))
        (set! ifloat 0.0)
        (dotimes (i (- N 1))
          (set! ifloat (+ ifloat 1.0))
          (set! tmp1 (* dtx (pref velocx (+ (+ i 1) (* (+ j 1) N)))))
          (set! tmp2 (* dty (pref velocy (+ (+ i 1) (* (+ j 1) N)))))
          (set! x (- ifloat tmp1))
          (set! y (- jfloat tmp2))