We provide a new random number generator (RNG) with this package, available both in Fortran90 and C++. It resides in src/f90/mod/rngmod.f90 and src/cxx/cxxsupport/planck_rng.h and supersedes the previous RNG (which is still available at src/f90/mod/ran_tools_dist.f90).
It produces double precision real numbers with and is based on a xorshift method described by Marsaglia in Journal of Statistical Software 2003, vol 8 (cf. https://www.cs.hku.hk/). It accepts up to four different seeds simultaneously, allowing each sequence to have a theoretical period of . A Gaussian deviate RNG is also provided. See the respective routines documentation for details on their usage. Please note that we have not extensively tested this generator — it did not represent the main drive of this project.
Version 3.82, 2022-07-28