summaryrefslogtreecommitdiff
path: root/fourier_integrator.cpp
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2025-04-03 17:54:00 -0300
committerJaron Kent-Dobias <jaron@kent-dobias.com>2025-04-03 17:54:00 -0300
commita446c8ae813ae5e1c33d377dabd85d3d362b1878 (patch)
tree757a8f74f1c8cb2be42e03a337a6bf26d4ae5f2b /fourier_integrator.cpp
parent0120bef06d4e537703ee7f0852ad31eceb43abf3 (diff)
downloadcode-a446c8ae813ae5e1c33d377dabd85d3d362b1878.tar.gz
code-a446c8ae813ae5e1c33d377dabd85d3d362b1878.tar.bz2
code-a446c8ae813ae5e1c33d377dabd85d3d362b1878.zip
Removed some more duplicate code.
Diffstat (limited to 'fourier_integrator.cpp')
-rw-r--r--fourier_integrator.cpp4
1 files changed, 1 insertions, 3 deletions
diff --git a/fourier_integrator.cpp b/fourier_integrator.cpp
index 86c3b24..2fe7e49 100644
--- a/fourier_integrator.cpp
+++ b/fourier_integrator.cpp
@@ -103,9 +103,7 @@ int main(int argc, char* argv[]) {
Ct = fft.fourier(C);
Rt = fft.fourier(R);
- auto [RddfCt, dfCt] = RddfCtdfCt(fft, C, R, p, s, λ);
-
- z = ((Γ₀ * std::conj(Rt[0]) + pow(y, 2) * (RddfCt[0] * Ct[0] + dfCt[0] * std::conj(Rt[0]))) / Ct[0]).real();
+ z = estimateZ(fft, C, Ct, R, Rt, p, s, λ, y);
}
while (y += Δy, y <= yₘₐₓ) {