diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-04 16:31:11 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-04 16:31:11 -0300 |
| commit | 61a2a8aaf7a4c02f75ac85e42841285aba7ef5a1 (patch) | |
| tree | 2131f41358582b2fb21cf49d3b95a522f78ea18a /fourier.cpp | |
| parent | 82c348698ccaa41b2f455ba017113cab03590016 (diff) | |
| download | code-61a2a8aaf7a4c02f75ac85e42841285aba7ef5a1.tar.gz code-61a2a8aaf7a4c02f75ac85e42841285aba7ef5a1.tar.bz2 code-61a2a8aaf7a4c02f75ac85e42841285aba7ef5a1.zip | |
Adopted new scaling of Γ₀ to keep z and y moderate
Diffstat (limited to 'fourier.cpp')
| -rw-r--r-- | fourier.cpp | 5 |
1 files changed, 3 insertions, 2 deletions
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<Complex>, std::vector<Complex>> RddfCtdfCt(FourierTransfo return {RddfCt, dfCt}; } -Real estimateZ(FourierTransform& 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 y) { +Real estimateZ(FourierTransform& 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 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(); } |