Function dipole_kernel_shifted 

pub(crate) fn dipole_kernel_shifted(
    grid: &Grid,
    bdir: (f64, f64, f64),
) -> Vec<f64>

fftshift(get_dipole_kernel_fourier(...)) in column-major layout.