Closes #3.

This feature must be used with care: if the internal structure of a
linear solver is already allocated, then the symbolic factorization is
automatically reused, even if the sparsity pattern are not the same !