Random.h

#pragma once
 
#include <cmath>
#include <cstdint>
#include <numbers>
#include <random>
 
namespace environment::field::aggregation
{
 
class Random
{
public:
   explicit Random(std::uint64_t seed)
       : m_engine(seed)
   {
   }
 
   double Uniform() { return static_cast<double>(m_engine() >> 11U) * (1.0 / 9007199254740992.0); }
 
   double Uniform(double low, double high) { return low + (high - low) * Uniform(); }
 
   double Normal()
   {
      const double u1 = 1.0 - Uniform(); // in (0, 1], so the logarithm is finite
      const double u2 = Uniform();
      return std::sqrt(-2.0 * std::log(u1)) * std::cos(2.0 * std::numbers::pi * u2);
   }
 
   double Normal(double mean, double sigma) { return mean + sigma * Normal(); }
 
   double LogNormal(double mu, double sigma) { return std::exp(Normal(mu, sigma)); }
 
private:
   std::mt19937_64 m_engine;
};
 
} // namespace environment::field::aggregation