Fill Γ under any conditions

1 files changed, 4 insertions, 1 deletions

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);