Cholesky kl method (the second KL method)

@karlnapf 
The same at here. I use log_det in MatrixOperations.cpp 
seperate PR and see my other comments. This will not be possible to merge if the class is in all PRs.

So send the matrix operations thing first