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.