summaryrefslogtreecommitdiff
path: root/fourier_integrator.cpp
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2025-04-05 11:25:41 -0300
committerJaron Kent-Dobias <jaron@kent-dobias.com>2025-04-05 11:25:41 -0300
commit685c562740aed1c62ec5ae27d5c138ce57e2dcee (patch)
treeb3823fcebc34b030d57f097f471119de0e393d2e /fourier_integrator.cpp
parent96875ddafbf1ff821a4a1fd628c2e4758096700c (diff)
downloadcode-685c562740aed1c62ec5ae27d5c138ce57e2dcee.tar.gz
code-685c562740aed1c62ec5ae27d5c138ce57e2dcee.tar.bz2
code-685c562740aed1c62ec5ae27d5c138ce57e2dcee.zip
Rescale τₘₐₓ with τ₀ to make approximations more consistent accross parameters
Diffstat (limited to 'fourier_integrator.cpp')
-rw-r--r--fourier_integrator.cpp4
1 files changed, 2 insertions, 2 deletions
diff --git a/fourier_integrator.cpp b/fourier_integrator.cpp
index 7669ca0..69c8b03 100644
--- a/fourier_integrator.cpp
+++ b/fourier_integrator.cpp
@@ -65,8 +65,8 @@ int main(int argc, char* argv[]) {
unsigned n = pow(2, log2n);
- Real Δτ = τₘₐₓ / M_PI / n;
- Real Δω = M_PI / τₘₐₓ;
+ Real Δτ = (1 + τ₀ / 2) * τₘₐₓ / M_PI / n;
+ Real Δω = M_PI / ((1 + τ₀ / 2) * τₘₐₓ);
Real z = 0.5;
Real Γ₀ = 1 + τ₀ / 2;