Function star_qsm 

pub fn star_qsm(
    field: &[f64],
    mask: &[u8],
    grid: &Grid,
    bdir: (f64, f64, f64),
    params: &StarParams,
    progress: impl FnMut(usize, usize),
) -> Vec<f64>

STAR-QSM as in STI Suite 3.0 (QSM_star).

See the module docs for the algorithm. field is the local field in ppm, bdir the B0 direction in voxel axes (normalised here). progress is called after every SpaRSA iteration with (iteration, 1000).

Returns

χ in ppm on the input grid, zero outside the mask.

Panics

If params fails StarParams::validate (the pipeline checks this first and returns InvalidConfig), or field / mask do not match grid.