/******************************************************** BCL Example Problems ==================== file xbdlvriis2iso.c ```````````````````` 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 <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "xprb.h"
#include "xprs.h"

#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 prob)
{
 XPRBctr ctr, CSupply[NSupp], CDemand[NCust];
 int s,c,i;
 XPRBvar x[NSupp][NCust];

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

 XPRSprob op;

/****VARIABLES****/
 for(s=0;s<NSupp;s++) for(c=0; x[s][c]=XPRBnewvar(prob,XPRB_PL,XPRBnewname("x_s%d_%d",s,c), ctr=XPRBnewctr(prob,"OBJ",XPRB_N); for(s=0;s<NSupp;s++) Objective: Minimize total cost Set objective function Satisfy demand of each customer CDemand[c]=XPRBnewctr(prob,"Demand", Keep within supply at CSupply[s]=XPRBnewctr(prob,"Supply",XPRB_L); if(IFVAN[s][c]!=0) sense to minimization XPRB_LP,="infeasible" Solve the LP-problem status: if (XPRBgetlpstat(prob)==XPRB_LP_OPTIMAL) Get value else IIS op=XPRBgetXPRSprob(prob); Retrieve Optimizer problem all Generate number independent Obtain variable and constraint names for later use in printout bytes required retrieving vnames=(char vindex=(char vindex[0]=vnames; for(i=1; vindex[i]=vindex[i-1]+strlen(vindex[i-1])+1; namelength=0;>namelength) namelength=strlen(vindex[i]);

                                  /* Retrieve constraint names */
  XPRSgetintattrib(op, XPRS_ORIGINALROWS, &nrow);
  XPRSgetnamelist(op, 1, NULL, 0, &len, 0, nrow-1);
  cnames = (char *)malloc(len * sizeof(char));
  cindex = (char **)malloc(nrow * sizeof(char *));
  XPRSgetnamelist(op, 1, cnames, len, NULL, 0, nrow-1);
  cindex[0]=cnames;
  for(i=1; i<nrow; cindex[i]=cindex[i-1]+strlen(cindex[i-1])+1; for(i=1; if>namelength) namelength=strlen(cindex[i]);


/**** Retrieve detailed IIS info (incl. isolations) ****/
  for(s=1;s<=numiis;s++)
  {
   XPRSgetiisdata(op, s, &numc, &numv, NULL, NULL, NULL, NULL,
                  NULL, NULL, NULL, NULL);

   XPRSiisisolations(op, s);     /* Find isolations */

   ciis = (int *)malloc(numc * sizeof(int));
   viis = (int *)malloc(numv * sizeof(int));
   ctrtype = (char *)malloc(numc*sizeof(char));
   bndtype = (char *)malloc(numv*sizeof(char));
   duals = (double *)malloc(numc*sizeof(double));
   rdcs = (double *)malloc(numv*sizeof(double));
   isolationrows = (char *)malloc(numc*sizeof(char));
   isolationbnds = (char *)malloc(numv*sizeof(char));

   XPRSgetiisdata(op, s, &numc, &numv, ciis, viis, ctrtype, bndtype,
                  duals, rdcs, isolationrows, isolationbnds);
   printf("IIS %d:  %d variables, %d constraints\n", s, numv, numc);
   printf("   %-*s  Type    Sense   Bound   Dual values   In iso. \n", namelength, "Name");
   if (numv>0)
   {                              /* Print all variables in the IIS */
    for(i=0;i<numv;i++) if (bndtype[i]== else vindex[viis[i]],="column" 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]);
     printf(" %-*s    %s  %c %10g  %10g       %s\n", namelength,
            cindex[ciis[i]], "row   ", ctrtype[i], rhs, duals[i],
            isotype[1+(int)isolationrows[i]]);
    }
    free(ciis);                   /* Free the array of constraints */
   }
  }

  free(vnames);
  free(cnames);
  free(vindex);
  free(cindex);
 }

}

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

    /**** Read data from files ****/
void readdata(void)
{
 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 prob;

 prob=XPRBnewprob("Delivery");   /* Initialize a new problem in BCL */
 readdata();                     /* Data input from file */
 moddelivery(prob);              /* Problem formulation & solving */

 return 0;
}