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

This somehow also fixes the problems with the higher order version of
the multilevel additive Schwarz preconditioner (this is tested now!)