Accuracy of preconditioned CG-type methods for least squares problems. (Q1416375)
From MaRDI portal
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Accuracy of preconditioned CG-type methods for least squares problems. |
scientific article |
Statements
Accuracy of preconditioned CG-type methods for least squares problems. (English)
0 references
14 December 2003
0 references
The conjugate gradient (CG) method with IMGS, an incomplete modified version of Gram-Schmidt orthogonalization to obtain an incomplete orthogonal factorization preconditioner, applied to the normal equations (PCGLS) is often used as the basic iterative method to solve the linear least squares problems. In this paper, a detailed analysis is given for understanding the effect of rounding errors on IMGS and determining the accuracy of computed solutions of PCGLS with IMGS for linear least squares problems in finite precision. It is shown that for a consistent system, the difference between the true residuals and the updated approximate residual vectors generated depends on the machine precision \(\varepsilon\), on the maximum growth in norm of the iterates over their initial values, the norm of the true solution, and the condition number of \(R\) which is affected by the drop set in incomplete Gram-Schmidt factorization. Similar results are obtained for the difference between the true and computed solution for inconsistent systems. Numerical tests are carried out to confirm the theoretical conclusions.
0 references
Least square problems
0 references
Incomplete modified Gram-Schmidt orthogonalization
0 references
conjugate gradient-type method
0 references
Rounding error analysis
0 references
Finite precision arithmetic
0 references
preconditioning
0 references
numerical examples
0 references
incomplete orthogonal factorization
0 references
incomplete Gram-Schmidt factorization
0 references
0 references
0 references
0 references