A hierarchical method for obtaining eigenvalue enclosures

DOI10.1090/S0025-5718-00-01238-2MaRDI QIDQ4501034zbMATH OpenFDO

Publication date 3 September 2000

Published in Mathematics of Computation (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/math/9803075

eigenvalues interval arithmetic Rayley-Ritz method spectral enclosures spectrum Sturm-Liouville operator Temple-Lehmann method

Mathematics Subject Classification ID

Eigenvalues, estimation of eigenvalues, upper and lower bounds of ordinary differential operators (34L15) Estimates of eigenvalues in context of PDEs (35P15) Numerical solution of eigenvalue problems involving ordinary differential equations (65L15) Finite element, Rayleigh-Ritz, Galerkin and collocation methods for ordinary differential equations (65L60) Error bounds for numerical methods for ordinary differential equations (65L70) Numerical methods for eigenvalue problems for boundary value problems involving PDEs (65N25)

Abstract: We propose a method for obtaining rigorous and accurate upper and lower bounds on the eigenvalues of ordinary and partial differential operators in bounded regions of Euclidean space. It uses a boundary condition homotopy method starting from a decomposition of the region into several simpler regions for which the eigenvalues can already be computed accurately. Several numerical examples are presented.

Review text

The author introduces a new method for obtaining guaranteed enclosures of the eigenvalues of a variety of selfadjoint differential and difference operators with discrete spectrum. First, the method is described for the Sturm-Liouville operator \(H\) acting in \(L^2(\alpha,\beta)\) of the form \[ Hf(x):=-{d\over{dx}}[a(x){{df}\over{dx}}]+V(x)f(x) \] subject to Neumann boundary conditions. The method is based upon decoupling the interval \((\alpha,\beta)\) into \(2^N\) subintervals \((\alpha_j,\alpha_{j+1})\). The subintervals should be such that the restrictions of H onto \(L^2(\alpha_j,\alpha_{j+1})\) are close to constant coefficient operators. The several stages algorithm is proposed for obtaining enclosures of eigenvalues of \(H\). If the algorithm does finish, then the eigenvalue enclosures obtained are correct. Extensions of the method to selfadjoint systems of Sturm-Liouville operators and to higher-order differential operators are also described. Moreover, the author shows how the method may be applied with modifications to discrete Laplacians on graphs and to the Laplacian in \(N\) dimensions in bounded regions with piecewise smooth boundaries.

Recommendations

Cites work

Cited in

(11)

This page was built for publication: A hierarchical method for obtaining eigenvalue enclosures

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4501034)