specialize raw float/raw double

Author: jube

include/gf2/core/Random.h

@@ -73,19 +73,25 @@ namespace gf {
     template<typename T>
     std::enable_if_t<std::is_floating_point_v<T>, T> compute_uniform_float()
     {
-      return static_cast<T>(compute_raw_float());
+      if constexpr (std::is_same_v<T, double>) {
+        return compute_raw_double();
+      } else if constexpr (std::is_same_v<T, float>) {
+        return compute_raw_float();
+      } else {
+        return static_cast<T>(compute_raw_float());
+      }
     }
 
     template<typename T>
     std::enable_if_t<std::is_floating_point_v<T>, T> compute_uniform_float(T max)
     {
-      return static_cast<T>(compute_raw_float() * max);
+      return compute_uniform_float<T>() * max;
     }
 
     template<typename T>
     std::enable_if_t<std::is_floating_point_v<T>, T> compute_uniform_float(T min, T max)
     {
-      return min + static_cast<T>(compute_raw_float() * (max - min));
+      return min + compute_uniform_float<T>() * (max - min);
     }
 
     template<typename T>
@@ -119,7 +125,8 @@ namespace gf {
     }
 
     uint64_t compute_raw_integer(uint64_t max);
-    double compute_raw_float();
+    double compute_raw_double();
+    float compute_raw_float();
 
   private:
     RandomEngine m_engine;

library/core/Random.cc

@@ -220,10 +220,16 @@ namespace gf {
     return r;
   }
 
-  double Random::compute_raw_float()
+  double Random::compute_raw_double()
   {
     const uint64_t value = m_engine();
     return double(value >> 11) * 0x1.0p-53;
   }
 
+  float Random::compute_raw_float()
+  {
+    const uint64_t value = m_engine();
+    return float(value >> 40) * 0x1.0p-24f;
+  }
+
 } // namespace gf