| ... | ... |
@@ -29,49 +29,50 @@ |
| 29 | 29 |
#endif |
| 30 | 30 |
|
| 31 | 31 |
#ifdef LEMON_HAVE_CBC |
| 32 | 32 |
#include <lemon/cbc.h> |
| 33 | 33 |
#endif |
| 34 | 34 |
|
| 35 | 35 |
|
| 36 | 36 |
using namespace lemon; |
| 37 | 37 |
|
| 38 | 38 |
void solveAndCheck(MipSolver& mip, MipSolver::ProblemType stat, |
| 39 | 39 |
double exp_opt) {
|
| 40 | 40 |
using std::string; |
| 41 | 41 |
|
| 42 | 42 |
mip.solve(); |
| 43 | 43 |
//int decimal,sign; |
| 44 | 44 |
std::ostringstream buf; |
| 45 | 45 |
buf << "Type should be: " << int(stat)<<" and it is "<<int(mip.type()); |
| 46 | 46 |
|
| 47 | 47 |
|
| 48 | 48 |
// itoa(stat,buf1, 10); |
| 49 | 49 |
check(mip.type()==stat, buf.str()); |
| 50 | 50 |
|
| 51 | 51 |
if (stat == MipSolver::OPTIMAL) {
|
| 52 | 52 |
std::ostringstream sbuf; |
| 53 |
|
|
| 53 |
sbuf << "Wrong optimal value ("<< mip.solValue()
|
|
| 54 |
<<" instead of " << exp_opt << ")"; |
|
| 54 | 55 |
check(std::abs(mip.solValue()-exp_opt) < 1e-3, sbuf.str()); |
| 55 | 56 |
//+ecvt(exp_opt,2) |
| 56 | 57 |
} |
| 57 | 58 |
} |
| 58 | 59 |
|
| 59 | 60 |
void aTest(MipSolver& mip) |
| 60 | 61 |
{
|
| 61 | 62 |
//The following example is very simple |
| 62 | 63 |
|
| 63 | 64 |
|
| 64 | 65 |
typedef MipSolver::Row Row; |
| 65 | 66 |
typedef MipSolver::Col Col; |
| 66 | 67 |
|
| 67 | 68 |
|
| 68 | 69 |
Col x1 = mip.addCol(); |
| 69 | 70 |
Col x2 = mip.addCol(); |
| 70 | 71 |
|
| 71 | 72 |
|
| 72 | 73 |
//Objective function |
| 73 | 74 |
mip.obj(x1); |
| 74 | 75 |
|
| 75 | 76 |
mip.max(); |
| 76 | 77 |
|
| 77 | 78 |
//Unconstrained optimization |
0 comments (0 inline)