pub fn compute_gs_update(
    x: &mut DVector<f64>,
    p_lt: &CscMatrix<f64>,
    residue: &mut DVector<f64>
)