From 4ff08f3f3e6b14e41dc389fe341e626ee4318da6 Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Thu, 7 Oct 2021 12:41:29 +0200
Subject: [PATCH] deparabolic13 in hyperpoint

---
 hyperpoint.cpp | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/hyperpoint.cpp b/hyperpoint.cpp
index b284a89f..788b8802 100644
--- a/hyperpoint.cpp
+++ b/hyperpoint.cpp
@@ -753,6 +753,16 @@ EX hyperpoint deparabolic10(const hyperpoint h) {
   return point31(x, h[1] * exp(x), h[2] * exp(x));
   }
 
+/** this looks correct */
+EX hyperpoint deparabolic13(hyperpoint h) {
+  h /= (1 + h[3]);
+  hyperpoint one = point3(1,0,0);
+  h -= one;
+  h /= sqhypot_d(3, h);
+  h[0] += .5;
+  return point3(log(2) + log(-h[0]), h[1] * 2, h[2] * 2);
+  }
+
 EX transmatrix spintoc(const hyperpoint& H, int t, int f) {
   transmatrix T = Id;
   ld R = hypot(H[f], H[t]);