From c106b20715cd64b9f1bd2eef6ef16efe9b50c1b2 Mon Sep 17 00:00:00 2001
From: James Willis <james.s.willis@durham.ac.uk>
Date: Fri, 26 May 2017 15:29:50 +0100
Subject: [PATCH] Add h_max when computing max_d instead of hi.
---
src/runner_doiact_vec.c | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/src/runner_doiact_vec.c b/src/runner_doiact_vec.c
index f50c46baaa..5e259044f0 100644
--- a/src/runner_doiact_vec.c
+++ b/src/runner_doiact_vec.c
@@ -305,7 +305,8 @@ __attribute__((always_inline)) INLINE static void populate_max_d_no_cache(
h = p->h;
d = sort_i[k].d + dx_max;
- max_di[k] = d + h * kernel_gamma - rshift;
+ //max_di[k] = d + h * kernel_gamma - rshift;
+ max_di[k] = d + hi_max;
/* If the particle is out of range set the index to
* the last active particle within range. */
@@ -331,7 +332,8 @@ __attribute__((always_inline)) INLINE static void populate_max_d_no_cache(
d = sort_j[k].d - dx_max;
/*TODO: don't think rshift should be taken off here, waiting on Pedro. */
- max_dj[k] = d - h * kernel_gamma - rshift;
+ //max_dj[k] = d - h * kernel_gamma - rshift;
+ max_dj[k] = d - hj_max;
/* If the particle is out of range set the index to
* the last active particle within range. */
--
GitLab