1 files changed, 58 insertions, 0 deletions

diff --git a/external/include/glm/gtx/vector_angle.inl b/external/include/glm/gtx/vector_angle.inl
new file mode 100644
index 0000000..05c3028
--- /dev/null
+++ b/external/include/glm/gtx/vector_angle.inl
@@ -0,0 +1,58 @@
+/// @ref gtx_vector_angle
+/// @file glm/gtx/vector_angle.inl
+
+namespace glm
+{
+	template <typename genType> 
+	GLM_FUNC_QUALIFIER genType angle
+	(
+		genType const & x,
+		genType const & y
+	)
+	{
+		GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'angle' only accept floating-point inputs");
+		return acos(clamp(dot(x, y), genType(-1), genType(1)));
+	}
+
+	template <typename T, precision P, template <typename, precision> class vecType> 
+	GLM_FUNC_QUALIFIER T angle
+	(
+		vecType<T, P> const & x,
+		vecType<T, P> const & y
+	)
+	{
+		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'angle' only accept floating-point inputs");
+		return acos(clamp(dot(x, y), T(-1), T(1)));
+	}
+
+	//! \todo epsilon is hard coded to 0.01
+	template <typename T, precision P>
+	GLM_FUNC_QUALIFIER T orientedAngle
+	(
+		tvec2<T, P> const & x,
+		tvec2<T, P> const & y
+	)
+	{
+		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'orientedAngle' only accept floating-point inputs");
+		T const Angle(acos(clamp(dot(x, y), T(-1), T(1))));
+
+		if(all(epsilonEqual(y, glm::rotate(x, Angle), T(0.0001))))
+			return Angle;
+		else
+			return -Angle;
+	}
+
+	template <typename T, precision P>
+	GLM_FUNC_QUALIFIER T orientedAngle
+	(
+		tvec3<T, P> const & x,
+		tvec3<T, P> const & y,
+		tvec3<T, P> const & ref
+	)
+	{
+		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'orientedAngle' only accept floating-point inputs");
+
+		T const Angle(acos(clamp(dot(x, y), T(-1), T(1))));
+		return mix(Angle, -Angle, dot(ref, cross(x, y)) < T(0));
+	}
+}//namespace glm