Add ranqd1 from Numerical Recipes in C, 2nd Edition.

Author: polytypic

internals/testing/ranqd1_test.cpp

@@ -0,0 +1,21 @@
+#include "dumpster_v1/ranqd1.hpp"
+
+#include "testing_v1/test.hpp"
+
+using namespace testing_v1;
+using namespace dumpster_v1;
+
+auto ranqd1_test = test([]() {
+  uint32_t seed = 0;
+  verify((seed = ranqd1(seed)) == 0x3c6ef35f);
+  verify((seed = ranqd1(seed)) == 0x47502932);
+  verify((seed = ranqd1(seed)) == 0xd1ccf6e9);
+  verify((seed = ranqd1(seed)) == 0xaaf95334);
+  verify((seed = ranqd1(seed)) == 0x6252e503);
+  verify((seed = ranqd1(seed)) == 0x9f2ec686);
+  verify((seed = ranqd1(seed)) == 0x57fe6c2d);
+  verify((seed = ranqd1(seed)) == 0xa3d95fa8);
+  verify((seed = ranqd1(seed)) == 0x81fdbee7);
+  verify((seed = ranqd1(seed)) == 0x94f0af1a);
+  verify((seed = ranqd1(seed)) == 0xcbf633b1);
+});

provides/include/dumpster_v1/ranqd1.hpp

@@ -0,0 +1,7 @@
+#pragma once
+
+#include "dumpster_v1/synopsis.hpp"
+
+inline uint32_t dumpster_v1::ranqd1(uint32_t seed) {
+  return seed * 1664525 + 1013904223;
+}

provides/include/dumpster_v1/synopsis.hpp

@@ -1,5 +1,7 @@
 #pragma once
 
+#include <cstdint>
+
 #include "dumpster_v1/private.hpp"
 
 namespace dumpster_v1 {
@@ -53,6 +55,11 @@ template <class Action> struct finally_t : Private::finally_t<Action> {
 template <class Action>
 finally_t<std::remove_cvref_t<Action>> finally(Action &&action);
 
+// ranqd1.hpp ==================================================================
+
+/// The `ranqd1` generator from Numerical Recipes in C, 2nd Edition.
+uint32_t ranqd1(uint32_t seed);
+
 // insertion_sort.hpp ==========================================================
 
 /// Sorts given sentinel terminated sequence to ascending order.  Intended for