" * where P is a permutation matrix, L is lower triangular with unit" . . . .