improve interpolate when it needs to be called locally many times

Currently, interpolate! and the subroutines that e.g. compute the moments of the interpolated functions for an H1Pk interpolation are only efficient when called once for all cells at once. When it needs to be called many times on smaller sets (like in the new periodic boundary functionality from PRs #31 and #34) this is not efficient and causes lots of allocations.