diff --git a/src/gravity_derivatives.h b/src/gravity_derivatives.h index fe42a4e0a808984f3dad163857abb280103f8e13..d317f33484dd73318a028b083318ebf495f8fd9b 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"