Skip to main content
Engineering LibreTexts

8.5: Evaluation of the Cooley-Tukey FFT Algorithms

  • Page ID
    2011
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    The evaluation of any FFT algorithm starts with a count of the real (or floating point) arithmetic. The Table 8.5.1 below gives the number of real multiplications and additions required to calculate a length-N FFT of complex data. Results of programs with one, two, three and five butterflies are given to show the improvement that can be expected from removing unnecessary multiplications and additions. Results of radices two, four, eight and sixteen for the Cooley-Tukey FFT as well as of the split-radix FFT are given to show the relative merits of the various structures. Comparisons of these data should be made with the table of counts for the PFA and WFTA programs in The Prime Factor and Winograd Fourier Transform Algorithms . All programs use the four-multiply-two-add complex multiply algorithm. A similar table can be developed for the three-multiply-three-add algorithm, but the relative results are the same.

    From the table it is seen that a greater improvement is obtained going from radix-2 to 4 than from 4 to 8 or 16. This is partly because length 2 and 4 butterflies have no multiplications while length 8, 16 and higher do. It is also seen that going from one to two butterflies gives more improvement than going from two to higher values. From an operation count point of view and from practical experience, a three butterfly radix-4 or a two butterfly radix-8 FFT is a good compromise. The radix-8 and 16 programs become long, especially with multiple butterflies, and they give a limited choice of transform length unless combined with some length 2 and 4 butterflies.

    N M1 M2 M3 M5 A1 A2 A3 A5
    2 4 0 0 0 6 4 4 4
    4 16 4 0 0 24 18 16 16
    8 48 20 8 4 72 58 52 52
    16 128 68 40 28 192 162 148 148
    32 320 196 136 108 480 418 388 388
    64 768 516 392 332 1152 1026 964 964
    128 1792 1284 1032 908 2688 2434 2308 2308
    256 4096 3076 2568 2316 6144 5634 5380 5380
    512 9216 7172 6152 5644 13824 12802 12292 12292
    1024 20480 16388 14344 13324 30720 28674 27652 27652
    2048 45056 36868 32776 30732 67584 63490 61444 61444
    4096 98304 81924 73736 69644 147456 139266 135172 135172
    4 12 0 0 0 22 16 16 16
    16 96 36 28 24 176 146 144 144
    64 576 324 284 264 1056 930 920 920
    256 3072 2052 1884 1800 5632 5122 5080 5080
    1024 15360 11268 10588 10248 28160 26114 25944 25944
    4096 73728 57348 54620 53256 135168 126978 126296 126296
    8 32 4 4 4 66 52 52 52
    64 512 260 252 248 1056 930 928 928
    512 6144 4100 4028 3992 12672 11650 11632 11632
    4096 65536 49156 48572 48280 135168 126978 126832 126832
    16 80 20 20 20 178 148 148 148
    256 2560 1540 1532 1528 5696 5186 5184 5184
    4096 61440 45060 44924 44856 136704 128514 128480 128480
    2 0 0 0 0 4 4 4 4
    4 8 0 0 0 20 16 16 16
    8 24 8 4 4 60 52 52 52
    16 72 32 28 24 164 144 144 144
    32 184 104 92 84 412 372 372 372
    64 456 288 268 248 996 912 912 912
    128 1080 744 700 660 2332 2164 2164 2164
    256 2504 1824 1740 1656 5348 5008 5008 5008
    512 5688 4328 4156 3988 12060 11380 11380 11380
    1024 12744 10016 9676 9336 26852 25488 25488 25488
    2048 28216 22760 22076 21396 59164 56436 56436 56436
    4096 61896 50976 49612 48248 129252 123792 123792 123792

    Table 8.5.1: Number of Real Multiplications and Additions for Complex Single Radix FFTs

    In Table 8.2.1 Mi and Ai refer to the number of real multiplications and real additions used by an FFT with i separately written butterflies. The first block has the counts for Radix-2, the second for Radix-4, the third for Radix-8, the fourth for Radix-16, and the last for the Split-Radix FFT. For the split-radix FFT, M3 and A3 refer to the two- butterfly-plus program and M5 and A5 refer to the three-butterfly program.

    The first evaluations of FFT algorithms were in terms of the number of real multiplications required as that was the slowest operation on the computer and, therefore, controlled the execution speed. Later with hardware arithmetic both the number of multiplications and additions became important. Modern systems have arithmetic speeds such that indexing and data transfer times become important factors. Morris has looked at some of these problems and has developed a procedure called autogen to write partially straight-line program code to significantly reduce overhead and speed up FFT run times. Some hardware, such as the TMS320 signal processing chip, has the multiply and add operations combined. Some machines have vector instructions or have parallel processors. Because the execution speed of an FFT depends not only on the algorithm, but also on the hardware architecture and compiler, experiments must be run on the system to be used.

    In many cases the unscrambler or bit-reverse-counter requires 10% of the execution time, therefore, if possible, it should be eliminated. In high-speed convolution where the convolution is done by multiplication of DFT's, a decimation-in-frequency FFT can be combined with a decimation-in-time inverse FFT to require no unscrambler. It is also possible for a radix-2 FFT to do the unscrambling inside the FFT but the structure is not very regular.

    Although there can be significant differences in the efficiencies of the various Cooley-Tukey and Split-Radix FFTs, the number of multiplications and additions for all of them is on the order of \(N\log N\). That is fundamental to the class of algorithms.

    Contributor

    • ContribEEBurrus

    This page titled 8.5: Evaluation of the Cooley-Tukey FFT Algorithms is shared under a CC BY license and was authored, remixed, and/or curated by C. Sidney Burrus.

    • Was this article helpful?