{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

Cholesky

Cholesky - LDL t factorization for spd matrices The...

This preview shows page 1. Sign up to view the full content.

The Cholesky Factorization Symmetric matrices are important because they are common in applications, have some very nice properties, and because the symmetry can be exploited by algorithms to save time and memory. For example, we know that if A = A t has an LU factorization, then A = LDL t can be computed in about n 3 / 3 flops. Because of small pivots, the LDL t algorithm is not backward stable for general symmetric matrices. However, there are matrices for which pivoting is never needed: for which the diagonal element is always the | largest | element in its column. The most important of these are the symmetric positive definite (spd) matrices. A R n × n is spd if it is symmetric and if for all x 6 = 0 R n , x t Ax > 0. All matrices of the form X t X are spd iff X has linearly independent columns. A matrix has an LDL t factorization with d ii > 0 iff it is spd. In that case we may write A = LDL t = LD 1 / 2 D 1 / 2 L t GG t . The Cholesky factorization , which computes A = GG t directly, is a simple and popular alternative to the
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: LDL t factorization for spd matrices. The classical algorithm appears easily by looking at the columns of A = GG t with lower triangular G . Assume we know the ﬁrst k-1 columns of G , and look at the k th column of A = GG t : a k ≡ Ae k = GG t e k = [ g 1 ,g 2 ,...,g n ] z. Here z t is the k th row of (lower triangular) G , so we can write a k = k X i =1 g ki g i , or g kk g k = a k-k-1 X i =1 g ki g i . In particular g 2 kk = a kk-∑ k-1 j =1 g 2 kj , so we take the positive root and solve for g k : g k = ( a k-k-1 X i =1 g ki g i ) /g kk . This method runs to completion (no zero or complex roots) iﬀ A is spd. It requires 1 3 n 3 + O( n 2 ) ﬂops, and no extra memory if the lower triangular of A is overwritten by that of G . With respect to rounding errors, the computed ˜ G satisﬁes ˜ G ˜ G t = A + δA, where k δA k ≤ 12 n 2 μ k A k ....
View Full Document

{[ snackBarMessage ]}

Ask a homework question - tutors are online