Implementation of a Phi-Divergence Two-Stage Linear Program with Recourse (PhiLP-2) for Tucson water data. This project originally started life as an implementation fo the Likelihood Robust Two-Stage Linear Program with Recourse, but was later generalized to include other phi-divergences as well.

Classes are as follows:

• @LPMODEL contains the structure of a two-stage stochastic linear program with recourse.
• @LRLP implements the modified Bender's Decomposition used to solve the PhiLP-2 model, formerly the Likelihood Robust model.
• @PhiDivergence is a class to contain the phi-divergences and their important properties.
• @Solution class is used to store the candidate and best solution as the modified Bender's Decomposition algorithm runs.

A few functions of interest:

• SolveLRLP.m implements the basic Bender's Decomposition loop using the LRLP class
• GatherLROData.m loops over a variety of parameter values to collect data on optimal solution of the PhiLP-2 problem.

=====

DR. LANSEY'S STUDENTS:

For the code to build an LP based on the excel spreadsheets, you want the LPModel class, contained in the folder @LPModel

The file OneStageAZ.m shows you the basic use of the LPModel class.

The doc folder contains some technical documentation on how LPModel constructs the water model.

I currently have LPModel set to never write to Excel (it doesn't work on Linux). To enable writing to excel, open @LPModel/LPModel.m, and set the variable writeToExcel = true. This is line 52 as of the time I'm writing this.

For reference, the files in LPModel:

• LPModel.m contains the basic information about LPModel--data and some small methods.
• BuildA.m builds the submatrices A, A_st and A_log. Alicia can tell you all about those.
• ExpandA.m creates the full constraint matrix out of A, A_st and A_log.
• BuildVectors.m Builds the cost, right hand side, and lower and upper bound vectors. It should work for both one and two stage problems. NOTE BuildVectors will only work for Time_Lag = 1 or Time_Lag = 12. This is based on what I see in sf.xlsx
• GenerateDiagonal.m is used by ExpandA.m to build the full constraint matrix. Don't modify this file.
• PrintConstraint.m prints the constraints in a human readable format, which should be helpful for debugging any changes to the code, or creating new spreadsheets.