Add 3D Yukawa kernel

Does this also need some changes to the local/multipole expansion code? There's a Y2D{Local,Multipole}Expansion for the 2D case.