Function tgv_qsm_with_progress 

pub fn tgv_qsm_with_progress<F>(
    phase: &[f32],
    mask: &[u8],
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f32,
    vsy: f32,
    vsz: f32,
    params: &TgvParams,
    b0_dir: (f32, f32, f32),
    progress: F,
) -> Vec<f32>
where
    F: Fn(usize, usize),

TGV-QSM with progress callback (optimized version)