Use PDELab::interpolate
The current code contains a hand-coded interpolation method that does interesting trickery to work for dimensions 1, 2, 3. As an alternative, this patch simply offloads the work to dune-pdelab, which does the right thing off-the-shelf.
Edited by Sander, Oliver