warning when inverting a singular matrix

2018-01-13 19:21:08 +01:00 · 2018-01-13 19:21:08 +01:00 · 4de6e9d3e4
parent 0f13eae87d
commit 4de6e9d3e4
1 changed files with 7 additions and 1 deletions
--- a/hyperpoint.cpp
+++ b/hyperpoint.cpp
@ -384,7 +384,13 @@ transmatrix inverse(const transmatrix& T) {
  
  ld d = det(T);
  transmatrix T2;
-  if(d == 0) return T2;
+  if(d == 0) { 
+    printf("Warning: inverting a singular matrix\n");
+    display(T);
+    display(T2);
+    T2 = Id; 
+    return T2; 
+    }
  
  for(int i=0; i<3; i++) 
  for(int j=0; j<3; j++)