rand-clcg4.h

Go to the documentation of this file.

#ifndef INC_clcg4_h
#define INC_clcg4_h
 
#include <stdio.h>
#include "ross-random-internal.h"
 
typedef int32_t * tw_seed;
 
struct tw_rng
{
        /*
         * equals a[i]^{m[i]-2} mod m[i]
         */
        long long       b[4];
 
        /*
         * a[j]^{2^w} et a[j]^{2^{v+w}}.
         */
        int32_t m[4];
        int32_t a[4];
        int32_t aw[4];
        int32_t avw[4];
 
        // the seed..
        int32_t seed[4];
};
 
enum SeedType
{
        InitialSeed, LastSeed, NewSeed
};
 
typedef enum SeedType SeedType;
 
struct tw_rng_stream
{
    unsigned long count;
        int32_t  Ig[4];
        int32_t  Lg[4];
        int32_t  Cg[4];
 
        //tw_rng        *rng;
 
#ifdef RAND_NORMAL
        double   tw_normal_u1;
        double   tw_normal_u2;
        int      tw_normal_flipflop;
#endif
};
 
extern tw_rng     *rng_init(int v, int w);
extern tw_rng   *rng_core_init(int v, int w);
extern void     rng_set_initial_seed();
extern void     rng_init_generator(tw_rng_stream * g, SeedType Where, tw_rng * the_rng);
extern void     rng_set_seed(tw_rng_stream * g, uint32_t * s, tw_rng * the_rng);
extern void     rng_get_state(tw_rng_stream * g, uint32_t * s);
extern void     rng_write_state(tw_rng_stream * g, FILE *f);
extern double   rng_gen_val(tw_rng_stream * g);
extern double   rng_gen_reverse_val(tw_rng_stream * g);
 
#endif