From 4410597e441f38320f095f1b12d03861181f73ce Mon Sep 17 00:00:00 2001
From: lisyarus <lisyarus@gmail.com>
Date: Wed, 23 Sep 2020 22:23:37 +0300
Subject: [PATCH] Add geom::ray primitive

---
 libs/geom/include/psemek/geom/ray.hpp | 59 +++++++++++++++++++++++++++
 1 file changed, 59 insertions(+)
 create mode 100644 libs/geom/include/psemek/geom/ray.hpp

diff --git a/libs/geom/include/psemek/geom/ray.hpp b/libs/geom/include/psemek/geom/ray.hpp
new file mode 100644
index 00000000..a26a743f
--- /dev/null
+++ b/libs/geom/include/psemek/geom/ray.hpp
@@ -0,0 +1,59 @@
+#pragma once
+
+#include <psemek/geom/vector.hpp>
+#include <psemek/geom/point.hpp>
+
+namespace psemek::geom
+{
+
+	template <typename T, std::size_t N>
+	struct ray
+	{
+		point<T, N> origin;
+		vector<T, N> direction;
+
+		auto operator()(T const & t) const
+		{
+			return origin + t * direction;
+		}
+
+		ray & operator += (vector<T, N> const & v)
+		{
+			origin += v;
+			return *this;
+		}
+
+		ray & operator -= (vector<T, N> const & v)
+		{
+			origin -= v;
+			return *this;
+		}
+	};
+
+	template <typename T, std::size_t N>
+	ray(point<T, N>, vector<T, N>) -> ray<T, N>;
+
+	template <typename T, std::size_t N>
+	ray<T, N> advance(ray<T, N> r, T const & t)
+	{
+		r.origin = r(t);
+		return r;
+	}
+
+	template <typename T, std::size_t N>
+	ray<T, N> operator + (ray<T, N> const & r, vector<T, N> const & v)
+	{
+		auto res = r;
+		res += v;
+		return res;
+	}
+
+	template <typename T, std::size_t N>
+	ray<T, N> operator - (ray<T, N> const & r, vector<T, N> const & v)
+	{
+		auto res = r;
+		res -= v;
+		return res;
+	}
+
+}