* general iterative solvers (pcg, multigrid, ...)
* patchwise assembly of stiffness matrices
* default quadrature rules based on polynomial degrees
* improve handling of boundary conditions