/******************************************************** Xpress-BCL C++ Example Problems =============================== file xbdlvriis2iso.cxx `````````````````````` Transportation problem (infeasible data). Retrieving and printing IIS. - Using Optimizer functions to retrieve detailed IIS information including isolation rows/bounds - (c) 2008 Fair Isaac Corporation author: S.Heipcke, Jan. 2008, rev. Mar. 2011 ********************************************************/ #include <iostream>
#include <cstdio>
#include <cstring>
#include "xprb_cpp.h"
#include "xprs.h"

using namespace std;
using namespace ::dashoptimization;

#define NSupp 10                  // Number of suppliers
#define NCust 7                   // Number of customers
#define MaxArcs 100               // Max. num. of non-zero cost values

#define VANFILE XPRBDATAPATH "/delivery/ifvan.dat"    // Van data file
#define COSTFILE XPRBDATAPATH "/delivery/cost.dat"    // Cost data file

/****DATA****/
// Supplier:      London  Luton  B'ham Bristl  Derby Stckpt   York
double SUPPLY[] = {140.0, 200.0,  50.0,  10.0, 400.0, 200.0,  20.0,
// Supplier:       Derby  Soton Scnthp
                   90.0,  30.0,  12};
// Customer:       London Livpol Doncst   York   Hull  Manchr Shffld
double DEMAND[] = {1230.3, 560.4, 117.1, 592.8, 310.0, 1247.0,  86};

double COST[NSupp][NCust];        // Cost per supplier-customer pair

double IFVAN[NSupp][NCust];       // Non-zero if route uses vans instead 
                                  // of lorries
double VANCAP=40;               // Capacity on routes that use vans

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

void modDelivery(XPRBprob &p)
{
 XPRBexpr lobj, lc;
 int s,c,i;
 XPRBvar x[NSupp][NCust];

 int numv, numc, numiis, len, ncol, nrow;
 char *vnames, *cnames;
 int *viis,*ciis;
 char **vindex,**cindex;
 double bnd, rhs;
 char *ctrtype = NULL;
 char *bndtype = NULL;
 double *duals = NULL;
 double *rdcs = NULL;
 char *isolationrows = NULL;
 char *isolationbnds = NULL;
 const char *isotype[] = {"N/A", "No ", "Yes"};

 XPRSprob op;
  
/****VARIABLES****/
 for(s=0;s<NSupp;s++) for(c=0; x[s][c]=p.newVar(XPRBnewname("x_s%d",s)); for(s=0;s<NSupp;s++) Objective: Minimize total cost lobj +=COST[s][c]*x[s][c]; Set objective function Satisfy demand of each customer lc=0;>= DEMAND[c]);
 }
 for(c=4; c<NCust; Satisfy demand of each customer lc=0; for(s=5;s<NSupp;s++) +=x[s][c];>= DEMAND[c]);
 }
        
 for(s=0;s<5;s++)                 // Keep within supply at each supplier
 {
  lc=0;
  for(c=0; c<4; c++)  lc+= x[s][c];
  p.newCtr("Supply", lc <= SUPPLY[s]);
 }          
 for(s=5;s<NSupp;s++) Keep within supply at each supplier lc=0; for(c=4; lc+=x[s][c]; <=SUPPLY[s]); for(s=0;s<NSupp;s++) if(IFVAN[s][c]!=0) Set objective sense to minimization Solve the LP-problem cout <<="LP status: " if (p.getLPStat()==XPRB_LP_OPTIMAL) Get value else Problem is infeasible op=p.getXPRSprob(); Retrieve Optimizer problem all IIS numiis=p.getNumIIS(); Generate get number of independent Alternatively: use functions generate Obtain variable and constraint names for later in printout bytes required retrieving vnames=new vindex=new vindex[0]=vnames; for(i=1; vindex[i]=vindex[i-1]+strlen(vindex[i-1])+1; cnames=new cindex=new cindex[0]=cnames; cindex[i]=cindex[i-1]+strlen(cindex[i-1])+1; detailed info Find isolations ciis=new viis=new ctrtype=new bndtype=new duals=new rdcs=new isolationrows=new isolationbnds=new s numv>0)
   {                              // Print all variables in the IIS
    for(i=0;i<numv;i++) if (bndtype[i]== else cout <<=" " bnd delete Free the array of variables>0)
   {                              // Print all constraints in the IIS
    for(i=0;i<numc;i++) 
    {
     XPRSgetrhs(op, &rhs, ciis[i], ciis[i]);
     cout << " " << cindex[ciis[i]]<< "    row     " << ctrtype[i];
     cout << "   " << rhs << "   " << duals[i] << "       ";
     cout << isotype[1+(int)isolationrows[i]] << endl;
    } 
    delete [] ciis;               // Free the array of constraints
   } 
  }

  delete [] vnames;
  delete [] cnames;
  delete [] vindex;
  delete [] cindex;
 }   
    
}

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

    /**** Read data from files ****/
void readData()
{
 FILE *datafile;
 int s,c;  

        // Initialize data tables to 0
 for(s=0;s<NSupp;s++)
  for(c=0; c<NCust; c++)
  {
   COST[s][c] = 0;
   IFVAN[s][c] = 0;
  }  
        // Read the demand data file
 datafile=fopen(COSTFILE,"r");
 for(s=0;s<NSupp;s++)
  XPRBreadarrlinecb(XPRB_FGETS, datafile, 99, "g,", COST[s], NCust);
 fclose(datafile);

        // Read the van data file
 datafile=fopen(VANFILE,"r");
 for(s=0;s<NSupp;s++)
  XPRBreadarrlinecb(XPRB_FGETS, datafile, 99, "g,", IFVAN[s], NCust);
 fclose(datafile);
}

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

int main(int argc, char **argv)
{
 XPRBprob p("Delivery");          // Initialize a new problem in BCL
 readData();                      // Data input from file
 modDelivery(p);                  // Problem formulation and solution

 return 0;
}