mt19937ar.c File Reference

random numbers (mersenne twister). More...

#include <stdio.h>
#include "p2.h"

Include dependency graph for mt19937ar.c:

Go to the source code of this file.

Macros
#define	N   624
#define	M   397
#define	MATRIX_A   0x9908b0dfUL /* constant vector a */
#define	UMASK   0x80000000UL /* most significant w-r bits */
#define	LMASK   0x7fffffffUL /* least significant r bits */
#define	MIXBITS(u, v)   ( ((u) & UMASK) | ((v) & LMASK) )
#define	TWIST(u, v)   ((MIXBITS(u,v) >> 1) ^ ((v)&1UL ? MATRIX_A : 0UL))

Functions
void	init_genrand (unsigned long s)
void	init_by_array (unsigned long init_key[], int key_length)
static void	next_state (void)
unsigned long	potion_rand_int (void)
	generates a random number on [0,0xffffffff]-interval More...
double	potion_rand_double (void)
	generates a random number on [0,1) with 53-bit resolution More...

Variables
static unsigned long	state [N]
static int	left = 1
static int	initf = 0
static unsigned long *	next

Detailed Description

random numbers (mersenne twister).

As Lobby (global long) or PNInteger (0-1.0)

A C-program for MT19937, with initialization improved 2002/2/10. Coded by Takuji Nishimura and Makoto Matsumoto. This is a faster version by taking Shawn Cokus's optimization, Matthe Bellew's simplification, Isaku Wada's real version.

Before using, initialize the state by using init_genrand(seed) or init_by_array(init_key, key_length).

Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, All rights reserved.

Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:

1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
3. The names of its contributors may not be used to endorse or promote products derived from this software without specific prior written permission.

THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

Any feedback is very welcome. http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space)

Definition in file mt19937ar.c.

Macro Definition Documentation

#define N   624

Definition at line 52 of file mt19937ar.c.

#define M   397

Definition at line 53 of file mt19937ar.c.

#define MATRIX_A   0x9908b0dfUL /* constant vector a */

Definition at line 54 of file mt19937ar.c.

#define UMASK   0x80000000UL /* most significant w-r bits */

Definition at line 55 of file mt19937ar.c.

#define LMASK   0x7fffffffUL /* least significant r bits */

Definition at line 56 of file mt19937ar.c.

#define MIXBITS	(		u,
			v
	)		( ((u) & UMASK) | ((v) & LMASK) )

Definition at line 57 of file mt19937ar.c.

#define TWIST	(		u,
			v
	)		((MIXBITS(u,v) >> 1) ^ ((v)&1UL ? MATRIX_A : 0UL))

Definition at line 58 of file mt19937ar.c.

Function Documentation

void init_genrand

(

unsigned long

s

)

Definition at line 66 of file mt19937ar.c.

void init_by_array	(	unsigned long	init_key[],
		int	key_length
	)

Definition at line 84 of file mt19937ar.c.

static void next_state ( void  )

static

Definition at line 109 of file mt19937ar.c.

unsigned long potion_rand_int

(

void

)

generates a random number on [0,0xffffffff]-interval

Definition at line 138 of file mt19937ar.c.

double potion_rand_double

(

void

)

generates a random number on [0,1) with 53-bit resolution

Definition at line 154 of file mt19937ar.c.

Variable Documentation

unsigned long state[N]

static

Definition at line 60 of file mt19937ar.c.

int left = 1

static

Definition at line 61 of file mt19937ar.c.

int initf = 0

static

Definition at line 62 of file mt19937ar.c.

unsigned long* next

static

Definition at line 63 of file mt19937ar.c.