pub fn compute_precond_residue(
    residue: &mut DVector<f64>,
    a: &CsrMatrix<f64>,
    z: &DVector<f64>,
    alpha: f64
)