rand_state.h 文件参考

#include "../lmmp.h"
#include "../impl/longlong.h"

rand_state.h 的引用(Include)关系图:

此图展示该文件直接或间接的被哪些文件引用了:

浏览源代码.

结构体
struct	pcg64_128_state
struct	pcg64_le_seq_t
struct	xoshiro256pp_state

宏定义
#define	INLINE_   static inline
#define	PCG128_DEFAULT_MULTIPLIER_HI   0x2360ED051FC65DA4ULL
#define	PCG128_DEFAULT_MULTIPLIER_LO   0x4385DF649FCCF645ULL
#define	PCG64_LE_INCREMENT   1442695040888963407ULL
#define	PCG64_LE_MULTIPLIER   6364136223846793005ULL
#define	PRIME64_0   0x9E3779B185EBCA87ULL
#define	PRIME64_1   0xC2B2AE3D27D4EB4FULL
#define	PRIME64_2   0x165667B19E3779F9ULL
#define	PRIME64_3   0x85EBCA77C2B2AE63ULL
#define	PRIME64_4   0x27D4EB2F165667C5ULL

函数
static mp_limb_t	lmmp_pcg64_128_random (pcg64_128_state *rng)
static void	lmmp_pcg64_128_srandom (pcg64_128_state *rng, mp_limb_t seed)
static mp_limb_t	lmmp_seed_generator (mp_limb_t seed)
	种子生成器
static mp_limb_t	lmmp_xoshiro256pp_random (xoshiro256pp_state *rng)
static void	lmmp_xoshiro256pp_srandom (xoshiro256pp_state *rng, mp_limb_t seed)
static void	pcg64_128_action (mp_limb_t state[2], const mp_limb_t inc[2])
static mp_limb_t	pcg64_le_action (mp_limb_t *restrict state)
static void	pcg64_le_seq_init (pcg64_le_seq_t *rng, mp_size_t i, mp_limb_t seed)
static void	pcg64_le_seq_next (mp_ptr restrict dst, mp_size_t n, pcg64_le_seq_t *rng)
static mp_limb_t	rotl (const mp_limb_t x, int k)

---

结构体说明

pcg64_128_state

struct pcg64_128_state

在文件 rand_state.h 第 25 行定义.

pcg64_128_state 的协作图:

成员变量
mp_limb_t	inc[2]
mp_limb_t	state[2]

pcg64_le_seq_t

struct pcg64_le_seq_t

在文件 rand_state.h 第 133 行定义.

pcg64_le_seq_t 的协作图:

成员变量
mp_size_t	k
mp_limb_t *restrict	state

xoshiro256pp_state

struct xoshiro256pp_state

在文件 rand_state.h 第 30 行定义.

xoshiro256pp_state 的协作图:

成员变量
mp_limb_t	s[4]

宏定义说明

INLINE_

#define INLINE_   static inline

在文件 rand_state.h 第 38 行定义.

PCG128_DEFAULT_MULTIPLIER_HI

#define PCG128_DEFAULT_MULTIPLIER_HI   0x2360ED051FC65DA4ULL

在文件 rand_state.h 第 34 行定义.

PCG128_DEFAULT_MULTIPLIER_LO

#define PCG128_DEFAULT_MULTIPLIER_LO   0x4385DF649FCCF645ULL

在文件 rand_state.h 第 35 行定义.

PCG64_LE_INCREMENT

#define PCG64_LE_INCREMENT   1442695040888963407ULL

在文件 rand_state.h 第 131 行定义.

PCG64_LE_MULTIPLIER

#define PCG64_LE_MULTIPLIER   6364136223846793005ULL

在文件 rand_state.h 第 130 行定义.

PRIME64_0

#define PRIME64_0   0x9E3779B185EBCA87ULL

PRIME64_1

#define PRIME64_1   0xC2B2AE3D27D4EB4FULL

PRIME64_2

#define PRIME64_2   0x165667B19E3779F9ULL

PRIME64_3

#define PRIME64_3   0x85EBCA77C2B2AE63ULL

PRIME64_4

#define PRIME64_4   0x27D4EB2F165667C5ULL

函数说明

lmmp_pcg64_128_random()

static mp_limb_t lmmp_pcg64_128_random ( pcg64_128_state *  rng )

inlinestatic

在文件 rand_state.h 第 77 行定义.

                                                              {
    lmmp_param_assert(rng != NULL);
    mp_limb_t oldstate[2] = {rng->state[0], rng->state[1]};
    pcg64_128_action(rng->state, rng->inc);
 
    // XSL-RR
    mp_limb_t xsl = ((oldstate[1]) ^ oldstate[0]);
    
    mp_byte_t rot = (mp_byte_t)(oldstate[1] >> 58); 
    return (xsl >> rot) | (xsl << ((-rot) & 63));
}

引用了 pcg64_128_state::inc, lmmp_param_assert, pcg64_128_action() , 以及 pcg64_128_state::state.

被这些函数引用 lmmp_seed_random_().

函数调用图:

这是这个函数的调用关系图:

lmmp_pcg64_128_srandom()

static void lmmp_pcg64_128_srandom ( pcg64_128_state *  rng, mp_limb_t  seed  )

inlinestatic

在文件 rand_state.h 第 62 行定义.

                                                                          {
    lmmp_param_assert(rng != NULL);
 
    rng->state[0] = lmmp_seed_generator(seed);
    rng->state[1] = lmmp_seed_generator(seed << 17);
    rng->inc[0] = lmmp_seed_generator(seed << 7); 
    rng->inc[0] |= 1ull;               
    rng->inc[1] = lmmp_seed_generator(seed << 21);
 
    // warm up
    pcg64_128_action(rng->state, rng->inc);
    pcg64_128_action(rng->state, rng->inc);
}

引用了 pcg64_128_state::inc, lmmp_param_assert, lmmp_seed_generator(), pcg64_128_action() , 以及 pcg64_128_state::state.

被这些函数引用 lmmp_seed_random_().

函数调用图:

这是这个函数的调用关系图:

lmmp_seed_generator()

static mp_limb_t lmmp_seed_generator ( mp_limb_t  seed )

inlinestatic

种子生成器

参数

seed

低熵种子

返回

高熵种子

在文件 rand_state.h 第 46 行定义.

                                                      {
    mp_limb_t z = (seed += 0x9e3779b97f4a7c15);
    z = (z ^ (z >> 30)) * 0xbf58476d1ce4e5b9;
    z = (z ^ (z >> 27)) * 0x94d049bb133111eb;
    return z ^ (z >> 31);
}

被这些函数引用 lmmp_global_rng_init_(), lmmp_pcg64_128_srandom(), lmmp_random_(), lmmp_strong_rng_extern_(), lmmp_strong_rng_init_(), lmmp_xoshiro256pp_srandom(), pcg64_le_seq_init() , 以及 pcg64_le_seq_next().

这是这个函数的调用关系图:

lmmp_xoshiro256pp_random()

static mp_limb_t lmmp_xoshiro256pp_random ( xoshiro256pp_state *  rng )

inlinestatic

在文件 rand_state.h 第 94 行定义.

                                                                    {
    lmmp_param_assert(rng != NULL);
    const mp_limb_t r = rotl(rng->s[0] + rng->s[3], 23) + rng->s[0];
    const mp_limb_t t = rng->s[1] << 17;
 
    rng->s[2] ^= rng->s[0];
    rng->s[3] ^= rng->s[1];
    rng->s[1] ^= rng->s[2];
    rng->s[0] ^= rng->s[3];
    rng->s[2] ^= t;
    rng->s[3] = rotl(rng->s[3], 45);
 
    return r;
}

引用了 lmmp_param_assert, rotl() , 以及 xoshiro256pp_state::s.

被这些函数引用 lmmp_seed_random_() , 以及 lmmp_xoshiro256pp_srandom().

函数调用图:

这是这个函数的调用关系图:

lmmp_xoshiro256pp_srandom()

static void lmmp_xoshiro256pp_srandom ( xoshiro256pp_state *  rng, mp_limb_t  seed  )

inlinestatic

在文件 rand_state.h 第 109 行定义.

                                                                                {
    lmmp_param_assert(rng != NULL);
 
    rng->s[0] = lmmp_seed_generator(seed);
    rng->s[1] = lmmp_seed_generator(seed << 17);
    rng->s[2] = lmmp_seed_generator(0x9e37b97f8a5a7c19ULL ^ seed);
    rng->s[3] = lmmp_seed_generator(0x8a07e6c7f6b9c92eULL ^ seed);
 
    // assert that the initial state is non-zero
    if (rng->s[0] == 0 && rng->s[1] == 0 && rng->s[2] == 0 && rng->s[3] == 0) {
        rng->s[0] = 0xa12ef3383a3d2eefULL + seed;
        rng->s[1] = 0xcdfabe82ecd412ceULL + seed;
        rng->s[2] = 0x90b5ec55c9235815ULL + seed;
        rng->s[3] = 0xcfb28093ca79a3a7ULL + seed;
    }
    // warm up
    lmmp_xoshiro256pp_random(rng);
    lmmp_xoshiro256pp_random(rng);
    lmmp_xoshiro256pp_random(rng);
}

引用了 lmmp_param_assert, lmmp_seed_generator(), lmmp_xoshiro256pp_random() , 以及 xoshiro256pp_state::s.

被这些函数引用 lmmp_seed_random_().

函数调用图:

这是这个函数的调用关系图:

pcg64_128_action()

static void pcg64_128_action ( mp_limb_t  state[2], const mp_limb_t  inc[2]  )

inlinestatic

在文件 rand_state.h 第 53 行定义.

                                                                          {
    // state = (state * PCG64_MULT + inc) mod 2^128
    mp_limb_t tmp[2] = {0, 0};
    _umul64to128_(state[0], PCG128_DEFAULT_MULTIPLIER_LO, tmp, tmp + 1);
    tmp[1] += state[1] * PCG128_DEFAULT_MULTIPLIER_LO;
    tmp[1] += state[0] * PCG128_DEFAULT_MULTIPLIER_HI;
    _u128add(state, tmp, inc);
}

引用了 _u128add, _umul64to128_(), PCG128_DEFAULT_MULTIPLIER_HI , 以及 PCG128_DEFAULT_MULTIPLIER_LO.

被这些函数引用 lmmp_pcg64_128_random() , 以及 lmmp_pcg64_128_srandom().

函数调用图:

这是这个函数的调用关系图:

pcg64_le_action()

static mp_limb_t pcg64_le_action ( mp_limb_t *restrict  state )

inlinestatic

在文件 rand_state.h 第 178 行定义.

                                                             {
    mp_limb_t old_state = *state;
    *state = old_state * PCG64_LE_MULTIPLIER + PCG64_LE_INCREMENT;
 
    // RXS-M-XS
    mp_limb_t x = old_state;
    int count = x >> 59;
    x ^= x >> (5 + count);
    x *= 12605985483714917081ULL;
    x ^= x >> 43;
 
    return x;
}

引用了 PCG64_LE_INCREMENT , 以及 PCG64_LE_MULTIPLIER.

被这些函数引用 pcg64_le_seq_next().

这是这个函数的调用关系图:

pcg64_le_seq_init()

static void pcg64_le_seq_init ( pcg64_le_seq_t *  rng, mp_size_t  i, mp_limb_t  seed  )

inlinestatic

在文件 rand_state.h 第 138 行定义.

                                                                                 {
    lmmp_param_assert(rng != NULL);
    lmmp_param_assert(rng->k > 0);
    lmmp_param_assert(rng->state != NULL);
 
#define PRIME64_0 0x9E3779B185EBCA87ULL
#define PRIME64_1 0xC2B2AE3D27D4EB4FULL
#define PRIME64_2 0x165667B19E3779F9ULL
#define PRIME64_3 0x85EBCA77C2B2AE63ULL
#define PRIME64_4 0x27D4EB2F165667C5ULL
 
    mp_limb_t s0, s1, s2, s3;
 
    for (; i + 3 < rng->k; i += 4) {
        s0 = rotl(seed + i + 0, 41);
        s1 = rotl(seed + i + 1, 29);
        s2 = rotl(seed + i + 2, 23);
        s3 = rotl(seed + i + 3, 7);
        s0 *= PRIME64_0;
        s1 *= PRIME64_1;
        s2 *= PRIME64_2;
        s3 *= PRIME64_3;
        rng->state[i + 0] = lmmp_seed_generator(s0 ^ rotl(s0, 17));
        rng->state[i + 1] = lmmp_seed_generator(s1 ^ rotl(s1, 21));
        rng->state[i + 2] = lmmp_seed_generator(s2 ^ rotl(s2, 13));
        rng->state[i + 3] = lmmp_seed_generator(s3 ^ rotl(s3, 33));
    }
    for (; i < rng->k; i++) {
        s0 = rotl(seed + i, 31);
        s0 *= PRIME64_4;
        rng->state[i] = lmmp_seed_generator(s0 ^ rotl(s0, 27));
    }
 
#undef PRIME64_0
#undef PRIME64_1
#undef PRIME64_2
#undef PRIME64_3
#undef PRIME64_4
}

引用了 pcg64_le_seq_t::k, lmmp_param_assert, lmmp_seed_generator(), PRIME64_0, PRIME64_1, PRIME64_2, PRIME64_3, PRIME64_4, rotl(), s1, s2, s3 , 以及 pcg64_le_seq_t::state.

被这些函数引用 lmmp_strong_rng_extern_() , 以及 lmmp_strong_rng_init_().

函数调用图:

这是这个函数的调用关系图:

pcg64_le_seq_next()

static void pcg64_le_seq_next ( mp_ptr restrict  dst, mp_size_t  n, pcg64_le_seq_t *  rng  )

inlinestatic

在文件 rand_state.h 第 192 行定义.

                                                                                      {
    lmmp_param_assert(dst != NULL);
    lmmp_param_assert(rng != NULL);
    lmmp_param_assert(n <= rng->k);
    mp_size_t i;
    mp_limb_t mixn = lmmp_seed_generator(n * 0xb9ce52b55c72d585ULL);
    for (i = 0; i + 3 < n; i += 4) {
        dst[i + 0] = pcg64_le_action(&rng->state[i + 0]) ^ mixn;
        dst[i + 1] = pcg64_le_action(&rng->state[i + 1]) ^ mixn;
        dst[i + 2] = pcg64_le_action(&rng->state[i + 2]) ^ mixn;
        dst[i + 3] = pcg64_le_action(&rng->state[i + 3]) ^ mixn;
    }
    for (; i < n; i++) {
        dst[i] = pcg64_le_action(&rng->state[i]) ^ mixn;
    }
    for (; i < rng->k; i++) {
        mp_limb_t old_state = rng->state[i];
        rng->state[i] = old_state * PCG64_LE_MULTIPLIER + PCG64_LE_INCREMENT;
    }
}

引用了 pcg64_le_seq_t::k, k, lmmp_param_assert, lmmp_seed_generator(), pcg64_le_action(), PCG64_LE_INCREMENT, PCG64_LE_MULTIPLIER , 以及 pcg64_le_seq_t::state.

被这些函数引用 lmmp_strong_random_().

函数调用图:

这是这个函数的调用关系图:

rotl()

static mp_limb_t rotl ( const mp_limb_t  x, int  k  )

inlinestatic

在文件 rand_state.h 第 89 行定义.

                                                 {
    const int shift = k & 63;
    return (x << shift) | (x >> (64 - shift));
}

引用了 k.

被这些函数引用 lmmp_random_(), lmmp_strong_rng_extern_(), lmmp_strong_rng_init_(), lmmp_xoshiro256pp_random() , 以及 pcg64_le_seq_init().

这是这个函数的调用关系图: