decomposition of a Hermitian, positive-definite matrix into the product of a lower triangular matrix and its conjugate transpose. (used for Monte-Carlo simulations)

where is a lower triangular matrix with real and positive diagonal entries, and is the conjugate transpose of .