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(-)
(limited to 'log-fourier_integrator.cpp')
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