Commit 35a653b7 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Same correction to dS/dt in vector mode for GADGET_SPH

parent 6d077c19
......@@ -644,7 +644,7 @@ __attribute__((always_inline)) INLINE static void runner_iact_vec_force(
pi[k]->force.v_sig = fmaxf(pi[k]->force.v_sig, v_sig.f[k]);
pj[k]->force.v_sig = fmaxf(pj[k]->force.v_sig, v_sig.f[k]);
pi[k]->entropy_dt += entropy_dt.f[k] * mj.f[k];
pj[k]->entropy_dt -= entropy_dt.f[k] * mi.f[k];
pj[k]->entropy_dt += entropy_dt.f[k] * mi.f[k];
}
#else
......
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