Function linear_ip::solve

pub fn solve<T, M, N>(
    A: &OMatrix<T, M, N>,
    b: &OVector<T, M>,
    c: &OVector<T, N>,
    eps: T,
    theta: T,
    gamma: T,
    max_it: usize
) -> (OVector<T, N>, OVector<T, M>)where
    T: Scalar + RealField + Copy,
    M: Dim + DimName + DimMin<M, Output = M> + DimSub<Const<1>>,
    N: Dim + DimName + DimMin<N, Output = N>,
    DefaultAllocator: Allocator<T, M> + Allocator<T, M, N> + Allocator<T, N, M> + Allocator<T, N, N> + Allocator<T, N> + Allocator<T, M, M> + Allocator<(T, usize), M> + Allocator<(usize, usize), M> + Allocator<T, <M as DimSub<Const<1>>>::Output>,

A is Dim(MxN), x,c,s are Dim(N), y,b are Dim(M) Returns (x (primal), y (dual))