Start from β = 0 again

1 files changed, 2 insertions, 3 deletions

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++) {