Line data    Source code

       1              : #include "Set.H"
       2              : #include <random>
       3              : namespace Set
       4              : {
       5              : }
       6              : 
       7              : namespace Util
       8              : {
       9       313960 : Set::Scalar Random()
      10              : {
      11       313960 :     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            0 : }
      21              : 
      22              : }