diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-20 00:58:53 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-20 00:58:53 -0300 |
| commit | 399301c9f9c5bb801af389d50ca7dc2adbd562b3 (patch) | |
| tree | 75dcffbe4b1886bb73e5c543131bf30b632dabbf /log-fourier_integrator.cpp | |
| parent | 90c013f8c5c8663f0f5c6aad8c01574db3e4668f (diff) | |
| download | code-399301c9f9c5bb801af389d50ca7dc2adbd562b3.tar.gz code-399301c9f9c5bb801af389d50ca7dc2adbd562b3.tar.bz2 code-399301c9f9c5bb801af389d50ca7dc2adbd562b3.zip | |
Diffstat (limited to 'log-fourier_integrator.cpp')
| -rw-r--r-- | log-fourier_integrator.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 51f5326..b21f275 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -76,10 +76,10 @@ int main(int argc, char* argv[]) { /* Start from the exact solution for β = 0 */ for (unsigned n = 0; n < N; n++) { - if (τ₀ > 0) { + if (τ₀ != 1) { Cₜ₋₁[n] = Γ₀ * (exp(-μ * fft.t(n)) - μ * τ₀ * exp(-fft.t(n) / τ₀)) / (μ - pow(μ, 3) * pow(τ₀, 2)); } else { - Cₜ₋₁[n] = Γ₀ * exp(-μ * fft.t(n)) / μ; + Cₜ₋₁[n] = Γ₀ * exp(-fft.t(n)) * (1 + fft.t(n)); } Rₜ₋₁[n] = exp(-μ * fft.t(n)); |