athos@1241: // -*- c++ -*-
athos@1241: #ifndef LEMON_LP_SOLVER_GLPK_H
athos@1241: #define LEMON_LP_SOLVER_GLPK_H
athos@1241: 
athos@1241: ///\ingroup misc
athos@1241: ///\file
athos@1241: 
athos@1241: // #include <stdio.h>
athos@1241: /* #include <stdlib.h> */
athos@1241: /* #include <iostream> */
athos@1241: /* #include <map> */
athos@1241: /* #include <limits> */
athos@1241: // #include <stdio>
athos@1241: //#include <stdlib>
athos@1241: extern "C" {
athos@1241: #include "glpk.h"
athos@1241: }
athos@1241: 
athos@1241: /* #include <iostream> */
athos@1241: /* #include <vector> */
athos@1241: /* #include <string> */
athos@1241: /* #include <list> */
athos@1241: /* #include <memory> */
athos@1241: /* #include <utility> */
athos@1241: 
athos@1241: //#include <lemon/invalid.h>
athos@1241: //#include <expression.h>
athos@1241: #include <lp_solver_base.h>
athos@1241: //#include <stp.h>
athos@1241: //#include <lemon/max_flow.h>
athos@1241: //#include <augmenting_flow.h>
athos@1241: //#include <iter_map.h>
athos@1241: 
athos@1241: using std::cout;
athos@1241: using std::cin;
athos@1241: using std::endl;
athos@1241: 
athos@1241: namespace lemon {
athos@1241:   
athos@1241:   
athos@1241:   template <typename Value>
athos@1241:   const Value LpSolverBase<Value>::INF=std::numeric_limits<Value>::infinity();
athos@1241: 
athos@1241: 
athos@1241:   /// \brief Wrapper for GLPK solver
athos@1241:   /// 
athos@1241:   /// This class implements a lemon wrapper for GLPK.
athos@1241:   class LpGlpk : public LpSolverBase<double> {
athos@1241:   public:
athos@1241:     typedef LpSolverBase<double> Parent;
athos@1241: 
athos@1241:   public:
athos@1241:     /// \e
athos@1241:     LPX* lp;
athos@1241: 
athos@1241:   public:
athos@1241:     /// \e
athos@1241:     LpGlpk() : Parent(), 
athos@1241: 			lp(lpx_create_prob()) {
athos@1241:       int_row_map.push_back(Row());
athos@1241:       int_col_map.push_back(Col());
athos@1241:       lpx_set_int_parm(lp, LPX_K_DUAL, 1);
athos@1241:     }
athos@1241:     /// \e
athos@1241:     ~LpGlpk() {
athos@1241:       lpx_delete_prob(lp);
athos@1241:     }
athos@1241: 
athos@1241:     //MATRIX INDEPEDENT MANIPULATING FUNCTIONS
athos@1241: 
athos@1241:     /// \e
athos@1241:     void setMinimize() { 
athos@1241:       lpx_set_obj_dir(lp, LPX_MIN);
athos@1241:     }
athos@1241:     /// \e
athos@1241:     void setMaximize() { 
athos@1241:       lpx_set_obj_dir(lp, LPX_MAX);
athos@1241:     }
athos@1241: 
athos@1241:     //LOW LEVEL INTERFACE, MATRIX MANIPULATING FUNCTIONS
athos@1241: 
athos@1241:   protected:
athos@1241:     /// \e
athos@1241:     int _addCol() { 
athos@1241:       int i=lpx_add_cols(lp, 1);
athos@1241:       _setColLowerBound(i, -INF);
athos@1241:       _setColUpperBound(i, INF);
athos@1241:       return i;
athos@1241:     }
athos@1241:     /// \e
athos@1241:     int _addRow() { 
athos@1241:       int i=lpx_add_rows(lp, 1);
athos@1241:       return i;
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual void _setRowCoeffs(int i, 
athos@1241: 			       const std::vector<std::pair<int, double> >& coeffs) {
athos@1241:       int mem_length=1+colNum();
athos@1241:       int* indices = new int[mem_length];
athos@1241:       double* doubles = new double[mem_length];
athos@1241:       int length=0;
athos@1241:       for (std::vector<std::pair<int, double> >::
athos@1241: 	     const_iterator it=coeffs.begin(); it!=coeffs.end(); ++it) {
athos@1241: 	++length;
athos@1241: 	indices[length]=it->first;
athos@1241: 	doubles[length]=it->second;
athos@1241:       }
athos@1241:       lpx_set_mat_row(lp, i, length, indices, doubles);
athos@1241:       delete [] indices;
athos@1241:       delete [] doubles;
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual void _getRowCoeffs(int i, 
athos@1241: 			       std::vector<std::pair<int, double> >& coeffs) {
athos@1241:       int mem_length=1+colNum();
athos@1241:       int* indices = new int[mem_length];
athos@1241:       double* doubles = new double[mem_length];
athos@1241:       int length=lpx_get_mat_row(lp, i, indices, doubles);
athos@1241:       for (int i=1; i<=length; ++i) {
athos@1241: 	coeffs.push_back(std::make_pair(indices[i], doubles[i]));
athos@1241:       }
athos@1241:       delete [] indices;
athos@1241:       delete [] doubles;
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual void _setColCoeffs(int i, 
athos@1241: 			       const std::vector<std::pair<int, double> >& coeffs) {
athos@1241:       int mem_length=1+rowNum();
athos@1241:       int* indices = new int[mem_length];
athos@1241:       double* doubles = new double[mem_length];
athos@1241:       int length=0;
athos@1241:       for (std::vector<std::pair<int, double> >::
athos@1241: 	     const_iterator it=coeffs.begin(); it!=coeffs.end(); ++it) {
athos@1241: 	++length;
athos@1241: 	indices[length]=it->first;
athos@1241: 	doubles[length]=it->second;
athos@1241:       }
athos@1241:       lpx_set_mat_col(lp, i, length, indices, doubles);
athos@1241:       delete [] indices;
athos@1241:       delete [] doubles;
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual void _getColCoeffs(int i, 
athos@1241: 			       std::vector<std::pair<int, double> >& coeffs) {
athos@1241:       int mem_length=1+rowNum();
athos@1241:       int* indices = new int[mem_length];
athos@1241:       double* doubles = new double[mem_length];
athos@1241:       int length=lpx_get_mat_col(lp, i, indices, doubles);
athos@1241:       for (int i=1; i<=length; ++i) {
athos@1241: 	coeffs.push_back(std::make_pair(indices[i], doubles[i]));
athos@1241:       }
athos@1241:       delete [] indices;
athos@1241:       delete [] doubles;
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual void _eraseCol(int i) {
athos@1241:       int cols[2];
athos@1241:       cols[1]=i;
athos@1241:       lpx_del_cols(lp, 1, cols);
athos@1241:     }
athos@1241:     virtual void _eraseRow(int i) {
athos@1241:       int rows[2];
athos@1241:       rows[1]=i;
athos@1241:       lpx_del_rows(lp, 1, rows);
athos@1241:     }
athos@1241:     void _setCoeff(int col, int row, double value) {
athos@1243:       ///FIXME Of course this is not efficient at all, but GLPK knows not more.
athos@1243:       int change_this;
athos@1243:       bool get_set_row;
athos@1243:       //The only thing we can do is optimize on whether working with a row 
athos@1243:       //or a coloumn
athos@1243:       int row_num = rowNum();
athos@1243:       int col_num = colNum();
athos@1243:       if (col_num<row_num){
athos@1243: 	//There are more rows than coloumns
athos@1243: 	get_set_row=true;
athos@1243: 	int mem_length=1+row_num;
athos@1243: 	int* indices = new int[mem_length];
athos@1243: 	double* doubles = new double[mem_length];
athos@1243: 	int length=lpx_get_mat_col(lp, i, indices, doubles);
athos@1243:       }else{
athos@1243: 	get_set_row=false;
athos@1243: 	int mem_length=1+col_num;
athos@1243: 	int* indices = new int[mem_length];
athos@1243: 	double* doubles = new double[mem_length];
athos@1243: 	int length=lpx_get_mat_row(lp, i, indices, doubles);
athos@1243:       }
athos@1243:       //Itten      
athos@1243: int* indices = new int[mem_length];
athos@1243:       double* doubles = new double[mem_length];
athos@1243:       int length=lpx_get_mat_col(lp, i, indices, doubles);
athos@1243:  
athos@1243:       delete [] indices;
athos@1243:       delete [] doubles;
athos@1243: 
athos@1241:     }
athos@1241:     double _getCoeff(int col, int row) {
athos@1241:       /// FIXME not yet implemented
athos@1241:       return 0.0;
athos@1241:     }
athos@1241:     virtual void _setColLowerBound(int i, double lo) {
athos@1241:       if (lo==INF) {
athos@1241: 	//FIXME error
athos@1241:       }
athos@1241:       int b=lpx_get_col_type(lp, i);
athos@1241:       double up=lpx_get_col_ub(lp, i);	
athos@1241:       if (lo==-INF) {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	case LPX_LO:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_FR, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_UP:
athos@1241: 	  break;
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_UP, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       } else {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	case LPX_LO:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_LO, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_UP:	  
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  if (lo==up) 
athos@1241: 	    lpx_set_col_bnds(lp, i, LPX_FX, lo, up);
athos@1241: 	  else 
athos@1241: 	    lpx_set_col_bnds(lp, i, LPX_DB, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       }
athos@1241:     }
athos@1241:     virtual double _getColLowerBound(int i) {
athos@1241:       int b=lpx_get_col_type(lp, i);
athos@1241:       switch (b) {
athos@1241:       case LPX_FR:
athos@1241: 	return -INF;
athos@1241:       case LPX_LO:
athos@1241: 	return lpx_get_col_lb(lp, i);
athos@1241:       case LPX_UP:
athos@1241: 	return -INF;
athos@1241:       case LPX_DB:
athos@1241:       case LPX_FX:
athos@1241: 	return lpx_get_col_lb(lp, i);
athos@1241:       default: ;
athos@1241: 	//FIXME error
athos@1241: 	return 0.0;
athos@1241:       }
athos@1241:     }
athos@1241:     virtual void _setColUpperBound(int i, double up) {
athos@1241:       if (up==-INF) {
athos@1241: 	//FIXME error
athos@1241:       }
athos@1241:       int b=lpx_get_col_type(lp, i);
athos@1241:       double lo=lpx_get_col_lb(lp, i);
athos@1241:       if (up==INF) {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	case LPX_LO:
athos@1241: 	  break;
athos@1241: 	case LPX_UP:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_FR, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_LO, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       } else {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_UP, lo, up);
athos@1241: 	case LPX_LO:
athos@1241: 	  if (lo==up) 
athos@1241: 	    lpx_set_col_bnds(lp, i, LPX_FX, lo, up);
athos@1241: 	  else
athos@1241: 	    lpx_set_col_bnds(lp, i, LPX_DB, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_UP:
athos@1241: 	  lpx_set_col_bnds(lp, i, LPX_UP, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  if (lo==up) 
athos@1241: 	    lpx_set_col_bnds(lp, i, LPX_FX, lo, up);
athos@1241: 	  else 
athos@1241: 	    lpx_set_col_bnds(lp, i, LPX_DB, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       }
athos@1241:     }
athos@1241:     virtual double _getColUpperBound(int i) {
athos@1241:       int b=lpx_get_col_type(lp, i);
athos@1241:       switch (b) {
athos@1241:       case LPX_FR:
athos@1241:       case LPX_LO:
athos@1241: 	return INF;
athos@1241:       case LPX_UP:
athos@1241:       case LPX_DB:
athos@1241:       case LPX_FX:
athos@1241: 	return lpx_get_col_ub(lp, i);
athos@1241:       default: ;
athos@1241: 	//FIXME error
athos@1241: 	return 0.0;
athos@1241:       }
athos@1241:     }
athos@1241:     virtual void _setRowLowerBound(int i, double lo) {
athos@1241:       if (lo==INF) {
athos@1241: 	//FIXME error
athos@1241:       }
athos@1241:       int b=lpx_get_row_type(lp, i);
athos@1241:       double up=lpx_get_row_ub(lp, i);	
athos@1241:       if (lo==-INF) {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	case LPX_LO:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_FR, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_UP:
athos@1241: 	  break;
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_UP, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       } else {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	case LPX_LO:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_LO, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_UP:	  
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  if (lo==up) 
athos@1241: 	    lpx_set_row_bnds(lp, i, LPX_FX, lo, up);
athos@1241: 	  else 
athos@1241: 	    lpx_set_row_bnds(lp, i, LPX_DB, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       }
athos@1241:     }
athos@1241:     virtual double _getRowLowerBound(int i) {
athos@1241:       int b=lpx_get_row_type(lp, i);
athos@1241:       switch (b) {
athos@1241:       case LPX_FR:
athos@1241: 	return -INF;
athos@1241:       case LPX_LO:
athos@1241: 	return lpx_get_row_lb(lp, i);
athos@1241:       case LPX_UP:
athos@1241: 	return -INF;
athos@1241:       case LPX_DB:
athos@1241:       case LPX_FX:
athos@1241: 	return lpx_get_row_lb(lp, i);
athos@1241:       default: ;
athos@1241: 	//FIXME error
athos@1241: 	return 0.0;
athos@1241:       }
athos@1241:     }
athos@1241:     virtual void _setRowUpperBound(int i, double up) {
athos@1241:       if (up==-INF) {
athos@1241: 	//FIXME error
athos@1241:       }
athos@1241:       int b=lpx_get_row_type(lp, i);
athos@1241:       double lo=lpx_get_row_lb(lp, i);
athos@1241:       if (up==INF) {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	case LPX_LO:
athos@1241: 	  break;
athos@1241: 	case LPX_UP:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_FR, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_LO, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       } else {
athos@1241: 	switch (b) {
athos@1241: 	case LPX_FR:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_UP, lo, up);
athos@1241: 	case LPX_LO:
athos@1241: 	  if (lo==up) 
athos@1241: 	    lpx_set_row_bnds(lp, i, LPX_FX, lo, up);
athos@1241: 	  else
athos@1241: 	    lpx_set_row_bnds(lp, i, LPX_DB, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_UP:
athos@1241: 	  lpx_set_row_bnds(lp, i, LPX_UP, lo, up);
athos@1241: 	  break;
athos@1241: 	case LPX_DB:
athos@1241: 	case LPX_FX:
athos@1241: 	  if (lo==up) 
athos@1241: 	    lpx_set_row_bnds(lp, i, LPX_FX, lo, up);
athos@1241: 	  else 
athos@1241: 	    lpx_set_row_bnds(lp, i, LPX_DB, lo, up);
athos@1241: 	  break;
athos@1241: 	default: ;
athos@1241: 	  //FIXME error
athos@1241: 	}
athos@1241:       }
athos@1241:     }
athos@1241:     virtual double _getRowUpperBound(int i) {
athos@1241:       int b=lpx_get_row_type(lp, i);
athos@1241:       switch (b) {
athos@1241:       case LPX_FR:
athos@1241:       case LPX_LO:
athos@1241: 	return INF;
athos@1241:       case LPX_UP:
athos@1241:       case LPX_DB:
athos@1241:       case LPX_FX:
athos@1241: 	return lpx_get_row_ub(lp, i);
athos@1241:       default: ;
athos@1241: 	//FIXME error
athos@1241: 	return 0.0;
athos@1241:       }
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual double _getObjCoeff(int i) { 
athos@1241:       return lpx_get_obj_coef(lp, i);
athos@1241:     }
athos@1241:     /// \e
athos@1241:     virtual void _setObjCoeff(int i, double obj_coef) { 
athos@1241:       lpx_set_obj_coef(lp, i, obj_coef);
athos@1241:     }
athos@1241:   public:
athos@1241:     /// \e
athos@1241:     void solveSimplex() { lpx_simplex(lp); }
athos@1241:     /// \e
athos@1241:     void solvePrimalSimplex() { lpx_simplex(lp); }
athos@1241:     /// \e
athos@1241:     void solveDualSimplex() { lpx_simplex(lp); }
athos@1241:   protected:
athos@1241:     virtual double _getPrimal(int i) {
athos@1241:       return lpx_get_col_prim(lp, i);
athos@1241:     }
athos@1241:   public:
athos@1241:     /// \e
athos@1241:     double getObjVal() { return lpx_get_obj_val(lp); }
athos@1241:     /// \e
athos@1241:     int rowNum() const { return lpx_get_num_rows(lp); }
athos@1241:     /// \e
athos@1241:     int colNum() const { return lpx_get_num_cols(lp); }
athos@1241:     /// \e
athos@1241:     int warmUp() { return lpx_warm_up(lp); }
athos@1241:     /// \e
athos@1241:     void printWarmUpStatus(int i) {
athos@1241:       switch (i) {
athos@1241:       case LPX_E_OK: cout << "LPX_E_OK" << endl; break;
athos@1241:       case LPX_E_EMPTY: cout << "LPX_E_EMPTY" << endl; break;	
athos@1241:       case LPX_E_BADB: cout << "LPX_E_BADB" << endl; break;
athos@1241:       case LPX_E_SING: cout << "LPX_E_SING" << endl; break;
athos@1241:       }
athos@1241:     }
athos@1241:     /// \e
athos@1241:     int getPrimalStatus() { return lpx_get_prim_stat(lp); }
athos@1241:     /// \e
athos@1241:     void printPrimalStatus(int i) {
athos@1241:       switch (i) {
athos@1241:       case LPX_P_UNDEF: cout << "LPX_P_UNDEF" << endl; break;
athos@1241:       case LPX_P_FEAS: cout << "LPX_P_FEAS" << endl; break;	
athos@1241:       case LPX_P_INFEAS: cout << "LPX_P_INFEAS" << endl; break;
athos@1241:       case LPX_P_NOFEAS: cout << "LPX_P_NOFEAS" << endl; break;
athos@1241:       }
athos@1241:     }
athos@1241:     /// \e
athos@1241:     int getDualStatus() { return lpx_get_dual_stat(lp); }
athos@1241:     /// \e
athos@1241:     void printDualStatus(int i) {
athos@1241:       switch (i) {
athos@1241:       case LPX_D_UNDEF: cout << "LPX_D_UNDEF" << endl; break;
athos@1241:       case LPX_D_FEAS: cout << "LPX_D_FEAS" << endl; break;	
athos@1241:       case LPX_D_INFEAS: cout << "LPX_D_INFEAS" << endl; break;
athos@1241:       case LPX_D_NOFEAS: cout << "LPX_D_NOFEAS" << endl; break;
athos@1241:       }
athos@1241:     }
athos@1241:     /// Returns the status of the slack variable assigned to row \c row.
athos@1241:     int getRowStat(const Row& row) { 
athos@1241:       return lpx_get_row_stat(lp, row_iter_map[row]); 
athos@1241:     }
athos@1241:     /// \e
athos@1241:     void printRowStatus(int i) {
athos@1241:       switch (i) {
athos@1241:       case LPX_BS: cout << "LPX_BS" << endl; break;
athos@1241:       case LPX_NL: cout << "LPX_NL" << endl; break;	
athos@1241:       case LPX_NU: cout << "LPX_NU" << endl; break;
athos@1241:       case LPX_NF: cout << "LPX_NF" << endl; break;
athos@1241:       case LPX_NS: cout << "LPX_NS" << endl; break;
athos@1241:       }
athos@1241:     }
athos@1241:     /// Returns the status of the variable assigned to column \c col.
athos@1241:     int getColStat(const Col& col) { 
athos@1241:       return lpx_get_col_stat(lp, col_iter_map[col]); 
athos@1241:     }
athos@1241:     /// \e
athos@1241:     void printColStatus(int i) {
athos@1241:       switch (i) {
athos@1241:       case LPX_BS: cout << "LPX_BS" << endl; break;
athos@1241:       case LPX_NL: cout << "LPX_NL" << endl; break;	
athos@1241:       case LPX_NU: cout << "LPX_NU" << endl; break;
athos@1241:       case LPX_NF: cout << "LPX_NF" << endl; break;
athos@1241:       case LPX_NS: cout << "LPX_NS" << endl; break;
athos@1241:       }
athos@1241:     }
athos@1241: 
athos@1241:     // MIP
athos@1241:     /// \e
athos@1241:     void solveBandB() { lpx_integer(lp); }
athos@1241:     /// \e
athos@1241:     void setLP() { lpx_set_class(lp, LPX_LP); }
athos@1241:     /// \e
athos@1241:     void setMIP() { lpx_set_class(lp, LPX_MIP); }
athos@1241:   protected:
athos@1241:     /// \e
athos@1241:     void _setColCont(int i) { lpx_set_col_kind(lp, i, LPX_CV); }
athos@1241:     /// \e
athos@1241:     void _setColInt(int i) { lpx_set_col_kind(lp, i, LPX_IV); }
athos@1241:     /// \e
athos@1241:     double _getMIPPrimal(int i) { return lpx_mip_col_val(lp, i); }
athos@1241:   };
athos@1241:   
athos@1241:   /// @}
athos@1241: 
athos@1241: } //namespace lemon
athos@1241: 
athos@1241: #endif //LEMON_LP_SOLVER_GLPK_H