Allocation free shapefunctions for GWP and Largange Spaces

[cmuenger]

• Rework of shape function evaluation, avoiding allocation of a workspace during the evaluation of higher-order polynomials. (preparation for compatibility with GPU)
• Unifying code for Lagrange global and local spaces into a single evaluation function for arbitrary order. (2D meshes/triangle only).

Breaking changes:

• Indexing/ordering of global and local Lagrange functions of order 1 and 2. Now consistent with higher-order spaces.
• By default, DoFs will be placed everywhere on the supplied mesh during the construction of the Lagrange space, if no restrictions are explicitly given by the user. No longer implicitly assumes the Dirichlet boundary condition on open meshes.