Commit 807f39bf authored by Volker Springel's avatar Volker Springel
Browse files

fixed small inconsistency in PRESSURE_ENTROPY_SPH when...

fixed small inconsistency in PRESSURE_ENTROPY_SPH when REUSE_HYDRO_ACCELERATIONS_FROM_PREVIOUS_STEP is not active
parent 4d022ab5
......@@ -557,6 +557,9 @@ void sim::hydro_force(int step_indicator)
if(step_indicator == SECOND_HALF_STEP)
{
Sp.SphP[target].EntropyPred = Sp.SphP[target].Entropy;
#ifdef PRESSURE_ENTROPY_SPH
Sp.SphP[target].EntropyToInvGammaPred = pow(Sp.SphP[target].EntropyPred, 1.0 / GAMMA);
#endif
Sp.SphP[target].set_thermodynamic_variables();
Sp.SphP[target].VelPred[0] += (Sp.SphP[target].HydroAccel[0] - Old[i].HydroAccel[0]) * dt_hydrokick;
......
Markdown is supported
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