From 1c8f0e7e37ac6f1a8aacca11d28c96bb4cccf52a Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Sun, 19 Aug 2018 22:54:11 +0200
Subject: [PATCH] inverse_error as a separate function (helps with debugging)

---
 hyperpoint.cpp | 16 +++++++++-------
 1 file changed, 9 insertions(+), 7 deletions(-)

diff --git a/hyperpoint.cpp b/hyperpoint.cpp
index 183a98e9..57b52c1a 100644
--- a/hyperpoint.cpp
+++ b/hyperpoint.cpp
@@ -471,18 +471,20 @@ ld det(const transmatrix& T) {
     det -= T[0][i] * T[1][(i+2)%3] * T[2][(i+1)%3];
   return det;
   }
-  
+
+void inverse_error(const transmatrix& T) {
+  printf("Warning: inverting a singular matrix\n");
+  display(T);
+  }
+
 transmatrix inverse(const transmatrix& T) {
   profile_start(7);
   
   ld d = det(T);
   transmatrix T2;
-  if(d == 0) { 
-    printf("Warning: inverting a singular matrix\n");
-    display(T);
-    display(T2);
-    T2 = Id; 
-    return T2; 
+  if(d == 0) {
+    inverse_error(T); 
+    return Id;
     }
   
   for(int i=0; i<3; i++)