Let local GFE functions store a basis

This is the basis used as interpolation weights (for geodesic
FE functions) or for interpolating in the embedding space
(for projection-based FE functions).  The basis is not actually
used yet, but it needs to be given to the class in the
constructor.