RandomNumberProvider.hpp

Go to the documentation of this file.

#ifndef SPATIALSEIR_INCLUDEFILES
#define SPATIALSEIR_INCLUDEFILES

#include<math.h>
#include<cstring>
#include<vector>
#endif

#ifndef RANDOM_NUMBER_PROVIDER_INC
#define RANDOM_NUMBER_PROVIDER_INC


#ifdef LSS_USE_BOOST

#include<boost/random/mersenne_twister.hpp>
#include<boost/random/uniform_real.hpp>
#include<boost/random/uniform_int_distribution.hpp>
#include<boost/random/gamma_distribution.hpp>
#include<boost/random/poisson_distribution.hpp>
#include<boost/random/normal_distribution.hpp>
#include<boost/random/binomial_distribution.hpp>
#include <boost/math/distributions/normal.hpp>



namespace SpatialSEIR
{
    class RandomNumberProvider
    {
        boost::random::mt19937 generator;
        boost::random::uniform_real_distribution<> unidist;
        boost::random::uniform_int_distribution<> unidist_int;
        boost::random::gamma_distribution<> gammadist;
        public:
            //Random Methods
            RandomNumberProvider(unsigned int seed);
            double uniform();   
            int uniform_int();
            int uniform_int(int a, int b);
            int poisson(int mu);
            double normal(double mu, double sd);
            double* uniform(int n);
            double* uniform(int n, double* output);
            double gamma();
            double gamma(double a);
            double* gamma(int n);
            double* gamma(int n, double* output);
            double beta(double a, double b);
            int binom(int n, double p);


            //Density Functions
            double dpois(int x, double mu);
            double dnorm(double x, double mu, double sd);
            double dbinom(int x, int n, double p);
            double dgamma(double x, double a, double b);

            // Utility Functions
            double choose(int n,int k);
            double choosePartial(int n,int k);

            double factorial(int k);

            //Data
            double* logFactorialMemo;
            int* maxFactorial;

            ~RandomNumberProvider();
    };
}
#endif
#ifndef LSS_USE_BOOST

#include<Rcpp.h>

namespace SpatialSEIR
{
    class RandomNumberProvider
    {
        public:
            //Random Methods
            RandomNumberProvider(unsigned int seed);
            double uniform();   
            int uniform_int();
            int uniform_int(int a, int b);
            int poisson(int mu);
            double normal(double mu, double sd);
            double* uniform(int n);
            double* uniform(int n, double* output);
            double gamma();
            double gamma(double a);
            double* gamma(int n);
            double* gamma(int n, double* output);
            double beta(double a, double b);
            int binom(int n, double p);

            //Density Functions
            double dpois(int x, double mu);
            double dnorm(double x, double mu, double sd);
            double dbinom(int x, int n, double p);
            double dgamma(double x, double a, double b);

            ~RandomNumberProvider();
    };
}


#endif
#endif