Add the implementation of memory-saving mode (for Cholesky only).

When the user has specified a matrix is SPD, he can also choose
whether to solve the problem in memory saving mode or not. This
is a little bit slower, though.