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