pub fn compute_jacobi_update(
    x: &mut DVector<f64>,
    p: &CsrMatrix<f64>,
    residue: &DVector<f64>,
    support: &mut DVector<f64>
)