diff options
| -rw-r--r-- | log-fourier_integrator.cpp | 5 |
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); |