/* -*- 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-GLPK mip solver interface.

#include <lemon/mip_glpk.h>

#if GLP_MAJOR_VERSION > 4 || (GLP_MAJOR_VERSION == 4 && GLP_MINOR_VERSION > 15)

#define LEMON_glp(func) (glp_##func)

#define LEMON_lpx(func) (lpx_##func)

#define LEMON_GLP(def) (GLP_##def)

#define LEMON_LPX(def) (LPX_##def)

#else

#define LEMON_glp(func) (lpx_##func)

#define LEMON_lpx(func) (lpx_##func)

#define LEMON_GLP(def) (LPX_##def)

#define LEMON_LPX(def) (LPX_##def)

#endif

namespace lemon {

  MipGlpk::MipGlpk() {

#if !(GLP_MAJOR_VERSION > 4 || (GLP_MAJOR_VERSION == 4 && GLP_MINOR_VERSION > 15))

    LEMON_lpx(set_class)(lp,LEMON_GLP(MIP));

#endif

  }

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

    switch (col_type){

      case INT:

	LEMON_glp(set_col_kind)(lp,i,LEMON_GLP(IV));

	break;

      case REAL:

	LEMON_glp(set_col_kind)(lp,i,LEMON_GLP(CV));

	break;

    default:;

        //FIXME problem

    }

  }

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

    switch (LEMON_glp(get_col_kind)(lp,i)){

    case LEMON_GLP(IV):

      return INT;//Or binary

    case LEMON_GLP(CV):

      return REAL;

    default:

      return REAL;//Error!

    }

  }

  LpGlpk::SolveExitStatus MipGlpk::_solve() {

    int result = LEMON_lpx(simplex)(lp);

    // hack: mip does not contain integer variable 

#if GLP_MAJOR_VERSION > 4 || (GLP_MAJOR_VERSION == 4 && GLP_MINOR_VERSION > 15)

    int tmp = -1;

    if (LEMON_glp(get_num_int(lp)) == 0) {

      tmp = LEMON_lpx(add_cols)(lp, 1);

      LEMON_glp(set_col_bnds)(lp, tmp, LEMON_GLP(FX), 0.0, 0.0);

      LEMON_glp(set_col_kind)(lp, tmp, LEMON_GLP(IV));

    }

#endif

    if (LEMON_lpx(get_status)(lp)==LEMON_LPX(OPT)) {

      //Maybe we could try the routine lpx_intopt(lp), a revised

      //version of lpx_integer

      result = LEMON_lpx(integer)(lp);

      switch (result){

      case LEMON_LPX(E_OK):

	solved = true;

	break;

      default:

	solved = false;

      }  

    } else {

      solved = false;

    }

#if GLP_MAJOR_VERSION > 4 || (GLP_MAJOR_VERSION == 4 && GLP_MINOR_VERSION > 15)

    if (tmp != -1) {

      int tmpa[2];

      tmpa[1] = tmp;

      LEMON_lpx(del_cols)(lp, 1, tmpa);

    }

#endif

    return solved ? SOLVED : UNSOLVED;

  }

  LpGlpk::SolutionStatus MipGlpk::_getMipStatus() const {

    if (LEMON_lpx(get_status)(lp)==LEMON_LPX(OPT)){

      //Meg kell nezni: ha az LP is infinite, akkor ez is, ha az is

      //infeasible, akkor ez is, de ez lehet maskepp is infeasible.

      int stat= LEMON_lpx(mip_status)(lp);

      switch (stat) {

      case LEMON_LPX(I_UNDEF)://Undefined (no solve has been run yet)

	return UNDEFINED;

      case LEMON_LPX(I_NOFEAS)://There is no feasible integral solution

	return INFEASIBLE;

	//     case LEMON_LPX(UNBND)://Unbounded

	//       return INFINITE;

      case LEMON_LPX(I_FEAS)://Feasible

	return FEASIBLE;

      case LEMON_LPX(I_OPT)://Feasible

	return OPTIMAL;

      default:

	return UNDEFINED; //to avoid gcc warning

      //FIXME error

      }

    }

    else 

      return UNDEFINED; //Maybe we could refine this: what does the LP

			//relaxation look like

  }  

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

    return LEMON_glp(mip_col_val)(lp,i);

  }

  MipGlpk::Value MipGlpk::_getPrimalValue() const {

    return LEMON_glp(mip_obj_val)(lp);

  }

} //END OF NAMESPACE LEMON