checkPolytopeStatus: checkPolytopeStatus check if the polytope has solution or not

In Irstea/cpgsR: Uniform sampling in convex polytope

Description

checkPolytopeStatus check if the polytope has solution or not

Usage

checkPolytopeStatus(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

print a message to tell if the polygon is ok or not

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)
X0 <- checkPolytopeStatus(A,b)

#this one is empty
C <- matrix(c(1,rep(0,n-1)),1)
v <- 3
X0 <- checkPolytopeStatus(A,b,C,v)

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