From 18b258ef2aeea50b95dbb5f6b6b1dd98784e9b86 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Thu, 23 Apr 2020 12:12:20 +0200
Subject: [PATCH] Improve comments in the MAC functions
---
src/multipole_accept.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/multipole_accept.h b/src/multipole_accept.h
index 272b470bd4..233633f342 100644
--- a/src/multipole_accept.h
+++ b/src/multipole_accept.h
@@ -93,14 +93,14 @@ __attribute__((nonnull, pure)) INLINE static int gravity_M2L_accept(
const int cond_2 = max_softening * max_softening < r2;
/* Condition 3: The contribution is accurate enough
- * (E_BA / r^(p+2) < eps a_min) */
+ * (E_BA / r^(p+2) < eps * a_min) */
const int cond_3 = E_BA_term < eps * min_a_grav * r_to_p_plus2;
return cond_1 && cond_2 && cond_3;
} else {
- /* Condition 1: We are in the converging part of the Taylor expansion */
+ /* Condition 1: We are obeying the purely geometric criterion */
const int cond_1 = rho_sum * rho_sum < theta_crit2 * r2;
/* Condition 2: We are not below softening */
@@ -169,7 +169,7 @@ __attribute__((nonnull, pure)) INLINE static int gravity_M2P_accept(
} else {
- /* Condition 1: We are in the converging part of the Taylor expansion */
+ /* Condition 1: We are obeying the purely geometric criterion */
const int cond_1 = (B->r_max) * (B->r_max) < theta_crit2 * r2;
/* Condition 2: We are not below softening */
--
GitLab