## 2008.56: Analysis of the Cholesky Decomposition of a Semi-definite Matrix

2008.56:
Nicholas J. Higham
(1990)
*Analysis of the Cholesky Decomposition of a Semi-definite Matrix.*
In:
M. G. Cox and S. J. Hammarling, (eds).
Reliable Numerical Computation.
Oxford University Press, Oxford, UK, pp. 161-185.
ISBN 0-19-853564-3

*This is the latest version of this eprint.*

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader 200 Kb |

## Abstract

Perturbation theory is developed for the Cholesky decomposition of an $n \times n$ symmetric positive semidefinite matrix $A$ of rank~$r$. The matrix $W=\All^{-1}\A{12}$ is found to play a key role in the perturbation bounds, where $\All$ and $\A{12}$ are $r \times r$ and $r \times (n-r)$ submatrices of $A$ respectively.

A backward error analysis is given; it shows that the computed Cholesky factors are the exact ones of a matrix whose distance from $A$ is bounded by $4r(r+1)\bigl(\norm{W}+1\bigr)^2u\norm{A}+O(u^2)$, where $u$ is the unit roundoff. For the complete pivoting strategy it is shown that $\norm{W}^2 \le {1 \over 3}(n-r)(4^r- 1)$, and empirical evidence that $\norm{W}$ is usually small is presented. The overall conclusion is that the Cholesky algorithm with complete pivoting is stable for semi-definite matrices.

Similar perturbation results are derived for the QR decomposition with column pivoting and for the LU decomposition with complete pivoting. The results give new insight into the reliability of these decompositions in rank estimation.

Item Type: | Book Section |
---|---|

Uncontrolled Keywords: | Cholesky decomposition, positive semi-definite matrix, perturbation theory, backward error analysis, QR factorization, rank estimation, LINPACK. |

Subjects: | MSC 2000 > 15 Linear and multilinear algebra; matrix theory MSC 2000 > 65 Numerical analysis |

MIMS number: | 2008.56 |

Deposited By: | Nick Higham |

Deposited On: | 19 November 2008 |

### Available Versions of this Item

- Analysis of the Cholesky Decomposition of a Semi-definite Matrix (deposited 19 November 2008)
**[Currently Displayed]**

Download Statistics: last 4 weeks

Repository Staff Only: edit this item