From 571ab90eb3f21f72785ce46055ba98f8746ea4fc Mon Sep 17 00:00:00 2001
From: Orestis Karapiperis <karapiperis@strw.leidenuniv.nl>
Date: Sat, 5 Apr 2025 11:57:51 +0200
Subject: [PATCH] Small fix not yet adressed from when moving the Balsara
 computation to where the switch is needed

---
 src/hydro/SPHENIX/hydro.h | 15 ++++-----------
 1 file changed, 4 insertions(+), 11 deletions(-)

diff --git a/src/hydro/SPHENIX/hydro.h b/src/hydro/SPHENIX/hydro.h
index dc6826c49e..283be043ed 100644
--- a/src/hydro/SPHENIX/hydro.h
+++ b/src/hydro/SPHENIX/hydro.h
@@ -838,9 +838,9 @@ __attribute__((always_inline)) INLINE static void hydro_prepare_force(
   const float pressure = hydro_get_comoving_pressure(p);
   const float pressure_including_floor =
       pressure_floor_get_comoving_pressure(p, pressure_floor, pressure, cosmo);
-  const float soundspeed_physical =
-      gas_soundspeed_from_pressure(p->rho, pressure_including_floor) *
-      cosmo->a_factor_sound_speed;
+  const float soundspeed =
+      gas_soundspeed_from_pressure(p->rho, pressure_including_floor);
+  const float soundspeed_physical = soundspeed * cosmo->a_factor_sound_speed;
   const float fac_B = cosmo->a_factor_Balsara_eps;
   
   const float sound_crossing_time_inverse =
@@ -853,14 +853,7 @@ __attribute__((always_inline)) INLINE static void hydro_prepare_force(
 
   /* Compute the norm of div v */
   const float abs_div_v = fabsf(p->viscosity.div_v);
-
-  /* Compute the sound speed  */
-  const float pressure = hydro_get_comoving_pressure(p);
-  const float pressure_including_floor =
-      pressure_floor_get_comoving_pressure(p, pressure_floor, pressure, cosmo);
-  const float soundspeed =
-      gas_soundspeed_from_pressure(p->rho, pressure_including_floor);
-
+  
   /* Get the squares of the quantities necessary for the Balsara-like switch */
   const float fac_B_2 = fac_B * fac_B;
   const float curl_v_2 = curl_v * curl_v;
-- 
GitLab