/********************************************************
Xpress-BCL Java Example Problems
================================
file xbcoco1.java
`````````````````
Coco Problem Phase 1.
Initial formulation: data, variables and constraints fixed.
(c) 2008 Fair Isaac Corporation
author: S.Heipcke, Jan. 2000, rev. Mar. 2011
********************************************************/
import com.dashoptimization.*;
public class xbcoco1 {
public static void main(String[] args) {
XPRBvar make11,make21,make12,make22;
try (XPRBprob p = new XPRBprob("Coco1")) { /* Initialize BCL and create a new problem */
/****VARIABLES****/
make11 = p.newVar("make11"); /* Amount of prod. 1 to make at factory 1 */
make21 = p.newVar("make21"); /* Amount of prod. 2 to make at factory 1 */
make12 = p.newVar("make12"); /* Amount of prod. 1 to make at factory 2 */
make22 = p.newVar("make22"); /* Amount of prod. 2 to make at factory 2 */
/****OBJECTIVE****/
/* Define & set objective function: maximize total profit */
p.setObj(make11.mul(50).add(make21.mul(125)).add(make12.mul(47))
.add(make22.mul(132)) );
/****CONSTRAINTS****/
p.newCtr("MxMake1", make11.add(make21).lEql(400) );
/* Capacity limit at factory 1 */
p.newCtr("MxMake2", make12.add(make22).lEql(500) );
/* Capacity limit at factory 2 */
p.newCtr("MxSell1", make11.add(make12).lEql(650) );
/* Limit on the amount of prod. 1 to be sold */
p.newCtr("MxSell2", make21.add(make22).lEql(600) );
/* Limit on the amount of prod. 2 to be sold */
/****SOLVING + OUTPUT****/
p.setSense(XPRB.MAXIM); /* Choose the sense of the optimization */
p.lpOptimize(""); /* Solve the LP-problem */
System.out.println("Objective: " + p.getObjVal()); /* Get objective value */
/* Print out the solution values for all variables */
System.out.println(make11.getName() +": " + make11.getSol());
System.out.println(make21.getName() +": " + make21.getSol());
System.out.println(make12.getName() +": " + make12.getSol());
System.out.println(make22.getName() +": " + make22.getSol());
}
}
}
|
/********************************************************
Xpress-BCL Java Example Problems
================================
file xbcoco2.java
`````````````````
Coco Problem Phase 2.
Use parameters, data tables and subscripted variables
to separate the model structure from the data.
Read data tables in from text data files.
(c) 2008 Fair Isaac Corporation
author: S.Heipcke, Jan. 2000, rev. Mar. 2011
********************************************************/
import java.io.*;
import com.dashoptimization.*;
public class xbcoco2 {
static final int NP = 2; /* Number of products (p) */
static final int NF = 2; /* factories (f) */
static final int NR = 2; /* raw materials (r) */
static final String REVFILE = System.getProperty("XPRBDATA") +
"/coco/rev.dat";
static final String CMAKEFILE = System.getProperty("XPRBDATA") +
"/coco/cmake.dat";
static final String CBUYFILE = System.getProperty("XPRBDATA") +
"/coco/cbuy.dat";
static final String REQFILE = System.getProperty("XPRBDATA") +
"/coco/req.dat";
static final String MXSELLFILE= System.getProperty("XPRBDATA") +
"/coco/maxsell.dat";
static final String MXMAKEFILE= System.getProperty("XPRBDATA") +
"/coco/mxmake.dat";
/****TABLES****/
static double[] REV; /* Unit selling price of product p */
static double[][] CMAK; /* Unit cost to make product p at factory f */
static double[] CBUY; /* Unit cost to buy raw material r */
static double[][] REQ; /* Requirement by unit of prod. p for raw mat. r */
static double[] MXSELL; /* Max. amount of p that can be sold */
static double[] MXMAKE; /* Max. amount factory f can make over all prod.s */
static double[][] PROFIT; /* Profit contribution of product p at factory f */
/***********************************************************************/
static void modCoco2() {
XPRBvar[][] make;
XPRBexpr lobj, lc;
int p,f;
try (XPRBprob pb = new XPRBprob("Coco2")) { /* Initialize BCL and create a new problem */
/****VARIABLES****/
make = new XPRBvar[NP][NF];
for(p=0;p<NP;p++) /* Amount of prod. p to make at factory f */
for(f=0;f<NF;f++)
make[p][f] = pb.newVar("make_"+ (p+1) +"f"+ (f+1));
/****OBJECTIVE****/
lobj = new XPRBexpr();
for(p=0;p<NP;p++) /* Objective: maximize total profit */
for(f=0;f<NF;f++) lobj.add(make[p][f].mul(PROFIT[p][f]));
pb.setObj(lobj);
/****CONSTRAINTS****/
for(p=0;p<NP;p++) { /* Limit on the amount of product p to be sold */
lc = new XPRBexpr();
for(f=0;f<NF;f++) lc .add(make[p][f]);
pb.newCtr("MxSell", lc.lEql(MXSELL[p]) );
}
for(f=0;f<NF;f++) { /* Capacity limit at factory f */
lc = new XPRBexpr();
for(p=0;p<NP;p++) lc.add(make[p][f]);
pb.newCtr("MxMake", lc.lEql(MXMAKE[f]) );
}
/****SOLVING + OUTPUT****/
pb.setSense(XPRB.MAXIM); /* Choose the sense of the optimization */
pb.lpOptimize(""); /* Solve the LP-problem */
System.out.println("Objective: " + pb.getObjVal()); /* Get objective value */
for(p=0;p<NP;p++) /* Print the solution values */
for(f=0;f<NF;f++)
System.out.print(make[p][f].getName()+ ":"+ make[p][f].getSol() +" ");
System.out.println();
}
}
/***********************************************************************/
/**** Initialize the stream tokenizer ****/
static StreamTokenizer initST(FileReader file) {
StreamTokenizer st=null;
st= new StreamTokenizer(file); /* Initialize the stream tokenizer */
st.commentChar('!'); /* Use the character '!' for comments */
st.eolIsSignificant(true); /* Return end-of-line character */
st.ordinaryChar(','); /* Use ',' as separator */
st.parseNumbers(); /* Read numbers as numbers (not strings)*/
return st;
}
/**** Read single line data file ****/
static void readOneLineFile(String filename, double[] data) throws IOException {
FileReader datafile=null;
StreamTokenizer st;
int i=0;
datafile = new FileReader(filename);
st = initST(datafile);
do {
do {
st.nextToken();
} while(st.ttype==st.TT_EOL); /* Skip empty lines */
while(st.ttype == st.TT_NUMBER) {
data[i++] = st.nval;
if(st.nextToken() != ',') break;
st.nextToken();
}
} while( st.ttype == st.TT_EOL );
datafile.close();
}
/**** Read multi line data file ****/
static void readMultiLineFile(String filename, double[][] data) throws IOException {
FileReader datafile=null;
StreamTokenizer st;
int i=0, j=0;
datafile = new FileReader(filename);
st = initST(datafile);
do {
do {
st.nextToken();
} while(st.ttype==st.TT_EOL); /* Skip empty lines */
while(st.ttype == st.TT_NUMBER) {
data[j][i++] = st.nval;
if(st.nextToken() != ',') { j++; i=0; break;}
st.nextToken();
}
} while( st.ttype == st.TT_EOL );
datafile.close();
}
/**** Read data from files ****/
static void readData() throws IOException {
int p,f,r;
REV = new double[NP];
CMAK = new double[NP][NF];
CBUY = new double[NR];
REQ = new double[NP][NR];
MXSELL = new double[NP];
MXMAKE = new double[NF];
PROFIT = new double[NP][NF];
/* Read the revenu data file */
readOneLineFile(REVFILE,REV);
/* Read the production cost data file */
readMultiLineFile(CMAKEFILE,CMAK);
/* Read the raw material cost data file */
readOneLineFile(CBUYFILE,CBUY);
/* Read the resource requirement data file */
readMultiLineFile(REQFILE,REQ);
/* Read the max. sales quantities data file */
readOneLineFile(MXSELLFILE,MXSELL);
/* Read the production capacities data file */
readOneLineFile(MXMAKEFILE,MXMAKE);
/* Calculate the table PROFIT */
for(p=0;p<NP;p++)
for(f=0;f<NF;f++) {
PROFIT[p][f] = REV[p] - CMAK[p][f];
for(r=0;r<NR;r++) PROFIT[p][f] -= (REQ[p][r]*CBUY[r]);
}
}
/***********************************************************************/
public static void main(String[] args) {
try {
readData(); /* Data input from file */
}
catch(IOException e) {
System.err.println(e.getMessage());
System.exit(1);
}
modCoco2(); /* Model and solve the problem */
}
}
|
/********************************************************
Xpress-BCL Java Example Problems
================================
file xbcoco3.java
`````````````````
Coco Problem Phase 3.
Introduce time periods and inventory.
(c) 2008 Fair Isaac Corporation
author: S.Heipcke, Jan. 2000, rev. Mar. 2011
********************************************************/
import java.io.*;
import com.dashoptimization.*;
public class xbcoco3 {
static final int NP = 2; /* Number of products (p) */
static final int NF = 2; /* factories (f) */
static final int NR = 2; /* raw materials (r) */
static final int NT = 4; /* time periods (t) */
static final String REVFILE = System.getProperty("XPRBDATA") +
"/coco/revt.dat";
static final String CMAKEFILE = System.getProperty("XPRBDATA") +
"/coco/cmake.dat";
static final String CBUYFILE = System.getProperty("XPRBDATA") +
"/coco/cbuyt.dat";
static final String REQFILE = System.getProperty("XPRBDATA") +
"/coco/req.dat";
static final String MXSELLFILE = System.getProperty("XPRBDATA") +
"/coco/maxsellt.dat";
static final String MXMAKEFILE = System.getProperty("XPRBDATA") +
"/coco/mxmake.dat";
static final String PSTOCK0FILE = System.getProperty("XPRBDATA") +
"/coco/pstock0.dat";
static final String RSTOCK0FILE = System.getProperty("XPRBDATA") +
"/coco/rstock0.dat";
/****TABLES****/
static double[][] REV; /* Unit selling price of product p in period t */
static double[][] CMAK; /* Unit cost to make product p at factory f */
static double[][] CBUY; /* Unit cost to buy raw material r in period t */
static double[][] REQ; /* Requirement by unit of prod. p for raw mat. r */
static double[][] MXSELL; /* Max. amount of p that can be sold in period t */
static double[] MXMAKE; /* Max. amount factory f can make over all prod.s */
static double[][] PSTOCK0; /* Initial product p stock level at factory f */
static double[][] RSTOCK0; /* Initial raw material r stock level at factory f*/
/****DATA****/
static final double CPSTOCK = 2.0; /* Unit cost to store any product p */
static final double CRSTOCK = 1.0; /* Unit cost to store any raw mat. r */
static final double MXRSTOCK = 300; /* Max. amount of r that can be stored
each f and t */
/***********************************************************************/
static void modCoco3() {
XPRBvar[][][] make, sell, pstock, buy, rstock;
XPRBexpr lobj, lc;
int p,f,r,t;
try (XPRBprob pb = new XPRBprob("Coco3")) { /* Initialize BCL and create a new problem */
/****VARIABLES****/
make = new XPRBvar[NP][NF][NT];
sell = new XPRBvar[NP][NF][NT];
pstock = new XPRBvar[NP][NF][NT+1];
buy = new XPRBvar[NR][NF][NT];
rstock = new XPRBvar[NR][NF][NT+1];
for(p=0;p<NP;p++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++) {
make[p][f][t] = pb.newVar("make_p" + (p+1) + "_f" + (f+1));
/* Amount of prod. p to make at factory f in period t */
sell[p][f][t] = pb.newVar("sell_p" + (p+1) + "_f" + (f+1));
/* Amount of prod. p sold from factory f in period t */
}
for(t=0;t<NT+1;t++)
pstock[p][f][t] = pb.newVar("pstock_p" + (p+1) + "_f" + (f+1));
/* Stock level of prod. p at factory f at start of period t */
}
for(r=0;r<NR;r++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++)
buy[r][f][t] = pb.newVar("buy_r" + (r+1) + "_f" + (f+1));
/* Amount of raw material r bought for factory f in period t */
for(t=0;t<NT+1;t++)
rstock[r][f][t] = pb.newVar("rstock_r" + (r+1) + "_f" + (f+1));
/* Stock level of raw mat. r at factory f at start of per. t */
}
/****OBJECTIVE****/
lobj = new XPRBexpr();
for(f=0;f<NF;f++) { /* Objective: maximize total profit */
for(p=0;p<NP;p++) {
for(t=0;t<NT;t++)
lobj.add(sell[p][f][t].mul(REV[p][t]))
.add(make[p][f][t].mul(-CMAK[p][f]));
for(t=1;t<NT+1;t++) lobj.add(pstock[p][f][t].mul(-CPSTOCK));
}
for(r=0;r<NR;r++) {
for(t=0;t<NT;t++) lobj.add(buy[r][f][t].mul(-CBUY[r][t]));
for(t=1;t<NT+1;t++) lobj.add(rstock[r][f][t].mul(-CRSTOCK));
}
}
pb.setObj(lobj); /* Set objective function */
/****CONSTRAINTS****/
for(p=0;p<NP;p++) /* Product stock balance */
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
pb.newCtr("PBal", pstock[p][f][t].add(make[p][f][t])
.eql(sell[p][f][t].add(pstock[p][f][t+1])) );
for(r=0;r<NR;r++) /* Raw material stock balance */
for(f=0;f<NF;f++)
for(t=0;t<NT;t++) {
lc = new XPRBexpr();
for(p=0;p<NP;p++) lc.add(make[p][f][t].mul(REQ[p][r]));
pb.newCtr("RBal", rstock[r][f][t].add(buy[r][f][t])
.eql(lc.add(rstock[r][f][t+1])) );
}
for(p=0;p<NP;p++)
for(t=0;t<NT;t++)
{ /* Limit on the amount of product p to be sold */
lc = new XPRBexpr();
for(f=0;f<NF;f++) lc.add(sell[p][f][t]);
pb.newCtr("MxSell", lc.lEql(MXSELL[p][t]) );
}
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
{ /* Capacity limit at factory f */
lc = new XPRBexpr();
for(p=0;p<NP;p++) lc.add(make[p][f][t]);
pb.newCtr("MxMake", lc.lEql(MXMAKE[f]) );
}
for(f=0;f<NF;f++)
for(t=1;t<NT+1;t++)
{ /* Raw material stock limit */
lc = new XPRBexpr();
for(r=0;r<NR;r++) lc.add(rstock[r][f][t]);
pb.newCtr("MxRStock", lc.lEql(MXRSTOCK) );
}
/****BOUNDS****/
for(p=0;p<NP;p++)
for(f=0;f<NF;f++)
pstock[p][f][0].fix(PSTOCK0[p][f]); /* Initial product levels */
for(r=0;r<NR;r++)
for(f=0;f<NF;f++)
rstock[r][f][0].fix(RSTOCK0[r][f]); /* Initial raw mat. levels */
/****SOLVING + OUTPUT****/
pb.setSense(XPRB.MAXIM); /* Choose the sense of the optimization */
pb.lpOptimize(""); /* Solve the LP-problem */
System.out.println("Objective: " + pb.getObjVal()); /* Get objective value */
/* Uncomment to print out the solution values */
/* for(p=0;p<NP;p++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++)
System.out.print(make[p][f][t].getName() + ":" + make[p][f][t].getSol()
+ " " + sell[p][f][t].getName() +":" + sell[p][f][t].getSol());
for(t=0;t<NT+1;t++)
System.out.print(pstock[p][f][t].getName() +":" + pstock[p][f][t].getSol() + " ");
System.out.println();
}
for(r=0;r<NR;r++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++)
System.out.print(buy[r][f][t].getName() +":" + buy[r][f][t].getSol()
+ " ");
for(t=0;t<NT+1;t++)
System.out.print(rstock[r][f][t].getName() +":" + rstock[r][f][t].getSol()
+ " ");
}
System.out.println();
*/
}
}
/***********************************************************************/
/**** Initialize the stream tokenizer ****/
static StreamTokenizer initST(FileReader file) {
StreamTokenizer st=null;
st= new StreamTokenizer(file); /* Initialize the stream tokenizer */
st.commentChar('!'); /* Use the character '!' for comments */
st.eolIsSignificant(true); /* Return end-of-line character */
st.ordinaryChar(','); /* Use ',' as separator */
st.parseNumbers(); /* Read numbers as numbers (not strings)*/
return st;
}
/**** Read single line data file ****/
static void readOneLineFile(String filename, double[] data) throws IOException {
FileReader datafile=null;
StreamTokenizer st;
int i=0;
datafile = new FileReader(filename);
st = initST(datafile);
do {
do {
st.nextToken();
} while(st.ttype==st.TT_EOL); /* Skip empty lines */
while(st.ttype == st.TT_NUMBER) {
data[i++] = st.nval;
if(st.nextToken() != ',') break;
st.nextToken();
}
} while( st.ttype == st.TT_EOL );
datafile.close();
}
/**** Read multi line data file ****/
static void readMultiLineFile(String filename, double[][] data) throws IOException {
FileReader datafile=null;
StreamTokenizer st;
int i=0, j=0;
datafile = new FileReader(filename);
st = initST(datafile);
do {
do {
st.nextToken();
} while(st.ttype==st.TT_EOL); /* Skip empty lines */
while(st.ttype == st.TT_NUMBER) {
data[j][i++] = st.nval;
if(st.nextToken() != ',') { j++; i=0; break;}
st.nextToken();
}
} while( st.ttype == st.TT_EOL );
datafile.close();
}
/**** Read data from files ****/
static void readData() throws IOException {
int p,r;
REV = new double[NP][NT];
CMAK = new double[NP][NF];
CBUY = new double[NR][NT];
REQ = new double[NP][NR];
MXSELL = new double[NP][NT];
MXMAKE = new double[NF];
PSTOCK0 = new double[NP][NF];
RSTOCK0 = new double[NR][NF];
/* Read the revenu data file */
readMultiLineFile(REVFILE,REV);
/* Read the production cost data file */
readMultiLineFile(CMAKEFILE,CMAK);
/* Read the raw material cost data file */
readMultiLineFile(CBUYFILE,CBUY);
/* Read the resource requirement data file */
readMultiLineFile(REQFILE,REQ);
/* Read the max. sales quantities data file */
readMultiLineFile(MXSELLFILE,MXSELL);
/* Read the production capacities data file */
readOneLineFile(MXMAKEFILE,MXMAKE);
/* Read the product stock data file */
readMultiLineFile(PSTOCK0FILE,PSTOCK0);
/* Read the raw material stock data file */
readMultiLineFile(RSTOCK0FILE,RSTOCK0);
}
/***********************************************************************/
public static void main(String[] args) {
try {
readData(); /* Data input from file */
}
catch(IOException e) {
System.err.println(e.getMessage());
System.exit(1);
}
modCoco3(); /* Model and solve the problem */
}
}
|
/********************************************************
Xpress-BCL Java Example Problems
================================
file xbcoco.java
````````````````
Complete Coco Problem.
Specify phase by PHASE parameter.
Data input in the model, not via data files.
(c) 2008 Fair Isaac Corporation
author: S.Heipcke, Jan. 2000, rev. Mar. 2011
********************************************************/
import com.dashoptimization.*;
public class xbcoco {
static final int PHASE = 5;
/* Phase = 3: Multi-period parameterised model; mines always open
* Phase = 4: Mines may open/closed freely; when closed save 20000 per month
* Phase = 5: Once closed always closed; larger saving */
static final int NP = 2; /* Number of products (p) */
static final int NF = 2; /* factories (f) */
static final int NR = 2; /* raw materials (r) */
static final int NT = 4; /* time periods (t) */
/****DATA****/
static final double[][] REV =
{{400, 380, 405, 350},
{410, 397, 412, 397}};
/* Unit selling price of prod. p in period t */
static final double[][] CMAK =
{{150, 153},
{ 75, 68}}; /* Unit cost to make product p at factory f */
static final double[][] CBUY =
{{100, 98, 97, 100},
{200, 195, 198, 200}};
/* Unit cost to buy raw material r in period t */
static final double[] COPEN = {50000, 63000};
/* Fixed cost of factory f being open for one period */
static final double CPSTOCK = 2.0; /* Unit cost to store any product p */
static final double CRSTOCK = 1.0; /* Unit cost to store any raw material r */
static final double[][] REQ =
{{1.0, 0.5},
{1.3, 0.4}}; /* Requirement by unit of prod. p for raw material r */
static final double[][] MXSELL =
{{650, 600, 500, 400},
{600, 500, 300, 250}};
/* Max. amount of p that can be sold in period t */
static final double[] MXMAKE = {400, 500};
/* Max. amount factory f can make over all products */
static final double MXRSTOCK = 300;
/* Max. amount of r that can be stored each f and t */
static final double[][] PSTOCK0 =
{{50, 100},
{50, 50}}; /* Initial product p stock level at factory f */
static final double[][] RSTOCK0 =
{{100, 150},
{ 50, 100}}; /* Initial raw material r stock level at factory f*/
/***********************************************************************/
public static void main(String[] args) {
XPRB bcl;
XPRBvar[][][] make, sell, pstock, buy, rstock;
XPRBvar[][] openm;
XPRBexpr lobj, lc;
int p,f,r,t;
try (XPRBprob pb = new XPRBprob("Coco")) { /* Initialize BCL and create a new problem */
/****VARIABLES****/
make = new XPRBvar[NP][NF][NT];
sell = new XPRBvar[NP][NF][NT];
pstock = new XPRBvar[NP][NF][NT+1];
buy = new XPRBvar[NR][NF][NT];
rstock = new XPRBvar[NR][NF][NT+1];
openm = new XPRBvar[NF][NT];
for(p=0;p<NP;p++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++) {
make[p][f][t] = pb.newVar("make_p" + (p+1) + "_f" + (f+1));
/* Amount of prod. p to make at factory f in period t */
sell[p][f][t] = pb.newVar("sell_p" + (p+1) + "_f" + (f+1));
/* Amount of prod. p sold from factory f in period t */
}
for(t=0;t<NT+1;t++)
pstock[p][f][t] = pb.newVar("pstock_p" + (p+1) + "_f" + (f+1));
/* Stock level of prod. p at factory f at start of period t */
}
for(r=0;r<NR;r++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++)
buy[r][f][t] = pb.newVar("buy_r" + (r+1) + "_f" + (f+1));
/* Amount of raw material r bought for factory f in period t */
for(t=0;t<NT+1;t++)
rstock[r][f][t] = pb.newVar("rstock_r" + (r+1) + "_f" + (f+1));
/* Stock level of raw mat. r at factory f at start of per. t */
}
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
openm[f][t] = pb.newVar("open_f" + (f+1), XPRB.BV);
/* 1 if factory f is open in period t, else 0 */
/****OBJECTIVE****/
lobj = new XPRBexpr();
for(f=0;f<NF;f++) { /* Objective: maximize total profit */
for(p=0;p<NP;p++) {
for(t=0;t<NT;t++)
lobj.add(sell[p][f][t].mul(REV[p][t]))
.add(make[p][f][t].mul(-CMAK[p][f]));
for(t=1;t<NT+1;t++) lobj.add(pstock[p][f][t].mul(-CPSTOCK));
}
if(PHASE==4)
for(t=0;t<NT;t++) lobj.add(openm[f][t].mul(20000-COPEN[f]));
else if(PHASE==5)
for(t=0;t<NT;t++) lobj.add(openm[f][t].mul(-COPEN[f]));
for(r=0;r<NR;r++) {
for(t=0;t<NT;t++) lobj.add(buy[r][f][t].mul(-CBUY[r][t]));
for(t=1;t<NT+1;t++) lobj.add(rstock[r][f][t].mul(-CRSTOCK));
}
}
pb.setObj(lobj); /* Set objective function */
/****CONSTRAINTS****/
for(p=0;p<NP;p++) /* Product stock balance */
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
pb.newCtr("PBal", pstock[p][f][t].add(make[p][f][t])
.eql(sell[p][f][t].add(pstock[p][f][t+1])) );
for(r=0;r<NR;r++) /* Raw material stock balance */
for(f=0;f<NF;f++)
for(t=0;t<NT;t++) {
lc = new XPRBexpr();
for(p=0;p<NP;p++) lc.add(make[p][f][t].mul(REQ[p][r]));
pb.newCtr("RBal", rstock[r][f][t].add(buy[r][f][t])
.eql(lc.add(rstock[r][f][t+1])) );
}
for(p=0;p<NP;p++)
for(t=0;t<NT;t++)
{ /* Limit on the amount of product p to be sold */
lc = new XPRBexpr();
for(f=0;f<NF;f++) lc.add(sell[p][f][t]);
pb.newCtr("MxSell", lc.lEql(MXSELL[p][t]) );
}
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
{ /* Capacity limit at factory f */
lc = new XPRBexpr();
for(p=0;p<NP;p++) lc.add(make[p][f][t]);
pb.newCtr("MxMake", lc.lEql(openm[f][t].mul(MXMAKE[f])) );
}
for(f=0;f<NF;f++)
for(t=1;t<NT+1;t++)
{ /* Raw material stock limit */
lc = new XPRBexpr();
for(r=0;r<NR;r++) lc.add(rstock[r][f][t]);
pb.newCtr("MxRStock", lc.lEql(MXRSTOCK) );
}
if(PHASE==5)
for(f=0;f<NF;f++)
for(t=0;t<NT-1;t++) /* Once closed, always closed */
pb.newCtr("Closed", openm[f][t+1].lEql(openm[f][t]) );
/****BOUNDS****/
for(p=0;p<NP;p++)
for(f=0;f<NF;f++)
pstock[p][f][0].fix(PSTOCK0[p][f]); /* Initial product levels */
for(r=0;r<NR;r++)
for(f=0;f<NF;f++)
rstock[r][f][0].fix(RSTOCK0[r][f]); /* Initial raw mat. levels */
if(PHASE<=3)
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
openm[f][t].fix(1);
/****SOLVING + OUTPUT****/
pb.setSense(XPRB.MAXIM); /* Choose the sense of the optimization */
pb.mipOptimize(""); /* Solve the MIP-problem */
System.out.println("Objective: " + pb.getObjVal()); /* Get objective value */
/* Uncomment to print out the solution values */
/* for(p=0;p<NP;p++)
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
System.out.print(make[p][f][t].getName() + ":" + make[p][f][t].getSol()
+ " " + sell[p][f][t].getName() +":" + sell[p][f][t].getSol());
System.out.println();
for(p=0;p<NP;p++)
for(f=0;f<NF;f++)
for(t=0;t<NT+1;t++)
System.out.print(pstock[p][f][t].getName() +":" + pstock[p][f][t].getSol() + " ");
System.out.println();
for(r=0;r<NR;r++)
for(f=0;f<NF;f++) {
for(t=0;t<NT;t++)
System.out.print(buy[r][f][t].getName() +":" + buy[r][f][t].getSol()
+ " ");
for(t=0;t<NT+1;t++)
System.out.print(rstock[r][f][t].getName() +":" + rstock[r][f][t].getSol()
+ " ");
}
for(f=0;f<NF;f++)
for(t=0;t<NT;t++)
System.out.print(openm[f][t].getName() +":" + openm[f][t].getSol() + " ");
System.out.println();
*/
}
}
}
|
