MinAres: An Iterative Solver for Symmetric Linear Systems

Abstract

We introduce an iterative solver named MINARES for symmetric linear systems Ax ≈ b, 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 MINRES, 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 \|xk - x\|, \|xk - x\|A and \|rk\| when A is positive definite.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…