From 501212c09ca114648d76a6aea2a35606968e7bee Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 10 May 2025 15:11:57 -0300
Subject: Tweak to μ
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
log-fourier_integrator.cpp | 6 ++++--
1 file 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])) {
--
cgit v1.2.3-70-g09d2