From 48e3caa0152e0606ce10157b9bd61431117baae6 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sun, 20 Apr 2025 00:32:46 -0300
Subject: Fixed serious mistake in log inverse transform

---
 log-fourier.cpp | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

(limited to 'log-fourier.cpp')

diff --git a/log-fourier.cpp b/log-fourier.cpp
index 3ffa3c6..7461a70 100644
--- a/log-fourier.cpp
+++ b/log-fourier.cpp
@@ -1,4 +1,5 @@
 #include "log-fourier.hpp"
+#include <complex>
 
 LogarithmicFourierTransform::LogarithmicFourierTransform(unsigned N, Real k, Real Δτ, unsigned pad) : N(N), pad(pad), k(k), Δτ(Δτ) {
   τₛ = -0.5 * N;
@@ -83,7 +84,11 @@ std::vector<Real> LogarithmicFourierTransform::inverse(const std::vector<Complex
   for (Real σ : σs) {
     for (unsigned n = 0; n < pad * N; n++) {
       if (n < N) {
-        a[n] = ĉ[n] * exp((1 - k) * ω(n));
+        if (σ < 0) {
+          a[n] = std::conj(ĉ[n]) * exp((1 - k) * ω(n));
+        } else {
+          a[n] = ĉ[n] * exp((1 - k) * ω(n));
+        }
       } else {
         a[n] = 0;
       }
-- 
cgit v1.2.3-70-g09d2