RhodeCode

namespace lemon {

  CbcMip::CbcMip() {

    _prob = new CoinModel();

    _prob->setProblemName("LEMON");

    _osi_solver = 0;

    _cbc_model = 0;

  }

  CbcMip::CbcMip(const CbcMip& other) {

    _prob = new CoinModel(*other._prob);

    _osi_solver = 0;

    _cbc_model = 0;

  }

  CbcMip::~CbcMip() {

    delete _prob;

    if (_osi_solver) delete _osi_solver;

    if (_cbc_model) delete _cbc_model;

  }

    _osi_solver->loadFromCoinModel(*_prob);

    if (_cbc_model) {

      delete _cbc_model;

    }

    _cbc_model= new CbcModel(*_osi_solver);

    _cbc_model->initialSolve();

    _cbc_model->solver()->setHintParam(OsiDoReducePrint, true, OsiHintTry);

    if (!_cbc_model->isInitialSolveAbandoned() &&

        _cbc_model->isInitialSolveProvenOptimal() &&

        !_cbc_model->isInitialSolveProvenPrimalInfeasible() &&

        !_cbc_model->isInitialSolveProvenDualInfeasible()) {

      _cbc_model = 0;

    }

    _prob = new CoinModel();

    rows.clear();

    cols.clear();

  }

  }

} //END OF NAMESPACE LEMON

    virtual SolveExitStatus _solve();

    virtual ProblemType _getType() const;

    virtual Value _getSol(int i) const;

    virtual Value _getSolValue() const;

    virtual void _clear();

  };

}

#endif

#include <lemon/clp.h>

#include <coin/ClpSimplex.hpp>

namespace lemon {

  ClpLp::ClpLp() {

    _prob = new ClpSimplex();

    _init_temporals();

  }

  ClpLp::ClpLp(const ClpLp& other) {

    _prob = new ClpSimplex(*other._prob);

    rows = other.rows;

    cols = other.cols;

    _init_temporals();

  }

  ClpLp::~ClpLp() {

    delete _prob;

    _clear_temporals();

  }

  void ClpLp::_init_temporals() {

    delete _prob;

    _prob = new ClpSimplex();

    rows.clear();

    cols.clear();

    _col_names_ref.clear();

    _clear_temporals();

  }

  }

} //END OF NAMESPACE LEMON

    virtual VarStatus _getColStatus(int i) const;

    virtual VarStatus _getRowStatus(int i) const;

    virtual ProblemType _getPrimalType() const;

    virtual ProblemType _getDualType() const;

    virtual void _clear();

  public:

    ///Solves LP with primal simplex method.

    SolveExitStatus solvePrimal();

    ///Solves LP with dual simplex method.

    SolveExitStatus solveDual();

    SolveExitStatus solveBarrier();

    ///Returns the constraint identifier understood by CLP.

    int clpRow(Row r) const { return rows(id(r)); }

    ///Returns the variable identifier understood by CLP.

    int clpCol(Col c) const { return cols(id(c)); }

  };

} //END OF NAMESPACE LEMON

#endif //LEMON_CLP_H

      delete _cnt;

      CPXcloseCPLEX(&_env);

    }

  }

  CplexBase::CplexBase() : LpBase() {

    int status;

    _prob = CPXcreateprob(cplexEnv(), &status, "Cplex problem");

  }

  CplexBase::CplexBase(const CplexEnv& env)

    : LpBase(), _env(env) {

    int status;

    _prob = CPXcreateprob(cplexEnv(), &status, "Cplex problem");

  }

  CplexBase::CplexBase(const CplexBase& cplex)

    : LpBase() {

    int status;

    _prob = CPXcloneprob(cplexEnv(), cplex._prob, &status);

    rows = cplex.rows;

    cols = cplex.cols;

  }

  CplexBase::~CplexBase() {

    CPXfreeprob(cplexEnv(),&_prob);

  }

  int CplexBase::_addCol() {

    int i = CPXgetnumcols(cplexEnv(), _prob);

  void CplexBase::_clear() {

    CPXfreeprob(cplexEnv(),&_prob);

    int status;

    _prob = CPXcreateprob(cplexEnv(), &status, "Cplex problem");

    rows.clear();

    cols.clear();

  }

  // CplexLp members

  CplexLp::CplexLp()

    : LpBase(), LpSolver(), CplexBase() {}

  CplexLp::CplexLp(const CplexEnv& env)

    : LpBase(), LpSolver(), CplexBase(env) {}

    } else {

      return UNSOLVED;

    }

#endif

  }

  CplexLp::SolveExitStatus CplexLp::_solve() {

    _clear_temporals();

    return convertStatus(CPXlpopt(cplexEnv(), _prob));

  }

  CplexLp::SolveExitStatus CplexLp::solvePrimal() {

    _clear_temporals();

    return convertStatus(CPXprimopt(cplexEnv(), _prob));

  }

  CplexLp::SolveExitStatus CplexLp::solveDual() {

    _clear_temporals();

    return convertStatus(CPXdualopt(cplexEnv(), _prob));

  }

  CplexLp::SolveExitStatus CplexLp::solveBarrier() {

    _clear_temporals();

    return convertStatus(CPXbaropt(cplexEnv(), _prob));

  }

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

    Value x;

    CPXgetx(cplexEnv(), _prob, &x, i, i);

    return x;

  }

  CplexLp::Value CplexLp::_getDualRay(int i) const {

    if (_dual_ray.empty()) {

    }

    return _dual_ray[i];

  }

  // This table lists the statuses, returned by the CPXgetstat()

  // routine, for solutions to LP problems or mixed integer problems. If

  // no solution exists, the return value is zero.

  // For Simplex, Barrier

  // 1          CPX_OPTIMAL

  //          Optimal solution found

  // 2          CPX_INFEASIBLE

  //          Aborted in barrier, primal and dual feasible

  // 18          CPX_ABORT_CROSSOVER

  //          Aborted in crossover

  // 19          CPX_INForUNBD

  //          Infeasible or unbounded

  // 20   CPX_PIVOT

  //       User pivot used

  //

  // ??case CPX_ABORT_DUAL_INFEAS

  // ??case CPX_ABORT_CROSSOVER

  // ??case CPX_INForUNBD

  // ??case CPX_PIVOT

  //Some more interesting stuff:

  // CPX_PARAM_PROBMETHOD  1062  int  LPMETHOD

      case CPX_STAT_INFEASIBLE:

        return INFEASIBLE;

      default:

        return UNDEFINED;

      }

#else

    statusSwitch(cplexEnv(),stat);

    //CPXgetstat(cplexEnv(), _prob);

    switch (stat) {

    case 0:

      return UNDEFINED; //Undefined

    case CPX_OPTIMAL://Optimal

      return OPTIMAL;

    case CPX_UNBOUNDED://Unbounded

      return INFEASIBLE;//In case of dual simplex

      //return UNBOUNDED;

      //       return FEASIBLE;

    default:

      return UNDEFINED; //Everything else comes here

      //FIXME error

    }

#endif

  }

  // CPX_STAT_ABORT_DUAL_OBJ_LIM

  // CPX_STAT_ABORT_IT_LIM

  // CPX_STAT_ABORT_OBJ_LIM

  // CPX_STAT_ABORT_PRIM_OBJ_LIM

  // CPX_STAT_ABORT_TIME_LIM

  // CPX_STAT_ABORT_USER

  // CPX_STAT_FEASIBLE_RELAXED

  // CPX_STAT_INFEASIBLE

      LEMON_ASSERT(false, "Invalid column type");

      return ColTypes();

    }

  }

  CplexMip::SolveExitStatus CplexMip::_solve() {

    int status;

    status = CPXmipopt (cplexEnv(), _prob);

    if (status==0)

      return SOLVED;

    else

      return UNSOLVED;

  }

    virtual void _setObjCoeff(int i, Value obj_coef);

    virtual Value _getObjCoeff(int i) const;

    virtual void _setSense(Sense sense);

    virtual Sense _getSense() const;

    virtual void _clear();

  public:

    /// Returns the used \c CplexEnv instance

    const CplexEnv& env() const { return _env; }

    ///

    const cpxenv* cplexEnv() const { return _env.cplexEnv(); }

    cpxlp* cplexLp() { return _prob; }

    const cpxlp* cplexLp() const { return _prob; }

  };

  /// \brief Interface for the CPLEX LP solver

  ///

  /// This class implements an interface for the CPLEX LP solver.

  ///\ingroup lp_group

namespace lemon {

  // GlpkBase members

  GlpkBase::GlpkBase() : LpBase() {

    lp = glp_create_prob();

    glp_create_index(lp);

  }

  GlpkBase::GlpkBase(const GlpkBase &other) : LpBase() {

    lp = glp_create_prob();

    glp_copy_prob(lp, other.lp, GLP_ON);

    glp_create_index(lp);

    rows = other.rows;

    cols = other.cols;

  }

  GlpkBase::~GlpkBase() {

    glp_delete_prob(lp);

  }

  int GlpkBase::_addCol() {

    int i = glp_add_cols(lp, 1);

    rows.clear();

    cols.clear();

  }

  void GlpkBase::freeEnv() {

    glp_free_env();

  }

  GlpkBase::FreeEnvHelper GlpkBase::freeEnvHelper;

  // GlpkLp members

  GlpkLp::GlpkLp()

    : LpBase(), LpSolver(), GlpkBase() {

    presolver(false);

  }

  GlpkLp::GlpkLp(const GlpkLp& other)

    : LpBase(other), LpSolver(other), GlpkBase(other) {

    presolver(false);

  }

  GlpkLp* GlpkLp::newSolver() const { return new GlpkLp; }

  GlpkLp* GlpkLp::cloneSolver() const { return new GlpkLp(*this); }

  const char* GlpkLp::_solverName() const { return "GlpkLp"; }

  }

  GlpkLp::SolveExitStatus GlpkLp::solvePrimal() {

    _clear_temporals();

    glp_smcp smcp;

    glp_init_smcp(&smcp);

    smcp.presolve = _presolve;

    // If the basis is not valid we get an error return value.

    // In this case we can try to create a new basis.

    switch (glp_simplex(lp, &smcp)) {

    case 0:

      break;

    case GLP_EBADB:

  }

  GlpkLp::SolveExitStatus GlpkLp::solveDual() {

    _clear_temporals();

    glp_smcp smcp;

    glp_init_smcp(&smcp);

    smcp.meth = GLP_DUAL;

    smcp.presolve = _presolve;

    // If the basis is not valid we get an error return value.

    // In this case we can try to create a new basis.

    switch (glp_simplex(lp, &smcp)) {

    case 0:

      break;

      return  GlpkLp::ProblemType();

    }

  }

  void GlpkLp::presolver(bool presolve) {

    _presolve = presolve;

  }

  // GlpkMip members

  GlpkMip::GlpkMip()

    : LpBase(), MipSolver(), GlpkBase() {

  }

  GlpkMip::GlpkMip(const GlpkMip& other)

    : LpBase(), MipSolver(), GlpkBase(other) {

  }

  void GlpkMip::_setColType(int i, GlpkMip::ColTypes col_type) {

    switch (col_type) {

    case INTEGER:

      glp_set_col_kind(lp, i, GLP_IV);

      break;

    case REAL:

    }

  }

  GlpkMip::SolveExitStatus GlpkMip::_solve() {

    glp_smcp smcp;

    glp_init_smcp(&smcp);

    smcp.meth = GLP_DUAL;

    // If the basis is not valid we get an error return value.

    // In this case we can try to create a new basis.

    switch (glp_simplex(lp, &smcp)) {

    case 0:

      break;

    case GLP_EBADB:

      return UNSOLVED;

    }

    if (glp_get_status(lp) != GLP_OPT) return SOLVED;

    glp_iocp iocp;

    glp_init_iocp(&iocp);

    if (glp_intopt(lp, &iocp) != 0) return UNSOLVED;

    return SOLVED;

  }

  GlpkMip::ProblemType GlpkMip::_getType() const {

    switch (glp_get_status(lp)) {

    return glp_mip_obj_val(lp);

  }

  GlpkMip* GlpkMip::newSolver() const { return new GlpkMip; }

  GlpkMip* GlpkMip::cloneSolver() const {return new GlpkMip(*this); }

  const char* GlpkMip::_solverName() const { return "GlpkMip"; }

} //END OF NAMESPACE LEMON

    virtual void _setObjCoeff(int i, Value obj_coef);

    virtual Value _getObjCoeff(int i) const;

    virtual void _setSense(Sense);

    virtual Sense _getSense() const;

    virtual void _clear();

  private:

    static void freeEnv();

    struct FreeEnvHelper {

      ~FreeEnvHelper() {

        freeEnv();

      }

    };

    static FreeEnvHelper freeEnvHelper;

  public:

    ///Pointer to the underlying GLPK data structure.

    LPX *lpx() {return lp;}

    ///Const pointer to the underlying GLPK data structure.

    const LPX *lpx() const {return lp;}

    ///Turns on or off the presolver

    ///Turns on (\c b is \c true) or off (\c b is \c false) the presolver

    ///

    ///The presolver is off by default.

    void presolver(bool presolve);

  };

  /// \brief Interface for the GLPK MIP solver

  ///

  /// This class implements an interface for the GLPK MIP solver.

  ///\ingroup lp_group

  class GlpkMip : public MipSolver, public GlpkBase {

  public:

    virtual ColTypes _getColType(int col) const;

    virtual void _setColType(int col, ColTypes col_type);

    virtual SolveExitStatus _solve();

    virtual ProblemType _getType() const;

    virtual Value _getSol(int i) const;

    virtual Value _getSolValue() const;

  };

} //END OF NAMESPACE LEMON

#endif //LEMON_GLPK_H

    ///Direction of the optimization

    enum Sense {

      /// Minimization

      MIN,

      /// Maximization

      MAX

    };

    ///The floating point type used by the solver

    typedef double Value;

    ///The infinity constant

    static const Value INF;

    ///The not a number constant

    static const Value NaN;

    friend class Col;

    virtual void _setSense(Sense) = 0;

    virtual Sense _getSense() const = 0;

    virtual void _clear() = 0;

    virtual const char* _solverName() const = 0;

    //Own protected stuff

    //Constant component of the objective function

    Value obj_const_comp;

    LpBase() : rows(), cols(), obj_const_comp(0) {}

  public:

    void max() { _setSense(MAX); }

    ///Set the sense to maximization

    void min() { _setSense(MIN); }

    ///Clears the problem

    void clear() { _clear(); }

    ///@}

  };

  /// Addition

  ///\relates LpBase::Expr

  ///

  void SkeletonSolverBase::_setSense(Sense) {}

  SkeletonSolverBase::Sense SkeletonSolverBase::_getSense() const

  { return MIN; }

  void SkeletonSolverBase::_clear() {

    row_num = col_num = 0;

  }

  LpSkeleton::SolveExitStatus LpSkeleton::_solve() { return SOLVED; }

  LpSkeleton::Value LpSkeleton::_getPrimal(int) const { return 0; }

  LpSkeleton::Value LpSkeleton::_getDual(int) const { return 0; }

  LpSkeleton::Value LpSkeleton::_getPrimalValue() const { return 0; }

  LpSkeleton::Value LpSkeleton::_getPrimalRay(int) const { return 0; }

  LpSkeleton::Value LpSkeleton::_getDualRay(int) const { return 0; }

    ///\e

    virtual void _setSense(Sense);

    ///\e

    virtual Sense _getSense() const;

    ///\e

    virtual void _clear();

  };

  /// \brief Skeleton class for an LP solver interface

  ///

  ///This class does nothing, but it can serve as a skeleton when

  ///implementing an interface to new solvers.

  ///\ingroup lp_group

 * purpose.

 *

 */

#include <iostream>

#include <lemon/soplex.h>

#include <soplex.h>

///\file

///\brief Implementation of the LEMON-SOPLEX lp solver interface.

namespace lemon {

  SoplexLp::SoplexLp() {

    soplex = new soplex::SoPlex;

  }

  SoplexLp::~SoplexLp() {

    delete soplex;

  }

  SoplexLp::SoplexLp(const SoplexLp& lp) {

    rows = lp.rows;

    (*static_cast<soplex::SPxLP*>(soplex)) = *(lp.soplex);

    _col_names = lp._col_names;

    _col_names_ref = lp._col_names_ref;

    _row_names = lp._row_names;

    _row_names_ref = lp._row_names_ref;

  }

  void SoplexLp::_clear_temporals() {

    _primal_values.clear();

    _dual_values.clear();

  }

  SoplexLp* SoplexLp::newSolver() const {

  SoplexLp::Value SoplexLp::_getObjCoeff(int i) const {

    return soplex->obj(i);

  }

  SoplexLp::SolveExitStatus SoplexLp::_solve() {

    _clear_temporals();

    soplex::SPxSolver::Status status = soplex->solve();

    switch (status) {

    case soplex::SPxSolver::OPTIMAL:

    case soplex::SPxSolver::INFEASIBLE:

    case soplex::SPxSolver::UNBOUNDED:

      return SOLVED;

    _col_names_ref.clear();

    _row_names.clear();

    _row_names_ref.clear();

    cols.clear();

    rows.clear();

    _clear_temporals();

  }