From 077b51f1eb359aae744aae655c4199e8dd458cb0 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Sun, 18 May 2025 18:35:11 -0300 Subject: Simplified initial C --- log-fourier_integrator.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'log-fourier_integrator.cpp') diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index f553569..30354a6 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -114,7 +114,7 @@ int main(int argc, char* argv[]) { if (τ₀ == 2) { Cₜ₋₁[n] = Γ₀ * std::exp(-fft.t(n) / 2) * (1 + fft.t(n) / 2); } else { - Cₜ₋₁[n] = Γ₀ * (std::exp(-μ₀ * fft.t(n)) - μ₀ * τ₀ * std::exp(-fft.t(n) / τ₀)) / (μ₀ - pow(μ₀, 3) * pow(τ₀, 2)); + Cₜ₋₁[n] = Γ₀ * (std::exp(-μ₀ * fft.t(n)) / μ₀ - τ₀ * std::exp(-fft.t(n) / τ₀)) / (1 - pow(μ₀ * τ₀, 2)); } } else { Cₜ₋₁[n] = Γ₀ * std::exp(-μ₀ * fft.t(n)) / μ₀; -- cgit v1.2.3-70-g09d2