LU Decomposition, why?
There are at least three main advantages with LU decomposition compared with standard Gaussian elimination:
- It is straightforward to compute the determinant of a matrix
- If we have to solve sets of linear equations with the same matrix but with different vectors \( \mathbf{y} \), the number of FLOPS is of the order \( n^3 \).
- The inverse is such an operation