diff --git a/src/feedback/EAGLE/feedback.c b/src/feedback/EAGLE/feedback.c index 8d81eae0ac8932fcae39c8456710992048f73a80..76370156e4cda9c7abab8f538bd9d2517be1ebf8 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 e6a3a9ec9340fd72eb707fab016c692256c234f2..63d5332f4d7eed815d7e32119d3ce6db67eb3af6 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;