Use numerical limits to query precision of types

Signed-off-by: henrykotze <henry@flycloudline.com>

Author: henrykotze

src/systems/lighter_than_air_dynamics/LighterThanAirDynamics.cc

@@ -45,6 +45,7 @@
 #include "gz/transport/Node.hh"
 
 #include <mutex>
+#include <limits>
 
 #include "LighterThanAirDynamics.hh"
 
@@ -418,7 +419,7 @@ void LighterThanAirDynamics::PreUpdate(
   double visc_normal_y_ = 0.0;
   double visc_normal_z_ = 0.0;
 
-  if(visc_normal_mag_ > __DBL_EPSILON__){
+  if(visc_normal_mag_ > std::numeric_limits<double>::epsilon()){
 
     visc_normal_y_ = vel_eps_v[1]/visc_normal_mag_;
     visc_normal_z_ = vel_eps_v[2]/visc_normal_mag_;