From e0126108e6a9e6b1b1a30ee830158004559445b8 Mon Sep 17 00:00:00 2001
From: Martin Reinecke <martin@mpa-garching.mpg.de>
Date: Thu, 25 Jul 2019 22:17:36 +0200
Subject: [PATCH] flip entire vectors, not their components

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

diff --git a/pocketfft_hdronly.h b/pocketfft_hdronly.h
index ea6538e..a27b8a1 100644
--- a/pocketfft_hdronly.h
+++ b/pocketfft_hdronly.h
@@ -3210,14 +3210,12 @@ template<typename T> POCKETFFT_NOINLINE void general_r(
             tdatav[i][j] = tin[it.iofs(j,i)];
         if ((!r2c) && forward)
           for (size_t i=2; i<len; i+=2)
-            for (size_t j=0; j<vlen; ++j)
-              tdatav[i][j] = -tdatav[i][j];
+            tdatav[i] = -tdatav[i];
         forward ? plan->forward (tdatav, fct)
                 : plan->backward(tdatav, fct);
         if (r2c && (!forward))
           for (size_t i=2; i<len; i+=2)
-            for (size_t j=0; j<vlen; ++j)
-              tdatav[i][j] = -tdatav[i][j];
+            tdatav[i] = -tdatav[i];
         for (size_t i=0; i<len; ++i)
           for (size_t j=0; j<vlen; ++j)
             out[it.oofs(j,i)] = tdatav[i][j];
-- 
GitLab