From 7a73ec6d83bb528f05da2f12686491a0ecd37a1e Mon Sep 17 00:00:00 2001 From: Matthieu Schaller <matthieu.schaller@durham.ac.uk> Date: Thu, 21 Jul 2016 12:44:13 +0100 Subject: [PATCH] Don't update the signal velocity of particle j in non-sym routines. --- src/hydro/Default/hydro_iact.h | 2 -- 1 file changed, 2 deletions(-) diff --git a/src/hydro/Default/hydro_iact.h b/src/hydro/Default/hydro_iact.h index e4062f6a64..2748bdca4e 100644 --- a/src/hydro/Default/hydro_iact.h +++ b/src/hydro/Default/hydro_iact.h @@ -750,7 +750,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_nonsym_force( /* Update the signal velocity. */ pi->force.v_sig = fmaxf(pi->force.v_sig, v_sig); - pj->force.v_sig = fmaxf(pj->force.v_sig, v_sig); } /** @@ -946,7 +945,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_nonsym_vec_force( pi[k]->force.u_dt += piu_dt.f[k]; pi[k]->h_dt -= pih_dt.f[k]; pi[k]->force.v_sig = vi_sig.f[k]; - pj[k]->force.v_sig = vj_sig.f[k]; for (j = 0; j < 3; j++) pi[k]->a_hydro[j] -= pia[j].f[k]; } -- GitLab