Volume 53, pp. 113-216, 2020.

Block generalized locally Toeplitz sequences: theory and applications in the multidimensional case

Giovanni Barbarino, Carlo Garoni, and Stefano Serra-Capizzano

Abstract

In computational mathematics, when dealing with a large linear discrete problem (e.g., a linear system) arising from the numerical discretization of a partial differential equation (PDE), knowledge of the spectral distribution of the associated matrix has proved to be useful information for designing/analyzing appropriate solvers–-especially, preconditioned Krylov and multigrid solvers–-for the considered problem. Actually, this spectral information is of interest also in itself as long as the eigenvalues of the aforementioned matrix represent physical quantities of interest, which is the case for several problems from engineering and applied sciences (e.g., the study of natural vibration frequencies in an elastic material). The theory of multilevel generalized locally Toeplitz (GLT) sequences is a powerful apparatus for computing the asymptotic spectral distribution of matrices $A_n$ arising from virtually any kind of numerical discretization of PDEs. Indeed, when the mesh-fineness parameter $n$ tends to infinity, these matrices $A_n$ give rise to a sequence $\{A_n{\}}_n$, which often turns out to be a multilevel GLT sequence or one of its “relatives”, i.e., a multilevel block GLT sequence or a (multilevel) reduced GLT sequence. In particular, multilevel block GLT sequences are encountered in the discretization of systems of PDEs as well as in the higher-order finite element or discontinuous Galerkin approximation of scalar/vectorial PDEs. In this work, we systematically develop the theory of multilevel block GLT sequences as an extension of the theories of (unilevel) GLT sequences [Garoni and Serra-Capizzano, Generalized Locally Toeplitz Sequences: Theory and Applications. Vol. I., Springer, Cham, 2017], multilevel GLT sequences [Garoni and Serra-Capizzano, Generalized Locally Toeplitz Sequences: Theory and Applications. Vol. II., Springer, Cham, 2018], and block GLT sequences [Barbarino, Garoni, and Serra-Capizzano, Electron. Trans. Numer. Anal., 53 (2020), pp. 28–112]. We also present several emblematic applications of this theory in the context of PDE discretizations.

Full Text (PDF) [1.3 MB], BibTeX , DOI: 10.1553/etna_vol53s113

Key words

asymptotic distribution of singular values and eigenvalues, multilevel block Toeplitz matrices, multilevel block generalized locally Toeplitz matrices, numerical discretization of partial differential equations, finite differences, finite elements, isogeometric analysis, discontinuous Galerkin methods, tensor products, B-splines

AMS subject classifications

15A18, 15B05, 47B06, 65N06, 65N30, 65N25, 15A60, 15A69, 65D07

Links to the cited ETNA articles

[8]	Vol. 53 (2020), pp. 28-112 Giovanni Barbarino, Carlo Garoni, and Stefano Serra-Capizzano: Block generalized locally Toeplitz sequences: theory and applications in the unidimensional case

ETNA articles which cite this article

Vol. 55 (2022), pp. 585-617 Giovanni Barbarino, Carlo Garoni, Mariarosa Mazza, and Stefano Serra-Capizzano: Rectangular GLT sequences

Vol. 58 (2023), pp. 136-163 Matthias Bolten, Sven-Erik Ekström, Isabella Furci, and Stefano Serra-Capizzano: A note on the spectral analysis of matrix sequences via GLT momentary symbols: from all-at-once solution of parabolic problems to distributed fractional order matrices

Vol. 59 (2023), pp. 157-178 Marco Donatelli, Paola Ferrari, and Silvia Gazzola: Symmetrization techniques in image deblurring