RNG.h

Go to the documentation of this file.

#ifndef RNG_H
#define RNG_H

//Used by the randomise function
#include <time.h>
#include <cmath>
#include "ExtendedMath.h"

class RNG
{
public:

        RNG(){random_seed_LCG=0; a=1103515245; c=12345;m=1024*1024*1024;type=1;p=607;q=273;random_seed_LFG.resize(p); seed_LFG();}

        // (note the call to seed_LFG is only required really if using that type of generator, but the other one is always required
        void set_RandomSeed(Mdouble new_seed){random_seed_LCG=new_seed; seed_LFG();}

        Mdouble get_RN(Mdouble min, Mdouble max);
        
        Mdouble test();
        
        void set_LCGParms(int new_a,int new_c,int new_m)
        {a=new_a;c=new_c;m=new_m;}
        
        void randomise() {set_RandomSeed(time(NULL));}
        
        void set_LFGParms(int new_p, int new_q);
        
        void set_RNtypeLCG(){type=0;}
        
        void set_RNtypeLFG(){type=1;}
        
private:

        Mdouble get_LCG(Mdouble min, Mdouble max);

        Mdouble get_LFG(Mdouble min, Mdouble max);
        
        void seed_LFG();

        unsigned long int random_seed_LCG;
        
        vector<Mdouble> random_seed_LFG;
        
        long int a,c,m;
        
        long int p,q;
        
        int type;
        
};


#endif