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

A=LLA = LL^{*}

where LL is a lower triangular matrix with real and positive diagonal entries, and LL^{*} is the conjugate transpose of LL.