pub fn get_precond_alpha_k(
    a: &CsrMatrix<f64>,
    residue: &DVector<f64>,
    z: &DVector<f64>
) -> f64