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 + DimMin<M, Output = M> + DimSub<Const<1>>,
N: Dim + 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>,
Expand description
A is Dim(MxN), x,c,s are Dim(N), y,b are Dim(M) Returns (x (primal), y (dual))