From 442b83dbe83f1a05956f871f97129123f36f75f1 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 19 Apr 2025 17:40:10 -0300
Subject: New way to regularize the factor

---
 log-fourier_integrator.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index dc4f48f..21b672c 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -172,7 +172,7 @@ int main(int argc, char* argv[]) {
           /* We found an upper bound */
           μ₂ = μ;
         }
-        μ *= sqrt(sqrt(Cₜ[0]));
+        μ *= sqrt(sqrt(std::tanh(Cₜ[0]-1)+1));
       } else {
         /* Once the bounds are set, we can use bisection */
         if (Cₜ[0] > 1) {
-- 
cgit v1.2.3-70-g09d2