This article concerns the spectral analysis of matrix-sequences which can be written as a non-Hermitian perturbation of a given Hermitian matrix-sequence. The main result reads as follows. Suppose that for every n there is a Hermitian matrix Xn of size n and that {Xn}n∼λf, that is, the matrix-sequence {Xn}n enjoys an asymptotic spectral distribution, in the Weyl sense, described by a Lebesgue measurable function f; if (Formula presented.) with ‖·‖2 being the Schatten 2 norm, then {Xn+Yn}n∼λf. In a previous article by Leonid Golinskii and the second author, a similar result was proved, but under the technical restrictive assumption that the involved matrix-sequences {Xn}n and {Yn}n are uniformly bounded in spectral norm. Nevertheless, the result had a remarkable impact in the analysis of both spectral distribution and clustering of matrix-sequences arising from various applications, including the numerical approximation of partial differential equations (PDEs) and the preconditioning of PDE discretization matrices. The new result considerably extends the spectral analysis tools provided by the former one, and in fact we are now allowed to analyze linear PDEs with (unbounded) variable coefficients, preconditioned matrix-sequences, and so forth. A few selected applications are considered, extensive numerical experiments are discussed, and a further conjecture is illustrated at the end of the article.

Non-Hermitian perturbations of Hermitian matrix-sequences and applications to the spectral analysis of the numerical approximation of partial differential equations

Serra Capizzano S.
2020-01-01

Abstract

This article concerns the spectral analysis of matrix-sequences which can be written as a non-Hermitian perturbation of a given Hermitian matrix-sequence. The main result reads as follows. Suppose that for every n there is a Hermitian matrix Xn of size n and that {Xn}n∼λf, that is, the matrix-sequence {Xn}n enjoys an asymptotic spectral distribution, in the Weyl sense, described by a Lebesgue measurable function f; if (Formula presented.) with ‖·‖2 being the Schatten 2 norm, then {Xn+Yn}n∼λf. In a previous article by Leonid Golinskii and the second author, a similar result was proved, but under the technical restrictive assumption that the involved matrix-sequences {Xn}n and {Yn}n are uniformly bounded in spectral norm. Nevertheless, the result had a remarkable impact in the analysis of both spectral distribution and clustering of matrix-sequences arising from various applications, including the numerical approximation of partial differential equations (PDEs) and the preconditioning of PDE discretization matrices. The new result considerably extends the spectral analysis tools provided by the former one, and in fact we are now allowed to analyze linear PDEs with (unbounded) variable coefficients, preconditioned matrix-sequences, and so forth. A few selected applications are considered, extensive numerical experiments are discussed, and a further conjecture is illustrated at the end of the article.
2020
approximation of PDEs; perturbation results; preconditioning; spectral distribution in the Weyl sense
Barbarino, G.; Serra Capizzano, S.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11383/2119503
 Attenzione

L'Ateneo sottopone a validazione solo i file PDF allegati

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 24
  • ???jsp.display-item.citation.isi??? 17
social impact