Commit 625aefaa authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Give a 10% tolerance to the check of valid valocities when drifting multipoles.

parent 552a9c16
......@@ -240,17 +240,17 @@ INLINE static void gravity_drift(struct gravity_tensors *m, double dt) {
m->CoM[2] += dz;
#ifdef SWIFT_DEBUG_CHECKS
if (m->m_pole.vel[0] > m->m_pole.max_delta_vel[0])
if (m->m_pole.vel[0] > m->m_pole.max_delta_vel[0] * 1.1)
error("Invalid maximal velocity");
if (m->m_pole.vel[0] < m->m_pole.min_delta_vel[0])
if (m->m_pole.vel[0] < m->m_pole.min_delta_vel[0] * 1.1)
error("Invalid minimal velocity");
if (m->m_pole.vel[1] > m->m_pole.max_delta_vel[1])
if (m->m_pole.vel[1] > m->m_pole.max_delta_vel[1] * 1.1)
error("Invalid maximal velocity");
if (m->m_pole.vel[1] < m->m_pole.min_delta_vel[1])
if (m->m_pole.vel[1] < m->m_pole.min_delta_vel[1] * 1.1)
error("Invalid minimal velocity");
if (m->m_pole.vel[2] > m->m_pole.max_delta_vel[2])
if (m->m_pole.vel[2] > m->m_pole.max_delta_vel[2] * 1.1)
error("Invalid maximal velocity");
if (m->m_pole.vel[2] < m->m_pole.min_delta_vel[2])
if (m->m_pole.vel[2] < m->m_pole.min_delta_vel[2] * 1.1)
error("Invalid minimal velocity");
#endif
......
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