Function compute_alpha 

pub fn compute_alpha(
    d: &DVector<f64>,
    residue: &DVector<f64>,
    a: &CsrMatrix<f64>,
    support: &mut DVector<f64>,
) -> f64