RhodeCode

    messageLevel(MESSAGE_NO_OUTPUT);

    presolver(false);

  }

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

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

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

  void GlpkLp::_clear_temporals() {

    _primal_ray.clear();

    _dual_ray.clear();

  }

  GlpkLp::SolveExitStatus GlpkLp::_solve() {

    return solvePrimal();

  }

  GlpkLp::SolveExitStatus GlpkLp::solvePrimal() {

    _clear_temporals();

    glp_smcp smcp;

    glp_init_smcp(&smcp);

    switch (_message_level) {

    case MESSAGE_NO_OUTPUT:

      smcp.msg_lev = GLP_MSG_OFF;

      break;

    case MESSAGE_ERROR_MESSAGE:

      smcp.msg_lev = GLP_MSG_ERR;

      break;

    case MESSAGE_NORMAL_OUTPUT:

      smcp.msg_lev = GLP_MSG_ON;

      break;

    case MESSAGE_FULL_OUTPUT:

      smcp.msg_lev = GLP_MSG_ALL;

      break;

    }

    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:

    case GLP_ESING:

    case GLP_ECOND:

      glp_adv_basis(lp, 0);

      if (glp_simplex(lp, &smcp) != 0) return UNSOLVED;

      break;

    default:

      return UNSOLVED;

    }

    return SOLVED;

  }

  GlpkLp::SolveExitStatus GlpkLp::solveDual() {

    _clear_temporals();

    glp_smcp smcp;

    glp_init_smcp(&smcp);

    switch (_message_level) {

    case MESSAGE_NO_OUTPUT:

      smcp.msg_lev = GLP_MSG_OFF;

      break;

    case MESSAGE_ERROR_MESSAGE:

      smcp.msg_lev = GLP_MSG_ERR;

      break;

    case MESSAGE_NORMAL_OUTPUT:

      smcp.msg_lev = GLP_MSG_ON;

      break;

    case MESSAGE_FULL_OUTPUT:

      smcp.msg_lev = GLP_MSG_ALL;

      break;

    }

    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;

    case GLP_EBADB:

    case GLP_ESING:

    case GLP_ECOND:

      glp_adv_basis(lp, 0);

      if (glp_simplex(lp, &smcp) != 0) return UNSOLVED;

      break;

    default:

      return UNSOLVED;

    }

    return SOLVED;

  }

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

    return glp_get_col_prim(lp, i);

  }

  GlpkLp::Value GlpkLp::_getDual(int i) const {

    return glp_get_row_dual(lp, i);

  }

  GlpkLp::Value GlpkLp::_getPrimalValue() const {

    return glp_get_obj_val(lp);

  }

  GlpkLp::VarStatus GlpkLp::_getColStatus(int i) const {

    switch (glp_get_col_stat(lp, i)) {

    case GLP_BS:

      return BASIC;

    case GLP_UP:

      return UPPER;

    case GLP_LO:

      return LOWER;

    case GLP_NF:

      return FREE;

    case GLP_NS:

      return FIXED;

    default:

      LEMON_ASSERT(false, "Wrong column status");

      return GlpkLp::VarStatus();

    }

  }

  GlpkLp::VarStatus GlpkLp::_getRowStatus(int i) const {

    switch (glp_get_row_stat(lp, i)) {

    case GLP_BS:

      return BASIC;

    case GLP_UP:

      return UPPER;

    case GLP_LO:

      return LOWER;

    case GLP_NF:

      return FREE;

    case INTEGER:

      glp_set_col_kind(lp, i, GLP_IV);

      break;

    case REAL:

      glp_set_col_kind(lp, i, GLP_CV);

      break;

    }

  }

  GlpkMip::ColTypes GlpkMip::_getColType(int i) const {

    switch (glp_get_col_kind(lp, i)) {

    case GLP_IV:

    case GLP_BV:

      return INTEGER;

    default:

      return REAL;

    }

  }

  GlpkMip::SolveExitStatus GlpkMip::_solve() {

    glp_smcp smcp;

    glp_init_smcp(&smcp);

    switch (_message_level) {

    case MESSAGE_NO_OUTPUT:

      smcp.msg_lev = GLP_MSG_OFF;

      break;

    case MESSAGE_ERROR_MESSAGE:

      smcp.msg_lev = GLP_MSG_ERR;

      break;

    case MESSAGE_NORMAL_OUTPUT:

      smcp.msg_lev = GLP_MSG_ON;

      break;

    case MESSAGE_FULL_OUTPUT:

      smcp.msg_lev = GLP_MSG_ALL;

      break;

    }

    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:

    case GLP_ESING:

    case GLP_ECOND:

      glp_adv_basis(lp, 0);

      if (glp_simplex(lp, &smcp) != 0) return UNSOLVED;

      break;

    default:

      return UNSOLVED;

    }

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

    glp_iocp iocp;

    glp_init_iocp(&iocp);

    switch (_message_level) {

    case MESSAGE_NO_OUTPUT:

      iocp.msg_lev = GLP_MSG_OFF;

      break;

    case MESSAGE_ERROR_MESSAGE:

      iocp.msg_lev = GLP_MSG_ERR;

      break;

    case MESSAGE_NORMAL_OUTPUT:

      iocp.msg_lev = GLP_MSG_ON;

      break;

    case MESSAGE_FULL_OUTPUT:

      iocp.msg_lev = GLP_MSG_ALL;

      break;

    }

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

    return SOLVED;

  }

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

    switch (glp_get_status(lp)) {

    case GLP_OPT:

      switch (glp_mip_status(lp)) {

      case GLP_UNDEF:

        return UNDEFINED;

      case GLP_NOFEAS:

        return INFEASIBLE;

      case GLP_FEAS:

        return FEASIBLE;

      case GLP_OPT:

        return OPTIMAL;

      default:

        LEMON_ASSERT(false, "Wrong problem type.");

        return GlpkMip::ProblemType();

      }

    case GLP_NOFEAS: