RandomGenerator.hh

[詳解]

 
//
// Copyright (C) 2011-2020 Yokokura, Yuki
// This program is free software;
// you can redistribute it and/or modify it under the terms of the FreeBSD License.
// For details, see the License.txt file.
 
#ifndef RANDOMGENERATOR
#define RANDOMGENERATOR
 
#include <random>
#include "Matrix.hh"
 
namespace ARCS {    // ARCS名前空間
    class RandomGenerator {
        public:
            RandomGenerator(const double MinOrMean, const double MaxOrStdDev)
                : RandomDevice(), MersenneTwister(RandomDevice()),
                  RandomInt(MinOrMean, MaxOrStdDev),
                  RandomDouble(MinOrMean, MaxOrStdDev),
                  GaussianRandom(MinOrMean, MaxOrStdDev)
            {
                
            }
 
            RandomGenerator(RandomGenerator&& r)
                : RandomDevice(), MersenneTwister(r.MersenneTwister),
                  RandomInt(r.RandomInt), RandomDouble(r.RandomDouble), GaussianRandom(r.GaussianRandom)
            {
                
            }
 
            ~RandomGenerator(){
                
            }
            
            int GetIntegerRandom(void){
                return RandomInt(MersenneTwister);
            }
 
            double GetDoubleRandom(void){
                return RandomDouble(MersenneTwister);
            }
            
            double GetGaussianRandom(void){
                return GaussianRandom(MersenneTwister);
            }
 
            void ResetSeed(void){
                RandomDevice.entropy();
                MersenneTwister.seed(RandomDevice());
            }
 
            template <size_t N, size_t M>
            void GetRandomMatrix(Matrix<N,M>& Y){
                for(size_t n = 1; n <= N; ++n){
                    for(size_t m = 1; m <= M; ++m){
                        Y.SetElement(n, m, GetDoubleRandom());  // 行列の各要素に乱数値を書き込む
                    }
                }
            }
 
            template <size_t N, size_t M>
            void GetGaussianRandomMatrix(Matrix<N,M>& Y){
                for(size_t n = 1; n <= N; ++n){
                    for(size_t m = 1; m <= M; ++m){
                        Y.SetElement(n, m, GetGaussianRandom());// 行列の各要素に乱数値を書き込む
                    }
                }
            }
            
        private:
            RandomGenerator(const RandomGenerator&) = delete;                   
            const RandomGenerator& operator=(const RandomGenerator&) = delete;  
            std::random_device RandomDevice;                
            std::mt19937 MersenneTwister;                   
            std::uniform_int_distribution<> RandomInt;      
            std::uniform_real_distribution<> RandomDouble;  
            std::normal_distribution<> GaussianRandom;      
    };
}
 
#endif