Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters

Abstract : This paper presents a posteriori error estimates for conforming numerical approximations of eigenvalue clusters of second-order self-adjoint elliptic linear operators with compact resolvent. Given a cluster of eigenvalues, we estimate the error in the sum of the eigenvalues, as well as the error in the eigenvectors represented through the density matrix, i.e., the orthogonal projector on the associated eigenspace. This allows us to deal with degenerate (multiple) eigenvalues within the framework. All the bounds are valid under the only assumption that the cluster is separated from the surrounding smaller and larger eigenvalues; we show how this assumption can be numerically checked. Our bounds are guaranteed and converge with the same speed as the exact errors. They can be turned into fully computable bounds as soon as an estimate on the dual norm of the residual is available, which is presented in two particular cases: the Laplace eigenvalue problem discretized with conforming finite elements, and a Schrödinger operator with periodic boundary conditions of the form −∆ + V discretized with planewaves. For these two cases, numerical illustrations are provided on a set of test problems.
Document type :
Preprints, Working Papers, ...
Complete list of metadatas

Cited literature [42 references]  Display  Hide  Download

https://hal.archives-ouvertes.fr/hal-02127954
Contributor : Geneviève Dusson <>
Submitted on : Monday, May 13, 2019 - 7:28:57 PM
Last modification on : Tuesday, September 17, 2019 - 1:24:41 AM

File

eig_clusters.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-02127954, version 1

Citation

Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Vohralík. Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters. 2019. ⟨hal-02127954⟩

Share

Metrics

Record views

124

Files downloads

121