From 20546133c3d14fc1024bd16febce0eaf5f19411e Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Tue, 30 May 2023 18:38:58 +0200
Subject: [PATCH] Berger sphere should now work correctly in subdivided

---
 reg3.cpp | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/reg3.cpp b/reg3.cpp
index 6aea3ba9..6e7f26cd 100644
--- a/reg3.cpp
+++ b/reg3.cpp
@@ -79,6 +79,16 @@ EX void make_vertices_only(vector<hyperpoint>& vo, const vector<vector<hyperpoin
   }
 EX }
 
+transmatrix sfix_rgpushxto0(hyperpoint h) {
+  if(stretch::applicable() && !fake::in()) return stretch::translate(h);
+  return rgpushxto0(h);
+  };
+
+transmatrix sfix_gpushxto0(hyperpoint h) {
+  if(stretch::applicable() && !fake::in()) return stretch::itranslate(h);
+  return gpushxto0(h);
+  };
+
 #if MAXMDIM >= 4
 void subcellshape::compute_sub() {
   hyperpoint gres = Hypc;
@@ -90,8 +100,8 @@ void subcellshape::compute_sub() {
     gres += res;
     }
   cellcenter = normalize(gres);
-  to_cellcenter = rgpushxto0(cellcenter);
-  from_cellcenter = gpushxto0(cellcenter);
+  to_cellcenter = sfix_rgpushxto0(cellcenter);
+  from_cellcenter = sfix_gpushxto0(cellcenter);
   compute_common();
   }