[physics] Detect rank_deficiency for meshes in make_incompressible()

tum-pbs · Jan 4, 2025 · 2d66768 · 2d66768
1 parent e181b87
commit 2d66768
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/phi/physics/fluid.py b/phi/physics/fluid.py
@@ -163,6 +163,8 @@ def make_incompressible(velocity: Field,
         div = field.where(field.is_finite(div), div, 0)
     if not input_velocity.extrapolation.is_flexible and all_active:
         solve = solve.with_preprocessing(_balance_divergence, active)
+        if solve.rank_deficiency is None and not div.is_grid:
+            solve = copy_with(solve, rank_deficiency=1)
     if solve.x0 is None:
         pressure_extrapolation = _pressure_extrapolation(input_velocity.extrapolation)
         solve = copy_with(solve, x0=Field(div.geometry, 0, pressure_extrapolation))  # convert=False