[PYTHON MATRIX-SIG] Uniform Random Number Generator extension class

[P. Dubois]

On the matrix extension site I have posted URNG.tar. It has been tested with
the 0.11 version of the matrix extension.

This class supplies two modules:

URNGmodule.c
   Uniform random number generator objects can be created with either
a seed initialized from the clock or a seed supplied by the creating call.
The random number generator objects contain methods to return either a single
random number or a vector of random numbers (i.e., a Matrix of shape (n,)).
The individual generators are completely independent and do not affect the
stream produced by other generators. The generator is a Cray ranf() 
compatible generator.

Ranf.py
   This is the interface intended for users. Here is the entire text:

import URNG
import Matrix
# x=CreateGenerator(seed) creates an random number generator stream
#   seed < 0  ==>  Use the default initial seed value.
#   seed = 0  ==>  Set a "random" value for the seed from the system clock.
#   seed > 0  ==>  Set seed directly (32 bits only).
#   x.ranf() samples from that stream.
#   x.sample(n) returns a vector from that stream.
#
# ranf() returns a stream of random numbers
# random_sample(n) returns a vector of length n filled with random numbers
# random_matrix(n,m,...) returns an arbitrarily shaped matrix filled with
#                        random numbers.

CreateGenerator = URNG.CreateGenerator
standard_generator = CreateGenerator(-1)

def ranf():
	"Returns a single random number from the standard generator."
	return standard_generator.ranf()

def random_sample(n):
	"""Returns a vector of length n of random numbers 
	from the standard generator."""
	return  standard_generator.sample(n)

def random_matrix(*n):
	"""Returns a matrix with shape given by n of random numbers from 
	the standard generator."""
	if not n:
		raise 'random_matrix cannot be an empty shape'
        m = 1
	for i in n:
		m = m * i
	x = Matrix.reshape(random_sample(m),n)
        return x







=================
MATRIX-SIG  - SIG on Matrix Math for Python

send messages to: matrix-sig@python.org
administrivia to: matrix-sig-request@python.org
=================