Spectral Analysis of Block Diagonally Preconditioned Multiple Saddle-Point Matrices with Inexact Schur Complements
Marco Pilotto, Luca Bergamaschi, Angeles Martinez
TL;DR
The paper develops a polynomial-based framework to bound the eigenvalues of block-diagonal Schur-complement preconditioned multiple saddle-point systems with inexact Schur complements, valid for an arbitrary number of blocks $N$. By linking the spectrum to zeros of a monic sequence of Chebyshev-like polynomials $U_k(\\lambda,\\boldsymbol{\\gamma}_E,\\boldsymbol{\\gamma}_R)$ and introducing the parameter rays $\\gamma_E^{(i)},\\gamma_R^{(i)}$, the authors derive interlacing properties, extremal-zero bounds, and explicit eigenvalue intervals. They further show how zeros move with the parameters, enabling evaluation of bounds at endpoint values, and provide concrete bounds for the double saddle-point case $N=2$. Numerical experiments on random and realistic Biot-type problems validate the theoretical bounds and demonstrate the effectiveness of inexact block-diagonal preconditioning in MINRES. The results offer rigorous, computable spectral estimates that inform preconditioner design for multiphysics systems with multiple saddle points.
Abstract
We derive eigenvalue bounds for symmetric block-tridiagonal multiple saddle-point systems preconditioned with block-diagonal Schur complement matrices. This analysis applies to an arbitrary number of blocks and accounts for the case where the Schur complements are approximated, generalizing the findings in [Bergamaschi et al., Linear Algebra and its Applications, 2026]. Numerical experiments are carried out to validate the proposed estimates.
