From c78ba29b0792c5129c6fd7f2c38d87310ae07d43 Mon Sep 17 00:00:00 2001
From: loikki <loic.hausammann@protonmail.ch>
Date: Wed, 14 Aug 2019 15:34:36 +0200
Subject: [PATCH] Add kernel_gamma to the pressure floor
---
src/pressure_floor/GEAR/pressure_floor.h | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/pressure_floor/GEAR/pressure_floor.h b/src/pressure_floor/GEAR/pressure_floor.h
index 873615cf7f..ca0cec5f59 100644
--- a/src/pressure_floor/GEAR/pressure_floor.h
+++ b/src/pressure_floor/GEAR/pressure_floor.h
@@ -66,11 +66,11 @@ struct pressure_floor_properties {
__attribute__((always_inline)) static INLINE float pressure_floor_get_physical_pressure(
const struct part* p, const float pressure_physical, const struct cosmology *cosmo) {
- const float h_phys = p->h * cosmo->a_inv;
+ const float H_phys = p->h * cosmo->a_inv * kernel_gamma;
const float rho = hydro_get_physical_density(p, cosmo);
/* Compute the pressure floor */
- float floor = h_phys * h_phys * rho * pressure_floor_props.constants -
+ float floor = H_phys * H_phys * rho * pressure_floor_props.constants -
p->pressure_floor_data.sigma2;
floor *= rho * hydro_one_over_gamma;
@@ -94,8 +94,8 @@ __attribute__((always_inline)) static INLINE float pressure_floor_get_comoving_p
const float rho = hydro_get_comoving_density(p);
/* Compute the pressure floor */
- float floor = p->h * p->h * rho * pressure_floor_props.constants -
- p->pressure_floor_data.sigma2 * cosmo->a * cosmo->a;
+ float floor = kernel_gamma * kernel_gamma * p->h * p->h * rho * pressure_floor_props.constants;
+ floor -= p->pressure_floor_data.sigma2 * cosmo->a * cosmo->a;
floor *= a_coef * rho * hydro_one_over_gamma;
return fmaxf(pressure_comoving, floor);
--
GitLab