From 11344d18210b582c547c697f86bfac97c360d6ca Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 19 Apr 2025 21:24:04 -0300
Subject: Fixed energy

---
 fourier.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/fourier.cpp b/fourier.cpp
index d637985..3821623 100644
--- a/fourier.cpp
+++ b/fourier.cpp
@@ -85,8 +85,8 @@ Real energy(const std::vector<Real>& C, const std::vector<Real>& R, unsigned p,
   Real e = 0;
 
   for (unsigned n = 0; n < C.size() / 4 -1; n++) {
-    Real h₂ₙ   =  Δτ;
-    Real h₂ₙ₊₁ = Δτ;
+    Real h₂ₙ   = M_PI * Δτ;
+    Real h₂ₙ₊₁ = M_PI * Δτ;
     Real f₂ₙ   = R[2*n]   * df(λ, p, s, C[2*n]);
     Real f₂ₙ₊₁ = R[2*n+1] * df(λ, p, s, C[2*n+1]);
     Real f₂ₙ₊₂ = R[2*n+2] * df(λ, p, s, C[2*n+2]);
-- 
cgit v1.2.3-70-g09d2