diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-05-10 15:11:57 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-05-10 15:11:57 -0300 |
| commit | 501212c09ca114648d76a6aea2a35606968e7bee (patch) | |
| tree | f369ecc28d22e6da42b464e735a416bc2010a80b /log-fourier_integrator.cpp | |
| parent | d93dc5dba5d9e00298038c3dac3b51aa554d2ca8 (diff) | |
| download | code-501212c09ca114648d76a6aea2a35606968e7bee.tar.gz code-501212c09ca114648d76a6aea2a35606968e7bee.tar.bz2 code-501212c09ca114648d76a6aea2a35606968e7bee.zip | |
Tweak to μ
Diffstat (limited to 'log-fourier_integrator.cpp')
| -rw-r--r-- | log-fourier_integrator.cpp | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 9db21f0..04dccfd 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -145,7 +145,9 @@ int main(int argc, char* argv[]) { std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); std::vector<Real> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); - μₜ *= pow(tanh(Cₜ₊₁[0]-1)+1, x); + Real C₀ = (Cₜ₊₁[0] + Cₜ₊₁[1]) / 2; + + μₜ *= pow(tanh(C₀-1)+1, x); ΔCₜ = 0; for (unsigned i = 0; i < N; i++) { @@ -175,7 +177,7 @@ int main(int argc, char* argv[]) { } std::cerr << "\x1b[2K" << "\r"; - std::cerr << β << " " << μₜ << " " << ΔCₜ << " " << γ << " " << Cₜ[0]; + std::cerr << β << " " << μₜ << " " << ΔCₜ << " " << γ << " " << C₀; } if (std::isnan(Cₜ[0])) { |