RandomLib
1.10

RandomLib is a C++ class which implements the Mersenne Twister random number generator, MT19937 and the SIMDoriented Fast Mersenne Twister random number generator, SFMT19937. For a description of MT19937 see
Makoto Matsumoto and Takuji Nishimura,
Mersenne Twister: A 623Dimensionally Equidistributed Uniform PseudoRandom Number Generator,
ACM TOMACS 8, 3–30 (1998).
For a description of SFMT19937 see
Mutsuo Saito,
An Application of Finite Field: Design and Implementation of 128bit InstructionBased Fast Pseudorandom Number Generator,
Master's Thesis, Dept. of Math., Hiroshima University (Feb. 2007),
Mutsuo Saito and Makoto Matsumoto,
SIMDoriented Fast Mersenne Twister: a 128bit Pseudorandom Number Generator,
accepted in the proceedings of MCQMC2006.
MT19937 and SFMT19937 are highquality random number generators with an exceptionally large period 2^{19937} − 1 or about 10^{6001}; it passes all current tests for randomness.
The emphasis in this implementation is on providing a reliable source of random numbers for scientific applications where there's a premium on
Random provides the basic functionality of the MT19937 and SFMT19937 random generator, converting the random data into various formats. It provides methods for returning random integers of various sizes (short int, int, long int, etc.), for returning random integers in the semiclosed interval [0,n) and the closed interval [m,n].
You can obtain uniform random reals at various precisions; these are defined by rounding a random number uniformly sampled in (0,1) and exactly rounding it (down, up, or nearest) to a subset of representable reals. Thus Float() is the result of rounding a random number in (0,1) down to the nearest representable double.
Boolean() returns true with probability 1/2. Prob(x) returns true with probability x. Prob(a, b) returns true with probability a/b.
Bits<n>() returns n bits of randomness in a bitset<n>.
In addition, Random provides facilities for setting seeds, for selecting a "random" seed, for saving and restoring its state, and for jumping the generator forwards and backwards.
NormalDistribution and ExponentialDistribution are classes which sample from the normal and exponential distributions. RandomSelect selects from an arbitrary discrete distribution.
Finally, RandomNumber provides support for infinite precision random numbers. These are used by ExactExponential and ExactNormal which sample exactly from the exponential and normal distributions. DiscreteNormal and DiscreteNormalAlt sample exactly from the discrete normal distribution (important in some applications in cryptography).
Both 32bit and 64bit versions of MT19937 and SFMT19937 are provided with the 32bit version of SFMT19937 being the default generator. See Selection of default generator for a comparison between the various generators and how to change the default generator.
RandomLib depends on no external libraries (other than the standard C++ library). However, boost serialization can optionally be used for saving and restoring the state of the generator. The example code for parallelization, RandomThread.cpp, uses OpenMP.
My interest in random number generators extends back through much of my professional career in plasma physics, chaos theory, and computational chemistry. I wrote a random number library for Fortran 77 and Fortran 90 which implemented one of Knuth's recommended random number generators (see http://w3.pppl.gov/ntcc/RNG/). With the current C++ random number library, I switched to a more robust underlying generator, SFMT19937, provided more flexible seeding options, and provided exact implementations for uniform real and integer distributions.
Undoubtedly, bugs lurk in this code and in the documentation. Please report any you find to charl. es@k arney .com