From 4abab5ad6c8465a7528ccdd5f49367da05f78bbd Mon Sep 17 00:00:00 2001
From: Vladimir Azarov <avm@intermediate-node.net>
Date: Tue, 1 Oct 2024 15:47:05 +0200
Subject: Initial version

---
 src/prng/drand48.c | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)
 create mode 100644 src/prng/drand48.c

(limited to 'src/prng/drand48.c')

diff --git a/src/prng/drand48.c b/src/prng/drand48.c
new file mode 100644
index 0000000..08283e2
--- /dev/null
+++ b/src/prng/drand48.c
@@ -0,0 +1,17 @@
+#include <stdlib.h>
+#include <inttypes.h>
+#include "rand48.h"
+
+double erand48(unsigned short s[3])
+{
+	union {
+		uint64_t u;
+		double f;
+	} x = { 0x3ff0000000000000ULL | __rand48_step(s, __seed48+3)<<4 };
+	return x.f - 1.0;
+}
+
+double drand48(void)
+{
+	return erand48(__seed48);
+}
-- 
cgit v1.2.3