From 150d96e74e778048e8e87437d54109d621583db8 Mon Sep 17 00:00:00 2001
From: lisyarus <lisyarus@gmail.com>
Date: Sun, 13 Sep 2020 11:03:06 +0300
Subject: [PATCH] Add point in box generation

---
 .../include/psemek/pcg/random/point_box.hpp   | 42 +++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100644 libs/pcg/include/psemek/pcg/random/point_box.hpp

diff --git a/libs/pcg/include/psemek/pcg/random/point_box.hpp b/libs/pcg/include/psemek/pcg/random/point_box.hpp
new file mode 100644
index 00000000..af90bd53
--- /dev/null
+++ b/libs/pcg/include/psemek/pcg/random/point_box.hpp
@@ -0,0 +1,42 @@
+#pragma once
+
+#include <psemek/geom/point.hpp>
+#include <psemek/geom/box.hpp>
+
+#include <random>
+
+namespace psemek::pcg
+{
+
+	template <typename T, std::size_t N>
+	struct box_point_distribution
+	{
+		std::uniform_real_distribution<T> d[N];
+
+		box_point_distribution(geom::box<T, N> const & b)
+		{
+			for (std::size_t i = 0; i < N; ++i)
+			{
+				d[i] = std::uniform_real_distribution<T>{b[i].min, b[i].max};
+			}
+		}
+
+		box_point_distribution()
+		{
+			for (std::size_t i = 0; i < N; ++i)
+			{
+				d[i] = std::uniform_real_distribution<T>{T{0}, T{1}};
+			}
+		}
+
+		template <typename RNG>
+		auto operator()(RNG && rng)
+		{
+			geom::point<T, N> p;
+			for (std::size_t i = 0; i < N; ++i)
+				p[i] = d[i](rng);
+			return p;
+		}
+	};
+
+}