support for solving lower triangular system by forward substitution.
Several checks for detection of singularity in triangular matrices.
basic implementation of solving LDU x = b problem.