Function ls_solver::pgr_solve::compute_z

pub fn compute_z(
    z: &mut DVector<f64>,
    p: &CsrMatrix<f64>,
    residue: &DVector<f64>
)