Commit 08c85d30 authored by Martin Reinecke's avatar Martin Reinecke

bug fix in dcst4

parent d250df55
...@@ -2549,9 +2549,9 @@ template<typename T0> class T_dcst4 ...@@ -2549,9 +2549,9 @@ template<typename T0> class T_dcst4
{ {
if ((N&1)==0) 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) for (size_t i=0; i<N/2; ++i)
C2[i] = tw[8*N-8*i-1]; C2[i] = conj(tw[8*i+1]);
} }
} }
......
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