mersenne_twister.h

/* mersenne_twister.h
 *
 * This is a version of the Mersenne Twister (aka mt19937) that is designed for
 * parallel use. Usage:
 *
 * struct mt19937p myMT19937;
 * sgenrand(seed,&myMT19937);
 * randomNumber = genrand(&myMT19937);
 *
 * Simply create different mt19937p structures for your different streams.
 * Disclaimer & Warning: This has only been lightly tested and has shown that it
 * generates identical sequences as mt19937.c! Use at your own risk!
 *
 * Massaged into parallel for my Austen McDonald <austen@cc.gatech.edu>
 * http://www.prism.gatech.edu/~gte363v/
 * November, 2000
 */

#ifndef MERSENNE_TWISTER_H
#define MERSENNE_TWISTER_H

#define MATRIX_A 0x9908b0df /* constant vector a */
#define N 624

struct mt19937p {
	unsigned long mt[N];
	int mti;
	unsigned long mag01[2];
};

/* initializing the array with a NONZERO seed */
void sgenrand(unsigned long seed, struct mt19937p* config);

/* double */ /* generating reals */
 unsigned long  /* for integer generation */
genrand(struct mt19937p* config);
#endif