diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-19 14:38:03 -0300 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-19 14:38:03 -0300 |
commit | 2aaf47af7a5a49033f54d6619ed9e4a448d69a0a (patch) | |
tree | 8726d52fefc171aa2065bfe34b0ef7d381f98f60 | |
parent | ae23cdf557611d213b6aa03597ac4f8dd139f98b (diff) | |
download | code-2aaf47af7a5a49033f54d6619ed9e4a448d69a0a.tar.gz code-2aaf47af7a5a49033f54d6619ed9e4a448d69a0a.tar.bz2 code-2aaf47af7a5a49033f54d6619ed9e4a448d69a0a.zip |
Changed notation
-rw-r--r-- | log-fourier_integrator.cpp | 14 |
1 files changed, 0 insertions, 14 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index bf83da5..ffbce1d 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -87,13 +87,10 @@ int main(int argc, char* argv[]) { Real ΔC = 100; while (ΔC > ε) { std::vector<Real> RddfC(N); -// std::vector<Real> dfC(N); for (unsigned n = 0; n < N; n++) { RddfC[n] = Rₜ[n] * ddf(λ, p, s, Cₜ[n]); -// dfC[n] = df(λ, p, s, Cₜ[n]); } std::vector<Complex> RddfCt = fft.fourier(RddfC, false); -// std::vector<Complex> dfCt = fft.fourier(dfC, true); std::vector<Complex> Ȓₜ₊₁(N); std::vector<Complex> Ĉₜ₊₁(N); @@ -104,17 +101,6 @@ int main(int argc, char* argv[]) { } std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); - /* - for (unsigned n = 0; n < N; n++) { - RddfC[n] = Rₜ₊₁[n] * ddf(λ, p, s, Cₜ[n]); - } - RddfCt = fft.fourier(RddfC, false); - - for (unsigned n = 0; n < N; 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> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); ΔC = 0; |