diff --git a/src/hydro/Planetary/hydro.h b/src/hydro/Planetary/hydro.h index 887693accfb28ae16fe8a8f044c1d72e850e80cf..957e96dcf391b9027016926a969b28366590664f 100644 --- a/src/hydro/Planetary/hydro.h +++ b/src/hydro/Planetary/hydro.h @@ -523,7 +523,7 @@ __attribute__((always_inline)) INLINE static void hydro_prepare_force( if (p->h == hydro_props->h_max) { rho_dh = 0.f; } - const float grad_h_term = + const float grad_h_term = 1.f / (1.f + hydro_dimension_inv * p->h * rho_dh * rho_inv); /* Compute the Balsara switch */