Commit 41436dcb authored by Loic Hausammann's avatar Loic Hausammann Committed by Loic Hausammann
Browse files

GEAR: remove sigma in SFR

parent 389d10d8
......@@ -103,8 +103,8 @@ pressure_floor_get_comoving_pressure(const struct part* p,
/* Compute the pressure floor */
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;
pressure_floor_props.constants * cosmo->a_inv;
floor -= p->pressure_floor_data.sigma2;
floor *= a_coef * rho * hydro_one_over_gamma;
return fmaxf(pressure_comoving, floor);
......
......@@ -75,11 +75,10 @@ INLINE static int star_formation_is_star_forming(
}
/* Get the required variables */
const float sigma2 = p->pressure_floor_data.sigma2 * cosmo->a * cosmo->a;
const float density = hydro_get_physical_density(p, cosmo);
const float n_jeans_2_3 = starform->n_jeans_2_3;
const float h = p->h * kernel_gamma * cosmo->a;
const float density = hydro_get_physical_density(p, cosmo);
// TODO use GRACKLE */
const float mu = hydro_props->mu_neutral;
......@@ -88,9 +87,8 @@ INLINE static int star_formation_is_star_forming(
const float coef =
M_PI_4 / (phys_const->const_newton_G * n_jeans_2_3 * h * h);
const float density_criterion =
coef * (hydro_gamma * phys_const->const_boltzmann_k * temperature /
(mu * phys_const->const_proton_mass) +
sigma2);
coef * hydro_gamma * phys_const->const_boltzmann_k * temperature /
(mu * phys_const->const_proton_mass);
/* Check the density criterion */
return density > density_criterion;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment