A modified Cholesky algorithm based on a symmetric indefinite factorization

Cheng, Sheung Hun and Higham, Nicholas J. (1998) A modified Cholesky algorithm based on a symmetric indefinite factorization. SIAM Journal On Matrix Analysis And Applications, 19 (4). pp. 1097-1110. ISSN 1095-7162

PDF 30289.pdf Download (356kB)

Abstract

Given a symmetric and not necessarily positive definite matrix A, a modified Cholesky algorithm computes a Cholesky factorization P(A+E)PT = RT R, where P is a permutation matrix and E is a perturbation chosen to make A+E positive definite. The aims include producing a small-normed E and making A+E reasonably well conditioned. Modified Cholesky factorizations are widely used in optimization. We propose a new modified Cholesky algorithm based on a symmetric indefinite factorization computed using a new pivoting strategy of Ashcraft, Grimes, and Lewis. We analyze the effectiveness of the algorithm, both in theory and practice, showing that the algorithm is competitive with the existing algorithms of Gill, Murray, and Wright and Schnabel and Eskow. Attractive features of the new algorithm include easy-to-interpret inequalities that explain the extent to which it satisfies its design goals, and the fact that it can be implemented in terms of existing software.

Item Type:	Article
Uncontrolled Keywords:	modified Cholesky factorization, optimization, Newton's method, symmetric indefinite factorization
Subjects:	MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User:	Ms Lucy van Russelt
Date Deposited:	03 Jul 2006
Last Modified:	20 Oct 2017 14:12
URI:	https://eprints.maths.manchester.ac.uk/id/eprint/342

Actions (login required)

View Item