Function ls_solver::cg_solve::compute_beta

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