Function sharp_padded 

fn sharp_padded(
    field: &[f64],
    mask: &[u8],
    grid: &Grid,
    threshold: f64,
    radius: f64,
) -> (Vec<f64>, Vec<u8>)

SHARP on QSM.jl’s padded box (steps 2-5 of the module docs).