Function apply_reference 

pub fn apply_reference(
    chi: &[f64],
    mask: &[u8],
    method: QsmReference,
) -> Vec<f64>

Apply QSM referencing to a susceptibility map.

• Mean: subtract the mean of masked voxels, zero outside mask
• None: zero outside mask only