Correct construction of lagrange function space ready for support of curvilinear meshes

Some parts of lagrange.jl use the skeleton function under the assumption that it provides the nodes of the mesh, whereas it actually would return the vertices including the control vertices. Currently, we only support linear elements, but it becomes a problem when curvilinear elements are used.

The implementation of the GWP function could serve as a template.