Evaluate a function on a set of vertices

variational formulation
2

See: Mapping from mesh.coordinates() to u.vector() - #2 by dokken