diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-20 20:28:14 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-20 20:28:14 -0300 |
| commit | 88d5c33164724434a41b9c5064aa786a3cb51029 (patch) | |
| tree | ae9ae84ab6eea4232a5f0da7753370a624b67b09 /log-fourier.cpp | |
| parent | a85abcca348e8e83b476eca92eef7c27631fd192 (diff) | |
| download | code-88d5c33164724434a41b9c5064aa786a3cb51029.tar.gz code-88d5c33164724434a41b9c5064aa786a3cb51029.tar.bz2 code-88d5c33164724434a41b9c5064aa786a3cb51029.zip | |
Change convention to fixed Γ₀, and changed behavior on NaN
Diffstat (limited to 'log-fourier.cpp')
| -rw-r--r-- | log-fourier.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/log-fourier.cpp b/log-fourier.cpp index f7d6f4b..47d5c5c 100644 --- a/log-fourier.cpp +++ b/log-fourier.cpp @@ -173,7 +173,7 @@ std::tuple<std::vector<Complex>, std::vector<Complex>> RddfCtdfCt(LogarithmicFou Real estimateZ(LogarithmicFourierTransform& fft, const std::vector<Real>& C, const std::vector<Complex>& Ct, const std::vector<Real>& R, const std::vector<Complex>& Rt, unsigned p, unsigned s, Real λ, Real τ₀, Real β) { auto [RddfCt, dfCt] = RddfCtdfCt(fft, C, R, p, s, λ); - Real Γ₀ = 1 + τ₀; + Real Γ₀ = 1.0; return ((2 * Γ₀ * std::conj(Rt[0]) + pow(β, 2) * (RddfCt[0] * Ct[0] + dfCt[0] * std::conj(Rt[0]))) / Ct[0]).real(); } |