Do not cache both K and its LU decomposition

[mancellin]

Since #235, the default setting of the code is to cache the matrices $S$ and $K$ and the LU decomposition of $K$. It means a 50% increase in the RAM usage, which can be an issue for larger problems.

It would make more sense to keep only the LU decomposition of $K$.