From 42265258f20e675c9540b22df3d48ceedb05d92f Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Wed, 7 Apr 2021 14:46:46 +0200
Subject: [PATCH] hyperpoint:: det function is now based on the largest
 coordinate, making it more precise

---
 hyperpoint.cpp | 12 +++++++-----
 1 file changed, 7 insertions(+), 5 deletions(-)

diff --git a/hyperpoint.cpp b/hyperpoint.cpp
index a4ed4378..a234d4d1 100644
--- a/hyperpoint.cpp
+++ b/hyperpoint.cpp
@@ -936,12 +936,14 @@ EX ld det(const transmatrix& T) {
     ld det = 1;
     transmatrix M = T;
     for(int a=0; a<MDIM; a++) {
+      int max_at = a;
       for(int b=a; b<MDIM; b++)
-        if(M[b][a]) {
-          if(b != a)
-            for(int c=a; c<MDIM; c++) tie(M[b][c], M[a][c]) = make_pair(-M[a][c], M[b][c]);
-          break;
-          }
+        if(abs(M[b][a]) > abs(M[max_at][a]))
+          max_at = b;
+      
+      if(max_at != a)
+        for(int c=a; c<MDIM; c++) tie(M[max_at][c], M[a][c]) = make_pair(-M[a][c], M[max_at][c]);
+
       if(!M[a][a]) return 0;
       for(int b=a+1; b<MDIM; b++) {
         ld co = -M[b][a] / M[a][a];