0
4
0
| ... | ... |
@@ -501,48 +501,52 @@ |
| 501 | 501 |
case MAX: |
| 502 | 502 |
glp_set_obj_dir(lp, GLP_MAX); |
| 503 | 503 |
break; |
| 504 | 504 |
} |
| 505 | 505 |
} |
| 506 | 506 |
|
| 507 | 507 |
GlpkBase::Sense GlpkBase::_getSense() const {
|
| 508 | 508 |
switch(glp_get_obj_dir(lp)) {
|
| 509 | 509 |
case GLP_MIN: |
| 510 | 510 |
return MIN; |
| 511 | 511 |
case GLP_MAX: |
| 512 | 512 |
return MAX; |
| 513 | 513 |
default: |
| 514 | 514 |
LEMON_ASSERT(false, "Wrong sense"); |
| 515 | 515 |
return GlpkBase::Sense(); |
| 516 | 516 |
} |
| 517 | 517 |
} |
| 518 | 518 |
|
| 519 | 519 |
void GlpkBase::_clear() {
|
| 520 | 520 |
glp_erase_prob(lp); |
| 521 | 521 |
rows.clear(); |
| 522 | 522 |
cols.clear(); |
| 523 | 523 |
} |
| 524 | 524 |
|
| 525 |
void GlpkBase::freeEnv() {
|
|
| 526 |
glp_free_env(); |
|
| 527 |
} |
|
| 528 |
|
|
| 525 | 529 |
// GlpkLp members |
| 526 | 530 |
|
| 527 | 531 |
GlpkLp::GlpkLp() |
| 528 | 532 |
: LpBase(), GlpkBase(), LpSolver() {
|
| 529 | 533 |
messageLevel(MESSAGE_NO_OUTPUT); |
| 530 | 534 |
} |
| 531 | 535 |
|
| 532 | 536 |
GlpkLp::GlpkLp(const GlpkLp& other) |
| 533 | 537 |
: LpBase(other), GlpkBase(other), LpSolver(other) {
|
| 534 | 538 |
messageLevel(MESSAGE_NO_OUTPUT); |
| 535 | 539 |
} |
| 536 | 540 |
|
| 537 | 541 |
GlpkLp* GlpkLp::_newSolver() const { return new GlpkLp; }
|
| 538 | 542 |
GlpkLp* GlpkLp::_cloneSolver() const { return new GlpkLp(*this); }
|
| 539 | 543 |
|
| 540 | 544 |
const char* GlpkLp::_solverName() const { return "GlpkLp"; }
|
| 541 | 545 |
|
| 542 | 546 |
void GlpkLp::_clear_temporals() {
|
| 543 | 547 |
_primal_ray.clear(); |
| 544 | 548 |
_dual_ray.clear(); |
| 545 | 549 |
} |
| 546 | 550 |
|
| 547 | 551 |
GlpkLp::SolveExitStatus GlpkLp::_solve() {
|
| 548 | 552 |
return solvePrimal(); |
| ... | ... |
@@ -81,48 +81,56 @@ |
| 81 | 81 |
virtual Value _getColLowerBound(int i) const; |
| 82 | 82 |
|
| 83 | 83 |
virtual void _setColUpperBound(int i, Value value); |
| 84 | 84 |
virtual Value _getColUpperBound(int i) const; |
| 85 | 85 |
|
| 86 | 86 |
virtual void _setRowLowerBound(int i, Value value); |
| 87 | 87 |
virtual Value _getRowLowerBound(int i) const; |
| 88 | 88 |
|
| 89 | 89 |
virtual void _setRowUpperBound(int i, Value value); |
| 90 | 90 |
virtual Value _getRowUpperBound(int i) const; |
| 91 | 91 |
|
| 92 | 92 |
virtual void _setObjCoeffs(ExprIterator b, ExprIterator e); |
| 93 | 93 |
virtual void _getObjCoeffs(InsertIterator b) const; |
| 94 | 94 |
|
| 95 | 95 |
virtual void _setObjCoeff(int i, Value obj_coef); |
| 96 | 96 |
virtual Value _getObjCoeff(int i) const; |
| 97 | 97 |
|
| 98 | 98 |
virtual void _setSense(Sense); |
| 99 | 99 |
virtual Sense _getSense() const; |
| 100 | 100 |
|
| 101 | 101 |
virtual void _clear(); |
| 102 | 102 |
|
| 103 | 103 |
public: |
| 104 | 104 |
|
| 105 |
/// \brief Deallocates the globally allocated memory of GLPK. |
|
| 106 |
|
|
| 107 |
/// Deallocates the globally allocated memory of GLPK. \note |
|
| 108 |
/// Usually, it do not have to be called, because the GLPK use |
|
| 109 |
/// only a small amount of global memory, and it is deallocated |
|
| 110 |
/// automatically at the end of program. |
|
| 111 |
static void freeEnv(); |
|
| 112 |
|
|
| 105 | 113 |
///Pointer to the underlying GLPK data structure. |
| 106 | 114 |
LPX *lpx() {return lp;}
|
| 107 | 115 |
///Const pointer to the underlying GLPK data structure. |
| 108 | 116 |
const LPX *lpx() const {return lp;}
|
| 109 | 117 |
|
| 110 | 118 |
///Returns the constraint identifier understood by GLPK. |
| 111 | 119 |
int lpxRow(Row r) const { return rows(id(r)); }
|
| 112 | 120 |
|
| 113 | 121 |
///Returns the variable identifier understood by GLPK. |
| 114 | 122 |
int lpxCol(Col c) const { return cols(id(c)); }
|
| 115 | 123 |
|
| 116 | 124 |
}; |
| 117 | 125 |
|
| 118 | 126 |
/// \brief Interface for the GLPK LP solver |
| 119 | 127 |
/// |
| 120 | 128 |
/// This class implements an interface for the GLPK LP solver. |
| 121 | 129 |
///\ingroup lp_group |
| 122 | 130 |
class GlpkLp : public GlpkBase, public LpSolver {
|
| 123 | 131 |
public: |
| 124 | 132 |
|
| 125 | 133 |
///\e |
| 126 | 134 |
GlpkLp(); |
| 127 | 135 |
///\e |
| 128 | 136 |
GlpkLp(const GlpkLp&); |
| ... | ... |
@@ -345,48 +345,49 @@ |
| 345 | 345 |
|
| 346 | 346 |
//Erase one constraint and return to maximization |
| 347 | 347 |
lp.erase(upright); |
| 348 | 348 |
lp.sense(lp.MAX); |
| 349 | 349 |
expected_opt=LpSolver::INF; |
| 350 | 350 |
solveAndCheck(lp, LpSolver::UNBOUNDED, expected_opt); |
| 351 | 351 |
|
| 352 | 352 |
//Infeasibilty |
| 353 | 353 |
lp.addRow(x1+x2 <=-2); |
| 354 | 354 |
solveAndCheck(lp, LpSolver::INFEASIBLE, expected_opt); |
| 355 | 355 |
|
| 356 | 356 |
} |
| 357 | 357 |
|
| 358 | 358 |
int main() |
| 359 | 359 |
{
|
| 360 | 360 |
LpSkeleton lp_skel; |
| 361 | 361 |
lpTest(lp_skel); |
| 362 | 362 |
|
| 363 | 363 |
#ifdef HAVE_GLPK |
| 364 | 364 |
{
|
| 365 | 365 |
GlpkLp lp_glpk1,lp_glpk2; |
| 366 | 366 |
lpTest(lp_glpk1); |
| 367 | 367 |
aTest(lp_glpk2); |
| 368 | 368 |
} |
| 369 |
GlpkLp::freeEnv(); |
|
| 369 | 370 |
#endif |
| 370 | 371 |
|
| 371 | 372 |
#ifdef HAVE_CPLEX |
| 372 | 373 |
try {
|
| 373 | 374 |
CplexLp lp_cplex1,lp_cplex2; |
| 374 | 375 |
lpTest(lp_cplex1); |
| 375 | 376 |
aTest(lp_cplex2); |
| 376 | 377 |
} catch (CplexEnv::LicenseError& error) {
|
| 377 | 378 |
#ifdef LEMON_FORCE_CPLEX_CHECK |
| 378 | 379 |
check(false, error.what()); |
| 379 | 380 |
#else |
| 380 | 381 |
std::cerr << error.what() << std::endl; |
| 381 | 382 |
std::cerr << "Cplex license check failed, lp check skipped" << std::endl; |
| 382 | 383 |
#endif |
| 383 | 384 |
} |
| 384 | 385 |
#endif |
| 385 | 386 |
|
| 386 | 387 |
#ifdef HAVE_SOPLEX |
| 387 | 388 |
{
|
| 388 | 389 |
SoplexLp lp_soplex1,lp_soplex2; |
| 389 | 390 |
lpTest(lp_soplex1); |
| 390 | 391 |
aTest(lp_soplex2); |
| 391 | 392 |
} |
| 392 | 393 |
#endif |
| ... | ... |
@@ -94,43 +94,44 @@ |
| 94 | 94 |
//Restrict x2 to integer |
| 95 | 95 |
mip.colType(x2,MipSolver::INTEGER); |
| 96 | 96 |
expected_opt=1.0/2.0; |
| 97 | 97 |
solveAndCheck(mip, MipSolver::OPTIMAL, expected_opt); |
| 98 | 98 |
|
| 99 | 99 |
|
| 100 | 100 |
//Restrict both to integer |
| 101 | 101 |
mip.colType(x1,MipSolver::INTEGER); |
| 102 | 102 |
expected_opt=0; |
| 103 | 103 |
solveAndCheck(mip, MipSolver::OPTIMAL, expected_opt); |
| 104 | 104 |
|
| 105 | 105 |
|
| 106 | 106 |
|
| 107 | 107 |
} |
| 108 | 108 |
|
| 109 | 109 |
|
| 110 | 110 |
int main() |
| 111 | 111 |
{
|
| 112 | 112 |
|
| 113 | 113 |
#ifdef HAVE_GLPK |
| 114 | 114 |
{
|
| 115 | 115 |
GlpkMip mip1; |
| 116 | 116 |
aTest(mip1); |
| 117 | 117 |
} |
| 118 |
GlpkLp::freeEnv(); |
|
| 118 | 119 |
#endif |
| 119 | 120 |
|
| 120 | 121 |
#ifdef HAVE_CPLEX |
| 121 | 122 |
try {
|
| 122 | 123 |
CplexMip mip2; |
| 123 | 124 |
aTest(mip2); |
| 124 | 125 |
} catch (CplexEnv::LicenseError& error) {
|
| 125 | 126 |
#ifdef LEMON_FORCE_CPLEX_CHECK |
| 126 | 127 |
check(false, error.what()); |
| 127 | 128 |
#else |
| 128 | 129 |
std::cerr << error.what() << std::endl; |
| 129 | 130 |
std::cerr << "Cplex license check failed, lp check skipped" << std::endl; |
| 130 | 131 |
#endif |
| 131 | 132 |
} |
| 132 | 133 |
#endif |
| 133 | 134 |
|
| 134 | 135 |
return 0; |
| 135 | 136 |
|
| 136 | 137 |
} |
0 comments (0 inline)