From c108da640f84c1e49c4dbd173b1c60d18a85dc6f Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Fri, 17 Mar 2017 10:11:41 +0800
Subject: [PATCH] Correct multipole order guard
---
src/multipole.h | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/multipole.h b/src/multipole.h
index a34af3bbd8..4c90b5cf51 100644
--- a/src/multipole.h
+++ b/src/multipole.h
@@ -606,7 +606,7 @@ INLINE static int gravity_multipole_equal(const struct gravity_tensors *ga,
return 0;
}
#endif
-#if SELF_GRAVITY_MULTIPOLE_ORDER > 2
+#if SELF_GRAVITY_MULTIPOLE_ORDER > 3
/* Check 4th order terms */
if (fabsf(ma->M_400 + mb->M_400) > 1e-6 * ma->M_000 &&
fabsf(ma->M_400 - mb->M_400) / fabsf(ma->M_400 + mb->M_400) > tolerance) {
--
GitLab