From 6e6839c130dec7b20021e8160e919cd0b1c31646 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Thu, 17 May 2018 10:58:16 +0100
Subject: [PATCH] Correct expression for the Ewald correction of the potential.
---
src/gravity.c | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/gravity.c b/src/gravity.c
index b21c104d06..5d572f98df 100644
--- a/src/gravity.c
+++ b/src/gravity.c
@@ -74,6 +74,7 @@ void gravity_exact_force_ewald_init(double boxSize) {
#ifdef SWIFT_GRAVITY_FORCE_CHECKS
+ const float boxSize_inv = 1.f / boxSize;
const float boxSize_inv2 = 1.f / (boxSize * boxSize);
int use_file = 0;
@@ -302,7 +303,7 @@ void gravity_exact_force_ewald_init(double boxSize) {
fewald_x[i][j][k] *= boxSize_inv2;
fewald_y[i][j][k] *= boxSize_inv2;
fewald_z[i][j][k] *= boxSize_inv2;
- potewald[i][j][k] *= boxSize_inv2;
+ potewald[i][j][k] *= boxSize_inv;
}
}
}
--
GitLab