/********************************************************
  Xpress-BCL C++ Example Problems
  ===============================

  file xbcoco2.cxx
  `````````````````
  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-2024 Fair Isaac Corporation
      author: S.Heipcke, Jan. 2000, rev. Mar. 2011
********************************************************/

#include <iostream>
#include <cstdio>
#include "xprb_cpp.h"

using namespace std;
using namespace ::dashoptimization;

#define NP 2            /* Number of products (p) */
#define NF 2            /*           factories (f) */
#define NR 2            /*           raw materials (r) */

#define REVFILE     XPRBDATAPATH "/coco/rev.dat"
#define CMAKEFILE   XPRBDATAPATH "/coco/cmake.dat"
#define CBUYFILE    XPRBDATAPATH "/coco/cbuy.dat"
#define REQFILE     XPRBDATAPATH "/coco/req.dat"
#define MXSELLFILE  XPRBDATAPATH "/coco/maxsell.dat"
#define MXMAKEFILE  XPRBDATAPATH "/coco/mxmake.dat"

/****TABLES****/
double REV[NP];         /* Unit selling price of product p */
double CMAK[NP][NF];    /* Unit cost to make product p at factory f */
double CBUY[NR];        /* Unit cost to buy raw material r */
double REQ[NP][NR];     /* Requirement by unit of prod. p for raw material r */
double MXSELL[NP];      /* Max. amount of p that can be sold */
double MXMAKE[NF];      /* Max. amount factory f can make over all products */
double PROFIT[NP][NF];  /* Profit contribution of product p at factory f */
 
/***********************************************************************/

void modCoco2(XPRBprob &pb)
{
 XPRBvar make[NP][NF];
 XPRBexpr lobj, lc;
 int p,f;
 
/****VARIABLES****/
 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(XPRBnewname("make_p%df%d",p+1,f+1));    

/****OBJECTIVE****/
 for(p=0;p<NP;p++)         /* Objective: maximize total profit */
  for(f=0;f<NF;f++)  lobj += PROFIT[p][f]* make[p][f]; 
 pb.setObj(pb.newCtr("OBJ",lobj)); 

/****CONSTRAINTS****/
   
 for(p=0;p<NP;p++)         /* Limit on the amount of product p to be sold */
 {
  lc=0;
  for(f=0;f<NF;f++)  lc += make[p][f];
  pb.newCtr("MxSell", lc <= MXSELL[p]);
 }   

 for(f=0;f<NF;f++)         /* Capacity limit at factory f */
 {
  lc=0;
  for(p=0;p<NP;p++) lc += make[p][f];
  pb.newCtr("MxMake", lc <= MXMAKE[f]);      
 }

/****SOLVING + OUTPUT****/
 pb.setSense(XPRB_MAXIM);  /* Choose the sense of the optimization */
 pb.lpOptimize("");        /* Solve the LP-problem */
 cout << "Objective: " << pb.getObjVal() << endl;  /* Get objective value */
    
 for(p=0;p<NP;p++)         /* Print the solution values */      
  for(f=0;f<NF;f++)   
   cout << make[p][f].getName() << ":" << make[p][f].getSol() << " ";
 cout << endl;
}

/***********************************************************************/

    /**** Read data from files ****/
void readData()
{
 FILE *datafile;
 int p,f,r;
        /* Read the revenu data file */
 datafile=fopen(REVFILE,"r");
 XPRBreadarrlinecb(XPRB_FGETS, datafile, 200, "g,", REV, NP);
 fclose(datafile);

        /* Read the production cost data file */
 datafile=fopen(CMAKEFILE,"r");
 for(p=0;p<NP;p++)
  XPRBreadarrlinecb(XPRB_FGETS, datafile, 200, "g,", CMAK[p], NF);
 fclose(datafile);

        /* Read the raw material cost data file */
 datafile=fopen(CBUYFILE,"r");
 XPRBreadarrlinecb(XPRB_FGETS, datafile, 200, "g,", CBUY, NR);
 fclose(datafile);

        /* Read the resource requirement data file */
 datafile=fopen(REQFILE,"r");
 for(p=0;p<NP;p++)
  XPRBreadarrlinecb(XPRB_FGETS, datafile, 200, "g,", REQ[p], NR);
 fclose(datafile);
  
        /* Read the max. sales quantities data file */
 datafile=fopen(MXSELLFILE,"r");
 XPRBreadarrlinecb(XPRB_FGETS, datafile, 200, "g,", MXSELL, NP);
 fclose(datafile);
  
        /* Read the production capacities data file */
 datafile=fopen(MXMAKEFILE,"r");
 XPRBreadarrlinecb(XPRB_FGETS, datafile, 200, "g,", MXMAKE, NF);
 fclose(datafile);
 
        /* 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]);
  }
}

/***********************************************************************/

int main(int argc, char **argv)
{
 XPRBprob pb("Coco2");  /* Initialize a new problem in BCL */
 readData();            /* Data input from file */
 modCoco2(pb);          /* Model and solve the problem */
  
 return 0;
}