From a5e659bad596ece23da895a21b91cd018a695596 Mon Sep 17 00:00:00 2001
From: notgiven688 <37874600+notgiven688@users.noreply.github.com>
Date: Fri, 7 Jun 2024 00:57:40 +0200
Subject: [PATCH] Optimized quaternion vector transformation (#153)

---
 src/Jitter2/LinearMath/JVector.cs | 50 +++++++++++++------------------
 1 file changed, 20 insertions(+), 30 deletions(-)

diff --git a/src/Jitter2/LinearMath/JVector.cs b/src/Jitter2/LinearMath/JVector.cs
index 71809a45..a26928db 100644
--- a/src/Jitter2/LinearMath/JVector.cs
+++ b/src/Jitter2/LinearMath/JVector.cs
@@ -256,22 +256,17 @@ public static void TransposedTransform(in JVector vector, in JMatrix matrix, out
     [MethodImpl(MethodImplOptions.AggressiveInlining)]
     public static void Transform(in JVector vector, in JQuaternion quaternion, out JVector result)
     {
-        float x = vector.X;
-        float y = vector.Y;
-        float z = vector.Z;
-        float qx = quaternion.X;
-        float qy = quaternion.Y;
-        float qz = quaternion.Z;
-        float qw = quaternion.W;
+        float numx = 2.0f * (quaternion.Y * vector.Z - quaternion.Z * vector.Y);
+        float numy = 2.0f * (quaternion.Z * vector.X - quaternion.X * vector.Z);
+        float numz = 2.0f * (quaternion.X * vector.Y - quaternion.Y * vector.X);
 
-        float num0 = qw * x + qy * z - qz * y;
-        float num1 = qw * y + qz * x - qx * z;
-        float num2 = qw * z + qx * y - qy * x;
-        float num3 = -qx * x - qy * y - qz * z;
+        float num00 = quaternion.Y * numz - quaternion.Z * numy;
+        float num11 = quaternion.Z * numx - quaternion.X * numz;
+        float num22 = quaternion.X * numy - quaternion.Y * numx;
 
-        result.X = num0 * qw - num3 * qx - num1 * qz + num2 * qy;
-        result.Y = num1 * qw - num3 * qy - num2 * qx + num0 * qz;
-        result.Z = num2 * qw - num3 * qz - num0 * qy + num1 * qx;
+        result.X = vector.X + quaternion.W * numx + num00;
+        result.Y = vector.Y + quaternion.W * numy + num11;
+        result.Z = vector.Z + quaternion.W * numz + num22;
     }
 
     /// <summary>
@@ -280,22 +275,17 @@ public static void Transform(in JVector vector, in JQuaternion quaternion, out J
     [MethodImpl(MethodImplOptions.AggressiveInlining)]
     public static void ConjugatedTransform(in JVector vector, in JQuaternion quaternion, out JVector result)
     {
-        float x = vector.X;
-        float y = vector.Y;
-        float z = vector.Z;
-        float qx = quaternion.X;
-        float qy = quaternion.Y;
-        float qz = quaternion.Z;
-        float qw = quaternion.W;
-
-        float num0 = qw * x - qy * z + qz * y;
-        float num1 = qw * y - qz * x + qx * z;
-        float num2 = qw * z - qx * y + qy * x;
-        float num3 = qx * x + qy * y + qz * z;
-
-        result.X = num0 * qw + num3 * qx + num1 * qz - num2 * qy;
-        result.Y = num1 * qw + num3 * qy + num2 * qx - num0 * qz;
-        result.Z = num2 * qw + num3 * qz + num0 * qy - num1 * qx;
+        float numx = 2.0f * (quaternion.Z * vector.Y - quaternion.Y * vector.Z);
+        float numy = 2.0f * (quaternion.X * vector.Z - quaternion.Z * vector.X);
+        float numz = 2.0f * (quaternion.Y * vector.X - quaternion.X * vector.Y);
+
+        float num00 = quaternion.Z * numy - quaternion.Y * numz;
+        float num11 = quaternion.X * numz - quaternion.Z * numx;
+        float num22 = quaternion.Y * numx - quaternion.X * numy;
+
+        result.X = vector.X + quaternion.W * numx + num00;
+        result.Y = vector.Y + quaternion.W * numy + num11;
+        result.Z = vector.Z + quaternion.W * numz + num22;
     }
 
     /// <summary>