diff options
Diffstat (limited to 'log-fourier_integrator.cpp')
| -rw-r--r-- | log-fourier_integrator.cpp | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 627f49e..0dcf49d 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -241,7 +241,11 @@ int main(int argc, char* argv[]) { logFourierSave(Cₜ, Rₜ, Ĉₜ, Ȓₜ, p, s, λ, τ₀, β, log2n, Δτ, logShift); - β += Δβ; + if (Ĉₜ[0].real() / Ĉₜ₋₁[0].real() > 2) { + Δβ *= 0.5; + } + + β = std::round(1e6 * (β + Δβ)) / 1e6; Cₜ₋₁ = Cₜ; Rₜ₋₁ = Rₜ; Ĉₜ₋₁ = Ĉₜ; |