Solving linear systems of the form \((A + \gamma UU^T)\, {\bf x} = {\bf b}\) by preconditioned iterative methods

We consider the iterative solution of large linear systems of equations in which the coefficient matrix is the sum of two terms, a sparse matrix \(A\) and a possibly dense, rank deficient matrix of the form \(\gamma UU^T\), where \(\gamma > 0\) is a parameter which in some applications may be tak...

Full description

Saved in:
Bibliographic Details
Published in:arXiv.org 2022-11
Main Authors: Benzi, Michele, Faccio, Chiara
Format: Article
Language:English
Subjects:
Fluid mechanics
Iterative methods
Iterative solution
Linear systems
Optimization
Preconditioning
Sparse matrices
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We consider the iterative solution of large linear systems of equations in which the coefficient matrix is the sum of two terms, a sparse matrix \(A\) and a possibly dense, rank deficient matrix of the form \(\gamma UU^T\), where \(\gamma > 0\) is a parameter which in some applications may be taken to be 1. The matrix \(A\) itself can be singular, but we assume that the symmetric part of \(A\) is positive semidefinite and that \(A+\gamma UU^T\) is nonsingular. Linear systems of this form arise frequently in fields like optimization, fluid mechanics, computational statistics, and others. We investigate preconditioning strategies based on an alternating splitting approach combined with the use of the Sherman-Morrison-Woodbury matrix identity. The potential of the proposed approach is demonstrated by means of numerical experiments on linear systems from different application areas.
ISSN:2331-8422