SDDiP

SDDiP is a version of SDDP for integer local or state variables. If a model includes integer variables, we cannot compute cut coefficients in the way we normally do in SDDP (by querying duals from an LP solver). In SDDiP, each stage problem is solved via its Lagrangian dual in the backward pass.

The implementation of SDDiP in this package comes from the paper by Zhou, J., Ahmed, S., Sun, X.A. (2016): Nested Decomposition of Multistage Stochastic Integer Programs with Binary State Variables.

Installation

Pkg.clone("https://github.com/lkapelevich/SDDiP.jl")

This package uses @odow's SDDP package, SDDP.jl:

Pkg.clone("https://github.com/odow/SDDP.jl")

Usage

Note there are some examples in the examples folder.

Calling the solver

Using SDDP.jl, a user should include a call to setSDDiPsolver! at the end of a stage problem definition:

setSDDiPsolver!(sp::JuMP.Model; method=Subgradient(0.0), pattern=Pattern(), MIPsolver=sp.solver, LPsolver=MIPsolver)

The keyword arguments are as follows.

method: An AbstractLagrangianMethod object from that defines how the Lagrangian will be solved. See the readme in src/Lagrangian.
pattern: Allows us to choose to add Lagrangian cuts, Benders cuts, or strengthened Benders cuts every so often. See Pattern below.
MIPsolver: The solver that will be used every time a MIP is solved.
LPsolver: The solver that is used ever time an LP is solved.

Declaring state variables

We can declare state variables using the @binarystate macro (like @state).

There are four required arguments. The first three come from @state:

A stage problem object,
State at the end of a stage (written like a JuMP variable),
State at the start of a stage.

The fourth argument is special for SDDiP:

Type of variable: binary (Bin), integer (Int), or continuous (Cont).

For example:

@binarystate(sp, 1 <= x[i=1:3] <= 10, x0 == ones(3)[i] , Int)

Although our state may be described by variables that are binary, integer, or continuous, only binary variables are ever registered as state variables in the SDDP model.

We will compute the binary expansion of integer and continuous variables (as described in Zhou, J., et. al. (2016)).

For continuous variables, we can also specify a precision for binary expansion as an optional fifth argument.

For example:

@binarystate(sp, 1 <= x[i=1:3] <= 10, x0 == ones(3)[i] , Cont, 0.01)

If a precision is not specified, the default precision is 0.1.

Pattern

The Pattern function has the definition:

Pattern(;benders=0, strengthened_benders=0, lagrangian=1)

with, for example,

Pattern(benders=0, strengthened_benders=1, lagrangian=4)