From 00ca38555f905506967180ede7494169102c60cd Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 19 Apr 2025 15:38:14 -0300
Subject: More feedback

---
 log-fourier_integrator.cpp | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

(limited to 'log-fourier_integrator.cpp')

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index eb72b49..5dc64ec 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -4,15 +4,18 @@
 #include <iostream>
 
 int main(int argc, char* argv[]) {
+  /* Model parameters */
   unsigned p = 2;
   unsigned s = 2;
   Real λ = 0.5;
   Real τ₀ = 0;
 
+  /* Log-Fourier parameters */
   unsigned log2n = 8;
   Real Δτ = 0.1;
   Real k = 0.1;
 
+  /* Iteration parameters */
   Real ε = 1e-13;
   Real γ = 1;
   Real βₘₐₓ = 0.7;
@@ -69,7 +72,7 @@ int main(int argc, char* argv[]) {
   std::vector<Complex> Ĉₜ(N);
   std::vector<Complex> Ȓₜ(N);
 
-  // start from the exact solution for β = 0
+  /* Start from the exact solution for β = 0 */
   for (unsigned n = 0; n < N; n++) {
     if (τ₀ > 0) {
       Cₜ[n] = Γ₀ * (exp(-μ * fft.t(n)) - μ * τ₀ * exp(-fft.t(n) / τ₀)) / (μ - pow(μ, 3) * pow(τ₀, 2));
@@ -135,7 +138,8 @@ int main(int argc, char* argv[]) {
         it = 0;
       }
 
-//    std::cerr << ΔC << std::endl;
+      std::cerr << β << " " << μ << " " << ΔC << " " << γ;
+      std::cerr << "\r";
     }
 
     /* Integrate the energy using Simpson's rule */
-- 
cgit v1.2.3-70-g09d2