From 3bf42c6ff3805a0d42bbc661794a95ff31bedc26 Mon Sep 17 00:00:00 2001
From: untodesu <kirill@untode.su>
Date: Sat, 15 Mar 2025 16:22:09 +0500
Subject: Add whatever I was working on for the last month

---
 deps/include/glm/ext/quaternion_trigonometric.inl | 37 +++++++++++++++++++++++
 1 file changed, 37 insertions(+)
 create mode 100644 deps/include/glm/ext/quaternion_trigonometric.inl

(limited to 'deps/include/glm/ext/quaternion_trigonometric.inl')

diff --git a/deps/include/glm/ext/quaternion_trigonometric.inl b/deps/include/glm/ext/quaternion_trigonometric.inl
new file mode 100644
index 0000000..896449a
--- /dev/null
+++ b/deps/include/glm/ext/quaternion_trigonometric.inl
@@ -0,0 +1,37 @@
+#include "scalar_constants.hpp"
+
+namespace glm
+{
+	template<typename T, qualifier Q>
+	GLM_FUNC_QUALIFIER T angle(qua<T, Q> const& x)
+	{
+		if (abs(x.w) > cos_one_over_two<T>())
+		{
+			T const a = asin(sqrt(x.x * x.x + x.y * x.y + x.z * x.z)) * static_cast<T>(2);
+			if(x.w < static_cast<T>(0))
+				return pi<T>() * static_cast<T>(2) - a;
+			return a;
+		}
+
+		return acos(x.w) * static_cast<T>(2);
+	}
+
+	template<typename T, qualifier Q>
+	GLM_FUNC_QUALIFIER vec<3, T, Q> axis(qua<T, Q> const& x)
+	{
+		T const tmp1 = static_cast<T>(1) - x.w * x.w;
+		if(tmp1 <= static_cast<T>(0))
+			return vec<3, T, Q>(0, 0, 1);
+		T const tmp2 = static_cast<T>(1) / sqrt(tmp1);
+		return vec<3, T, Q>(x.x * tmp2, x.y * tmp2, x.z * tmp2);
+	}
+
+	template<typename T, qualifier Q>
+	GLM_FUNC_QUALIFIER qua<T, Q> angleAxis(T const& angle, vec<3, T, Q> const& v)
+	{
+		T const a(angle);
+		T const s = glm::sin(a * static_cast<T>(0.5));
+
+		return qua<T, Q>(glm::cos(a * static_cast<T>(0.5)), v * s);
+	}
+}//namespace glm
-- 
cgit