From 3caf700868eba2fb6f1e09e9d9329c0597b44692 Mon Sep 17 00:00:00 2001 From: Zeno Rogue Date: Wed, 25 Mar 2020 13:43:35 +0100 Subject: [PATCH] solv:: fixed another potential problem --- nonisotropic.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/nonisotropic.cpp b/nonisotropic.cpp index 35496472..bb9a4925 100644 --- a/nonisotropic.cpp +++ b/nonisotropic.cpp @@ -224,7 +224,7 @@ EX namespace sn { EX hyperpoint table_to_azeq(hyperpoint x) { // Poincare to azimuthal equidistant ld hr = sqhypot_d(3, x); - if(hr == 0) return point3(0,0,0); + if(hr < 1e-5) return x * 2; if(hr >= 1) return x * 60; ld hz = (1 + hr) / (1 - hr); ld d = (hz+1) * acosh(hz) / sinh(acosh(hz));