A fast FFT-based discrete Legendre transform

An O(N(logN)2/ loglogN) algorithm for computing the discrete Legendre transform and its inverse is described. The algorithm combines a recently developed fast transform for converting between Legendre and Chebyshev coefficients with a Taylor series expansion for Chebyshev polynomials about equallyspaced points in the frequency domain. Both components are based on the FFT, and as an intermediate step we obtain an O(N logN) algorithm for evaluating a degree N−1 Chebyshev expansion at an N-point Legendre grid. Numerical results are given to demonstrate performance and accuracy.