From 00b8e8743c504fb5169977c4f92629cf4c6f9c7b Mon Sep 17 00:00:00 2001
From: James McClung <james@themcclungs.net>
Date: Mon, 3 Nov 2025 09:59:06 -0500
Subject: [PATCH 1/2] rng: rename GeneratorPtr

---
 src/include/rng.hxx | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/include/rng.hxx b/src/include/rng.hxx
index c8381fa72..7a30cc248 100644
--- a/src/include/rng.hxx
+++ b/src/include/rng.hxx
@@ -15,8 +15,8 @@ namespace rng
 namespace detail
 {
 using Generator = std::default_random_engine;
-using ProcessGeneratorPtr = std::shared_ptr<Generator>;
-ProcessGeneratorPtr process_generator;
+using GeneratorPtr = std::shared_ptr<Generator>;
+GeneratorPtr process_generator;
 
 int get_process_seed()
 {
@@ -25,7 +25,7 @@ int get_process_seed()
   return 10000 * (rank + 1);
 }
 
-ProcessGeneratorPtr get_process_generator()
+GeneratorPtr get_process_generator()
 {
   if (!process_generator) {
     process_generator = std::make_shared<Generator>(get_process_seed());
@@ -50,7 +50,7 @@ struct Uniform
   Real get() { return dist(*gen); }
 
 private:
-  detail::ProcessGeneratorPtr gen;
+  detail::GeneratorPtr gen;
   std::uniform_real_distribution<Real> dist;
 };
 
@@ -74,7 +74,7 @@ struct Normal
   }
 
 private:
-  detail::ProcessGeneratorPtr gen;
+  detail::GeneratorPtr gen;
   std::normal_distribution<Real> dist;
 };
 

From 680d0a62a44ad2d58af357079f9af6028b2eac9e Mon Sep 17 00:00:00 2001
From: James McClung <james@themcclungs.net>
Date: Mon, 3 Nov 2025 10:48:49 -0500
Subject: [PATCH 2/2] rng: +constant generator support

enables manual seeding
---
 src/include/rng.hxx | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/src/include/rng.hxx b/src/include/rng.hxx
index 7a30cc248..9b8c7236e 100644
--- a/src/include/rng.hxx
+++ b/src/include/rng.hxx
@@ -33,6 +33,12 @@ GeneratorPtr get_process_generator()
   return process_generator;
 }
 
+GeneratorPtr get_constant_generator(int seed)
+{
+  // FIXME this doesn't need to be a shared pointer
+  return std::make_shared<Generator>(seed);
+}
+
 } // namespace detail
 
 // ======================================================================
@@ -41,6 +47,10 @@ GeneratorPtr get_process_generator()
 template <typename Real>
 struct Uniform
 {
+  Uniform(Real min, Real max, int seed)
+    : dist(min, max), gen(detail::get_constant_generator(seed))
+  {}
+
   Uniform(Real min, Real max)
     : dist(min, max), gen(detail::get_process_generator())
   {}
@@ -60,9 +70,14 @@ private:
 template <typename Real>
 struct Normal
 {
+  Normal(Real mean, Real stdev, int seed)
+    : dist(mean, stdev), gen(detail::get_constant_generator(seed))
+  {}
+
   Normal(Real mean, Real stdev)
     : dist(mean, stdev), gen(detail::get_process_generator())
   {}
+
   Normal() : Normal(0, 1) {}
 
   Real get() { return dist(*gen); }