Commit c5c2b512 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Adjust gravity derivative tolerances

parent 8e4b0d26
......@@ -931,7 +931,7 @@ int main() {
clocks_set_cpufreq(cpufreq);
/* Relative tolerance */
const double tol = 1e-4;
double tol = 1e-4;
/* Get some randomness going */
const int seed = time(NULL);
......@@ -981,6 +981,8 @@ int main() {
#endif
#if SELF_GRAVITY_MULTIPOLE_ORDER > 2
tol *= 2.;
/* 3rd order terms */
test(pot.D_300, D_300(dx, dy, dz, r_inv), tol, min, "D_300");
test(pot.D_030, D_030(dx, dy, dz, r_inv), tol, min, "D_030");
......@@ -1015,6 +1017,8 @@ int main() {
#if SELF_GRAVITY_MULTIPOLE_ORDER > 4
tol *= 2.;
/* 5th order terms */
test(pot.D_500, D_500(dx, dy, dz, r_inv), tol, min, "D_500");
test(pot.D_050, D_050(dx, dy, dz, r_inv), tol, min, "D_050");
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment