summaryrefslogtreecommitdiff
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
parent96875ddafbf1ff821a4a1fd628c2e4758096700c (diff)
downloadcode-685c562740aed1c62ec5ae27d5c138ce57e2dcee.tar.gz
code-685c562740aed1c62ec5ae27d5c138ce57e2dcee.tar.bz2
code-685c562740aed1c62ec5ae27d5c138ce57e2dcee.zip
Rescale τₘₐₓ with τ₀ to make approximations more consistent accross parameters
-rw-r--r--fourier_integrator.cpp4
-rw-r--r--get_energy.cpp4
2 files changed, 4 insertions, 4 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;
diff --git a/get_energy.cpp b/get_energy.cpp
index 8c610c5..c33c04e 100644
--- a/get_energy.cpp
+++ b/get_energy.cpp
@@ -50,8 +50,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 y = y₀;