Legendre Series Analysis and Computation via Composed Abel-Fourier Transform

Legendre coefficients of an integrable function f(x) are proved to coincide with the Fourier coefficients with a nonnegative index of a suitable Abel-type transform of the function itself. The numerical computation of N Legendre coefficients can thus be carried out efficiently in O(NlogN) operations by means of a single fast Fourier transform of the Abel-type transform of f(x). Symmetries associated with the Abel-type transform are exploited to further reduce the computational complexity. The dual problem of calculating the sum of Legendre expansions at a prescribed set of points is also considered. We prove that a Legendre series can be written as the Abel transform of a suitable Fourier series. This fact allows us to state an efficient algorithm for the evaluation of Legendre expansions. Finally, some numerical tests are illustrated to exemplify and confirm the theoretical results.

Automated summary

This paper proves that the Legendre coefficients of an integrable function f(x) are equal to the Fourier coefficients with a nonnegative index of a suitable Abel-type transform of the function itself. The computation of N Legendre coefficients can be efficiently carried out using a single fast Fourier transform of the Abel-type transform of f(x), with a complexity of O(NlogN) operations. The symmetries associated with the Abel-type transform are utilized to further reduce the computational complexity. The paper also discusses the dual problem of calculating the sum of Legendre expansions at a prescribed set of points, and presents an efficient algorithm based on the fact that a Legendre series can be written as the Abel transform of a suitable Fourier series. Numerical tests are provided to illustrate and confirm the theoretical results.