add the multiLaplace method

@karlnapf 
take a look at this.
could you investigate a big (using debug messages) whether things are computed multiple times?
We had the situation with @votjakovr where the GP inference was effectively computed twice since some minor parameters changed and then update re-computed all matrix inversion.