pub fn compute_rel_err(solution: &DVector<f64>, x: &DVector<f64>) -> f64