If h is h_max then set rho_dh to 0, prevent insanely high grad_h_term
Compare changes
+ 7
− 9
@@ -518,17 +518,15 @@ __attribute__((always_inline)) INLINE static void hydro_prepare_force(
@@ -518,17 +518,15 @@ __attribute__((always_inline)) INLINE static void hydro_prepare_force(