R/linear.constraint.iq.R

In rhochreiter/scenportopt: Scenario-based Portfolio Optimization

# linear constraint: A(range) * factors <= b
linear.constraint.iq <- function(constraints.linear, range, b, factors=NULL, leq=TRUE) {
  add_a <- rep(0, constraints.linear$n)
  if (is.null(factors)) { add_a[range] <- 1 } else { add_a[range] <- factors }
  if (!leq) {
   add_a <- -1 * add_a
   sum <- -1 * sum
  }
  if (!is.null(constraints.linear$A)) { constraints.linear$A <- rbind(constraints.linear$A, add_a) } else { constraints.linear$A <- add_a }
  if (!is.null(constraints.linear$b)) { constraints.linear$b <- c(constraints.linear$b, b) } else { constraints.linear$b <- b }
  return(constraints.linear)
}