Function dipole_apply 

fn dipole_apply(v: &[f64], dk: &[f64], grid: &Grid) -> Vec<f64>

D(v) = real(ifft3(fft3(v)·dk)) with dk in fftfreq (unshifted) layout.