Commit 5836cc16 authored by Martin Reinecke's avatar Martin Reinecke
Browse files


parent 5331805d
......@@ -2225,7 +2225,7 @@ template<typename T> void pocketfft_general_hartley(const shape_t &shape,
plan->forward((T *)tdata, fct);
// Hartley order
data_out[it_out.offset()] = tdata[0];
size_t i=1, i1=1, i2=len()-1;
size_t i=1, i1=1, i2=len-1;
for (i=1; i<len-1; i+=2, ++i1, --i2)
data_out[it_out.offset()+i1*s_o] = tdata[i]+tdata[i+1];
Supports Markdown
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