Added function Matrix_CholeskyDouble from Christian's version of Matrix_utils.f90.