Retour

G-2023-40

MINARES: An iterative solver for symmetric linear systems

, et

référence BibTeX

We introduce an iterative solver named MINARES for symmetric linear systems Axb, where A is possibly singular. MINARES is based on the symmetric Lanczos process, like MINRES and MINRES-QLP, but it minimizes Ark in each Krylov subspace rather than rk, where rk is the current residual vector. When A is symmetric, MINARES minimizes the same quantity Ark as LSMR, but in more relevant Krylov subspaces, and it requires only one matrix-vector product Av per iteration, whereas LSMR would need two. Our numerical experiments with MINRES-QLP and LSMR show that MINARES is a pertinent alternative on consistent symmetric systems and the most suitable Krylov method for inconsistent symmetric systems. We derive properties of MINARES from an equivalent solver named CAR that is to MINARES as CR is to MINARES, is not based on the Lanczos process, and minimizes Ark in the same Krylov subspace as MINARES. We establish that MINARES and CAR generate monotonic xkx, xkxA and rk when A is positive definite.

, 18 pages

Axes de recherche

Application de recherche

Document

G2340.pdf (860 Ko)