Function ls_solver::cg_solve::compute_alpha

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