From f61ae18527265b36391dd03ca8af857b2c11f86d Mon Sep 17 00:00:00 2001
From: Bert Vandenbroucke <bert.vandenbroucke@ugent.be>
Date: Wed, 8 Mar 2017 10:17:53 +0000
Subject: [PATCH] Added a debug check for negative mass/energy to Gizmo
kick_extra. Fixed a typo in space.c.
---
src/hydro/Gizmo/hydro.h | 12 ++++++++++++
src/space.c | 2 +-
2 files changed, 13 insertions(+), 1 deletion(-)
diff --git a/src/hydro/Gizmo/hydro.h b/src/hydro/Gizmo/hydro.h
index ba8a676ee4..46d532a900 100644
--- a/src/hydro/Gizmo/hydro.h
+++ b/src/hydro/Gizmo/hydro.h
@@ -484,6 +484,18 @@ __attribute__((always_inline)) INLINE static void hydro_kick_extra(
p->conserved.energy += p->conserved.flux.energy;
#endif
+#ifdef SWIFT_DEBUG_CHECKS
+ if (p->conserved.mass < 0.) {
+ error(
+ "Negative mass after conserved variables update (mass: %g, dmass: %g)!",
+ p->conserved.mass, p->conserved.flux.mass);
+ }
+
+ if (p->conserved.energy < 0.) {
+ error("Negative energy after conserved variables update!");
+ }
+#endif
+
/* Add gravity. We only do this if we have gravity activated. */
if (p->gpart) {
/* Retrieve the current value of the gravitational acceleration from the
diff --git a/src/space.c b/src/space.c
index 625fe944c4..a64d68cc61 100644
--- a/src/space.c
+++ b/src/space.c
@@ -523,7 +523,7 @@ void space_rebuild(struct space *s, int verbose) {
/* Be verbose about this. */
#ifdef SWIFT_DEBUG_CHECKS
- if (s->e->nodeID == 0 || verbose) message("re)building space");
+ if (s->e->nodeID == 0 || verbose) message("(re)building space");
fflush(stdout);
#endif
--
GitLab