The fundamental idea of finite element methods is to introduce a variational formulation of the considered boundary value problem for a partial differential equation and then to replace the infinite-dimensional spaces by finite-dimensional spaces. To this end, triangulations of the domain are considered. Finite element functions are usually piecewise polynomials, i.e., the restriction to a mesh cell is either a polynomial or a function defined by mapping a polynomial from a reference cell.