new file mode 100644

@@ -0,0 +1,36 @@

+#ifndef THEFT_MT_H

+#define THEFT_MT_H

+

+#include <stdint.h>

+

+/* Wrapper for Mersenne Twister.

+ * See copyright and license in theft_mt.c, more details at:

+ *     http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html

+ *

+ * The code has been modified to store internal state in heap/stack

+ * allocated memory, rather than statically allocated memory, to allow

+ * multiple instances running in the same address space. */

+

+#define THEFT_MT_PARAM_N 312

+

+struct theft_mt {

+    uint64_t mt[THEFT_MT_PARAM_N]; /* the array for the state vector  */

+    int16_t mti;

+};

+

+/* Heap-allocate a mersenne twister struct. */

+struct theft_mt *theft_mt_init(uint64_t seed);

+

+/* Free a heap-allocated mersenne twister struct. */

+void theft_mt_free(struct theft_mt *mt);

+

+/* Reset a mersenne twister struct, possibly stack-allocated. */

+void theft_mt_reset(struct theft_mt *mt, uint64_t seed);

+

+/* Get a 64-bit random number. */

+uint64_t theft_mt_random(struct theft_mt *mt);

+

+/* Generate a random number on [0,1]-real-interval. */

+double theft_mt_random_double(struct theft_mt *mt);

+

+#endif