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.