Merged requested to merge feature/assemble-matrix-only-once into master
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
Make a first attempt to avoid the unnecessary calculation by
adding a corresponding control flag.