From b9e5a1f2448aeaca8aa807f3baccbe66f5dc49d2 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 7 May 2025 22:47:06 -0300
Subject: Average first two values of C to fix μ
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 log-fourier_integrator.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index 55f358a..017def7 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -140,7 +140,8 @@ 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++) {
-- 
cgit v1.2.3-70-g09d2