From 4525068dfc6374813b9b55c0e968ac7b77425ace Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 14 May 2025 09:53:21 -0300
Subject: Fill Γ under any conditions
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 log-fourier_integrator.cpp | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index 8a2fe5c..c2a6568 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -96,6 +96,10 @@ int main(int argc, char* argv[]) {
   std::vector<Complex> Ȓₜ₋₁(N);
   std::vector<Real> Γ(N);
 
+  for (unsigned n = 0; n < N; n++) { 
+    Γ[n] = Γ₀ / (1 + std::pow(τ₀ * fft.ν(n), 2));
+  }
+
   if (!loadData) {
     /* Start from the exact solution for β = 0 */
     for (unsigned n = 0; n < N; n++) {
@@ -112,7 +116,6 @@ int main(int argc, char* argv[]) {
 
       Ĉₜ₋₁[n] = 2 * Γ₀ / (pow(μ₀, 2) + pow(fft.ν(n), 2)) / (1 + pow(τ₀ * fft.ν(n), 2));
       Ȓₜ₋₁[n] = (Real)1.0 / (μ₀ + II * fft.ν(n));
-      Γ[n] = Γ₀ / (1 + std::pow(τ₀ * fft.ν(n), 2));
     }
   } else {
     logFourierLoad(Cₜ₋₁, Rₜ₋₁, Ĉₜ₋₁, Ȓₜ₋₁, p, s, λ, τ₀, β₀, log2n, Δτ, logShift);
-- 
cgit v1.2.3-70-g09d2