diff options
| -rw-r--r-- | log-fourier_integrator.cpp | 4 | 
1 files changed, 2 insertions, 2 deletions
| diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 28dd6e7..3049bd5 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -148,7 +148,7 @@ int main(int argc, char* argv[]) {            ΔC₀ = ΔC;          } -        std::cerr << β << " " << μ << " " << ΔC << " " << γ; +        std::cerr << β << " " << μ << " " << ΔC << " " << γ << " " << Cₜ[0];          std::cerr << "\r";        }        if (std::isnan(Cₜ[0])) { @@ -172,7 +172,7 @@ int main(int argc, char* argv[]) {            /* We found an upper bound */            μ₂ = μ;          } -        μ *= Cₜ[0]; +        μ *= sqrt(Cₜ[0]);        } else {          /* Once the bounds are set, we can use bisection */          if (Cₜ[0] > 1) { | 
