From d847056cb286ceca4ae646f7f23b8d79696b398a Mon Sep 17 00:00:00 2001
From: Bert Vandenbroucke <bert.vandenbroucke@gmail.com>
Date: Tue, 10 Apr 2018 14:19:20 +0100
Subject: [PATCH] Added debug check on particle drift distance.
---
src/cell.c | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/src/cell.c b/src/cell.c
index ea036c9708..a1f57edcca 100644
--- a/src/cell.c
+++ b/src/cell.c
@@ -2405,6 +2405,15 @@ void cell_drift_part(struct cell *c, const struct engine *e, int force) {
drift_part(p, xp, dt_drift, dt_kick_hydro, dt_kick_grav, dt_therm,
ti_old_part, ti_current);
+#ifdef SWIFT_DEBUG_CHECKS
+ /* Make sure the particle does not drift by more than a box length. */
+ if (fabsf(xp->v_full[0] * dt_drift) > e->s->dim[0] ||
+ fabsf(xp->v_full[1] * dt_drift) > e->s->dim[1] ||
+ fabsf(xp->v_full[2] * dt_drift) > e->s->dim[2]) {
+ error("Particle drifts by more than a box length!");
+ }
+#endif
+
/* Limit h to within the allowed range */
p->h = min(p->h, hydro_h_max);
--
GitLab