From a9ce3df5c59515c12d4fcdf3b7ad43cefdf122b6 Mon Sep 17 00:00:00 2001
From: Tyge Løvset <tylo@norceresearch.no>
Date: Mon, 21 Dec 2020 14:28:56 +0100
Subject: Added some benchmarks. Fix typo in cdeq.h

---
 benchmarks/crand_benchmark.cpp | 155 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 155 insertions(+)
 create mode 100644 benchmarks/crand_benchmark.cpp

(limited to 'benchmarks/crand_benchmark.cpp')

diff --git a/benchmarks/crand_benchmark.cpp b/benchmarks/crand_benchmark.cpp
new file mode 100644
index 00000000..78099905
--- /dev/null
+++ b/benchmarks/crand_benchmark.cpp
@@ -0,0 +1,155 @@
+#include <cstdint>
+#include <iostream>
+#include <ctime>
+#include <stc/crand.h>
+
+static inline uint64_t rotl64(const uint64_t x, const int k)
+	{ return (x << k) | (x >> (64 - k)); }
+
+static uint64_t splitmix64_x = 87213627321ull; /* The state can be seeded with any value. */
+
+uint64_t splitmix64(void) {
+	uint64_t z = (splitmix64_x += 0x9e3779b97f4a7c15);
+	z = (z ^ (z >> 30)) * 0xbf58476d1ce4e5b9;
+	z = (z ^ (z >> 27)) * 0x94d049bb133111eb;
+	return z ^ (z >> 31);
+}
+
+static void init_state(uint64_t *rng, uint64_t seed) {
+    splitmix64_x = seed;
+    for (int i=0; i<4; ++i) rng[i] = splitmix64();
+}
+
+/* jsf64 */
+
+static inline uint64_t jsf64(uint64_t *s) {
+    uint64_t e = s[0] - rotl64(s[1], 7);
+    s[0] = s[1] ^ rotl64(s[2], 13);
+    s[1] = s[2] + rotl64(s[3], 37);
+    s[2] = s[3] + e;
+    s[3] = e + s[0];
+    return s[3];
+}
+
+/* xoshiro256**  */
+
+static inline uint64_t xoshiro256starstar(uint64_t* s) {
+    const uint64_t result = rotl64(s[1] * 5, 7) * 9;
+    const uint64_t t = s[1] << 17;
+    s[2] ^= s[0];
+    s[3] ^= s[1];
+    s[1] ^= s[2];
+    s[0] ^= s[3];
+    s[2] ^= t;
+    s[3] = rotl64(s[3], 45);
+    return result;
+}
+
+
+inline unsigned long long wyhash64(uint64_t* s) {
+  *s += 0x60bee2bee120fc15ull;
+  __uint128_t tmp = (__uint128_t) *s * 0xa3b195354a39b70dull;
+  unsigned long long m1 = (tmp >> 64) ^ tmp;
+  tmp = (__uint128_t)m1 * 0x1b03738712fad5c9ull;
+  return (tmp >> 64) ^ tmp;
+}
+
+inline unsigned long long lehmer64(uint64_t* s) {
+  *(__uint128_t *)s *= 0xda942042e4dd58b5ull;
+  return *(__uint128_t *)s >> 64;
+}
+
+using namespace std;
+
+int main(void)
+{
+  enum {N = 524288000};
+  uint64_t* recipient = new uint64_t[N];
+  static crand_t rng;
+  init_state(rng.state, 12345123);
+
+  clock_t beg, end;
+  for (size_t ti = 0; ti < 4; ti++) {
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      recipient[i] = wyhash64(rng.state);
+    end = clock();
+    cerr << "ROUND " << ti+1 << endl
+         << "wyhash64:\t"
+         << (float(end - beg) / CLOCKS_PER_SEC)
+         << " s" << endl;
+    cout << "bogus:" << recipient[312] << endl;
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      recipient[i] = crand_next(&rng);
+    end = clock();
+    cerr << "stc crand:\t"
+         << (float(end - beg) / CLOCKS_PER_SEC)
+         << " s" << endl;
+    cout << "bogus:" << recipient[312] << endl;
+
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      recipient[i] = xoshiro256starstar(rng.state);
+    end = clock();
+    cerr << "xoshiro256**:\t"
+         << (float(end - beg) / CLOCKS_PER_SEC)
+         << " s" << endl;
+    cout << "bogus:" << recipient[312] << endl;
+
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      recipient[i] = lehmer64(rng.state);
+    end = clock();
+    cerr << "lehmer64:\t"
+         << ((float) end - beg) / CLOCKS_PER_SEC
+         << " s" << endl;
+    cout << "bogus:" << recipient[312] << endl;
+
+
+
+    cout << endl
+         << "Next we do random number computations only, doing no work."
+         << endl;
+    uint64_t s = 0;
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      s += wyhash64(rng.state);
+    end = clock();
+    cerr << "wyhash64:\t"
+         << ((float) end - beg) / CLOCKS_PER_SEC
+         << " s" << endl;
+    cout << "bogus:" << s << endl;
+
+    s = 0;
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      s += crand_next(&rng);
+    end = clock();
+    cerr << "stc crand:\t"
+         << ((float) end - beg) / CLOCKS_PER_SEC
+         << " s" << endl;
+    cout << "bogus:" << s << endl;
+
+    s = 0;
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      s += xoshiro256starstar(rng.state);
+    end = clock();
+    cerr << "xoshiro256**:\t"
+         << ((float) end - beg) / CLOCKS_PER_SEC
+         << " s" << endl;
+    cout << "bogus:" << s << endl;
+
+    beg = clock();
+    for (size_t i = 0; i < N; i++)
+      s += lehmer64(rng.state);
+    end = clock();
+    cerr << "lehmer64:\t"
+         << ((float) end - beg) / CLOCKS_PER_SEC
+         << " s" << endl;
+    cout << "bogus:" << s << endl << endl;
+  }
+  delete[] recipient;
+  return 0;
+}
-- 
cgit v1.2.3