pith. sign in

arxiv: 2201.01314 · v1 · pith:TQGEONUKnew · submitted 2022-01-04 · 🧮 math.NA · cs.NA· math-ph· math.MP· math.SP

SpecSolve: Spectral methods for spectral measures

classification 🧮 math.NA cs.NAmath-phmath.MPmath.SP
keywords spectralmethodsoperatorsspecsolvemeasuresadaptivecontinuousgeneral
0
0 comments X
read the original abstract

Self-adjoint operators on infinite-dimensional spaces with continuous spectra are abundant but do not possess a basis of eigenfunctions. Rather, diagonalization is achieved through spectral measures. The SpecSolve package [SIAM Rev., 63(3) (2021), pp. 489--524] computes spectral measures of general (self-adjoint) differential and integral operators by combining state-of-the-art adaptive spectral methods with an efficient resolvent-based strategy. The algorithm achieves arbitrarily high orders of convergence in terms of a smoothing parameter, allowing computation of both discrete and continuous spectral components. This article extends SpecSolve to two important classes of operators: singular integro-differential operators and general operator pencils. Essential computational steps are performed with off-the-shelf spectral methods, including spectral methods on the real line, the ultraspherical spectral method, Chebyshev and Fourier spectral methods, and the ($hp$-adaptive and sparse) ultraspherical spectral element method. This collection illustrates the power and flexibility of SpecSolve's "discretization-oblivious" paradigm.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.