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