pub fn sparse_direct_solve<F>( a: &CsrMatrix<F>, b: &[F], _symmetric: bool, _positive_definite: bool, ) -> SparseResult<Vec<F>>where F: Float + NumAssign + Sum + 'static + Debug,
Solve a sparse linear system using direct methods