diff --git a/src/Evolution/Systems/GrMhd/ValenciaDivClean/FixConservatives.cpp b/src/Evolution/Systems/GrMhd/ValenciaDivClean/FixConservatives.cpp index 81d24c9fd02d4..d3ea940d996e7 100644 --- 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.0 * tau_tilde; + 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