In this paper we describe a method for computing the Discrete Fourier Transform (DFT) of a sequence of n elements over a finite field GF ( p m ) {\text {GF}}({p^m}) with a number of bit operations O ( n m log ⁡ ( n m ) ⋅ P ( q ) ) O(nm\,\log (nm) \cdot P(q)) where P ( q ) P(q) is the number of bit operations required to multiply two q-bit integers and q ≅ 2 log 2 n + 4 log 2 m + 4 log 2 p q \cong 2\,{\log _2}n + 4\,{\log _2}m + 4\,{\log _2}p . This method is uniformly applicable to all instances and its order of complexity is not inferior to that of methods whose success depends upon the existence of certain primes. Our algorithm is a combination of known and novel techniques. In particular, the finite-field DFT is at first converted into a finite field convolution; the latter is then implemented as a two-dimensional Fourier transform over the complex field. The key feature of the method is the fusion of these two basic operations into a single integrated procedure centered on the Fast Fourier Transform algorithm.