Add lower bound of zero for TildeTau in FixConservative for zero B

jyoo1042 · Aug 15, 2024 · ed74557 · ed74557
1 parent 82e0ac0
commit ed74557
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/Evolution/Systems/GrMhd/ValenciaDivClean/FixConservatives.cpp b/src/Evolution/Systems/GrMhd/ValenciaDivClean/FixConservatives.cpp
@@ -327,6 +327,14 @@ bool FixConservatives::operator()(
       }
     }
 
+    else {
+      if (const auto tilde_tau_mask = tau_tilde > 0;
+          simd::any(tilde_tau_mask)) {
+        needed_fixing = true;
+        tau_tilde = simd::select(tilde_tau_mask, 0.0 * tau_tilde, tau_tilde);
+      }
+    }
+
     // Decrease momentum density if necessary
     auto s_tilde_squared = load(get(tilde_s_squared));
     // Equation B.24 of Foucart