/*********************************************************************** Xpress Optimizer Examples ========================= file mipsolenum.c ````````````````` Find the 10 best solutions with the MIP solution enumerator We take the power generation problem stored in hpw15.mat which seeks to optimise the operating pattern of a group of electricity generators. We run the MIP solution enumerator on the problem using the default setup obtaining the best 10 solutions. The best 10 solutions are stored to a MIP solution pool. The solutions' objectives and solution values are printed to screen. (c) 2017 Fair Isaac Corporation ***********************************************************************/ #include <stdio.h>
#include <stdlib.h>
#include "xprs.h"                     /* Optimizer header file */
#include "xprs_mse_defaulthandler.h"  /* 'Default' enumerated mip solution handler */

void errormsg(const char *sSubName,int nLineNo,int nErrCode);


void main(void) {
  int nReturn;
  XPRSprob prob;
  XPRSmipsolpool msp;
  XPRSmipsolenum mse;
  int nMaxSols, iPresolveOps, iMIPPresolve;
  int i, j, nSols, nCols, iSolutionId, iSolutionIdStatus;
  double dObj, dSol;
  const char *sProblem = "hpw15";

  nReturn = XPRSinit("");
  if (nReturn!=0 && nReturn!=16) {
    char errmsgbuf[2048];
    XPRSgetlicerrmsg(errmsgbuf,2048);
    printf("We were unable to license the Xpress-optimizer: %s\n", errmsgbuf);
    exit(nReturn);
  }


  nReturn = XPRS_mse_create(&mse);
  if (nReturn != 0) errormsg("XPRS_mse_create",__LINE__,nReturn);

  nReturn = XPRScreateprob(&prob);
  if (nReturn != 0) errormsg("XPRScreateprob",__LINE__,nReturn);

  nReturn = XPRS_msp_create(&msp);
  if (nReturn != 0) errormsg("XPRS_msp_create",__LINE__,nReturn);

  nReturn = XPRSreadprob(prob, sProblem, "");
  if (nReturn != 0) errormsg("XPRSreadprob",__LINE__,nReturn);

  nReturn = XPRSgetintattrib(prob, XPRS_COLS, &nCols);
  if (nReturn != 0) errormsg("XPRSgetintattrib",__LINE__,nReturn);

  /* Avoid duplicated solutions from heuristics. */
  nReturn = XPRS_msp_setintcontrol(msp, XPRS_MSP_DUPLICATESOLUTIONSPOLICY, 3);
  if (nReturn != 0) errormsg("XPRS_msp_setintcontrol",__LINE__,nReturn);

  /* Disable dual reductions to prevent dominated solutions from being */
  /* presolved away. */
  nReturn = XPRSsetintcontrol(prob, XPRS_MIPDUALREDUCTIONS, 2);
  if (nReturn != 0) errormsg("XPRSsetintcontrol",__LINE__,nReturn);

  /* Run the enumeration */
  nMaxSols = 10;
  nReturn = XPRS_mse_minim(mse, prob, msp, XPRS_mse_defaulthandler, NULL, &nMaxSols);
  if (nReturn != 0) errormsg("XPRS_mse_minim",__LINE__,nReturn);

  /* Print out the solutions found */
  nReturn = XPRS_mse_getintattrib(mse, XPRS_MSE_SOLUTIONS, &nSols);
  if (nReturn != 0) errormsg("XPRS_mse_getintattrib",__LINE__,nReturn);
  for(i = 1; i <= nSols; i++) {
    nReturn = XPRS_mse_getsollist(mse, XPRS_MSE_METRIC_MIPOBJECT, i, i, &iSolutionId, NULL, NULL);
    if (nReturn!=0) errormsg("XPRS_mse_getsollist",__LINE__,nReturn);
    nReturn = XPRS_mse_getsolmetric(mse, iSolutionId, &iSolutionIdStatus, XPRS_MSE_METRIC_MIPOBJECT, &dObj);
    if (nReturn!=0) errormsg("XPRS_mse_getsolmetric",__LINE__,nReturn);
    printf("--------\n%3i = %12,5f\n", i, dObj);

    for(j = 0; j < nCols; j++) {
      nReturn = XPRS_msp_getsol(msp, iSolutionId, &iSolutionIdStatus, &dSol, j, j, NULL);
      if (nReturn!=0) errormsg("XPRS_msp_getsol",__LINE__,nReturn);
      printf("%3i = %12,5f\n", j, dSol);
    }
  }

  XPRS_msp_destroy(msp);
  XPRSdestroyprob(prob);
  XPRS_mse_destroy(mse);
  XPRSfree();
}


/**********************************************************************************\
* Name:         errormsg                                                           *
* Purpose:      Display error information about failed subroutines.                *
* Arguments:    const char *sSubName   Subroutine name                             *
*               int nLineNo            Line number                                 *
*               int nErrCode           Error code                                  *
* Return Value: None                                                               *
\**********************************************************************************/
void errormsg(const char *sSubName,int nLineNo,int nErrCode)
{
   int nErrNo;   /* Optimizer error number */

   /* Print error message */
   printf("The subroutine %s has failed on line %d\n",sSubName,nLineNo);

   /* Append the error code, if it exists */
   if (nErrCode!=-1) printf("with error code %d\n",nErrCode);

   /* Free memory, close files and exit */
   XPRSfree();
   exit(nErrCode);
}