From 243ba096e8f087ca199ec167c950deb49aec4c45 Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Thu, 4 Feb 2021 15:04:29 +0100
Subject: [PATCH] sphere:: fixmatrix is now used so that numeric errors do not
 accumulate

---
 sphere.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/sphere.cpp b/sphere.cpp
index 67d31abd..3a6b1c1c 100644
--- a/sphere.cpp
+++ b/sphere.cpp
@@ -181,6 +181,8 @@ struct hrmap_spherical : hrmap_standard {
       transmatrix T = inverse_shift(gmatrix0[c1], gmatrix0[c2]);
       if(elliptic && T[LDIM][LDIM] < 0)
         T = centralsym * T;
+
+      fixmatrix(T);
       return T;
       }
     else {