utils_distributions.hpp

Go to the documentation of this file.

#pragma once
 
#include "distribution.hpp"
#include <random>
 
namespace netdem {
 
class RandomEngine {
public:
  RandomEngine(const RandomEngine &) = delete;
 
  RandomEngine &operator=(const RandomEngine &) = delete;
 
  static RandomEngine &GetInstance() {
    static RandomEngine instance;
    return instance;
  }
 
  double GetUniformDistribution(double bound_min, double bound_max) {
    std::random_device rd;
    std::mt19937 mt_eng = std::mt19937(rd());
    std::uniform_real_distribution<double> real_dist(bound_min, bound_max);
 
    return real_dist(mt_eng);
  }
 
  VecXT<double> GetUniformDistribution(double bound_min, double bound_max,
                                       int num) {
    std::random_device rd;
    std::mt19937 mt_eng = std::mt19937(rd());
    std::uniform_real_distribution<double> real_dist(bound_min, bound_max);
 
    VecXT<double> num_list(num, 0);
 
    for (int i = 0; i < num; i++) {
      num_list[i] = real_dist(mt_eng);
    }
 
    return num_list;
  }
 
private:
  RandomEngine() {}
};
 
} // namespace netdem