User-agent:
Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.9.1.9) Gecko/20100423 Thunderbird/3.0.4
At the moment, a linear solver is not associated with a particular
matrix. This makes it awkward and dangerous to re-use preconditioners,
factorisations, etc. I suggest that we change the interface of the
linear solver classes to accept a matrix at construction, and that the
solvers maintain a (smart) pointer to the matrix. Opinions?