RANLIB
General Random Number Generators (RNG's)
RANLIB
is a C library which
produces random samples from Probability Density Functions (PDF's),
including Beta, Chisquare Exponential, F, Gamma, Multivariate normal,
Noncentral chisquare, Noncentral F, Univariate normal, random permutations,
Real uniform, Binomial, Negative Binomial, Multinomial, Poisson
and Integer uniform,
by Barry Brown and James Lovato.
Note that this C version of RANLIB was NOT created by simply
running the original FORTRAN77 source code through the f2c program!
If you have ever looked at the resultant source code from such a process,
you know it was never meant for human eyes.
RANLIB relies on streams of uniform random numbers generated
by a lower level package called RNGLIB. A copy of RNGLIB
must be available in order for RANLIB to executed.
The RNGLIB routines provide 32 virtual random number generators.
Each generator can provide 1,048,576 blocks of numbers, and each block
is of length 1,073,741,824. Any generator can be set to the beginning
or end of the current block or to its starting value. Packaging is
provided so that if these capabilities are not needed, a single
generator with period 2.3 X 10^18 is seen.
The routines, and the probability density functions they sample, include:

GENBET, Beta distribution;

GENCHI, ChiSquare distribution;

GENEXP, Exponential distribution;

GENF, F distribution;

GENGAM, Gamma distribution;

GENMN, multivariate normal distribution;

GENMUL, multinomial distribution;

GENNCH, noncentral ChiSquare distribution;

GENNF, noncentral F distribution;

GENNOR, normal distribution;

GENUNF, uniform distribution on [0,1];

IGNBIN, binomial distribution;

IGNLGI, uniform distribution on integers between 1 and 2147483562;

IGNNBN, negative binomial distribution.

IGNPOI, Poisson distribution.

IGNUIN, uniform distribution on integers in a given range.
Licensing:
The computer code and data files made available on this web page
are distributed under
the GNU LGPL license.
Languages:
RANLIB is available in
a C version and
a C++ version and
a FORTRAN90 version and
a MATLAB version and
a Python version.
Related Data and Programs:
ASA183,
a C library which
implements a random number generator (RNG),
by Wichman and Hill.
This is a C version of Applied Statistics Algorithm 183.
LOG_NORMAL_TRUNCATED_AB,
a C library which
returns quantities associated with the log normal Probability
Distribution Function (PDF) truncated to the interval [A,B].
NORMAL,
a C library which
computes elements of a
sequence of pseudorandom normally distributed values.
RAND48,
a C program which
demonstrates the use of the rand48 family of random number generators
(RNG's) available in the C/C++ standard library.
RANDLC,
a C library which
implements a random number generator (RNG)
used by the NAS Benchmark programs.
RNGLIB,
a C library which
implements a random number generator (RNG) with splitting facilities,
allowing multiple independent streams to be computed,
by L'Ecuyer and Cote.
UNIFORM,
a C library which
computes elements of a pseudorandom uniform sequence.
VAN_DER_CORPUT,
a C library which
computes elements of a 1D van der Corput Quasi Monte Carlo (QMC) sequence
using a simple interface.
WALKER_SAMPLE,
a C library which
efficiently samples a discrete probability vector using
Walker sampling.
Author:
Original FORTRAN77 version by Barry Brown, James Lovato.
C version by John Burkardt.
Reference:

Joachim Ahrens, Ulrich Dieter,
Computer Methods for Sampling From the
Exponential and Normal Distributions,
Communications of the ACM,
Volume 15, Number 10, October 1972, pages 873882.

Joachim Ahrens, Ulrich Dieter,
Generating Gamma Variates by a Modified Rejection Technique,
Communications of the ACM,
Volume 25, Number 1, January 1982, pages 4754.

Joachim Ahrens, Ulrich Dieter,
Computer Generation of Poisson Deviates
From Modified Normal Distributions,
ACM Transactions on Mathematical Software,
Volume 8, Number 2, June 1982, pages 163179.

Joachim Ahrens, Ulrich Dieter,
Computer Methods for Sampling from Gamma, Beta, Poisson and
Binomial Distributions,
Computing,
Volume 12, Number 3, September 1974, pages 223246.

Joachim Ahrens, Ulrich Dieter,
Extensions of Forsythe's Method for Random
Sampling from the Normal Distribution,
Mathematics of Computation,
Volume 27, Number 124, October 1973, page 927937.

Russell Cheng,
Generating Beta Variates with Nonintegral Shape Parameters,
Communications of the ACM,
Volume 21, Number 4, April 1978, pages 317322.

Luc Devroye,
NonUniform Random Variate Generation,
Springer, 1986,
ISBN: 0387963057,
LC: QA274.D48.

Voratas Kachitvichyanukul, Bruce Schmeiser,
Binomial Random Variate Generation,
Communications of the ACM,
Volume 31, Number 2, February 1988, page 216222.

Pierre LEcuyer, Serge Cote,
Implementing a Random Number Package with Splitting Facilities,
ACM Transactions on Mathematical Software,
Volume 17, Number 1, March 1991, pages 98111.
Source Code:
Examples and Tests:
List of Routines:

CH_CAP capitalizes a single character.

COVAR computes the covariance of two vectors.

GENBET generates a beta random deviate.

GENCHI generates a ChiSquare random deviate.

GENEXP generates an exponential random deviate.

GENF generates an F random deviate.

GENGAM generates a Gamma random deviate.

GENMN generates a multivariate normal deviate.

GENMUL generates a multinomial random deviate.

GENNCH generates a noncentral ChiSquare random deviate.

GENNF generates a noncentral F random deviate.

GENNOR generates a normal random deviate.

GENPRM generates and applies a random permutation to an array.

GENUNF generates a uniform random deviate.

I4_MAX returns the maximum of two I4's.

I4_MIN returns the smaller of two I4's.

IGNBIN generates a binomial random deviate.

IGNNBN generates a negative binomial random deviate.

IGNPOI generates a Poisson random deviate.

IGNUIN generates a random integer in a given range.

LENNOB counts the length of a string, ignoring trailing blanks.

PHRTST converts a phrase to a pair of random number generator seeds.

PRCOMP prints covariance information.

R4_MAX returns the maximum of two R4's.

R4_MIN returns the minimum of two R4's.

S_EQI reports whether two strings are equal, ignoring case.

SDOT forms the dot product of two vectors.

SETCOV sets a covariance matrix from variance and common correlation.

SETGMN sets data for the generation of multivariate normal deviates.

SEXPO evaluates the standard exponential distribution.

SGAMMA returns a deviate from the standard Gamma distribution.

SNORM returns a deviate from the standard normal distribution.

SPOFA factors a real symmetric positive definite matrix.

STATS computes statistics for a given array.

TRSTAT returns the mean and variance for distributions.
You can go up one level to
the C source codes.
Last revised on 02 April 2013.