How to do LU decomposition with unit lower triangular matrix L, in case a decomposed matrix has zeros on diagonal? This is obviously possible for positive defined matrix.
For example suppose instead of 4 we have 0 here:
then u11 = 0 and determinant of the right part = 0.
What is limitations of LU decomposition?