pub fn init_support_vectors(
    method: &Method,
    a: &CsrMatrix<f64>,
    x: &DVector<f64>,
    b: &DVector<f64>
) -> (DVector<f64>, DVector<f64>)