pub fn solve<M: MatrixTrait + LinearAlgebra<M>>(A: &M, b: &M, sk: SolveKind) -> M {
    A.solve_mat(b, sk)
}
Data structure for Complete Pivoting LU decomposition

# Usage
```rust
use peroxide::fuga::*;

let a
src/traits/matrix.rs