From 0d664b462d02e651c4adf569007589b462e87462 Mon Sep 17 00:00:00 2001
From: John Harwell <john.r.harwell@gmail.com>
Date: Mon, 3 Mar 2025 13:21:56 -0600
Subject: [PATCH] Fix {normal, extreme_value} distributions to work with
 quantities

---
 src/core/include/mp-units/random.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/core/include/mp-units/random.h b/src/core/include/mp-units/random.h
index 5da89966..dd0db72b 100644
--- a/src/core/include/mp-units/random.h
+++ b/src/core/include/mp-units/random.h
@@ -289,7 +289,7 @@ struct extreme_value_distribution : public std::extreme_value_distribution<typen
   template<typename Generator>
   Q operator()(Generator& g)
   {
-    return Q(base::operator()(g));
+    return base::operator()(g) * Q::reference;
   }
 
   [[nodiscard]] Q a() const { return base::a() * Q::reference; }
@@ -313,7 +313,7 @@ struct normal_distribution : public std::normal_distribution<typename Q::rep> {
   template<typename Generator>
   Q operator()(Generator& g)
   {
-    return Q(base::operator()(g));
+    return base::operator()(g) * Q::reference;
   }
 
   [[nodiscard]] Q mean() const { return base::mean() * Q::reference; }