runif.h

Go to the documentation of this file.

 
#ifndef RUNIF_H
#define RUNIF_H
 
#include "OACommonDefines.h"
#define SEEDOK 1
 
#define SEEDBAD 0
 
#define SEED_VECTOR_LENGTH (97+1)
 
namespace oacpp 
{
    struct SeedSet
    {
        SeedSet() {};
        SeedSet(int one, int two, int three, int four) { is = one; js = two; ks = three; ls = four; };
        int is;
        int js;
        int ks;
        int ls;
    };
 
    class RUnif
    {
    public:
        RUnif();
        RUnif(int is, int js, int ks, int ls);
        
        ~RUnif() {};
        
        void seed(int is, int js, int ks, int ls );
        
        void seed(SeedSet & seedSet);
        
        SeedSet getSeedSet();
 
        void runif(std::vector<double> & x, int n);
 
        static int mod(int a, int b);
 
    private:
        static int seedok(int is, int js, int ks, int ls );
        
        void ranums(std::vector<double> & x, int n);
 
        int m_jent, m_i, m_j, m_k, m_l, ip, jp;
        std::vector<double> u;
        double c, cd, cm;
    };
}
 
#endif