From 11b1e4c4451f5db9b237d06c60b7d23d40852e24 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 10 May 2025 10:58:39 -0300
Subject: Fixed energy agian

---
 log-fourier.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'log-fourier.cpp')

diff --git a/log-fourier.cpp b/log-fourier.cpp
index b4f035b..8bd7e51 100644
--- a/log-fourier.cpp
+++ b/log-fourier.cpp
@@ -192,7 +192,7 @@ Real energy(const LogarithmicFourierTransform& fft, std::vector<Real>&  C, const
   for (unsigned n = 0; n < C.size(); n++) {
     if (C[n] > 1) n₀ = n % 2 == 0 ? n / 2 : (n + 1) / 2;
   }
-  Real E = fft.t(n₀) * df(λ, p, s, 1);
+  Real E = fft.t(2*n₀) * df(λ, p, s, 1);
   for (unsigned n = n₀; n < C.size()/2-1; n++) {
     Real R₂ₙ   = R[2*n];
     Real R₂ₙ₊₁ = R[2*n+1];
-- 
cgit v1.2.3-70-g09d2