diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-19 18:12:00 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-19 18:12:00 -0300 |
| commit | 84e41519a0d98c450c0d44194ca7c722f9d2d997 (patch) | |
| tree | 1559ac770560930066a3065bb173d9488eac3777 | |
| parent | 0e6ed4e2dc16d341c7107adfe7b6f3a619772704 (diff) | |
| download | code-84e41519a0d98c450c0d44194ca7c722f9d2d997.tar.gz code-84e41519a0d98c450c0d44194ca7c722f9d2d997.tar.bz2 code-84e41519a0d98c450c0d44194ca7c722f9d2d997.zip | |
Back to other iterator
| -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 c9ff482..bc53f99 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -114,8 +114,8 @@ int main(int argc, char* argv[]) { for (unsigned n = 0; n < N; n++) { Ȓₜ₊₁[n] = (1.0 + pow(β, 2) * RddfCt[n] * Ȓₜ[n]) / (μ + 1i * fft.ν(n)); - Ĉₜ₊₁[n] = - 2 * Γ₀ * Ȓₜ₊₁[n].imag() / (1 + pow(τ₀ * fft.ν(n), 2)) / fft.ν(n); -// Ĉₜ₊₁[n] = (2 * Γ₀ * std::conj(Ȓₜ[n]) / (1 + pow(τ₀ * fft.ν(n), 2)) + pow(β, 2) * (RddfCt[n] * Ĉₜ[n] + dfCt[n] * std::conj(Ȓₜ[n]))) / (μ + 1i * fft.ν(n)); +// Ĉₜ₊₁[n] = - 2 * Γ₀ * Ȓₜ₊₁[n].imag() / (1 + pow(τ₀ * fft.ν(n), 2)) / fft.ν(n); + Ĉₜ₊₁[n] = (2 * Γ₀ * std::conj(Ȓₜ[n]) / (1 + pow(τ₀ * fft.ν(n), 2)) + pow(β, 2) * (RddfCt[n] * Ĉₜ[n] + dfCt[n] * std::conj(Ȓₜ[n]))) / (μ + 1i * fft.ν(n)); } std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); |