Fixed serious mistake in log inverse transform

1 files changed, 6 insertions, 1 deletions

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;
       }