## ----knitr, echo=FALSE, results="hide"-----------------------------------
library("knitr")
opts_chunk$set(tidy=FALSE,dev="png",fig.show="hide",
               fig.width=4,fig.height=4.5,
               message=FALSE)

## ----style-knitr, eval=TRUE, echo=FALSE, results="asis"------------------
BiocStyle::latex()

## ------------------------------------------------------------------------
library("lpsymphony")

## ------------------------------------------------------------------------
obj <- c(3, 1, 3)

## ------------------------------------------------------------------------
mat <- matrix(c(-1, 0, 1, 2, 4, -3, 1, -3, 2), nrow = 3)
# mat <- matrix(c(-1,2,1,0,4,-3,1,-3,2), nrow = 3)

## ------------------------------------------------------------------------
dir <- c("<=", "<=", "<=")

## ------------------------------------------------------------------------
rhs <- c(4, 2, 3)

## ------------------------------------------------------------------------
max <- TRUE

## ------------------------------------------------------------------------
types <- c("I", "C", "I")

## ------------------------------------------------------------------------
lpsymphony_solve_LP(obj, mat, dir, rhs, types = types, max = max)

## ------------------------------------------------------------------------
bounds <- list(lower = list(ind = c(1,2,3), val = c(1,1,1)),
               upper = list(ind = c(1,2,3), val = c(5,5,5)))

## ------------------------------------------------------------------------
types <- c("C", "C", "C")
lpsymphony_solve_LP(obj, mat, dir, rhs, types = types, max = max, bounds = bounds)

## ------------------------------------------------------------------------
types <- c("I", "I", "I")
lpsymphony_solve_LP(obj, mat, dir, rhs, types = types, max = max, bounds = bounds)

## ------------------------------------------------------------------------
obj_dual <- rhs # rhs fom the Quick Start example

## ------------------------------------------------------------------------
mat_dual <- t(mat) # mat from the Quick Start example

## ------------------------------------------------------------------------
rhs_dual <- obj

## ------------------------------------------------------------------------
dir_dual <- c(">=", ">=", ">=")

## ------------------------------------------------------------------------
max_dual <- FALSE

## ------------------------------------------------------------------------
lpsymphony_solve_LP(obj_dual, mat_dual, dir_dual, rhs_dual, types = c("I","C","I"), max = max_dual)