Module 545 in TOMSCMFFT and RMFFT: Fortran subprograms for computing real and complex fast Fourier transforms of a one-dimensional or multidimensional data set. The programs are designed to minimize I/O for the case where the data lie on mass storage. (See D. Fraser, ACM TOMS 5 (1979) pp. 500-517.).