prng: make rand_r have 2^32 period instead of 2^31

this is a minor fix to increase the period of the obsolete rand_r a bit.
an include header in __rand48_step.c is fixed as well.

diff --git a/src/prng/__rand48_step.c b/src/prng/__rand48_step.c
index 755b4f2f3903e33357684d33d27f8b2e4de3d2ae..105f810ff017d51e8979b049bae0167ec612a497 100644 (file)
--- a/src/prng/__rand48_step.c
+++ b/src/prng/__rand48_step.c
@@ -1,5 +1,5 @@
 #include <stdlib.h>
-#include <inttypes.h>
+#include <stdint.h>
 
 uint64_t __rand48_step(unsigned short *xi, unsigned short *lc)
 {

diff --git a/src/prng/rand_r.c b/src/prng/rand_r.c
index e96cfba99cdb9c8cb4b2c0d834389b864e857e23..ef26dbde6b673e27ebdf65189680559733cb01f0 100644 (file)
--- a/src/prng/rand_r.c
+++ b/src/prng/rand_r.c
@@ -2,5 +2,5 @@
 
 int rand_r(unsigned *seed)
 {
-       return (*seed = *seed * 1103515245 + 12345) & 0x7fffffff;
+       return (*seed = *seed * 1103515245 + 12345)/2;
 }