From ab0491e9280917e79ebaa149084e30d9bab4d35e Mon Sep 17 00:00:00 2001
From: Bert Vandenbroucke <bert.vandenbroucke@ugent.be>
Date: Fri, 3 Mar 2017 12:33:40 +0000
Subject: [PATCH] Made maximal signal velocity estimate for isothermal gas a
bit more conservative.
---
src/hydro/Gizmo/hydro_iact.h | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/hydro/Gizmo/hydro_iact.h b/src/hydro/Gizmo/hydro_iact.h
index 0593408915..bdcea49e4f 100644
--- a/src/hydro/Gizmo/hydro_iact.h
+++ b/src/hydro/Gizmo/hydro_iact.h
@@ -233,7 +233,9 @@ __attribute__((always_inline)) INLINE static void runner_iact_fluxes_common(
/* calculate the maximal signal velocity */
if (Wi[0] > 0.0f && Wj[0] > 0.0f) {
#ifdef EOS_ISOTHERMAL_GAS
- vmax = 2. * const_isothermal_soundspeed;
+ /* we use a value that is slightly higher than necessary, since the correct
+ value does not always work */
+ vmax = 2.5 * const_isothermal_soundspeed;
#else
vmax =
sqrtf(hydro_gamma * Wi[4] / Wi[0]) + sqrtf(hydro_gamma * Wj[4] / Wj[0]);
--
GitLab