Commit c7378f2c authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Time integration mistake ?

parent 091a434f
......@@ -218,8 +218,8 @@ __attribute__((always_inline)) INLINE static void hydro_kick_extra(
p->entropy *= 0.5f;
/* Do not 'overcool' when timestep increases */
if (p->entropy + 0.5f * p->entropy_dt * dt < 0.5f * p->entropy)
p->entropy_dt = -0.5f * p->entropy / dt;
if (p->entropy + p->entropy_dt * half_dt < 0.5f * p->entropy)
p->entropy_dt = -0.5f * p->entropy / half_dt;
}
/**
......
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