Function dipole_forward 

pub(crate) fn dipole_forward(
    y: &[f64],
    kernel_shifted: &[f64],
    grid: &Grid,
) -> Vec<f64>

Dipole forward D(y) = real(ifft3(fft3(y) · kernel_shifted)).