LL^t and LDL^t variant. Added foward and backward substitution to solve linear systems with triangular matrix.