From 8d9b898eceb1470bd03b5e209b70c64f2025daaf Mon Sep 17 00:00:00 2001
From: Pedro Gonnet <gonnet@google.com>
Date: Thu, 6 Jul 2017 22:35:53 +0200
Subject: [PATCH] be safe and check the error relative to the stored old
dx_max_sort.
---
src/runner_doiact.h | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/runner_doiact.h b/src/runner_doiact.h
index 865d2e78cc..2ef1ccbad6 100644
--- a/src/runner_doiact.h
+++ b/src/runner_doiact.h
@@ -918,7 +918,7 @@ void DOPAIR1(struct runner *r, struct cell *ci, struct cell *cj, const int sid,
p->x[1] * runner_shift[sid][1] +
p->x[2] * runner_shift[sid][2];
if (fabsf(d - sort_i[pid].d) - ci->dx_max_sort >
- 1.0e-6 * max(fabsf(d), ci->dx_max_sort))
+ 1.0e-6 * max(fabsf(d), ci->dx_max_sort_old))
error(
"particle shift diff exceeds dx_max_sort in cell ci. ci->nodeID=%d "
"cj->nodeID=%d d=%e sort_i[pid].d=%e ci->dx_max_sort=%e "
@@ -932,7 +932,7 @@ void DOPAIR1(struct runner *r, struct cell *ci, struct cell *cj, const int sid,
p->x[1] * runner_shift[sid][1] +
p->x[2] * runner_shift[sid][2];
if (fabsf(d - sort_j[pjd].d) - cj->dx_max_sort >
- 1.0e-6 * max(fabsf(d), cj->dx_max_sort))
+ 1.0e-6 * max(fabsf(d), cj->dx_max_sort_old))
error(
"particle shift diff exceeds dx_max_sort in cell cj. cj->nodeID=%d "
"ci->nodeID=%d d=%e sort_j[pjd].d=%e cj->dx_max_sort=%e "
--
GitLab