Function run_iqsm_plus 

pub fn run_iqsm_plus(
    phases: &[&[f64]],
    magnitudes: &[&[f64]],
    mask: &[u8],
    metadata: &ScanMetadata,
    reference: QsmReference,
    tile: Option<(usize, usize)>,
) -> Result<Vec<f64>, PipelineError>

Run iQSM+ single-step reconstruction from wrapped phase (orientation-adaptive variant; the B0 direction is a genuine network input). See run_iqsm; the echoes are combined with (magnitude·TE)² weights, as in both of the authors’ implementations.