Add cholesky based version

- fixed small error in formula for K+sigma^2*I
- minor code beautification for output