The computation of the spectra of highly oscillatory fredholm integral operators

Hermann BRUNNER*, Arieh Iserles, Syvert P. Nørsett

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

6 Citations (Scopus)

Abstract

We are concerned with the computation of spectra of highly oscillatory Fredholm problems, in particular with the Fox-Li operator whereω>> 1. Our main tool is the finite section method: an eigenfunction is expanded in an orthonormal basis of the underlying space, resulting in an algebraic eigenvalue problem. We consider two competing bases: a basis of Legendre polynomials and a basis consisting of modified Fourier functions (cosines and shifted sines), and derive detailed asymptotic estimates of the rate of decay of the coefficients. Although the Legendre basis enjoys in principle much faster convergence, this does not lead to much smaller matrices. Since the computation of Legendre coefficients is expensive, while modified Fourier coefficients can be computed efficiently with FFT, we deduce that modified Fourier expansions, implemented in a manner that takes advantage of their structure, present a considerably more effective tool for the computation of highly oscillatory Fredholm spectra.

Original languageEnglish
Pages (from-to)467-519
Number of pages53
JournalJournal of Integral Equations and Applications
Volume23
Issue number4
DOIs
Publication statusPublished - Dec 2011

Scopus Subject Areas

  • Numerical Analysis
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'The computation of the spectra of highly oscillatory fredholm integral operators'. Together they form a unique fingerprint.

Cite this