From 807f39bf1e757fa42c757ce7f96a938ecd8f0368 Mon Sep 17 00:00:00 2001
From: Volker Springel <vspringel@mpa-garching.mpg.de>
Date: Fri, 23 Apr 2021 13:00:32 +0200
Subject: [PATCH] fixed small inconsistency in PRESSURE_ENTROPY_SPH when
 REUSE_HYDRO_ACCELERATIONS_FROM_PREVIOUS_STEP is not active

---
 src/time_integration/kicks.cc | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/time_integration/kicks.cc b/src/time_integration/kicks.cc
index bc2697b..1d552ae 100644
--- a/src/time_integration/kicks.cc
+++ b/src/time_integration/kicks.cc
@@ -557,6 +557,9 @@ void sim::hydro_force(int step_indicator)
       if(step_indicator == SECOND_HALF_STEP)
         {
           Sp.SphP[target].EntropyPred = Sp.SphP[target].Entropy;
+#ifdef PRESSURE_ENTROPY_SPH
+          Sp.SphP[target].EntropyToInvGammaPred = pow(Sp.SphP[target].EntropyPred, 1.0 / GAMMA);
+#endif
           Sp.SphP[target].set_thermodynamic_variables();
 
           Sp.SphP[target].VelPred[0] += (Sp.SphP[target].HydroAccel[0] - Old[i].HydroAccel[0]) * dt_hydrokick;
-- 
GitLab