From 95fedf825a3c2020c5fb94f965e5038f7006802f Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Mon, 3 Apr 2017 05:22:05 +0100
Subject: [PATCH] Updated tolerance requirements for the multipole construction
accuracy test.
---
src/cell.c | 2 +-
src/multipole.h | 1 +
2 files changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/cell.c b/src/cell.c
index 6a10eaf8cd..deadc4d0c7 100644
--- a/src/cell.c
+++ b/src/cell.c
@@ -1157,7 +1157,7 @@ void cell_check_multipole(struct cell *c, void *data) {
/* Now compare the multipole expansion */
if (!gravity_multipole_equal(&ma, c->multipole, tolerance)) {
- message("Multipoles are not equal at depth=%d!", c->depth);
+ message("Multipoles are not equal at depth=%d! tol=%f", c->depth, tolerance);
message("Correct answer:");
gravity_multipole_print(&ma.m_pole);
message("Recursive multipole:");
diff --git a/src/multipole.h b/src/multipole.h
index 14e76dbd68..4268ff6ccc 100644
--- a/src/multipole.h
+++ b/src/multipole.h
@@ -572,6 +572,7 @@ INLINE static int gravity_multipole_equal(const struct gravity_tensors *ga,
return 0;
}
#endif
+ tolerance *= 10.;
#if SELF_GRAVITY_MULTIPOLE_ORDER > 2
/* Manhattan Norm of 3rd order terms */
const float order3_norm =
--
GitLab