Commit e0126108 authored by Martin Reinecke's avatar Martin Reinecke

flip entire vectors, not their components

parent d72f8ade
......@@ -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];
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment