Skip to content
Snippets Groups Projects
Commit cc470f71 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Make use of the Plummer conversion constant for the time-step as well.

parent d3331ea8
No related branches found
No related tags found
1 merge request!516More cosmology work
...@@ -106,7 +106,6 @@ gravity_compute_timestep_self(const struct gpart* const gp, ...@@ -106,7 +106,6 @@ gravity_compute_timestep_self(const struct gpart* const gp,
const float epsilon = gravity_get_softening(gp, grav_props); const float epsilon = gravity_get_softening(gp, grav_props);
/* Note that 0.66666667 = 2. (from Gadget) / 3. (Plummer softening) */
const float dt = sqrtf(2. * kernel_gravity_softening_plummer_equivalent_inv * const float dt = sqrtf(2. * kernel_gravity_softening_plummer_equivalent_inv *
cosmo->a * grav_props->eta * epsilon * ac_inv); cosmo->a * grav_props->eta * epsilon * ac_inv);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment