Function pdf_default 

pub fn pdf_default(
    field: &[f64],
    mask: &[u8],
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f64,
    vsy: f64,
    vsz: f64,
) -> Vec<f64>

PDF with default parameters (adaptive max_iter matching QSM.jl: ceil(sqrt(nxnynz)))