From d78a80777e53339aae3de04caf93aa1307ee69c7 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 19 Apr 2025 16:57:25 -0300
Subject: Start from β = 0 again
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 log-fourier_integrator.cpp | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

(limited to 'log-fourier_integrator.cpp')

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index 110b9e4..ebd161b 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -87,7 +87,6 @@ int main(int argc, char* argv[]) {
 
   Real β = 0;
   while (β < βₘₐₓ) {
-    β += Δβ;
     Real μ₁ = 0;
     Real μ₂ = 0;
     while (true) {
@@ -160,11 +159,10 @@ int main(int argc, char* argv[]) {
         /* Once the bounds are set, we can use bisection */
         if (Cₜ[0] > 1) {
           μ₁ = μ;
-          μ = μ₁ + (μ₂ - μ₁) / 2;
         } else {
           μ₂ = μ;
-          μ = μ₁ + (μ₂ - μ₁) / 2;
         }
+        μ = (μ₁ + μ₂) / 2;
       }
     }
 
@@ -185,6 +183,7 @@ int main(int argc, char* argv[]) {
     E *= β;
 
     std::cerr << β << " " << μ << " " << Ĉₜ[0].real() << " " << E << " " << γ << std::endl;
+    β += Δβ;
   }
 
   for (unsigned i = 0; i < N; i++) {
-- 
cgit v1.2.3-70-g09d2