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

Don't update the signal velocity of particle j in non-sym routines.

parent a0bc969c
No related branches found
No related tags found
No related merge requests found
...@@ -750,7 +750,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_nonsym_force( ...@@ -750,7 +750,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_nonsym_force(
/* Update the signal velocity. */ /* Update the signal velocity. */
pi->force.v_sig = fmaxf(pi->force.v_sig, v_sig); 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( ...@@ -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]->force.u_dt += piu_dt.f[k];
pi[k]->h_dt -= pih_dt.f[k]; pi[k]->h_dt -= pih_dt.f[k];
pi[k]->force.v_sig = vi_sig.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]; for (j = 0; j < 3; j++) pi[k]->a_hydro[j] -= pia[j].f[k];
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment