diff --git a/src/cooling/const_du/cooling.h b/src/cooling/const_du/cooling.h
index e7f4f27db79712e6bff54559cf6781fb3aa4a438..448af9c3765e3bb6d4cbf4cc94e245a3976d5314 100644
--- a/src/cooling/const_du/cooling.h
+++ b/src/cooling/const_du/cooling.h
@@ -103,7 +103,7 @@ __attribute__((always_inline)) INLINE static float cooling_timestep(
 
   const float cooling_rate = cooling->cooling_rate;
   const float internal_energy = hydro_get_internal_energy(p, 0);
-  return cooling->cooling_tstep_mult * internal_energy / abs(cooling_rate);
+  return cooling->cooling_tstep_mult * internal_energy / fabsf(cooling_rate);
 }
 
 /**
diff --git a/src/cooling/const_lambda/cooling.h b/src/cooling/const_lambda/cooling.h
index ffa1b9945b5311916995754d8ff385fba524225f..1caf49bff1696efd2a201b5953732c437ffd67a9 100644
--- a/src/cooling/const_lambda/cooling.h
+++ b/src/cooling/const_lambda/cooling.h
@@ -117,7 +117,7 @@ __attribute__((always_inline)) INLINE static float cooling_timestep(
   /* Get current internal energy (dt=0) */
   const float u = hydro_get_internal_energy(p, 0.f);
 
-  return u / abs(du_dt);
+  return u / fabsf(du_dt);
 }
 
 /**