6#ifndef THEORETICA_PRNG_H
7#define THEORETICA_PRNG_H
11#include "../core/error.h"
30 std::vector<uint64_t> param;
38 const std::vector<uint64_t>& s) : f(p), x(
seed), param(s) {}
61 return x = f(x, param);
75 for (uint64_t i = 0; i < n; ++i)
81 inline uint64_t
last()
const {
98 inline void set_param(
const std::vector<uint64_t>& v) {
130 {48271, 0, ((uint64_t) 1 << 31) - 1}
177 uint64_t p1 = 2549536629329, uint64_t p2 = 136137137) {
191 uint64_t
seed, uint64_t offset = 765872292751861) {
197 offset = 765872292751861;
211 template<
typename Vector>
217 for (
unsigned int i = 0; i < rounds; ++i) {
220 const size_t index1 = g() % v.size();
221 const size_t index2 = g() % v.size();
224 const auto x1 = v[index1];
225 v[index1] = v[index2];
237 template<
typename Vector>
A pseudorandom number generator.
Definition prng.h:19
static PRNG middlesquare(uint64_t seed, uint64_t offset=765872292751861)
Returns a Middle-square generator.
Definition prng.h:190
uint64_t last() const
Return the last generated number.
Definition prng.h:81
static PRNG linear_congruential(uint64_t seed=1)
Returns a standard linear congruential generator.
Definition prng.h:123
void set_param(const std::vector< uint64_t > &v)
Set the generator's parameters.
Definition prng.h:98
void seed(uint64_t seed)
Seed the PRNG.
Definition prng.h:54
uint64_t operator()()
Generate a pseudorandom number.
Definition prng.h:67
PRNG(pseudorandom_function p, uint64_t seed)
Construct a PRNG with the given generating algorithm and seed.
Definition prng.h:43
static PRNG wyrand(uint64_t seed=1, uint64_t p1=2549536629329, uint64_t p2=136137137)
Returns a Wyrand generator.
Definition prng.h:176
void discard(uint64_t n)
Discard n numbers from the generator.
Definition prng.h:74
pseudorandom_function get_function() const
Get the generating function.
Definition prng.h:92
static PRNG xoshiro(uint64_t seed=1)
Returns a Xoshiro256++ generator.
Definition prng.h:149
std::vector< uint64_t > get_param() const
Get the generator's parameters.
Definition prng.h:108
static PRNG xoshiro(const std::vector< uint64_t > &p)
Returns a Xoshiro256++ generator.
Definition prng.h:138
static PRNG splitmix64(uint64_t seed=1)
Returns a Splitmix64 generator.
Definition prng.h:166
PRNG & operator>>(uint64_t &n)
Stream the next generated number.
Definition prng.h:114
void set_function(pseudorandom_function p)
Set the generating function.
Definition prng.h:87
void set_param(unsigned int i, uint64_t value)
Set a specific parameter by index.
Definition prng.h:103
PRNG(uint64_t seed)
Construct a PRNG with the default generator and the given seed.
Definition prng.h:48
uint64_t next()
Generate a pseudorandom number.
Definition prng.h:60
PRNG(pseudorandom_function p, uint64_t seed, const std::vector< uint64_t > &s)
Construct a PRNG with the given generating algorithm p, seed x and parameters s.
Definition prng.h:36
#define TH_MATH_ERROR(F_NAME, VALUE, EXCEPTION)
TH_MATH_ERROR is a macro which throws exceptions or modifies errno (depending on which compiling opti...
Definition error.h:225
Main namespace of the library which contains all functions and objects.
Definition algebra.h:27
uint64_t randgen_xoshiro(uint64_t &a, uint64_t &b, uint64_t &c, uint64_t &d)
Generate a pseudorandom natural number using the xoshiro256++ pseudorandom number generation algorith...
Definition pseudorandom.h:86
uint64_t randgen_splitmix64(uint64_t x)
Generate a pseudorandom natural number using the SplitMix64 pseudorandom number generation algorithm.
Definition pseudorandom.h:130
uint64_t(*)(uint64_t, std::vector< uint64_t > &) pseudorandom_function
A function pointer which wraps a pseudorandom generator, taking as input the previous generated value...
Definition pseudorandom.h:23
uint64_t randgen_middlesquare(uint64_t seed, uint64_t offset=765872292751861)
Generate a pseudorandom natural number using the middle-square pseudorandom number generation algorit...
Definition pseudorandom.h:194
uint64_t randgen_congruential(uint64_t x, uint64_t a=48271, uint64_t c=0, uint64_t m=((uint64_t) 1<< 31) - 1)
Generate a pseudorandom natural number using the congruential pseudorandom number generation algorith...
Definition pseudorandom.h:39
uint64_t randgen_wyrand(uint64_t &seed, uint64_t p1, uint64_t p2)
Generate a pseudorandom natural number using the Wyrand pseudorandom number generation,...
Definition pseudorandom.h:160
dual2 square(dual2 x)
Return the square of a second order dual number.
Definition dual2_functions.h:23
Vector & shuffle(Vector &v, PRNG &g, unsigned int rounds)
Shuffle a set by exchanging random couples of elements.
Definition prng.h:212
Pseudorandom number generation algorithms.