diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-05-08 17:32:18 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-05-08 17:32:18 -0300 |
| commit | e151b804071d69a41beef04a73c12c42b12bd775 (patch) | |
| tree | 01d989abdc354ab28e9a4a551052039972d2cc42 /log-fourier.cpp | |
| parent | 18c33edc2fdf6abc9f8f36ea67b256d4a885493a (diff) | |
| download | code-e151b804071d69a41beef04a73c12c42b12bd775.tar.gz code-e151b804071d69a41beef04a73c12c42b12bd775.tar.bz2 code-e151b804071d69a41beef04a73c12c42b12bd775.zip | |
Revert "Made log-Fourier padding symmetric, and began writing regular integrator"
This reverts commit 5fd9866479ec50051d2c9eeb4e217e9376e6f9b4.
Diffstat (limited to 'log-fourier.cpp')
| -rw-r--r-- | log-fourier.cpp | 8 |
1 files changed, 0 insertions, 8 deletions
diff --git a/log-fourier.cpp b/log-fourier.cpp index 9d1f2cb..9ae70be 100644 --- a/log-fourier.cpp +++ b/log-fourier.cpp @@ -64,8 +64,6 @@ std::vector<Complex> LogarithmicFourierTransform::fourier(const std::vector<Real for (unsigned n = 0; n < pad*N; n++) { if (n < N) { a[n] = c[n] * exp((1 - k) * τ(n)); - } else if (n >= (pad - 1) * N) { - a[n] = c[pad*N-n-1] * exp((1 - k) * τ(pad*N-n-1)); } else { a[n] = 0; } @@ -94,12 +92,6 @@ std::vector<Real> LogarithmicFourierTransform::inverse(const std::vector<Complex } else { a[n] = ĉ[n] * std::exp((1 - k) * ω(n)); } - } else if (n >= (pad - 1) * N) { - if (σ < 0) { - a[n] = ĉ[pad*N-n-1] * std::exp((1 - k) * ω(pad*N-n-1)); - } else { - a[n] = std::conj(ĉ[pad*N-n-1]) * std::exp((1 - k) * ω(pad*N-n-1)); - } } else { a[n] = 0; } |