Improved stability of cholesky decomposition functions

Main improvement is to not use np.isclose() in the sparse_cholesky function. Its tolerance is not high enough!