r565:7ab97e2a0c33
Wed, 01 Apr 2009 22:54:00
[Not Reviewed]
0 comments (0 inline)
0
2
0
50
2
7
1
| ... | ... |
@@ -530,17 +530,19 @@ |
| 530 | 530 |
|
| 531 | 531 |
// GlpkLp members |
| 532 | 532 |
|
| 533 | 533 |
GlpkLp::GlpkLp() |
| 534 | 534 |
: LpBase(), LpSolver(), GlpkBase() {
|
| 535 | 535 |
messageLevel(MESSAGE_NO_OUTPUT); |
| 536 |
presolver(false); |
|
| 536 | 537 |
} |
| 537 | 538 |
|
| 538 | 539 |
GlpkLp::GlpkLp(const GlpkLp& other) |
| 539 | 540 |
: LpBase(other), LpSolver(other), GlpkBase(other) {
|
| 540 | 541 |
messageLevel(MESSAGE_NO_OUTPUT); |
| 542 |
presolver(false); |
|
| 541 | 543 |
} |
| 542 | 544 |
|
| 543 | 545 |
GlpkLp* GlpkLp::newSolver() const { return new GlpkLp; }
|
| 544 | 546 |
GlpkLp* GlpkLp::cloneSolver() const { return new GlpkLp(*this); }
|
| 545 | 547 |
|
| 546 | 548 |
const char* GlpkLp::_solverName() const { return "GlpkLp"; }
|
| ... | ... |
@@ -571,14 +573,30 @@ |
| 571 | 573 |
smcp.msg_lev = GLP_MSG_ON; |
| 572 | 574 |
break; |
| 573 | 575 |
case MESSAGE_FULL_OUTPUT: |
| 574 | 576 |
smcp.msg_lev = GLP_MSG_ALL; |
| 575 | 577 |
break; |
| 576 | 578 |
} |
| 579 |
smcp.presolve = _presolve; |
|
| 577 | 580 |
|
| 581 |
// If the basis is not valid we get an error return value. |
|
| 582 |
// In this case we can try to create a new basis. |
|
| 583 |
switch (glp_simplex(lp, &smcp)) {
|
|
| 584 |
case 0: |
|
| 585 |
break; |
|
| 586 |
case GLP_EBADB: |
|
| 587 |
case GLP_ESING: |
|
| 588 |
case GLP_ECOND: |
|
| 589 |
lpx_set_int_parm(lp, LPX_K_MSGLEV, smcp.msg_lev); |
|
| 590 |
glp_adv_basis(lp, 0); |
|
| 578 | 591 |
if (glp_simplex(lp, &smcp) != 0) return UNSOLVED; |
| 592 |
break; |
|
| 593 |
default: |
|
| 594 |
return UNSOLVED; |
|
| 595 |
} |
|
| 596 |
|
|
| 579 | 597 |
return SOLVED; |
| 580 | 598 |
} |
| 581 | 599 |
|
| 582 | 600 |
GlpkLp::SolveExitStatus GlpkLp::solveDual() {
|
| 583 | 601 |
_clear_temporals(); |
| 584 | 602 |
|
| ... | ... |
@@ -597,14 +615,29 @@ |
| 597 | 615 |
break; |
| 598 | 616 |
case MESSAGE_FULL_OUTPUT: |
| 599 | 617 |
smcp.msg_lev = GLP_MSG_ALL; |
| 600 | 618 |
break; |
| 601 | 619 |
} |
| 602 | 620 |
smcp.meth = GLP_DUAL; |
| 621 |
smcp.presolve = _presolve; |
|
| 603 | 622 |
|
| 623 |
// If the basis is not valid we get an error return value. |
|
| 624 |
// In this case we can try to create a new basis. |
|
| 625 |
switch (glp_simplex(lp, &smcp)) {
|
|
| 626 |
case 0: |
|
| 627 |
break; |
|
| 628 |
case GLP_EBADB: |
|
| 629 |
case GLP_ESING: |
|
| 630 |
case GLP_ECOND: |
|
| 631 |
lpx_set_int_parm(lp, LPX_K_MSGLEV, smcp.msg_lev); |
|
| 632 |
glp_adv_basis(lp, 0); |
|
| 604 | 633 |
if (glp_simplex(lp, &smcp) != 0) return UNSOLVED; |
| 634 |
break; |
|
| 635 |
default: |
|
| 636 |
return UNSOLVED; |
|
| 637 |
} |
|
| 605 | 638 |
return SOLVED; |
| 606 | 639 |
} |
| 607 | 640 |
|
| 608 | 641 |
GlpkLp::Value GlpkLp::_getPrimal(int i) const {
|
| 609 | 642 |
return glp_get_col_prim(lp, i); |
| 610 | 643 |
} |
| ... | ... |
@@ -816,14 +849,14 @@ |
| 816 | 849 |
default: |
| 817 | 850 |
LEMON_ASSERT(false, "Wrong primal type"); |
| 818 | 851 |
return GlpkLp::ProblemType(); |
| 819 | 852 |
} |
| 820 | 853 |
} |
| 821 | 854 |
|
| 822 |
void GlpkLp::presolver(bool b) {
|
|
| 823 |
lpx_set_int_parm(lp, LPX_K_PRESOL, b ? 1 : 0); |
|
| 855 |
void GlpkLp::presolver(bool presolve) {
|
|
| 856 |
_presolve = presolve; |
|
| 824 | 857 |
} |
| 825 | 858 |
|
| 826 | 859 |
void GlpkLp::messageLevel(MessageLevel m) {
|
| 827 | 860 |
_message_level = m; |
| 828 | 861 |
} |
| 829 | 862 |
|
| ... | ... |
@@ -878,13 +911,28 @@ |
| 878 | 911 |
case MESSAGE_FULL_OUTPUT: |
| 879 | 912 |
smcp.msg_lev = GLP_MSG_ALL; |
| 880 | 913 |
break; |
| 881 | 914 |
} |
| 882 | 915 |
smcp.meth = GLP_DUAL; |
| 883 | 916 |
|
| 917 |
// If the basis is not valid we get an error return value. |
|
| 918 |
// In this case we can try to create a new basis. |
|
| 919 |
switch (glp_simplex(lp, &smcp)) {
|
|
| 920 |
case 0: |
|
| 921 |
break; |
|
| 922 |
case GLP_EBADB: |
|
| 923 |
case GLP_ESING: |
|
| 924 |
case GLP_ECOND: |
|
| 925 |
lpx_set_int_parm(lp, LPX_K_MSGLEV, smcp.msg_lev); |
|
| 926 |
glp_adv_basis(lp, 0); |
|
| 884 | 927 |
if (glp_simplex(lp, &smcp) != 0) return UNSOLVED; |
| 928 |
break; |
|
| 929 |
default: |
|
| 930 |
return UNSOLVED; |
|
| 931 |
} |
|
| 932 |
|
|
| 885 | 933 |
if (glp_get_status(lp) != GLP_OPT) return SOLVED; |
| 886 | 934 |
|
| 887 | 935 |
glp_iocp iocp; |
| 888 | 936 |
glp_init_iocp(&iocp); |
| 889 | 937 |
|
| 890 | 938 |
switch (_message_level) {
|
| ... | ... |
@@ -175,18 +175,24 @@ |
| 175 | 175 |
///Solve with primal simplex |
| 176 | 176 |
SolveExitStatus solvePrimal(); |
| 177 | 177 |
|
| 178 | 178 |
///Solve with dual simplex |
| 179 | 179 |
SolveExitStatus solveDual(); |
| 180 | 180 |
|
| 181 |
private: |
|
| 182 |
|
|
| 183 |
bool _presolve; |
|
| 184 |
|
|
| 185 |
public: |
|
| 186 |
|
|
| 181 | 187 |
///Turns on or off the presolver |
| 182 | 188 |
|
| 183 | 189 |
///Turns on (\c b is \c true) or off (\c b is \c false) the presolver |
| 184 | 190 |
/// |
| 185 | 191 |
///The presolver is off by default. |
| 186 |
void presolver(bool |
|
| 192 |
void presolver(bool presolve); |
|
| 187 | 193 |
|
| 188 | 194 |
///Enum for \c messageLevel() parameter |
| 189 | 195 |
enum MessageLevel {
|
| 190 | 196 |
/// no output (default value) |
| 191 | 197 |
MESSAGE_NO_OUTPUT = 0, |
| 192 | 198 |
/// error messages only |
0 comments (0 inline)