1 files changed, 36 insertions, 0 deletions
diff --git a/deps/theft/theft_mt.h b/deps/theft/theft_mt.h
new file mode 100644
index 0000000..65d7098
--- /dev/null
+++ b/deps/theft/theft_mt.h
@@ -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