[physics] Fix MacCormack advection

tum-pbs · Nov 25, 2024 · e0c0b7d · e0c0b7d
1 parent d7f7c82
commit e0c0b7d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/phi/physics/advect.py b/phi/physics/advect.py
@@ -200,7 +200,7 @@ def mac_cormack(field: Field,
     Returns:
         Advected field of type `type(field)`
     """
-    v0 = resample(velocity, to=field).values
+    v0 = sample(velocity, field.geometry, at=field.sampled_at, boundary=field.boundary)
     points_bwd = integrator(field, velocity, -dt, v0=v0)
     points_fwd = integrator(field, velocity, dt, v0=v0)
     # --- forward+backward semi-Lagrangian advection ---