From a3e6692ffc582ec4f022694d8a3009ee6c79cfeb Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Wed, 24 Apr 2019 13:08:23 +0200
Subject: [PATCH] Added more checks that the quantities used in the feedback
make sense.
---
src/feedback/EAGLE/feedback.c | 4 ++++
src/feedback/EAGLE/feedback_iact.h | 5 +++++
2 files changed, 9 insertions(+)
diff --git a/src/feedback/EAGLE/feedback.c b/src/feedback/EAGLE/feedback.c
index 8d81eae0ac..76370156e4 100644
--- a/src/feedback/EAGLE/feedback.c
+++ b/src/feedback/EAGLE/feedback.c
@@ -690,6 +690,10 @@ void compute_stellar_evolution(const struct feedback_props* feedback_props,
TIMER_TIC;
+#ifdef SWIFT_DEBUG_CHECKS
+ if (age < 0.f) error("Negative age for a star.");
+#endif
+
/* Allocate temporary array for calculating imf weights */
float stellar_yields[eagle_feedback_N_imf_bins];
diff --git a/src/feedback/EAGLE/feedback_iact.h b/src/feedback/EAGLE/feedback_iact.h
index e6a3a9ec93..63d5332f4d 100644
--- a/src/feedback/EAGLE/feedback_iact.h
+++ b/src/feedback/EAGLE/feedback_iact.h
@@ -114,6 +114,11 @@ runner_iact_nonsym_feedback_apply(const float r2, const float *dx,
Omega_frac = 0.f;
}
+#ifdef SWIFT_DEBUG_CHECKS
+ if (Omega_frac < 0. || Omega_frac > 1.)
+ error("Invalid fraction of material to dsitribute.");
+#endif
+
/* Update particle mass */
const double current_mass = hydro_get_mass(pj);
const double delta_mass = si->feedback_data.to_distribute.mass * Omega_frac;
--
GitLab