From bbabd47ad02b53e9257704235cf8c67383efcf79 Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Tue, 30 Jul 2019 13:02:24 +0200
Subject: [PATCH] solv:: fixmatrix

---
 hyperpoint.cpp |  2 +-
 sol.cpp        | 10 ++++++++++
 system.cpp     |  2 ++
 3 files changed, 13 insertions(+), 1 deletion(-)

diff --git a/hyperpoint.cpp b/hyperpoint.cpp
index e5f158b4..a56a2a7a 100644
--- a/hyperpoint.cpp
+++ b/hyperpoint.cpp
@@ -513,7 +513,7 @@ transmatrix rgpushxto0(const hyperpoint& H) {
 // (without using this, imprecision could accumulate)
 
 void fixmatrix(transmatrix& T) {
-  if(sol) ;
+  if(sol) ; // fixed inside solmul
   else if(euclid) {
     for(int x=0; x<GDIM; x++) for(int y=0; y<=x; y++) {
       ld dp = 0;
diff --git a/sol.cpp b/sol.cpp
index 5ea7dc83..23f819ad 100644
--- a/sol.cpp
+++ b/sol.cpp
@@ -315,6 +315,16 @@ namespace solv {
     transmatrix push_back = eupush( tC0(inverse(V)) );
     transmatrix space_to_view = V * push_back;
     
+    for(int x=0; x<3; x++) for(int y=0; y<=x; y++) {
+      auto& T = space_to_view;
+      ld dp = 0;
+      for(int z=0; z<3; z++) dp += T[z][x] * T[z][y];
+      
+      if(y == x) dp = 1 - sqrt(1/dp);
+      
+      for(int z=0; z<3; z++) T[z][x] -= dp * T[z][y];
+      }
+    
     transmatrix view_to_space = inverse(space_to_view);
     using namespace hyperpoint_vec;
     hyperpoint shift = /* inverse(V) * T * V * C0; */ view_to_space * T * C0;
diff --git a/system.cpp b/system.cpp
index 829e3f59..ebf56055 100644
--- a/system.cpp
+++ b/system.cpp
@@ -72,6 +72,8 @@ void welcomeMessage() {
     addMessage(XLAT("Good luck in the elliptic plane!"));
   else if(sphere)
     addMessage(XLAT("Welcome to Spherogue!"));
+  else if(sol)
+    addMessage(XLAT("Welcome to SolvRogue!"));
   else if(PURE && geometry == gNormal && !cheater)
     addMessage(XLAT("Welcome to the Heptagonal Mode!"));
   else if(cheater || autocheat)