From 047ce8432c172f7d0862323f93e9dcc19e2f9e06 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Thu, 12 Apr 2018 21:05:02 +0200
Subject: [PATCH] Missing factor of \pi in the argument of the Green function
for the top-level mesh potential calculation.
---
src/runner_doiact_fft.c | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/runner_doiact_fft.c b/src/runner_doiact_fft.c
index 1bc3846531..9bb5f799de 100644
--- a/src/runner_doiact_fft.c
+++ b/src/runner_doiact_fft.c
@@ -222,7 +222,8 @@ void runner_do_grav_fft(struct runner* r, int timer) {
/* Some common factors */
const double green_fac = -1. / (M_PI * box_size);
- const double a_smooth2 = 4. * M_PI * a_smooth * a_smooth / ((double)(N * N));
+ const double a_smooth2 =
+ 4. * M_PI * M_PI * a_smooth * a_smooth / ((double)(N * N));
const double k_fac = M_PI / (double)N;
/* Now de-convolve the CIC kernel and apply the Green function */
--
GitLab