Quick description
The LU factorization is perhaps the most common method for solving linear systems of equations. It is equivalent to Gaussian elimination, but preserves information so that new systems with the same matrix can be solved efficiently. This factorization has the form

where is lower triangular (that is,
if
), and
is upper triangular (that is
if
).
Usually it is used in a modified form where
is a permutation matrix (that is, all entries are zero or one, and only one entry in each row and in each column is non-zero).
Prerequisites
Linear algebra.
Example 1
To compute the LU factorization use recursion, which is essentially the algorithmic version of induction.
The case is easy:
.
If
![A = \left[\begin{array}{cc} A_1&b\\ c^T&d \end{array}\right]](/images/tex/ee50b1b0cbabfd37a197a5597ebaf7f1.png)
then we can recursively compute the LU factorization of .
Then an LU factorization of
can be given by
![A = \left[\begin{array}{cc} L_1 & \\ m^T&1\end{array}\right] \left[\begin{array}{cc} U_1 & v\\ &w\end{array}\right]](/images/tex/2f515730510324901e1e3df7e7ed65dc.png)
where ,
, and
, which can all be solved provided
is non-singular.
General discussion
This method can be modified to give a constructive proof of the existence of the LU factorization provided every principal (north-west) square submatrix of is non-singular.
This recursive construction creates an
matrix with 1's on the diagonal.
It can be easily modified to create a
matrix with 1's on the diagonal.
A different recursive construction can be given where we use the decomposition
![A = \left[\begin{array}{cc} \alpha&\beta^T\\ \gamma&A_1 \end{array}\right].](/images/tex/a71cf7585afb9d24c69b84639da0ccb7.png)
This gives equivalent computations, but done in a different order.
Comments
Title
Thu, 14/05/2009 - 10:49 — VickyCould this article perhaps have a more informative title? If you don't already know what LU factorization is, then the title gives you no clue where this tool might be helpful!