English

Comparison Theorems for Splittings of M-matrices in (block) Hessenberg Form

Numerical Analysis 2021-11-18 v2 Numerical Analysis

Abstract

Some variants of the (block) Gauss--Seidel iteration for the solution of linear systems with MM-matrices in (block) Hessenberg form are discussed. Comparison results for the asymptotic convergence rate of some regular splittings are derived: in particular, we prove that for a lower-Hessenberg M-matrix ρ(PGS)ρ(PS)ρ(PAGS)\rho(P_{GS})\geq \rho(P_S)\geq \rho(P_{AGS}), where PGS,PS,PAGSP_{GS}, P_S, P_{AGS} are the iteration matrices of the Gauss--Seidel, staircase, and anti-Gauss--Seidel method. This is a result that does not seem to follow from classical comparison results, as these splittings are not directly comparable. It is shown that the concept of stair partitioning provides a powerful tool for the design of new variants that are suited for parallel computation.

Keywords

Cite

@article{arxiv.2106.10492,
  title  = {Comparison Theorems for Splittings of M-matrices in (block) Hessenberg Form},
  author = {Luca Gemignani and Federico Poloni},
  journal= {arXiv preprint arXiv:2106.10492},
  year   = {2021}
}