Volume 54, pp. 370-391, 2021.

A two-level iterative scheme for general sparse linear systems based on approximate skew-symmetrizers

Murat Manguoğlu and Volker Mehrmann

Abstract

We propose a two-level iterative scheme for solving general sparse linear systems. The proposed scheme consists of a sparse preconditioner that increases the norm of the skew-symmetric part relative to the rest and makes the main diagonal of the coefficient matrix as close to the identity as possible so that the preconditioned system is as close to a shifted skew-symmetric matrix as possible. The preconditioned system is then solved via a particular Minimal Residual Method for Shifted Skew-Symmetric Systems (MRS). This leads to a two-level (inner and outer) iterative scheme where the MRS has short-term recurrences and satisfies an optimality condition. A preconditioner for the inner system is designed via a skew-symmetry-preserving deflation strategy based on the skew-Lanczos process. We demonstrate the robustness of the proposed scheme on sparse matrices from various applications.

Full Text (PDF) [979 KB], BibTeX

Key words

symmetrizer, skew-symmetrizer, Krylov subspace method, shifted skew-symmetric system, skew-Lanczos method

AMS subject classifications

65F08, 65F10, 65F50

Links to the cited ETNA articles

[28]	Vol. 23 (2006), pp. 158-179 Olaf Schenk and Klaus Gärtner: On fast factorization pivoting methods for sparse symmetric indefinite systems

ETNA articles which cite this article

Vol. 55 (2022), pp. 455-468 Chen Greif: Structured shifts for skew-symmetric matrices

< Back