/********************************************************
  BCL Example Problems
  ====================

  file xbworkrng.c
  ````````````````
  Workshop planning example.
  Test ranges and number printing format.

  (c) 2008-2024 Fair Isaac Corporation
      author: S.Heipcke, 2003, rev. Mar. 2011
********************************************************/

#include <stdio.h>
#include "xprb.h"

#define NProd 2                  /* Number of products */
#define NShop 3                  /* Number of workshops */
#define WMAX 40                  /* Maximum weekly working time */

/****DATA****/
int DUR[][NShop] =  {{5, 9, 7},  /* Duration of product p on shop s */
                    {10, 2, 5}};
int RES[]   = {10, 8};           /* Man hours per unit */
int PRICE[] = {108, 84};         /* Selling price per unit */

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

int main(int argc, char **argv)
{
 int p,s;
 XPRBctr c[NShop],cobj;
 XPRBvar x[NProd];                   /* Amount of product p */
 XPRBprob prob;

 prob=XPRBnewprob("Workshop");       /* Initialize a new problem in BCL */

/****VARIABLES****/
 for(p=0;p<NProd;p++)
  x[p] = XPRBnewvar(prob,XPRB_PL,"x",0,XPRB_INFINITY);

/****OBJECTIVE****/
 cobj = XPRBnewctr(prob,"OBJ",XPRB_N);  /* Objective: Maximize Benefit */
 for(p=0;p<NProd;p++)
  XPRBaddterm(cobj,x[p],PRICE[p]-5*RES[p]);
 XPRBsetobj(prob,cobj);              /* Select objective function */

/****CONSTRAINTS****/
 for(s=0;s<NShop;s++)                /* Limit on weekly working hours */
 {
  c[s]=XPRBnewctr(prob,"ResMax",XPRB_L);
  for(p=0;p<NProd;p++)
   XPRBaddterm(c[s],x[p],DUR[p][s]);
  XPRBaddterm(c[s],NULL,WMAX);
 }

/****SOLVING + OUTPUT****/
 XPRBsetrealfmt(prob, "%4.2e");
 for(p=0;p<NProd;p++)
 { XPRBprintvar(x[p]); printf(" "); }
 printf("\n");

 XPRBsetsense(prob,XPRB_MAXIM);      /* Set obj. sense to maximization */
 XPRBlpoptimize(prob,"");            /* Solve the LP-problem */
 printf("Objective: %g\n",XPRBgetobjval(prob));  /* Get objective value */

 XPRBsetrealfmt(prob, "%g, ");
 for(p=0;p<NProd;p++)                /* Print the solution values */
  XPRBprintvar(x[p]);
 printf("\n");

 XPRBsetrealfmt(prob, "%8.4f");
 XPRBprintprob(prob);
                                 /* Row ranges */
 printf("Ctr: Lower activity, Upper activity, Unit cost DN, Unit cost UP\n");
 for(s=0;s<NShop;s++)
  printf("%s: %g, %g, %g, %g\n", XPRBgetctrname(c[s]),
          XPRBgetctrrng(c[s],XPRB_LOACT), XPRBgetctrrng(c[s],XPRB_UPACT),
          XPRBgetctrrng(c[s],XPRB_UDN), XPRBgetctrrng(c[s],XPRB_UUP));

                                 /* Column ranges */
 printf("Var: Lower activity, Upper activity, Unit cost DN, Unit cost UP, Lower profit, Upper profit\n");
 for(p=0;p<NProd;p++)
  printf("%s: %g, %g, %g, %g, %g, %g\n", XPRBgetvarname(x[p]),
          XPRBgetvarrng(x[p],XPRB_LOACT), XPRBgetvarrng(x[p],XPRB_UPACT),
          XPRBgetvarrng(x[p],XPRB_UDN), XPRBgetvarrng(x[p],XPRB_UUP),
          XPRBgetvarrng(x[p],XPRB_LCOST), XPRBgetvarrng(x[p],XPRB_UCOST));

 return 0;
}