From 752c9294955123ff94eccfadb464cf3819533f7d Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Sat, 19 Apr 2025 18:22:58 -0300 Subject: Tweaked iteration --- log-fourier_integrator.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index cc42cba..add6ead 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -115,7 +115,7 @@ 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 * Γ₀ * std::conj(Ȓₜ₊₁[n]) / (1 + pow(τ₀ * fft.ν(n), 2)) + pow(β, 2) * (RddfCt[n] * Ĉₜ[n] + dfCt[n] * std::conj(Ȓₜ₊₁[n]))) / (μ + 1i * fft.ν(n)); } std::vector Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); -- cgit v1.2.3-70-g09d2