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_integrator.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_integrator.cpp')
| -rw-r--r-- | fourier_integrator.cpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/fourier_integrator.cpp b/fourier_integrator.cpp index d83fb9d..24a7f1d 100644 --- a/fourier_integrator.cpp +++ b/fourier_integrator.cpp @@ -17,7 +17,7 @@ int main(int argc, char* argv[]) { unsigned maxIterations = 1000; Real ε = 1e-14; - Real γ = 0; + Real γ = 1; bool loadData = false; @@ -68,8 +68,8 @@ int main(int argc, char* argv[]) { Real Δτ = τₘₐₓ / M_PI / n; Real Δω = M_PI / τₘₐₓ; - Real z = (-1+sqrt(1+2*τ₀)) / (2 * τ₀); - Real Γ₀ = 1; + Real z = 0.5; + Real Γ₀ = 1 + τ₀ / 2; std::vector<Real> C(2 * n); std::vector<Real> R(2 * n); @@ -103,7 +103,7 @@ int main(int argc, char* argv[]) { Ct = fft.fourier(C); Rt = fft.fourier(R); - z = estimateZ(fft, C, Ct, R, Rt, p, s, λ, y); + z = estimateZ(fft, C, Ct, R, Rt, p, s, λ, τ₀, y); } while (y += Δy, y <= yₘₐₓ) { |