diff --git a/src/feedback/EAGLE/feedback.c b/src/feedback/EAGLE/feedback.c index d041969e808ac95a4f7f4de06a4d2cb8ccaf41d7..18b7d7ceb955c3bfbb0a4df356856c60d94a0d46 100644 --- a/src/feedback/EAGLE/feedback.c +++ b/src/feedback/EAGLE/feedback.c @@ -719,7 +719,7 @@ void compute_stellar_evolution(const struct feedback_props* feedback_props, const float ngb_gas_mass = sp->feedback_data.to_collect.ngb_mass; /* Check if there are neighbours, otherwise exit */ - if (ngb_gas_mass == 0.f) { + if (ngb_gas_mass == 0.f || sp->density.wcount * pow_dimension(sp->h) < 1e-4) { feedback_reset_feedback(sp, feedback_props); return; }