From daa25bf89cb7adcd7136abe93e1d81e1ede6ce27 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 19 Apr 2025 17:37:39 -0300
Subject: More sqrts

---
 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 3049bd5..dc4f48f 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(Cₜ[0]);
+        μ *= sqrt(sqrt(Cₜ[0]));
       } else {
         /* Once the bounds are set, we can use bisection */
         if (Cₜ[0] > 1) {
-- 
cgit v1.2.3-70-g09d2