assemble matrix only once in the linear pde solver

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.