Function smv 

fn smv(
    x: &[f64],
    k: &[Complex64],
    fft: &mut Fft3dWorkspace,
    buf: &mut [Complex64],
)

SMV.m with a precomputed kernel: ifftn(fftn(x).*K), left in buf