This Article 
 Bibliographic References 
 Add to: 
A VLSI Constant Geometry Architecture for the Fast Hartley and Fourier Transforms
January 1992 (vol. 3 no. 1)
pp. 58-70
An application-specific architecture for the parallel calculation of the decimation in time and radix 2 fast Hartley (FHT) and Fourier (FFT) transforms is presented. A real sequence with N=2/sup n/ data items is considered as input. The system calculates the FHT and the FFT in n and n+1 stages. respectively. The modular and regular parallel architecture is based on a constant geometry algorithm using butterflies of four data items and the perfect unshuffle permutation. With this permutation, the mapping of the algorithm in VLSI technology is simplified and the communications among processors are minimized. Organization of the processor memory based on first-in, first-out (FIFO) queues facilitates a systolic data flow and permits the implementation in a direct way of the complex data movements and address sequences of the transforms. This is accomplished by means of simple multiplexing operations, using hardwired control. The total calculation time is (Nlog/sub 2/N)/4Q cycles for the FHT and N(1+log/sub 2/N)/4Q cycles for the FFT, where Q is the number of processors (Q= 2/sup q/, Q>or=N/4).

[1] J. I. Agbinya, "Fast interpolation algorithm using fast Hartley transform."Proc. IEEE, vol. 75, no. 4, pp. 523-524, 1987.
[2] R. D. Albon, G. E. Floyd, and J. E. Coles, "A mask programmable DSP array," inProc. IEEE Int. Conf Custom Integrated Circuits, IEEE Press, New York, 1989, pp. 2011-2013.
[3] R.N. Bracewell, "The discrete Hartley transform,"J. Opt. Soc. Amer., vol. 73, no. 12, pp. 1832-1835, 1983.
[4] R.N. Bracewell, "The fast Hartley transform,"Proc. IEEE, vol. 72, no. 8, pp. 1010-1018, 1984.
[5] R.N. Bracewell, O. Buneman, H. Hao, and J. Villasenor, "Fast two-dimensional Hartley transform,"Proc. IEEE, vol. 74, no. 9, pp. 1282-1283, 1986.
[6] R.N. Bracewell,The Fast Hartley Transform. New York: Oxford University Press, 1986.
[7] R.N. Bracewell, H. Bartelt, A. W. Lohmann, and N. Streibl, "Optical synthesis of the Hartley transform,"Appl. Opt., vol. 24, no. 8, pp. 1401-1402, 1985.
[8] G. E. J. Bold, "A comparison of the time involved in computing fast Hartley and fast Fourier transforms,"Proc. IEEE, vol. 73, no. 12, pp. 1863-1864, 1985.
[9] S. Boussakta and A. G. J. Holt, "Calculation of the discréte Hartley transform via the Fermat number transform using a VLSI chip,"IEE Proc., pt. G, vol. 135, no. 3, pp. 101-103, 1988.
[10] S. Boussakta and A. G. J. Holt, "Fast multidimensional discrete Hartley transform using Fermat number transform,"IEE Proc., pt. G, vol. 135, no. 6, pp. 253-257, 1988.
[11] S. Boussakta and A. G. J. Holt, "Proposed prime number Hartley transform implementation using transversal filter structures,"Electron. Lett., vol. 24, no. 15, pp. 926-928, 1988.
[12] S. Boussakta and A. G. J. Holt, "Prime-factor Hartley and Hartley-like transform calculation using transversal filter-type structures,"IEE Proc., pt. G, vol. 136, no. 5, pp. 269-277, 1989.
[13] O. Buneman, "In-situ bit reversed ordering for Hartley transforms,"IEEE Trans. Acoust., Speech, Signal Processing, vol. 37, no. 4, pp. 577-580, 1989.
[14] F. Catthoor and H. J. de Man, "Application-specific architectural methodologies for high throughput digital signal and image processing,"IEEE Trans. Acoust., Speech, Signal Processing, vol. 38, no. 2, pp. 339-349, 1990.
[15] C. Chakrabarti and J. JaJa, "Optimal systolic designs for computing DHT and DCT," inVLSI Signal Processing, Ill, Monterey, CA (New York, IEEE Press), 1988, pp. 411-422.
[16] J. W. Cooley and J. W. Tukey, "An algorithm for the machine calculation of complex Fourier series,"Math. Comput., vol. 19, pp. 297-301, 1965.
[17] D. M. Cottrell, R. A. Lilly, J. A. Davis, and T. Day, "Optical correlator performance of binary phase-only filters using Fourier and Hartley transforms,"Appl. Opt., vol. 26, no. 18, pp. 3755-3761, 1987.
[18] D. Remshaw and P. Denyer,VLSI Signal Processing: A Bit-Serial Approach, Addison Wesley, Reading, Mass., 1985.
[19] A. Drygajlo, "An application of the fast Hartley transform to digital spectral filtering," inProc. European Conf. Circuits Theory and Design 87, Paris, 1987, pp. 263-268.
[20] P. Duhamel and M. Vetterli, "Improved Fourier and Hartley transform algorithms: Application to cyclic convolution of real data,"IEEE Trans. Acoust., Speech, Signal Processing, vol. ASSP-35, no. 6, pp. 818-824, 1987.
[21] R.D. Fellman, "Design issues and an architecture for the monolithic implementation of a parallel digital signal processor,"IEEE Trans. Acoust., Speech, Signal Processing, vol. 38, no. 5, pp. 839-852, 1990.
[22] R.D. Fellman, R.T. Kaneshiro, and K. Konstantinides, "Design and evaluation of architecture for a digital signal processor for instrumentation applications,"IEEE Trans. Acoust., Speech, Signal Processing, vol. 38, no. 3, pp. 537-546, 1990.
[23] S. Gomez, S. Gonzalez, D. D. Hsu, and A. E. H. Kuo, "An application-specific FFT processor,"Electron. Eng., vol. 60, no. 738, pp. 99-106, 1988.
[24] H. Hao and R.N. Bracewell, "A three-dimensional DFT algorithm using the fast Hartley transform,"Proc. IEEE, vol. 75, no. 2, pp. 264-266, 1987.
[25] R. V. L. Hartley, "A more symmetrical Fourier analysis applied to transmission problems,"Proc. IRE, vol. 30, pp. 142-150, 1942.
[26] B. Holland and J. Mather, "Monolithic frequency domain processing with 450 MFLOPS throughput,"Electron. Eng., vol. 61, no. 752, pp. 29-36, 1989.
[27] B. Holland, J. Mather, and S. Brightfield, "Fourier chip tackles transforms,"ESD: The Electronic Syst. Design Mag., pp. 30-36, July 1989.
[28] H.S. Hou, "The Fast Hartley Transform Algorithm,"IEEE Trans. Computers, Vol. C-36, No. 2, Feb. 1987, pp. 147-156.
[29] H.T. Kung, "Why systolic architectures?,"IEEE Comput. Mag., vol. 15, no. 1, pp. 37-46, 1982.
[30] H. T. Kung and C. E. Leiserson, "Systolic arrays for VLSI," inSparse Matrix Proc. 1978, I. S. Duff and G. W. Steward, Eds., Soc. Indust. Appl. Math., 1979, pp. 256-282.
[31] S.Y. Kung,VLSI Array Processors, Prentice Hall, Englewood Cliffs, N.J. 1988.
[32] E. A. Lee, "Programmable DSP architectures: part I,"IEEE ASSP Mag., pp. 4-19, Oct. 1988.
[33] E. A. Lee, "Programmable DSP architectures: part II,"IEEE ASSP Mag., pp. 4-14, Jan. 1989.
[34] Y. Li and G. Eichmann, "Coherent optical generation of Hartley transform of real images,"Opt. Commun., vol. 56, no. 3, pp. 150-154, 1985.
[35] M. Marchesi, G. Orlandi, and F. Piazza, "A systolic circuit for fast Hartley transform," inProc. IEEE Conf. Circuits Syst., Helsinki, 1988, pp. 2685-2688.
[36] B. C. McKinney and F. El-Guibaly, "VLSI design of an FFT processor network,"Integration, VLSI J., vol. 8, no. 3, pp. 301-320, 1989.
[37] R. P. Millane, "Image reconstruction from the Hartley transform intensity."Opt. Commun., vol. 60. no. 5. pp. 269-274. 1986.
[38] M. k. O'Neill, "Faster than fast Fourier,"Byte, no. 4, pp. 293-300, 1988.
[39] C. H. Paik and M.D. Fox, "Fast Hartley transform for image processing,"IEEE Trans. Med Imaging, vol. 7, no. 2, pp. 149-153, 1988.
[40] L. R. Rabiner and B. Gold,Theory and Application of Digital Signal Processing. Englewood Cliffs, NJ: Prentice-Hall, 1975.
[41] H.V. Sorensen, D.L. Jones, C. S. Burrus, and M.T. Heideman, "On computing the discrete Hartley transform,"IEEE Trans. Acoust., Speech, Signal Processing, vol. ASSP-33, no. 4, pp. 1231-1238, 1985.
[42] H.V. Sorensen et al., "Real-Valued Fast Fourier Transform Algorithms,"IEEE Trans. ASSP, Vol. ASSP-35, No. 6, June 1987, pp. 849-863. (Corrections appear inIEEE Trans. Acoustics, Speech and Signal Processing, Vol. ASSP-35, No. 9, Sept. 1987, p. 1353.)
[43] M. C. Steckner and D. J. Drost, "Fast cepstrum analysis using the Hartley transform,"IEEE Trans. Acoust., Speech, Signal Processing, vol. ASSP- 37, no. 8, pp. 1300-1302, 1989.
[44] P. N. Swarztrauber, "Multiprocessor FFT's,"Parallel Comput., vol. 5, pp. 197-210, 1987.
[45] B. V. K. Vijaya Kumar, "Geometric moments computed from the Hartley transform,"Opt. Eng., vol. 25, no. 7, pp. 1327-1336, 1986.
[46] E. L. Zapata,F. Arguello, F. F. Rivera, and J. D. Bruguera, "Multidimensional fast Hartley transform into SIMD hypercubes,"Microprocessing and Microprogramming, vol. 29, no. 2, pp. 121-134, 1990.

Index Terms:
Index Termsfast Hartley transform; fast Fourier transform; FIFO queues; VLSI constant geometryarchitecture; application-specific architecture; parallel calculation; parallel architecture;constant geometry algorithm; butterflies; perfect unshuffle permutation; processormemory; systolic data flow; multiplexing operations; hardwired control; computationalcomplexity; fast Fourier transforms; parallel algorithms; parallel architectures; VLSI
E.L. Zapata, F. Argüello, "A VLSI Constant Geometry Architecture for the Fast Hartley and Fourier Transforms," IEEE Transactions on Parallel and Distributed Systems, vol. 3, no. 1, pp. 58-70, Jan. 1992, doi:10.1109/71.113082
Usage of this product signifies your acceptance of the Terms of Use.