From ebd38828e38760954affea2c5c9073497c2926a7 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 19 Apr 2025 14:54:41 -0300
Subject: Playing with step size contral

---
 log-fourier_integrator.cpp | 18 +++++++++++++++++-
 1 file changed, 17 insertions(+), 1 deletion(-)

(limited to 'log-fourier_integrator.cpp')

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index ffbce1d..eb72b49 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -85,6 +85,8 @@ int main(int argc, char* argv[]) {
   Real β = 0;
   while (β < βₘₐₓ) {
     Real ΔC = 100;
+    Real ΔC₀ = 100;
+    unsigned it = 0;
     while (ΔC > ε) {
       std::vector<Real> RddfC(N);
       for (unsigned n = 0; n < N; n++) {
@@ -119,6 +121,20 @@ int main(int argc, char* argv[]) {
 
       μ *= Cₜ[0];
 
+      if (ΔC < ΔC₀) {
+        ΔC₀ = ΔC;
+        it = 0;
+        γ = std::min(1.01 * γ, 1.0);
+      } else {
+        it++;
+      }
+
+      if (it > 1000) {
+        std::cerr << "Decreased!" << std::endl;
+        γ = std::max(0.5 * γ, 1e-6);
+        it = 0;
+      }
+
 //    std::cerr << ΔC << std::endl;
     }
 
@@ -138,7 +154,7 @@ int main(int argc, char* argv[]) {
     }
     E *= β;
 
-    std::cerr << β << " " << μ << " " << Ĉₜ[0].real() << " " << E << std::endl;
+    std::cerr << β << " " << μ << " " << Ĉₜ[0].real() << " " << E << " " << γ << std::endl;
     β += Δβ;
   }
 
-- 
cgit v1.2.3-70-g09d2