×

Modified Gram-Schmidt-based methods for block downdating the Cholesky factorization. (English) Zbl 1213.65050

Let \(A\) be an \(n\times n\) symmetric positive definite matrix and \(A = R^T R\) its Cholesky factorization, that is, \(R\) is upper triangular and has positive diagonal entries. Consider the modified matrix \(\tilde A = A - X^T X = R^T R - X^T X\), where \(X\) is an \(t\times n\) matrix, and assume that \(\tilde A\) remains positive definite. The paper is concerned with computing the Cholesky factorization of \(\tilde A\) by appropriately downdating the Cholesky factorization of \(A\) rather than computing the factorization of \(\tilde A\) explicitly from scratch. A hyperbolic modified Gram-Schmidt method is proposed for this purpose. Several numerical experiments with random matrices and an academic example illustrate the accuracy properties of the newly proposed method.

MSC:

65F05 Direct numerical methods for linear systems and matrix inversion
15B48 Positive matrices and their generalizations; cones of matrices
15B52 Random matrices (algebraic aspects)

Software:

LINPACK
Full Text: DOI

References:

[1] Bojanczyk, A. W.; Brent, R. P.; Van Dooren, P.; De Hoog, F. R., A note on downdating the Cholesky factorization, SIAM J. Sci. Stat. Comput., 8, 210-221 (1987) · Zbl 0617.65013
[2] Bjorck, A.; Park, H.; Elden, L., Accurate downdating of least squares solutions, SIAM J. Matrix Anal. Appl., 15, 549-568 (1994) · Zbl 0811.65034
[3] Barlow, J. L.; Smoktunowicz, A.; Erbay, H., Improved Gram-Schmidt type downdating methods, BIT, 45, 259-285 (2005) · Zbl 1087.65550
[4] G. Stewart, On the stability of sequential updates and downdates, Report TR-3238, Department of Computer Science, University of Maryland, College Park, MD 20742, March 1994.; G. Stewart, On the stability of sequential updates and downdates, Report TR-3238, Department of Computer Science, University of Maryland, College Park, MD 20742, March 1994.
[5] Yoo, K.; Park, H., Accurate downdating of a modified Gram-Schmidt QR decomposition, BIT, 36, 1, 166-181 (1996) · Zbl 0848.65014
[6] Eldén, L.; Park, H., Block downdating of least squares solutions, SIAM J. Matrix Anal. Appl., 15, 3, 1018-1034 (1994) · Zbl 0808.65041
[7] Rader, C. M.; Steinhardt, A. O., Hyperbolic Householder transformations, SIAM J. Matrix Anal. Appl., 9, 269-290 (1988) · Zbl 0647.65027
[8] Bojanczyk, A. W.; Steinhardt, A. O., Stabilized hyperbolic Householder transformations, IEEE Trans. Acoust. Speech Signal Processing, ASSP-37, 1286-1288 (1989) · Zbl 0697.93064
[9] Bojanczyk, A. W.; Steinhardt, A., Stability analysis of a Householder-based algorithm for downdating the Cholesky factorization, SIAM J. Sci. Stat. Comput., 12, 6, 1255-1265 (1991) · Zbl 0737.65013
[10] Björck, A.; Paige, C. C., Loss and recapture in the modified Gram-Schmidt algorithm, SIAM J. Matrix Anal. Appl., 13, 176-190 (1992) · Zbl 0747.65026
[11] Björck, A., Numerical Methods for Least Squares Problems (1996), SIAM: SIAM Philadelphia, PA, USA · Zbl 0847.65023
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.