minor fix
add the matrix operations class
the first KL method
Merge pull request #2295 from yorkerlin/fullDiagonalKLMethod

Full diagonal kl method (the first KL method)