Changes

tr = (currentPhaseShiftReal * real[off]) - (currentPhaseShiftImag * imag[off]);

ti = (currentPhaseShiftReal * imag[off]) + (currentPhaseShiftImag * real[off]);

real[off] = real[i] - tr;

imag[off] = imag[i] - ti;

real[i] += tr;

imag[i] += ti;

i += halfSize << 1;

}