Constructor computes the Cholesky decomposition of sigma in place without memory allocation. Furthermore it is assumed to be a symmetric matrix, but only the lower/left half is actually accessed.
Compute Cholesky decomposition in place. Only accesses lower/left half of the matrix. Returns false if the matrix is not positive definite.
Multivariate normal distribution. Beta version (has not properly tested).