Implement 'setup' for all lagrange element transfer operators

Indead of limiting this to P1 elements, we can now build transfer matrices for
lagrange elements of any order. This is achived by a compile time parameter "k".

For back compatibility, it has a default value k = 1.

To determine the size of the matrix, we need a global basis from dune-functions.