From c62be882352bbb6c164d0c6406784357c046a895 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Thu, 3 Apr 2025 18:33:16 -0300 Subject: Don't change γ every iteration MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- fourier_integrator.cpp | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) (limited to 'fourier_integrator.cpp') diff --git a/fourier_integrator.cpp b/fourier_integrator.cpp index 9e41b5f..d83fb9d 100644 --- a/fourier_integrator.cpp +++ b/fourier_integrator.cpp @@ -146,14 +146,16 @@ int main(int argc, char* argv[]) { z *= Cnew[0]; - if (ΔC < ΔCprev) { - γ = std::min(1.1 * γ, 1.0); - } else { - γ /= 2; + if (it % maxIterations == 0) { + if (ΔC < ΔCprev) { + γ = std::min(1.1 * γ, 1.0); + } else { + γ /= 2; + } + + ΔCprev = ΔC; } - ΔCprev = ΔC; - std::cerr << it << " " << p << " " << s << " " << τ₀ << " " << y << " " << sqrt(2 * ΔC / C.size()) << " " << γ << std::endl; } -- cgit v1.2.3-70-g09d2