Draft: Add moment dofs for all variations of PETSCDUALSPACELAGRANGE

With trimmed polynomial evaluation, we can now implement dualspaces for all the femtable elements that are modal (integrated moments) as opposed to nodal.

• The dual spaces have the same symmetries as before, because the moments are integrated against nodal polynomials. This means that the dual spaces are not nested (e.g. the moment dofs for third order elements are not contained in the moment dofs for fourth order elements)
• The method in dspacelagrange.c of building up dual spaces by kronecker products and pushing forward boundary dual spaces still works without a lot of changes, but those dof building routines are slow for dense matrices because the make heavy use of MatGetRow().