From 61a2a8aaf7a4c02f75ac85e42841285aba7ef5a1 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Fri, 4 Apr 2025 16:31:11 -0300 Subject: Adopted new scaling of Γ₀ to keep z and y moderate MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- fourier.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) (limited to 'fourier.cpp') diff --git a/fourier.cpp b/fourier.cpp index 8943521..d2f9ad8 100644 --- a/fourier.cpp +++ b/fourier.cpp @@ -85,8 +85,9 @@ std::tuple, std::vector> RddfCtdfCt(FourierTransfo return {RddfCt, dfCt}; } -Real estimateZ(FourierTransform& fft, const std::vector& C, const std::vector& Ct, const std::vector& R, const std::vector& Rt, unsigned p, unsigned s, Real λ, Real y) { +Real estimateZ(FourierTransform& fft, const std::vector& C, const std::vector& Ct, const std::vector& R, const std::vector& Rt, unsigned p, unsigned s, Real λ, Real τ₀, Real y) { auto [RddfCt, dfCt] = RddfCtdfCt(fft, C, R, p, s, λ); + Real Γ₀ = 1 + τ₀ / 2; - return ((std::conj(Rt[0]) + pow(y, 2) * (RddfCt[0] * Ct[0] + dfCt[0] * std::conj(Rt[0]))) / Ct[0]).real(); + return ((Γ₀ * std::conj(Rt[0]) + pow(y, 2) * (RddfCt[0] * Ct[0] + dfCt[0] * std::conj(Rt[0]))) / Ct[0]).real(); } -- cgit v1.2.3-70-g09d2