findingIncompatibleConstr: findingIncompatibleConstraints help to dectect constraints...

In Irstea/cpgsR: Uniform sampling in convex polytope

Description

findingIncompatibleConstraints help to dectect constraints that leads to an empty polytope

Usage

findingIncompatibleConstr(A, b, C = NULL, v = NULL)

Arguments

A	the matrix of inequality A.x<=b
b	the vector A.x<=b
C	the matrix of equality C.x=v (default NULL for no equality)
v	the vector of equality C.x=v (default NULL for no equality

Value

a list of vector. To be fitted, constraints corresponding to the first elements of each vector of the list should be relaxed. Those constraints appear to be incompatible with following constraints in the Vector. For example, list(c(A,B),c(D,E, F)) means that A and B should be relaxed and that A is incompatible with B and that D is incompatible with both E and F

Examples

n <- 20
A1 <- -diag(n)
b1 <- as.matrix(rep(0,n))
A2 <- diag(n)
b2 <- as.matrix(rep(1,n))
A <- rbind(A1,A2)
b <- rbind(b1,b2)
C <- matrix(c(1,rep(0,n-1)),1)
v <- 3
X0 <- findingIncompatibleConstr(A,b,C,v)

Irstea/cpgsR documentation built on Jan. 25, 2020, 5:36 p.m.