Draft: Lobatto shape functions

Summary

An (adaptive) Lobatto basis is a set of basis functions of local Lobatto shape functions that allows to have different polynomial orders on each element. It even allows to have different orders on sub-entities and thus opens the possibility for conforming bases with p-adaptivity.

Note, for a stable discretization you need to fulfill a minimum-rule that essentially means that the polynomial degree on sub-entities shared by some entities is a most the minimum of the polynomial degrees of the entities.

The implementation is based on the local Lobatto[Cube]LocalFiniteElement implementing the functions from

"Higher-Order Finite Element Methods", P. Soling, K, Segeth, I. Dolezel, 2004, Chapman & Hall/CRC