Function dipole_kernel_default 

pub fn dipole_kernel_default(
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f64,
    vsy: f64,
    vsz: f64,
) -> Vec<f64>

Generate dipole kernel with default B direction (0, 0, 1)