How to solve a linear system using LDL^T factorization

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 2K views
stunner5000pt
Messages
1,447
Reaction score
5
Use the [itex]LDL^T[/itex] factorization to solve the following linear system
[tex]\left(\begin{array}{cccc|1}4&1&-1&0&7\\1&3&-1&0&8\\-1&-1&5&2&-4\\0&0&2&4&6\end{array}\right)[/tex]
now i know ihow to get a matrix in the form LDL^T. But i was wondering how one would go about solving from there?



Please help!
 
Last edited:
Physics news on Phys.org
stunner5000pt said:
Use the [itex]LDL^T[/itex] factorization to solve the following linear system
[tex]\left(\begin{array}{cccc|1}4&1&-1&0&7\\1&3&-1&0&8\\-1&-1&5&2&-4\\0&0&2&4&6\end{array}\right)[/tex]
now i know ihow to get a matrix in the form LDL^T. But i was wondering how one would go about solving from there?
Please help!
That should be straight forward- the Cholesky decomposition is supposed to be the hard part! L here is a lower triangular matrix, LT is upper triangular, and D is diagonal, so going from LDLTX= A to DLTX= B is just a matter of "back substitution", starting from the value you get immediately in the last row and working up.
Since D is diagonal, going from DLTX= B to LTX= C is just dividing by the diagonal elements. Finally, since LT is upper triangular, going from LTX= C to X= D is again back substitution, this time working from the top row down.