Premium
Multigrid methods for block‐Toeplitz linear systems: convergence analysis and applications
Author(s) -
Donatelli Marco,
Ferrari Paola,
Furci Isabella,
SerraCapizzano Stefano,
Sesana Debora
Publication year - 2021
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.2356
Subject(s) - toeplitz matrix , transfer operator , mathematics , multigrid method , rate of convergence , block (permutation group theory) , scalar (mathematics) , linear system , algorithm , mathematical optimization , computer science , partial differential equation , mathematical analysis , pure mathematics , combinatorics , geometry , computer network , channel (broadcasting)
Abstract In the past decades, multigrid methods for linear systems having multilevel Toeplitz coefficient matrices with scalar entries have been widely studied. On the other hand, only few papers have investigated the case of block entries, where the entries are small generic matrices of fixed size instead of scalars. In that case the efforts of the researchers have been mainly devoted to specific applications, focusing on algorithmic proposals, but with marginal theoretical results. In this paper, we propose a general two‐grid convergence analysis, proving an optimal convergence rate independent of the matrix size, in the case of positive definite block‐Toeplitz matrices with generic blocks. In particular, the proof of the approximation property is not a straightforward generalization of the scalar case and, in fact, we have to require a specific commutativity condition on the block symbol of the grid transfer operator. According to the analysis, we define a class of grid transfer operators satisfying the previous theoretical conditions and we propose a strategy to ensure fast multigrid convergence even for more than two grids. Among the numerous applications that lead to the block‐Toeplitz structure, high‐order Lagrangian finite element methods and staggered discontinuous Galerkin methods are considered in the numerical results, confirming the effectiveness of our proposal and the correctness of the proposed theoretical analysis.