generalize 'interpolate_current' to work for 1D to 3D

In principle we could write it in a *completely* general way, but this
would involve numerically calculating the surface elements and
integrating the smoothing function. Given that the only useful cases are
1D, 2D, and 3D it makes sense to just hard-code the prefactors for these
cases and be done with it.