P-adaptivity

It might be interesting to explore the implementation of adapting the order of the basis functions for each element individually, as sometimes seen in computational mechanics. This will require some changes to the dof distribution algorithm and the dof handler should also track the non-conformities here in some way.