feature/order2C1space

Implement lowest order C^1 space (order=2). In addition to some other minor changes this adds a laplace constraint to the value projection because otherwise we have fewer dofs than the dimension of P^3: int laplace Pi^0 v = int_boundary v.n

Builds on div-free space and requires dune-fem feature/improvePythonQuadPoints branch.