Multigrid methods for block-Toeplitz linear systems: convergence analysis and applications

by   Marco Donatelli, et al.

In the past decades, multigrid methods for linear systems having multilevel Toeplitz coefficient matrices with scalar entries have been largely studied. On the other hand, only few papers have investigated the case of block entries, where the entries are small generic matrices instead of scalars. In that case the efforts of the researchers have been mainly devoted to specific applications, focusing on algorithmic proposals but with very 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 has 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. Furthermore, we define a class of grid transfer operators satisfying the previous theoretical conditions and we propose a strategy to insure 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.


page 1

page 2

page 3

page 4


A symbol based analysis for multigrid methods for Block-Circulant and Block-Toeplitz Systems

In the literature, there exist several studies on symbol-based multigrid...

Symbol Based Convergence Analysis in Block Multigrid Methods with applications for Stokes problems

The main focus of this paper is the study of efficient multigrid methods...

Convergence of Discontinuous Galerkin Schemes for the Euler Equations via Dissipative Weak Solutions

In this paper, we present convergence analysis of high-order finite elem...

Quasi-Perron-Frobenius property of a class of saddle point matrices

The saddle point matrices arising from many scientific computing fields ...

Admissible and attainable convergence behavior of block Arnoldi and GMRES

It is well-established that any non-increasing convergence curve is poss...

A matrix-theoretic spectral analysis of incompressible Navier-Stokes staggered DG approximation and related solvers

The incompressible Navier-Stokes equations are solved in a channel, usin...

Unconstrained Proximal Operator: the Optimal Parameter for the Douglas-Rachford Type Primal-Dual Methods

In this work, we propose an alternative parametrized form of the proxima...

Please sign up or login with your details

Forgot password? Click here to reset