xoroshiro.h

Go to the documentation of this file.

 
#pragma once
 
#include <stdint.h>
 
#define rotl(x, k) (((x) << (k)) | ((x) >> (64 - (k))))
 
#define random_u64(rng_s)                                                                                              \
    __extension__({                                                                                                \
        const uint64_t __res = rotl((rng_s)[1] * 5, 7) * 9;                                                    \
        const uint64_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);                                                                     \
                                                                                                                       \
        __res;                                                                                                 \
    })