Skip to content
Snippets Groups Projects
Commit efda478c authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Corrected typo in 4th order derivatives.

parent d53edf05
No related branches found
No related tags found
1 merge request!420Gravity speedup
Hide whitespace changes
Inline Side-by-side
src/gravity_derivatives.h
+ 3
− 3
View file @ efda478c
...@@ -179,9 +179,9 @@ __attribute__((always_inline)) INLINE static void compute_potential_derivatives( ...@@ -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_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_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_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_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_inv5; 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_inv5; pot->D_112 = r_z2 * r_x * r_z * r_inv9 + r_x * r_y * r_inv7;
#endif #endif
#if SELF_GRAVITY_MULTIPOLE_ORDER > 5 #if SELF_GRAVITY_MULTIPOLE_ORDER > 5
#error "Missing implementation for order >5" #error "Missing implementation for order >5"
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment