Mathematics – Numerical Analysis
Scientific paper
2010-03-22
Mathematics
Numerical Analysis
26 pages, 6 figures
Scientific paper
CG, SYMMLQ, and MINRES are Krylov subspace methods for solving symmetric systems of linear equations. When these methods are applied to an incompatible system (that is, a singular symmetric least-squares problem), CG could break down and SYMMLQ's solution could explode, while MINRES would give a least-squares solution but not necessarily the minimum-length (pseudoinverse) solution. This understanding motivates us to design a MINRES-like algorithm to compute minimum-length solutions to singular symmetric systems. MINRES uses QR factors of the tridiagonal matrix from the Lanczos process (where R is upper-tridiagonal). MINRES-QLP uses a QLP decomposition (where rotations on the right reduce R to lower-tridiagonal form). On ill-conditioned systems (singular or not), MINRES-QLP can give more accurate solutions than MINRES. We derive preconditioned MINRES-QLP, new stopping rules, and better estimates of the solution and residual norms, the matrix norm, and the condition number.
Choi Sou-Cheng T.
Paige Christopher C.
Saunders Michael A.
No associations
LandOfFree
MINRES-QLP: a Krylov subspace method for indefinite or singular symmetric systems does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with MINRES-QLP: a Krylov subspace method for indefinite or singular symmetric systems, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and MINRES-QLP: a Krylov subspace method for indefinite or singular symmetric systems will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-205879