madness::Random Class Reference

A random number generator (portable, vectorized, and thread-safe) More...

#include <ran.h>

Inheritance diagram for madness::Random:

Collaboration diagram for madness::Random:

Public Member Functions
	Random (unsigned int seed=5461)
virtual	~Random ()
double	get ()
void	getbytes (int n, unsigned char *MADNESS_RESTRICT v)
	Returns vector of random bytes in [0,256) More...
RandomState	getstate () const
	Returns full state of the generator. More...
template<typename T >
void	getv (int n, T *MADNESS_RESTRICT v)
	Returns a vector of uniform doubles in [0,1) More...
void	setstate (const RandomState &s)
	Restores state of the generator. More...
void	setstate (unsigned int seed)
	Sets state of the generator from integer. More...

Static Public Member Functions
static void	test ()
	Test the generator. More...

Private Member Functions
void	generate ()
unsigned int	simple ()
Private Member Functions inherited from madness::Mutex
	Mutex (int junk=0)
	Make and initialize a mutex ... initial state is unlocked. More...
virtual	~Mutex ()
void	lock () const
	Acquire the mutex waiting if necessary. More...
pthread_mutex_t *	ptr () const
	Return a pointer to the pthread mutex for use by a condition variable. More...
bool	try_lock () const
	Try to acquire the mutex ... return true on success, false on failure. More...
void	unlock () const
	Free a mutex owned by this thread. More...

Private Attributes
const double	beta
int	cur
const int	r
const int	s
unsigned int	simple_state
double *const	u

Detailed Description

A random number generator (portable, vectorized, and thread-safe)

Following Brent 1992, we use a 48-bit generalized Fibonacci generator

u[n] = alpha*u[n-r] + beta*u[n-s] mod m

with alpha=1, beta=7, r=1279, s=861, m=2^48. Double precision numbers are used to perform exact integer arithmetic. 48-bit because we have 52 bits of mantissa, alpha+1 is 3 bits and 1 bit spare.

The period is nominally 2^m (2^r - 1) / 2 but if p is the period, X[n] and X[n+p/2k] differ by at most k bits (0 < k < 48) so usage should be limited to the first 2^r-1 entries (about 10^385 values).

Each instance provides a separate stream, but it is up to the user to partition the sequence by selecting distinct seeds or other means.

The streams are thread safe.

A default stream is provided as madness::default_random_generator.

Constructor & Destructor Documentation

Random()

madness::Random::Random

(

unsigned int

seed = 5461

)

References setstate().

~Random()

madness::Random::~Random ( )

virtual

References u.

Member Function Documentation

generate()

void madness::Random::generate ( )

private

References beta, cur, r, s, and u.

Referenced by get(), getbytes(), getv(), and setstate().

get()

double madness::Random::get ( )

inline

References cur, generate(), r, and u.

Referenced by madness::RandomValue< double >(), and madness::RandomValue< float >().

getbytes()

void madness::Random::getbytes	(	int	n,
		unsigned char *MADNESS_RESTRICT	v
	)

Returns vector of random bytes in [0,256)

References cur, generate(), r, u, and v.

getstate()

RandomState madness::Random::getstate

(

)

const

Returns full state of the generator.

References cur, r, s, and u.

getv()

template<typename T >

void madness::Random::getv ( int  n, T *MADNESS_RESTRICT  v  )

inline

Returns a vector of uniform doubles in [0,1)

References cur, generate(), r, T(), u, and v.

Referenced by madness::RandomVector< double >(), madness::RandomVector< double_complex >(), madness::RandomVector< float >(), and madness::RandomVector< float_complex >().

setstate() [1/2]

void madness::Random::setstate

(

const RandomState &

s

)

Restores state of the generator.

References cur, r, s, and u.

Referenced by Random(), madness::Localizer::localize_boys(), madness::Localizer::localize_new(), main(), test(), test_Kcommutator(), test_lowrank_function(), and test_math().

setstate() [2/2]

void madness::Random::setstate

(

unsigned int

seed

)

Sets state of the generator from integer.

References generate(), lo, r, simple(), simple_state, and u.

simple()

unsigned int madness::Random::simple ( )

private

References simple_state.

Referenced by setstate().

test()

void madness::Random::test ( )

static

Test the generator.

References b, madness::f, and r.

Member Data Documentation

beta

const double madness::Random::beta

private

Referenced by generate().

cur

int madness::Random::cur

private

Referenced by generate(), get(), getbytes(), getstate(), getv(), and setstate().

r

const int madness::Random::r

private

Referenced by generate(), get(), getbytes(), getstate(), getv(), setstate(), and test().

s

const int madness::Random::s

private

Referenced by generate(), getstate(), and setstate().

simple_state

unsigned int madness::Random::simple_state

private

Referenced by setstate(), and simple().

u

double* const madness::Random::u

private

Referenced by ~Random(), generate(), get(), getbytes(), getstate(), getv(), and setstate().

---

The documentation for this class was generated from the following files:

• ran.h
• ran.cc