public static void Transform( uint bits, ComplexFloat[] source, ComplexFloat[] target, Direction direction) { FPGA.Const <uint> length = GeneratorTools.ArrayLength(bits); float[] cosMap = GeneratorTools.CosArray(length, direction); FPGA.Config.NoSync(cosMap, target); var eK = new ComplexFloat(); var oK = new ComplexFloat(); var rotated = new ComplexFloat(); var tmp = new ComplexFloat(); for (uint i = 0; i < source.Length; i++) { uint j = FPGA.Runtime.Reverse(i, bits); tmp = source[i]; target[j] = tmp; } uint m = 1; uint groupSize = length; for (uint i = 0; i < bits; i++) { FPGA.Config.SetInclusiveRange(0, bits, i); groupSize = groupSize >> 1; for (uint group = 0; group < m; group++) { FPGA.Config.SetInclusiveRange(0, bits, i); uint arg = groupSize * group; for (uint idx = group; idx < length; idx += m * 2) { FPGA.Config.SetInclusiveRange(0, bits, i); eK = target[idx]; oK = target[idx + m]; FTTools.Rotate(cosMap, bits, ref oK, ref rotated, arg); tmp.Re = eK.Re + rotated.Re; tmp.Im = eK.Im + rotated.Im; target[idx] = tmp; tmp.Re = eK.Re - rotated.Re; tmp.Im = eK.Im - rotated.Im; target[idx + m] = tmp; } } m = m << 1; } FTTools.CopyAndNormalize(bits, target, target, direction, ref tmp); }