Parametrize Lagrange bases with LocalFiniteElement Cache
Currently there are multiple implementations of Lagrange local-functions and correspondingly cache implementations to generate these local finite-element types for a given GeometryType. The LagrangeBasis, on the other hand, does not support to choose a flexible local finite-element type but just allows to switch between compiletime and runtime polynomial order.
This is very restrictive. If I have another Lagrange local basis, e.g. other positions of the lagrange nodes/Lobatto nodes, or just a different order of the nodes - I have no way to pass this implementation type.
I have done a proof-of-concept implementation to change the parametrization of LagrangeBasis from a simple polynomial degree to a LocalFunctionCache type. See !248
Any comments on this idea?