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