On fixed-point, Krylov, and 2× 2 block preconditioners for nonsymmetric problems

by   Ben S. Southworth, et al.

The solution of matrices with 2× 2 block structure arises in numerous areas of computational mathematics, such as PDE discretizations based on mixed-finite element methods, constrained optimization problems, or the implicit or steady state treatment of any system of PDEs with multiple dependent variables. Often, these systems are solved iteratively using Krylov methods and some form of block preconditioner. Under the assumption that one diagonal block is inverted exactly, this paper proves a direct equivalence between convergence of 2×2 block preconditioned Krylov or fixed-point iterations to a given tolerance, with convergence of the underlying preconditioned Schur-complement problem. In particular, results indicate that an effective Schur-complement preconditioner is a necessary and sufficient condition for rapid convergence of 2× 2 block-preconditioned GMRES, for arbitrary relative-residual stopping tolerances. A number of corollaries and related results give new insight into block preconditioning, such as the fact that approximate block-LDU or symmetric block-triangular preconditioners offer minimal reduction in iteration over block-triangular preconditioners, despite the additional computational cost. Theoretical results are verified numerically on a nonsymmetric steady linearized Navier-Stokes discretization, which also demonstrate that theory based on the assumption of an exact inverse of one diagonal block extends well to the more practical setting on inexact inverses.


page 1

page 2

page 3

page 4


A note on 2× 2 block-diagonal preconditioning

For 2x2 block matrices, it is well-known that block-triangular or block-...

A new block diagonal preconditioner for a class of 3× 3 block saddle point problems

We study the performance of a new block preconditioner for a class of 3×...

A Mass Conserving Mixed hp-FEM Scheme for Stokes Flow. Part III: Implementation and Preconditioning

This is the third part in a series on a mass conserving, high order, mix...

Uniform block diagonal preconditioners for divergence-conforming HDG Methods for the generalized Stokes problem and linear elasticity

We propose a uniform block diagonal preconditioner for the condensed H(d...

Schur complement based preconditioners for twofold and block tridiagonal saddle point problems

In this paper, two types of Schur complement based preconditioners are s...

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

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

GPMR: An Iterative Method for Unsymmetric Partitioned Linear Systems

We introduce an iterative method named GPMR for solving 2x2 block unsymm...

Please sign up or login with your details

Forgot password? Click here to reset