Function chi_sep_medi 

pub fn chi_sep_medi(
    local_field: &[f64],
    r2prime: &[f64],
    magnitude: &[f64],
    mask: &[u8],
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f64,
    vsy: f64,
    vsz: f64,
    bdir: (f64, f64, f64),
    cf: f64,
    lambda_para: f64,
    lambda_dia: f64,
    lambda_cpl: f64,
    dr_pos: f64,
    dr_neg: f64,
    percentage: f64,
    cg_tol: f64,
    cg_max_iter: usize,
    max_iter: usize,
    tol: f64,
) -> (Vec<f64>, Vec<f64>, Vec<f64>)

Chi-separation using MEDI-based coupled optimization.

Arguments

• local_field - Local field map in Hz
• r2prime - R2’ map in Hz
• magnitude - Magnitude image for edge weighting
• mask - Binary brain mask, 1 = brain
• cf - Central frequency in Hz (e.g. 123.2e6 for 3T)
• dr_pos - Paramagnetic relaxivity in Hz/ppm (default: 114.0)
• dr_neg - Diamagnetic relaxivity in Hz/ppm (default: 30.0)

Returns

(chi_pos, chi_neg, chi_total) — susceptibility maps in ppm