From 5e6634afbefd63f79124b7fbec46bf38e5f501c4 Mon Sep 17 00:00:00 2001 From: Matthieu Schaller <matthieu.schaller@durham.ac.uk> Date: Mon, 24 Jul 2017 11:02:41 +0100 Subject: [PATCH] Re-instated long-range truncation for periodic gravity. --- src/kernel_long_gravity.h | 13 ++++++------- 1 file changed, 6 insertions(+), 7 deletions(-) diff --git a/src/kernel_long_gravity.h b/src/kernel_long_gravity.h index 7b1c598464..6952681999 100644 --- a/src/kernel_long_gravity.h +++ b/src/kernel_long_gravity.h @@ -37,15 +37,14 @@ __attribute__((always_inline)) INLINE static void kernel_long_grav_eval( float u, float *const W) { - /* const float arg1 = u * 0.5f; */ - /* const float arg2 = u * one_over_sqrt_pi; */ - /* const float arg3 = -arg1 * arg1; */ + const float arg1 = u * 0.5f; + const float arg2 = u * one_over_sqrt_pi; + const float arg3 = -arg1 * arg1; - /* const float term1 = erfcf(arg1); */ - /* const float term2 = arg2 * expf(arg3); */ + const float term1 = erfcf(arg1); + const float term2 = arg2 * expf(arg3); - /* *W = term1 + term2; */ - *W = 1.f; + *W = term1 + term2; } #endif // SWIFT_KERNEL_LONG_GRAVITY_H -- GitLab