So far, the LinearPDESolver
(re)assembles the matrix in its
solve
function unconditionally. This might be unnecessary if the
matrix did not change compared to the previous call to solve
.
Make a first attempt to avoid the unnecessary calculation by
adding a corresponding control flag.