FFTPACK: A package of Fortran subprograms for the Fast Fourier Transform of periodic and other symmetric sequences.