COIN-OR::LEMON - Graph Library

source: lemon/test/mip_test.cc @ 1352:f85ee41c84bc

Last change on this file since 1352:f85ee41c84bc was 1300:62dba6c90f35, checked in by Alpar Juttner <alpar@…>, 11 years ago

Fix default LP/MIP setting (#479)

File size: 3.2 KB
RevLine 
[481]1/* -*- mode: C++; indent-tabs-mode: nil; -*-
2 *
3 * This file is a part of LEMON, a generic C++ optimization library.
4 *
[598]5 * Copyright (C) 2003-2009
[481]6 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
7 * (Egervary Research Group on Combinatorial Optimization, EGRES).
8 *
9 * Permission to use, modify and distribute this software is granted
10 * provided that this copyright notice appears in all copies. For
11 * precise terms see the accompanying LICENSE file.
12 *
13 * This software is provided "AS IS" with no warranty of any kind,
14 * express or implied, and with no claim as to its suitability for any
15 * purpose.
16 *
17 */
18
19#include "test_tools.h"
20
21#include <lemon/config.h>
22
[674]23#ifdef LEMON_HAVE_CPLEX
[484]24#include <lemon/cplex.h>
[481]25#endif
26
[674]27#ifdef LEMON_HAVE_GLPK
[484]28#include <lemon/glpk.h>
[481]29#endif
30
[674]31#ifdef LEMON_HAVE_CBC
[614]32#include <lemon/cbc.h>
33#endif
34
[1300]35#ifdef LEMON_HAVE_MIP
36#include <lemon/lp.h>
37#endif
38
[481]39
40using namespace lemon;
41
[482]42void solveAndCheck(MipSolver& mip, MipSolver::ProblemType stat,
[481]43                   double exp_opt) {
44  using std::string;
45
[482]46  mip.solve();
[481]47  //int decimal,sign;
48  std::ostringstream buf;
[482]49  buf << "Type should be: " << int(stat)<<" and it is "<<int(mip.type());
[481]50
51
52  //  itoa(stat,buf1, 10);
[482]53  check(mip.type()==stat, buf.str());
[481]54
[482]55  if (stat ==  MipSolver::OPTIMAL) {
[481]56    std::ostringstream sbuf;
[795]57    sbuf << "Wrong optimal value ("<< mip.solValue()
58         <<" instead of " << exp_opt << ")";
[482]59    check(std::abs(mip.solValue()-exp_opt) < 1e-3, sbuf.str());
[481]60    //+ecvt(exp_opt,2)
61  }
62}
63
[482]64void aTest(MipSolver& mip)
[481]65{
[614]66  //The following example is very simple
[481]67
68
[482]69  typedef MipSolver::Row Row;
70  typedef MipSolver::Col Col;
[481]71
72
73  Col x1 = mip.addCol();
74  Col x2 = mip.addCol();
75
76
77  //Objective function
78  mip.obj(x1);
79
80  mip.max();
81
82  //Unconstrained optimization
83  mip.solve();
84  //Check it out!
85
86  //Constraints
[614]87  mip.addRow(2 * x1 + x2 <= 2);
88  Row y2 = mip.addRow(x1 - 2 * x2 <= 0);
[481]89
90  //Nonnegativity of the variable x1
91  mip.colLowerBound(x1, 0);
92
[614]93
[481]94  //Maximization of x1
95  //over the triangle with vertices (0,0),(4/5,2/5),(0,2)
96  double expected_opt=4.0/5.0;
[482]97  solveAndCheck(mip, MipSolver::OPTIMAL, expected_opt);
[481]98
[614]99
[481]100  //Restrict x2 to integer
[482]101  mip.colType(x2,MipSolver::INTEGER);
[481]102  expected_opt=1.0/2.0;
[482]103  solveAndCheck(mip, MipSolver::OPTIMAL, expected_opt);
[481]104
105
106  //Restrict both to integer
[482]107  mip.colType(x1,MipSolver::INTEGER);
[481]108  expected_opt=0;
[482]109  solveAndCheck(mip, MipSolver::OPTIMAL, expected_opt);
[481]110
[614]111  //Erase a variable
112  mip.erase(x2);
113  mip.rowUpperBound(y2, 8);
114  expected_opt=1;
115  solveAndCheck(mip, MipSolver::OPTIMAL, expected_opt);
[481]116
117}
118
[614]119
[587]120template<class MIP>
121void cloneTest()
122{
[598]123
[587]124  MIP* mip = new MIP();
125  MIP* mipnew = mip->newSolver();
126  MIP* mipclone = mip->cloneSolver();
127  delete mip;
128  delete mipnew;
129  delete mipclone;
130}
[481]131
132int main()
133{
134
[1300]135#ifdef LEMON_HAVE_MIP
136  {
137    Mip mip1;
138    aTest(mip1);
139    cloneTest<Mip>();
140  }
141#endif
142
[674]143#ifdef LEMON_HAVE_GLPK
[482]144  {
[485]145    GlpkMip mip1;
[482]146    aTest(mip1);
[587]147    cloneTest<GlpkMip>();
[482]148  }
[481]149#endif
150
[674]151#ifdef LEMON_HAVE_CPLEX
[482]152  try {
[485]153    CplexMip mip2;
[482]154    aTest(mip2);
[598]155    cloneTest<CplexMip>();
[482]156  } catch (CplexEnv::LicenseError& error) {
157    check(false, error.what());
158  }
[481]159#endif
160
[674]161#ifdef LEMON_HAVE_CBC
[614]162  {
163    CbcMip mip1;
164    aTest(mip1);
165    cloneTest<CbcMip>();
166  }
167#endif
168
[481]169  return 0;
170
171}
Note: See TracBrowser for help on using the repository browser.