From 070f08a866da7319b6ddc03b1022d413340f0796 Mon Sep 17 00:00:00 2001
From: loikki <loic.hausammann@protonmail.ch>
Date: Thu, 17 Oct 2019 16:54:02 +0200
Subject: [PATCH] Fix stars interaction
---
src/runner_doiact_functions_stars.h | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/src/runner_doiact_functions_stars.h b/src/runner_doiact_functions_stars.h
index d452fba01b..34b1280420 100644
--- a/src/runner_doiact_functions_stars.h
+++ b/src/runner_doiact_functions_stars.h
@@ -298,9 +298,8 @@ void DO_SYM_PAIR1_STARS(struct runner *r, struct cell *ci, struct cell *cj,
struct xpart *restrict xparts_j = cj->hydro.xparts;
#endif
const double dj_min = sort_j[0].d;
- const float dx_max_rshift =
- (ci->stars.dx_max_sort + cj->hydro.dx_max_sort) - rshift;
const float dx_max = (ci->stars.dx_max_sort + cj->hydro.dx_max_sort);
+ const float dx_max_rshift = dx_max - rshift;
/* Loop over the sparts in ci. */
for (int pid = count_i - 1;
@@ -428,9 +427,8 @@ void DO_SYM_PAIR1_STARS(struct runner *r, struct cell *ci, struct cell *cj,
struct xpart *restrict xparts_i = ci->hydro.xparts;
struct spart *restrict sparts_j = cj->stars.parts;
const double di_max = sort_i[count_i - 1].d - rshift;
- const float dx_max_rshift =
- (ci->hydro.dx_max_sort + cj->stars.dx_max_sort) + rshift;
const float dx_max = (ci->hydro.dx_max_sort + cj->stars.dx_max_sort);
+ const float dx_max_rshift = dx_max - rshift;
/* Loop over the parts in cj. */
for (int pjd = 0; pjd < count_j && sort_j[pjd].d - hj_max - dx_max < di_max;
--
GitLab