ZIGGURAT_INLINE
Ziggurat Method for Exponential or Normal Random Numbers
ZIGGURAT_INLINE
is a C library which
rapidly generates random variates from the uniform, normal or
exponential distributions.
The uniform numbers are generated directly. The ziggurat method is used
to compute the normal and exponential values.
In this version, the underlying generators are implemented "inline",
invoking a function call only in exceptional cases. This results in
very fast execution.
Licensing:
The computer code and data files made available on this web page
are distributed under
the GNU LGPL license.
Languages:
ZIGGURAT_INLINE is available in
a C version and
a C++ version.
Related Data and Programs:
NORMAL,
a C library which
computes elements of a sequence of pseudorandom normally distributed values.
UNIFORM,
a C library which
computes elements of uniform pseudorandom sequence.
ZIGGURAT,
a C program which
generates points from a uniform, normal or exponential distribution, using
the ziggurat method.
ZIGGURAT_OPENMP,
a C program which
demonstrates how the ZIGGURAT library can be used to generate random numbers
in an OpenMP parallel program.
Reference:

Philip Leong, Guanglie Zhang, DongU Lee, Wayne Luk, John Villasenor,
A comment on the implementation of the ziggurat method,
Journal of Statistical Software,
Volume 12, Number 7, February 2005.

George Marsaglia, Wai Wan Tsang,
The Ziggurat Method for Generating Random Variables,
Journal of Statistical Software,
Volume 5, Number 8, October 2000, seven pages.
Source Code:
Examples and Tests:
List of Routines:

CONG_SEEDED evaluates the CONG generator for uint32_t's.

CONG_VALUE evaluates the CONG generator for uint32_t's.

CPU_TIME reports the total CPU time for a program.

EFIX generates variates when rejection occurs in the exponential code.

KISS_SEEDED evaluates the KISS random number generator.

KISS_VALUE evaluates the KISS generator for uint32_t's.

MWC_SEEDED evaluates the MWC multiplywithcarry random number generator.

MWC_VALUE evaluates the MWC generator for uint32_t's.

NFIX generates variates when rejection occurs in the normal code.

R4_EXP_SETUP sets data needed by R4_EXP.

R4_EXP_VALUE returns an exponentially distributed float.

R4_NOR_SETUP sets data needed by R4_NOR.

R4_NOR_VALUE returns a normally distributed float.

R4_UNI_VALUE returns a uniformly distributed float.

SHR3_SEEDED evaluates the SHR3 generator for uint32_t's.

SHR3_VALUE evaluates the SHR3 generator for uint32_t's.

TIMESTAMP prints the current YMDHMS date as a time stamp.

ZIGGET gets the seeds.

ZIGSET sets the seeds and creates the tables for the Ziggurat method.
You can go up one level to
the C source codes.
Last revised on 16 October 2013.