- const Double_t t1_real = REAL(data,stride,i) + REAL(data,stride,j);
- const Double_t t1_imag = IMAG(data,stride,i) + IMAG(data,stride,j);
- const Double_t t2_real = REAL(data,stride,i) - REAL(data,stride,j);
- const Double_t t2_imag = IMAG(data,stride,i) - IMAG(data,stride,j);
-
- REAL(data,stride,i) = t1_real;
- IMAG(data,stride,i) = t1_imag;
- REAL(data,stride,j) = w_real*t2_real - w_imag * t2_imag;
- IMAG(data,stride,j) = w_real*t2_imag + w_imag * t2_real;
+ const Double_t kt1real = REAL(data,stride,ki) + REAL(data,stride,kj);
+ const Double_t kt1imag = IMAG(data,stride,ki) + IMAG(data,stride,kj);
+ const Double_t kt2real = REAL(data,stride,ki) - REAL(data,stride,kj);
+ const Double_t kt2imag = IMAG(data,stride,ki) - IMAG(data,stride,kj);
+
+ REAL(data,stride,ki) = kt1real;
+ IMAG(data,stride,ki) = kt1imag;
+ REAL(data,stride,kj) = wreal*kt2real - wimag * kt2imag;
+ IMAG(data,stride,kj) = wreal*kt2imag + wimag * kt2real;