/******************************************************** 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 XPRB bcl; static XPRBprob pb; /***********************************************************************/ static void modCoco2() { XPRBvar[][] make; XPRBexpr lobj, lc; int p,f; bcl = new XPRB(); /* Initialize BCL */ pb = bcl.newProb("Coco2"); /* Create a new problem in BCL */ /****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 */
 }
}