Suppose we have a covariance matrix $\Sigma$. I know that the Cholesky decomposition $A^T A$ can be found from the LDL decomposition using $$ \Sigma = LDL^T = (LD^{\frac 1 2})(LD^{ \frac 1 2 })^T = A^TA $$.

I was wondering if you can go the other way -- find the LDL decomposition from the Cholesky, I couldn't see how to do it!


1 Answer 1


Let $S=D^{\frac12}$, so that $S$ is a diagonal matrix with diagonal $s$.

Note that when you compute $A=LS$, multiplying $L$ by $S$ multiplies each column of $L$ by the corresponding element of $s$, and that $L_{ii}=1$. So you need to find the $s_i$ that would give $A_{ii}/s_i=1$, and then divide every element of the $i$th column of $A^\top$ (which is the $i$th row of $A$) by $s_i$ to give the $i$th column of $L$.

Then the diagonal of $D$ is trivial to compute from $S$.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.