Cholesky decomposition
comp_chol(C, M)matrix to be filled out with the cholesky of M
symmetric positive definite matrix for which the upper triangle of the cholesky is to be computed
small positive constant to be added to the diagonal for numerical stability
set C equal to the matrix square root of M