Optimize memory usage for stiff time integrators #184
Description
Most stiff (ode) time integrators need compatible matrices for M,J and W. Usually these have the same type and the same sparsity pattern. Hence we can share data which is currently duplicated. Hence, we should add constructor for all assembled operators to construct matrices this way.