diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-27 09:44:11 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-27 09:44:11 -0300 |
| commit | bbe904371abfa4c8b2cc84855c7f4013dffc9322 (patch) | |
| tree | a8812136dad52175e4566ee4afb1ce32fce3a509 /log-fourier_integrator.cpp | |
| parent | e30dce3c53b7ce12f24f6e777032ac32baa0f2d6 (diff) | |
| download | code-bbe904371abfa4c8b2cc84855c7f4013dffc9322.tar.gz code-bbe904371abfa4c8b2cc84855c7f4013dffc9322.tar.bz2 code-bbe904371abfa4c8b2cc84855c7f4013dffc9322.zip | |
Control μ damping with γ
Diffstat (limited to 'log-fourier_integrator.cpp')
| -rw-r--r-- | log-fourier_integrator.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index c30d01d..52a5d02 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -132,7 +132,7 @@ int main(int argc, char* argv[]) { std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); std::vector<Real> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); - μₜ *= pow(tanh(Cₜ₊₁[0]-1)+1, 0.05); + μₜ *= (γ * pow(tanh(Cₜ₊₁[0]-1)+1, 0.5) + (1 - γ));; ΔCₜ = 0; for (unsigned i = 0; i < N; i++) { |