#include <string.h>
#include <assert.h>
#include <stdio.h>
#include <stdint.h>
Go to the source code of this file.
Data Structures | |
struct | MTGP64_PARAMS_FAST_T |
MTGP64 parameters. More... | |
struct | MTGP64_STATUS_FAST_T |
MTGP64 internal state array. More... | |
struct | MTGP64_FAST_T |
Typedefs | |
typedef struct MTGP64_PARAMS_FAST_T | mtgp64_params_fast_t |
typedef struct MTGP64_STATUS_FAST_T | mtgp64_status_fast_t |
typedef struct MTGP64_FAST_T | mtgp64_fast_t |
Functions | |
int | mtgp64_init (mtgp64_fast_t *mtgp64, const mtgp64_params_fast_t *para, uint64_t seed) |
This function allocates and initializes the internal state array with a 64-bit integer seed. | |
void | mtgp64_init_state (uint64_t array[], const mtgp64_params_fast_t *para, uint64_t seed) |
This function initializes the internal state array with a 64-bit integer seed. | |
int | mtgp64_init_by_array (mtgp64_fast_t *mtgp64, const mtgp64_params_fast_t *para, uint64_t *array, int length) |
This function allocates and initializes the internal state array with a 64-bit integer array. | |
int | mtgp64_init_by_str (mtgp64_fast_t *mtgp64, const mtgp64_params_fast_t *para, char *str) |
This function allocates and initializes the internal state array with a character array. | |
void | mtgp64_free (mtgp64_fast_t *mtgp64) |
This releases the memory allocated by mtgp64_init(), mtgp64_init_by_array(), mtgp64_init_by_str(). | |
void | mtgp64_print_idstring (const mtgp64_fast_t *mtgp64, FILE *fp) |
This function prints the Mersenne exponent and SHA1 of characteristic polynomial of generators state transition function. | |
static void | mtgp64_do_recursion (uint64_t *r, uint64_t x1, uint64_t x2, uint64_t y, int sh1, int sh2, uint64_t mask, uint64_t tbl[16]) |
This is a recursion formula of the generator. | |
static void | mtgp64_next_state (mtgp64_fast_t *mtgp64) |
The state transition function. | |
static uint64_t | mtgp64_genrand_uint64 (mtgp64_fast_t *mtgp64) |
This function generates and returns 64-bit unsigned integer. | |
static double | mtgp64_genrand_close1_open2 (mtgp64_fast_t *mtgp64) |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range [1, 2). | |
static double | mtgp64_genrand_close_open (mtgp64_fast_t *mtgp64) |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range [0, 1). | |
static double | mtgp64_genrand_open_close (mtgp64_fast_t *mtgp64) |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range (0, 1]. | |
static double | mtgp64_genrand_open_open (mtgp64_fast_t *mtgp64) |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range (0, 1). | |
static uint64_t | mtgp64_temper (const uint64_t tmp_tbl[16], uint64_t r, uint64_t t) |
The tempering function. | |
static double | mtgp64_temper_double (const uint64_t dbl_tmp_tbl[16], uint64_t r, uint64_t t) |
The tempering and converting function. | |
static double | mtgp64_temper_double_open (const uint64_t dbl_tmp_tbl[16], uint64_t r, uint64_t t) |
The tempering and converting function for generating floating point number f (0 < f < 1). | |
Variables | |
mtgp64_params_fast_t | mtgp64_params_fast_23209 [128] |
parameter constants tables for MEXP=23209. | |
mtgp64_params_fast_t | mtgp64_params_fast_44497 [128] |
parameter constants tables for MEXP=44497. | |
mtgp64_params_fast_t | mtgp64_params_fast_110503 [128] |
parameter constants tables for MEXP=110503. |
Copyright (C) 2009 Mutsuo Saito, Makoto Matsumoto and Hiroshima University. All rights reserved.
The new BSD License is applied to this software, see LICENSE.txt
typedef struct MTGP64_FAST_T mtgp64_fast_t |
typedef struct MTGP64_PARAMS_FAST_T mtgp64_params_fast_t |
typedef struct MTGP64_STATUS_FAST_T mtgp64_status_fast_t |
static void mtgp64_do_recursion | ( | uint64_t * | r, | |
uint64_t | x1, | |||
uint64_t | x2, | |||
uint64_t | y, | |||
int | sh1, | |||
int | sh2, | |||
uint64_t | mask, | |||
uint64_t | tbl[16] | |||
) | [inline, static] |
This is a recursion formula of the generator.
MTGP64 is a 64-bit generator, but using 32-bit operations to fit to graphic processors.
[out] | r | output |
[in] | x1 | the farthest part of state array. |
[in] | x2 | the second farthest part of state array. |
[in] | y | a part of state array. |
[in] | sh1 | the shift parameter 1. |
[in] | sh2 | the shift parameter 2. |
[in] | mask | the bit mask parameter. |
[in] | tbl | the matrix parameter. |
Referenced by mtgp64_next_state().
void mtgp64_free | ( | mtgp64_fast_t * | mtgp64 | ) |
This releases the memory allocated by mtgp64_init(), mtgp64_init_by_array(), mtgp64_init_by_str().
[in,out] | mtgp64 | MTGP all in one structure. |
References MTGP64_FAST_T::status.
static double mtgp64_genrand_close1_open2 | ( | mtgp64_fast_t * | mtgp64 | ) | [inline, static] |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range [1, 2).
mtgp64_init(), mtgp64_init_by_array() or mtgp64_init_by_str() must be called before this function.
[in,out] | mtgp64 | MTGP all in one structure. |
References MTGP64_STATUS_FAST_T::array, MTGP64_PARAMS_FAST_T::dbl_tmp_tbl, MTGP64_STATUS_FAST_T::idx, MTGP64_STATUS_FAST_T::large_mask, MTGP64_STATUS_FAST_T::large_size, mtgp64_next_state(), mtgp64_temper_double(), MTGP64_FAST_T::params, MTGP64_PARAMS_FAST_T::pos, MTGP64_STATUS_FAST_T::size, and MTGP64_FAST_T::status.
Referenced by mtgp64_genrand_close_open(), and mtgp64_genrand_open_close().
static double mtgp64_genrand_close_open | ( | mtgp64_fast_t * | mtgp64 | ) | [inline, static] |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range [0, 1).
mtgp64_init(), mtgp64_init_by_array() or mtgp64_init_by_str() must be called before this function.
[in,out] | mtgp64 | MTGP all in one structure. |
References mtgp64_genrand_close1_open2().
static double mtgp64_genrand_open_close | ( | mtgp64_fast_t * | mtgp64 | ) | [inline, static] |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range (0, 1].
mtgp64_init(), mtgp64_init_by_array() or mtgp64_init_by_str() must be called before this function.
[in,out] | mtgp64 | MTGP all in one structure. |
References mtgp64_genrand_close1_open2().
static double mtgp64_genrand_open_open | ( | mtgp64_fast_t * | mtgp64 | ) | [inline, static] |
This function generates and returns double precision pseudorandom number which distributes uniformly in the range (0, 1).
mtgp64_init(), mtgp64_init_by_array() or mtgp64_init_by_str() must be called before this function.
[in,out] | mtgp64 | MTGP all in one structure. |
References MTGP64_STATUS_FAST_T::array, MTGP64_PARAMS_FAST_T::dbl_tmp_tbl, MTGP64_STATUS_FAST_T::idx, MTGP64_STATUS_FAST_T::large_mask, MTGP64_STATUS_FAST_T::large_size, mtgp64_next_state(), mtgp64_temper_double_open(), MTGP64_FAST_T::params, MTGP64_PARAMS_FAST_T::pos, MTGP64_STATUS_FAST_T::size, and MTGP64_FAST_T::status.
static uint64_t mtgp64_genrand_uint64 | ( | mtgp64_fast_t * | mtgp64 | ) | [inline, static] |
This function generates and returns 64-bit unsigned integer.
mtgp64_init(), mtgp64_init_by_array() or mtgp64_init_by_str() must be called before this function.
[in,out] | mtgp64 | MTGP all in one structure. |
References MTGP64_STATUS_FAST_T::array, MTGP64_STATUS_FAST_T::idx, MTGP64_STATUS_FAST_T::large_mask, MTGP64_STATUS_FAST_T::large_size, mtgp64_next_state(), mtgp64_temper(), MTGP64_FAST_T::params, MTGP64_PARAMS_FAST_T::pos, MTGP64_STATUS_FAST_T::size, MTGP64_FAST_T::status, and MTGP64_PARAMS_FAST_T::tmp_tbl.
int mtgp64_init | ( | mtgp64_fast_t * | mtgp64, | |
const mtgp64_params_fast_t * | para, | |||
uint64_t | seed | |||
) |
This function allocates and initializes the internal state array with a 64-bit integer seed.
The allocated memory should be freed by calling mtgp64_free(). para should be one of the elements in the parameter table (mtgp64-param-ref.c).
[out] | mtgp64 | MTGP structure. |
[in] | para | parameter structure |
[in] | seed | a 64-bit integer used as the seed. |
References alloc_state(), MTGP64_STATUS_FAST_T::array, MTGP64_STATUS_FAST_T::idx, mtgp64_init_state(), MTGP64_STATUS_FAST_T::size, and MTGP64_FAST_T::status.
int mtgp64_init_by_array | ( | mtgp64_fast_t * | mtgp64, | |
const mtgp64_params_fast_t * | para, | |||
uint64_t * | array, | |||
int | length | |||
) |
This function allocates and initializes the internal state array with a 64-bit integer array.
The allocated memory should be freed by calling mtgp64_free(). para should be one of the elements in the parameter table (mtgp64-param-ref.c).
[out] | mtgp64 | MTGP structure. |
[in] | para | parameter structure |
[in] | array | a 64-bit integer array used as a seed. |
[in] | length | length of the array. |
References alloc_state(), MTGP64_STATUS_FAST_T::array, MTGP64_STATUS_FAST_T::idx, ini_func1(), ini_func2(), MTGP64_PARAMS_FAST_T::mexp, non_zero, MTGP64_FAST_T::params, MTGP64_STATUS_FAST_T::size, MTGP64_FAST_T::status, and MTGP64_PARAMS_FAST_T::tbl.
int mtgp64_init_by_str | ( | mtgp64_fast_t * | mtgp64, | |
const mtgp64_params_fast_t * | para, | |||
char * | array | |||
) |
This function allocates and initializes the internal state array with a character array.
The allocated memory should be freed by calling mtgp64_free(). para should be one of the elements in the parameter table (mtgp64-param-ref.c). This is the same algorithm with mtgp64_init_by_array(), but hope to be more useful.
[out] | mtgp64 | MTGP structure. |
[in] | para | parameter structure |
[in] | array | a character array used as a seed. (terminated by zero.) |
References alloc_state(), MTGP64_STATUS_FAST_T::array, MTGP64_STATUS_FAST_T::idx, ini_func1(), ini_func2(), MTGP64_PARAMS_FAST_T::mexp, non_zero, MTGP64_FAST_T::params, MTGP64_STATUS_FAST_T::size, MTGP64_FAST_T::status, and MTGP64_PARAMS_FAST_T::tbl.
void mtgp64_init_state | ( | uint64_t | array[], | |
const mtgp64_params_fast_t * | para, | |||
uint64_t | seed | |||
) |
This function initializes the internal state array with a 64-bit integer seed.
The allocated memory should be freed by calling mtgp64_free(). para should be one of the elements in the parameter table (mtgp64-param-ref.c).
This function is call by cuda program, because cuda program uses another structure and another allocation method.
[out] | array | MTGP internal status vector. |
[in] | para | parameter structure |
[in] | seed | a 64-bit integer used as the seed. |
References MTGP64_PARAMS_FAST_T::mexp, and MTGP64_PARAMS_FAST_T::tbl.
Referenced by make_kernel_data(), and mtgp64_init().
static void mtgp64_next_state | ( | mtgp64_fast_t * | mtgp64 | ) | [inline, static] |
The state transition function.
[in,out] | mtgp64 | the all in one structure |
References MTGP64_STATUS_FAST_T::array, MTGP64_STATUS_FAST_T::idx, MTGP64_STATUS_FAST_T::large_mask, MTGP64_STATUS_FAST_T::large_size, MTGP64_PARAMS_FAST_T::mask, mtgp64_do_recursion(), MTGP64_FAST_T::params, MTGP64_PARAMS_FAST_T::pos, MTGP64_PARAMS_FAST_T::sh1, MTGP64_PARAMS_FAST_T::sh2, MTGP64_STATUS_FAST_T::size, MTGP64_FAST_T::status, and MTGP64_PARAMS_FAST_T::tbl.
Referenced by mtgp64_genrand_close1_open2(), mtgp64_genrand_open_open(), and mtgp64_genrand_uint64().
void mtgp64_print_idstring | ( | const mtgp64_fast_t * | mtgp64, | |
FILE * | fp | |||
) |
This function prints the Mersenne exponent and SHA1 of characteristic polynomial of generators state transition function.
[in] | mtgp64 | MTGP all in one structure. |
[in,out] | fp | FILE pointer. |
References MTGP64_PARAMS_FAST_T::mexp, MTGP64_FAST_T::params, and MTGP64_PARAMS_FAST_T::poly_sha1.
static uint64_t mtgp64_temper | ( | const uint64_t | tmp_tbl[16], | |
uint64_t | r, | |||
uint64_t | t | |||
) | [inline, static] |
The tempering function.
[in] | tmp_tbl | the pre-computed tempering table. |
[in] | r | the value to be tempered. |
[in] | t | the tempering helper value. |
Referenced by mtgp64_genrand_uint64().
static double mtgp64_temper_double | ( | const uint64_t | dbl_tmp_tbl[16], | |
uint64_t | r, | |||
uint64_t | t | |||
) | [inline, static] |
The tempering and converting function.
[in] | dbl_tmp_tbl | the pre-computed tempering table. |
[in] | r | the value to be tempered. |
[in] | t | the tempering helper value. |
Referenced by mtgp64_genrand_close1_open2().
static double mtgp64_temper_double_open | ( | const uint64_t | dbl_tmp_tbl[16], | |
uint64_t | r, | |||
uint64_t | t | |||
) | [inline, static] |
The tempering and converting function for generating floating point number f (0 < f < 1).
[in] | dbl_tmp_tbl | the pre-computed tempering table. |
[in] | r | the value to be tempered. |
[in] | t | the tempering helper value. |
Referenced by mtgp64_genrand_open_open().
parameter constants tables for MEXP=110503.
size - pos > 1024.
parameter constants tables for MEXP=23209.
size - pos > 256.