Extract the nodal values of a function on the boundary of the mesh

See for instance: Question on Dolfin poisson equation for source measurement on borders - #26 by dokken