pub fn spsolve<F>(a: &CsrMatrix<F>, b: &[F]) -> SparseResult<Vec<F>>where F: Float + NumAssign + Sum + 'static + Debug,
Solve a sparse linear system Ax = b