Line data    Source code

       1             : #include "Set.H"
       2             : #include <random>
       3             : namespace Set
       4             : {
       5             : }
       6             : 
       7             : namespace Util
       8             : {
       9        3814 : Set::Scalar Random()
      10             : {
      11        3814 :     return ((Set::Scalar) rand()) / ((Set::Scalar) RAND_MAX);
      12             : }
      13           0 : Set::Scalar Gaussian(amrex::Real mean,amrex::Real std_deviation)
      14             : {
      15           0 :     std::random_device randomness_device{};
      16           0 :     std::mt19937 pseudorandom_generator{randomness_device()};
      17           0 :     std::normal_distribution<double> distribution{mean, std_deviation};
      18           0 :     auto sample = distribution(pseudorandom_generator);
      19           0 :     return sample;
      20             : }
      21             : 
      22             : }