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

Updated tolerance requirements for the multipole construction accuracy test.

parent d6a3001f
No related branches found
No related tags found
1 merge request!328Gravity multi dt
......@@ -1157,7 +1157,7 @@ void cell_check_multipole(struct cell *c, void *data) {
/* Now compare the multipole expansion */
if (!gravity_multipole_equal(&ma, c->multipole, tolerance)) {
message("Multipoles are not equal at depth=%d!", c->depth);
message("Multipoles are not equal at depth=%d! tol=%f", c->depth, tolerance);
message("Correct answer:");
gravity_multipole_print(&ma.m_pole);
message("Recursive multipole:");
......
......@@ -572,6 +572,7 @@ INLINE static int gravity_multipole_equal(const struct gravity_tensors *ga,
return 0;
}
#endif
tolerance *= 10.;
#if SELF_GRAVITY_MULTIPOLE_ORDER > 2
/* Manhattan Norm of 3rd order terms */
const float order3_norm =
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment