From 4409032a9dac162e3126bffac3dea482975d8556 Mon Sep 17 00:00:00 2001
From: Christophe Riccio <g.truc.creation@gmail.com>
Date: Sat, 27 Nov 2010 13:24:28 +0000
Subject: [PATCH] Optimized inverse

---
 glm/core/type_mat3x3.inl | 13 ++++++++++++-
 glm/core/type_mat4x4.inl | 23 ++++++++++++++++++++++-
 2 files changed, 34 insertions(+), 2 deletions(-)

diff --git a/glm/core/type_mat3x3.inl b/glm/core/type_mat3x3.inl
index 6eddeac2..7766dd4d 100644
--- a/glm/core/type_mat3x3.inl
+++ b/glm/core/type_mat3x3.inl
@@ -377,7 +377,7 @@ namespace detail
         T S20 = value[2][0];
         T S21 = value[2][1];
         T S22 = value[2][2];
-
+/*
         tmat3x3<T> Inverse(
             + (S11 * S22 - S21 * S12),
             - (S10 * S22 - S20 * S12),
@@ -388,6 +388,17 @@ namespace detail
             + (S01 * S12 - S11 * S02),
             - (S00 * S12 - S10 * S02),
             + (S00 * S11 - S10 * S01));
+*/
+        tmat3x3<T> Inverse(
+            S11 * S22 - S21 * S12,
+            S12 * S20 - S22 * S10,
+            S10 * S21 - S20 * S11,
+            S02 * S21 - S01 * S22,
+            S00 * S22 - S02 * S20,
+            S01 * S20 - S00 * S21,
+            S12 * S01 - S11 * S02,
+            S10 * S02 - S12 * S00,
+            S11 * S00 - S10 * S01);
 
         T Determinant = S00 * (S11 * S22 - S21 * S12)
                       - S10 * (S01 * S22 - S21 * S02)
diff --git a/glm/core/type_mat4x4.inl b/glm/core/type_mat4x4.inl
index d42e8219..153a2c26 100644
--- a/glm/core/type_mat4x4.inl
+++ b/glm/core/type_mat4x4.inl
@@ -413,7 +413,7 @@ namespace detail
         value_type SubFactor16 = this->value[1][0] * this->value[2][3] - this->value[2][0] * this->value[1][3];
         value_type SubFactor17 = this->value[1][0] * this->value[2][2] - this->value[2][0] * this->value[1][2];
         value_type SubFactor18 = this->value[1][0] * this->value[2][1] - this->value[2][0] * this->value[1][1];
-
+/*
         tmat4x4<T> Inverse(
             + (this->value[1][1] * SubFactor00 - this->value[1][2] * SubFactor01 + this->value[1][3] * SubFactor02),
             - (this->value[1][0] * SubFactor00 - this->value[1][2] * SubFactor03 + this->value[1][3] * SubFactor04),
@@ -434,6 +434,27 @@ namespace detail
             + (this->value[0][0] * SubFactor13 - this->value[0][2] * SubFactor16 + this->value[0][3] * SubFactor17),
             - (this->value[0][0] * SubFactor14 - this->value[0][1] * SubFactor16 + this->value[0][3] * SubFactor18),
             + (this->value[0][0] * SubFactor15 - this->value[0][1] * SubFactor17 + this->value[0][2] * SubFactor18));
+*/
+        tmat4x4<T> Inverse(
+            + this->value[1][1] * SubFactor00 - this->value[1][2] * SubFactor01 + this->value[1][3] * SubFactor02,
+            - this->value[1][0] * SubFactor00 + this->value[1][2] * SubFactor03 - this->value[1][3] * SubFactor04,
+            + this->value[1][0] * SubFactor01 - this->value[1][1] * SubFactor03 + this->value[1][3] * SubFactor05,
+            - this->value[1][0] * SubFactor02 + this->value[1][1] * SubFactor04 - this->value[1][2] * SubFactor05,
+
+            - this->value[0][1] * SubFactor00 + this->value[0][2] * SubFactor01 - this->value[0][3] * SubFactor02,
+            + this->value[0][0] * SubFactor00 - this->value[0][2] * SubFactor03 + this->value[0][3] * SubFactor04,
+            - this->value[0][0] * SubFactor01 + this->value[0][1] * SubFactor03 - this->value[0][3] * SubFactor05,
+            + this->value[0][0] * SubFactor02 - this->value[0][1] * SubFactor04 + this->value[0][2] * SubFactor05,
+
+            + this->value[0][1] * SubFactor06 - this->value[0][2] * SubFactor07 + this->value[0][3] * SubFactor08,
+            - this->value[0][0] * SubFactor06 + this->value[0][2] * SubFactor09 - this->value[0][3] * SubFactor10,
+            + this->value[0][0] * SubFactor11 - this->value[0][1] * SubFactor09 + this->value[0][3] * SubFactor12,
+            - this->value[0][0] * SubFactor08 + this->value[0][1] * SubFactor10 - this->value[0][2] * SubFactor12,
+
+            - this->value[0][1] * SubFactor13 + this->value[0][2] * SubFactor14 - this->value[0][3] * SubFactor15,
+            + this->value[0][0] * SubFactor13 - this->value[0][2] * SubFactor16 + this->value[0][3] * SubFactor17,
+            - this->value[0][0] * SubFactor14 + this->value[0][1] * SubFactor16 - this->value[0][3] * SubFactor18,
+            + this->value[0][0] * SubFactor15 - this->value[0][1] * SubFactor17 + this->value[0][2] * SubFactor18);
 
         value_type Determinant = 
 			+ this->value[0][0] * Inverse[0][0]