X-Git-Url: https://git.sesse.net/?p=stockfish;a=blobdiff_plain;f=src%2Frkiss.h;h=d017849568ac67aa70495f1d14dbf558c042da8f;hp=0863c9cf7520ce279c4f1666f97e17d4ae52f290;hb=976270916bbfb9c4aa1f2b8ee2897b582a238e25;hpb=76506bd3d10214aa3f3d12dad60e8da729fc318c diff --git a/src/rkiss.h b/src/rkiss.h index 0863c9cf..d0178495 100644 --- a/src/rkiss.h +++ b/src/rkiss.h @@ -20,20 +20,6 @@ available under the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. - - ** A small "keep it simple and stupid" RNG with some fancy merits: - ** - ** Quite platform independent - ** Passes ALL dieharder tests! Here *nix sys-rand() e.g. fails miserably:-) - ** ~12 times faster than my *nix sys-rand() - ** ~4 times faster than SSE2-version of Mersenne twister - ** Average cycle length: ~2^126 - ** 64 bit seed - ** Return doubles with a full 53 bit mantissa - ** Thread safe - ** - ** (c) Heinz van Saanen - */ #if !defined(RKISS_H_INCLUDED) @@ -41,12 +27,26 @@ #include "types.h" +/// RKISS is our pseudo random number generator (PRNG) used to compute hash keys. +/// George Marsaglia invented the RNG-Kiss-family in the early 90's. This is a +/// specific version that Heinz van Saanen derived from some public domain code +/// by Bob Jenkins. Following the feature list, as tested by Heinz. +/// +/// - Quite platform independent +/// - Passes ALL dieharder tests! Here *nix sys-rand() e.g. fails miserably:-) +/// - ~12 times faster than my *nix sys-rand() +/// - ~4 times faster than SSE2-version of Mersenne twister +/// - Average cycle length: ~2^126 +/// - 64 bit seed +/// - Return doubles with a full 53 bit mantissa +/// - Thread safe + class RKISS { // Keep variables always together struct S { uint64_t a, b, c, d; } s; - uint64_t rot(uint64_t x, uint64_t k) const { + uint64_t rotate(uint64_t x, uint64_t k) const { return (x << k) | (x >> (64 - k)); } @@ -54,9 +54,9 @@ class RKISS { uint64_t rand64() { const uint64_t - e = s.a - rot(s.b, 7); - s.a = s.b ^ rot(s.c, 13); - s.b = s.c + rot(s.d, 37); + e = s.a - rotate(s.b, 7); + s.a = s.b ^ rotate(s.c, 13); + s.b = s.c + rotate(s.d, 37); s.c = s.d + e; return s.d = e + s.a; }