diff --git a/src/hydro/Gizmo/hydro_iact.h b/src/hydro/Gizmo/hydro_iact.h index 0593408915766f0316382fc200de7b51bd3da7c4..bdcea49e4fdc2a9ce17c43b99dab55d61f4095c8 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]);