From efda478ceb8e013c4804442944105e9c85abf83b Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Wed, 16 Aug 2017 10:00:31 +0100
Subject: [PATCH] Corrected typo in 4th order derivatives.

---
 src/gravity_derivatives.h | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/gravity_derivatives.h b/src/gravity_derivatives.h
index fe42a4e0a8..d317f33484 100644
--- a/src/gravity_derivatives.h
+++ b/src/gravity_derivatives.h
@@ -179,9 +179,9 @@ __attribute__((always_inline)) INLINE static void compute_potential_derivatives(
   pot->D_220 = r_x2 * r_y2 * r_inv9 + r_x2 * r_inv7 + r_y2 * r_inv7 + r_inv5;
   pot->D_202 = r_x2 * r_z2 * r_inv9 + r_x2 * r_inv7 + r_z2 * r_inv7 + r_inv5;
   pot->D_022 = r_y2 * r_z2 * r_inv9 + r_y2 * r_inv7 + r_z2 * r_inv7 + r_inv5;
-  pot->D_211 = r_x2 * r_y * r_z * r_inv9 + r_y * r_z * r_inv5;
-  pot->D_121 = r_y2 * r_x * r_z * r_inv9 + r_x * r_z * r_inv5;
-  pot->D_112 = r_z2 * r_x * r_z * r_inv9 + r_x * r_y * r_inv5;
+  pot->D_211 = r_x2 * r_y * r_z * r_inv9 + r_y * r_z * r_inv7;
+  pot->D_121 = r_y2 * r_x * r_z * r_inv9 + r_x * r_z * r_inv7;
+  pot->D_112 = r_z2 * r_x * r_z * r_inv9 + r_x * r_y * r_inv7;
 #endif
 #if SELF_GRAVITY_MULTIPOLE_ORDER > 5
 #error "Missing implementation for order >5"
-- 
GitLab