1 // -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; indent-tabs-mode: nil; -*- 2 // 3 // rnorm.h: Rcpp R/C++ interface class library -- 4 // 5 // Copyright (C) 2010 - 2016 Douglas Bates, Dirk Eddelbuettel and Romain Francois 6 // 7 // This file is part of Rcpp. 8 // 9 // Rcpp is free software: you can redistribute it and/or modify it 10 // under the terms of the GNU General Public License as published by 11 // the Free Software Foundation, either version 2 of the License, or 12 // (at your option) any later version. 13 // 14 // Rcpp is distributed in the hope that it will be useful, but 15 // WITHOUT ANY WARRANTY; without even the implied warranty of 16 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 17 // GNU General Public License for more details. 18 // 19 // You should have received a copy of the GNU General Public License 20 // along with Rcpp. If not, see <http://www.gnu.org/licenses/>. 21 22 #ifndef Rcpp__stats__random_rnorm_h 23 #define Rcpp__stats__random_rnorm_h 24 25 namespace Rcpp { 26 namespace stats { 27 28 class NormGenerator : public Generator<double> { 29 public: 30 31 NormGenerator( double mean_ = 0.0 , double sd_ = 1.0 ) : mean(mean_)32 mean(mean_), sd(sd_) {} 33 operator()34 inline double operator()() const { 35 return mean + sd * ::norm_rand() ; 36 } 37 38 private: 39 double mean ; 40 double sd ; 41 }; 42 43 class NormGenerator__sd1 : public Generator<double> { 44 public: 45 mean(mean_)46 NormGenerator__sd1( double mean_ = 0.0 ) : mean(mean_) {} 47 operator()48 inline double operator()() const { 49 return mean + ::norm_rand() ; 50 } 51 52 private: 53 double mean ; 54 }; 55 56 class NormGenerator__mean0 : public Generator<double> { 57 public: 58 sd(sd_)59 NormGenerator__mean0( double sd_ = 1.0 ) : sd(sd_) {} 60 operator()61 inline double operator()() const { 62 return sd * ::norm_rand() ; 63 } 64 65 private: 66 double sd ; 67 }; 68 69 class NormGenerator__mean0__sd1 : public Generator<double> { 70 public: 71 NormGenerator__mean0__sd1()72 NormGenerator__mean0__sd1( ) {} 73 operator()74 inline double operator()() const { 75 return ::norm_rand() ; 76 } 77 78 }; 79 80 } // stats 81 } // Rcpp 82 83 #endif 84