Function compute_dipole_stencil 

pub fn compute_dipole_stencil(
    res: (f32, f32, f32),
    b0_dir: (f32, f32, f32),
) -> [[[f32; 3]; 3]; 3]

Compute dipole stencil (27-point spatial kernel)