From 08c85d30e1f1ee529da149a5533f597bfa02f5e4 Mon Sep 17 00:00:00 2001
From: Martin Reinecke <martin@mpa-garching.mpg.de>
Date: Mon, 21 Oct 2019 11:27:02 +0200
Subject: [PATCH] bug fix in dcst4

---
 pocketfft_hdronly.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/pocketfft_hdronly.h b/pocketfft_hdronly.h
index 0b1fe15..8c59c48 100644
--- a/pocketfft_hdronly.h
+++ b/pocketfft_hdronly.h
@@ -2549,9 +2549,9 @@ template<typename T0> class T_dcst4
       {
       if ((N&1)==0)
         {
-        sincos_2pibyn<T0> tw(8*N);
+        sincos_2pibyn<T0> tw(16*N);
         for (size_t i=0; i<N/2; ++i)
-          C2[i] = tw[8*N-8*i-1];
+          C2[i] = conj(tw[8*i+1]);
         }
       }
 
-- 
GitLab