rand_gauss

This routine returns a number out of a pseudo-random normal deviate (ie, its probability distribution is a Gaussian of mean 0 and variance 1).

Location in HEALPix directory tree: src/f90/mod/rngmod.f90 

FORMAT

var=rand_gauss( rng_handle )

ARGUMENTS

name & dimensionality	kind	in/out	description
rng_handle	planck_rng	INOUT	structure of type planck_rng containing on all information necessary to continue same random sequence.
var	DP	OUT	number belonging to a pseudo-random normal deviate.

EXAMPLE:

use healpix_types
use rngmod
type(planck_rng) :: rng_handle
real(dp) :: gauss

call rand_init(rng_handle, 12345, 6789012)
gauss = rand_gauss(rng_handle)

initiates a random sequence with the pair of seeds (12345, 6789012), and generates one number out of the normal deviate.

RELATED ROUTINES

This section lists the routines related to rand_gauss

planck_rng

derived type describing RNG state

rand_uni

function which returns a random uniform deviate.

rand_init

subroutine to initiate a random number sequence.