Function run_iqsm_plus 

pub fn run_iqsm_plus(
    phases: &[&[f64]],
    magnitudes: &[&[f64]],
    mask: &[u8],
    metadata: &ScanMetadata,
    reference: QsmReference,
) -> 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.