Full diagonal kl method (the first KL method)

@karlnapf 
take a look at this.
Currently I use the log_det in MatrixOperations.cpp.
 I will use some @lambday's work once other things in the implementation are fine.
Please provide full details in the doxygen for people who read documentation (also in the base class)