/* -*- C++ -*-

  *

  * This file is a part of LEMON, a generic C++ optimization library

  *

  * Copyright (C) 2003-2007

  * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

  * (Egervary Research Group on Combinatorial Optimization, EGRES).

  *

  * Permission to use, modify and distribute this software is granted

  * provided that this copyright notice appears in all copies. For

  * precise terms see the accompanying LICENSE file.

  *

  * This software is provided "AS IS" with no warranty of any kind,

  * express or implied, and with no claim as to its suitability for any

  * purpose.

  *

  */

 ///\file

 ///\brief Implementation of the LEMON-CPLEX mip solver interface.

 #include <lemon/mip_cplex.h>

 namespace lemon {

   MipCplex::MipCplex() {

     //This is unnecessary: setting integrality constraints on

     //variables will set this, too 

     ///\todo The constant CPXPROB_MIP is

     ///called CPXPROB_MILP in later versions

 #if CPX_VERSION < 800

     CPXchgprobtype( env,  lp, CPXPROB_MIP);

 #else

     CPXchgprobtype( env,  lp, CPXPROB_MILP);

 #endif

   }

   void MipCplex::_colType(int i, MipCplex::ColTypes col_type){

     // Note If a variable is to be changed to binary, a call to CPXchgbds

     // should also be made to change the bounds to 0 and 1.

     int indices[1];

     indices[0]=i;

     char ctype[1];

     switch (col_type){

       case INT:

 	ctype[0]=CPX_INTEGER;//'I'

 	break;

       case REAL:

 	ctype[0]=CPX_CONTINUOUS	;//'C'

 	break;

     default:;

         //FIXME problem

     }

     CPXchgctype (env, lp, 1, indices, ctype);

   }

   MipCplex::ColTypes MipCplex::_colType(int i) const {

     char ctype[1];

     CPXgetctype (env, lp, ctype, i, i);

     switch (ctype[0]){

     case CPX_INTEGER:

       return INT;

     case CPX_CONTINUOUS:

       return REAL;

     default:

       return REAL;//Error!

     }

   }

   LpCplex::SolveExitStatus MipCplex::_solve(){

     status = CPXmipopt (env, lp);

     if (status==0)

       return SOLVED;

     else

       return UNSOLVED;

   }

   LpCplex::SolutionStatus MipCplex::_getMipStatus() const {

     int stat = CPXgetstat(env, lp);

     //Fortunately, MIP statuses did not change for cplex 8.0

     switch (stat)

     {

       case CPXMIP_OPTIMAL:

         // Optimal integer solution has been found.

       case CPXMIP_OPTIMAL_TOL:

         // Optimal soluton with the tolerance defined by epgap or epagap has 

         // been found.

         return OPTIMAL;

 	//This also exists in later issues

 	//    case CPXMIP_UNBOUNDED:

         //return INFINITE;

       case CPXMIP_INFEASIBLE:

         return INFEASIBLE;

       default:

         return UNDEFINED;

     }

     //Unboundedness not treated well: the following is from cplex 9.0 doc

     // About Unboundedness

     // The treatment of models that are unbounded involves a few

     // subtleties. Specifically, a declaration of unboundedness means that

     // ILOG CPLEX has determined that the model has an unbounded

     // ray. Given any feasible solution x with objective z, a multiple of

     // the unbounded ray can be added to x to give a feasible solution

     // with objective z-1 (or z+1 for maximization models). Thus, if a

     // feasible solution exists, then the optimal objective is

     // unbounded. Note that ILOG CPLEX has not necessarily concluded that

     // a feasible solution exists. Users can call the routine CPXsolninfo

     // to determine whether ILOG CPLEX has also concluded that the model

     // has a feasible solution.

   }  

   MipCplex::Value MipCplex::_getPrimal(int i) const {

     Value x;

     CPXgetmipx(env, lp, &x, i, i);

     return x;

   }

   MipCplex::Value MipCplex::_getPrimalValue() const {

     Value objval;

     CPXgetmipobjval(env, lp, &objval);

     return objval;

   }

 } //END OF NAMESPACE LEMON