|
alpar@461
|
1 |
/* -*- mode: C++; indent-tabs-mode: nil; -*-
|
|
alpar@461
|
2 |
*
|
|
alpar@461
|
3 |
* This file is a part of LEMON, a generic C++ optimization library.
|
|
alpar@461
|
4 |
*
|
|
alpar@1092
|
5 |
* Copyright (C) 2003-2013
|
|
alpar@461
|
6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
|
|
alpar@461
|
7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES).
|
|
alpar@461
|
8 |
*
|
|
alpar@461
|
9 |
* Permission to use, modify and distribute this software is granted
|
|
alpar@461
|
10 |
* provided that this copyright notice appears in all copies. For
|
|
alpar@461
|
11 |
* precise terms see the accompanying LICENSE file.
|
|
alpar@461
|
12 |
*
|
|
alpar@461
|
13 |
* This software is provided "AS IS" with no warranty of any kind,
|
|
alpar@461
|
14 |
* express or implied, and with no claim as to its suitability for any
|
|
alpar@461
|
15 |
* purpose.
|
|
alpar@461
|
16 |
*
|
|
alpar@461
|
17 |
*/
|
|
alpar@461
|
18 |
|
|
alpar@461
|
19 |
#ifndef LEMON_CLP_H
|
|
alpar@461
|
20 |
#define LEMON_CLP_H
|
|
alpar@461
|
21 |
|
|
alpar@461
|
22 |
///\file
|
|
alpar@461
|
23 |
///\brief Header of the LEMON-CLP lp solver interface.
|
|
alpar@461
|
24 |
|
|
alpar@461
|
25 |
#include <vector>
|
|
alpar@461
|
26 |
#include <string>
|
|
alpar@461
|
27 |
|
|
alpar@461
|
28 |
#include <lemon/lp_base.h>
|
|
alpar@461
|
29 |
|
|
alpar@461
|
30 |
class ClpSimplex;
|
|
alpar@461
|
31 |
|
|
alpar@461
|
32 |
namespace lemon {
|
|
alpar@461
|
33 |
|
|
alpar@461
|
34 |
/// \ingroup lp_group
|
|
alpar@461
|
35 |
///
|
|
alpar@461
|
36 |
/// \brief Interface for the CLP solver
|
|
alpar@461
|
37 |
///
|
|
alpar@461
|
38 |
/// This class implements an interface for the Clp LP solver. The
|
|
alpar@461
|
39 |
/// Clp library is an object oriented lp solver library developed at
|
|
alpar@461
|
40 |
/// the IBM. The CLP is part of the COIN-OR package and it can be
|
|
alpar@461
|
41 |
/// used with Common Public License.
|
|
alpar@462
|
42 |
class ClpLp : public LpSolver {
|
|
alpar@461
|
43 |
protected:
|
|
alpar@461
|
44 |
|
|
alpar@461
|
45 |
ClpSimplex* _prob;
|
|
alpar@461
|
46 |
|
|
alpar@461
|
47 |
std::map<std::string, int> _col_names_ref;
|
|
alpar@461
|
48 |
std::map<std::string, int> _row_names_ref;
|
|
alpar@461
|
49 |
|
|
alpar@461
|
50 |
public:
|
|
alpar@461
|
51 |
|
|
alpar@461
|
52 |
/// \e
|
|
alpar@462
|
53 |
ClpLp();
|
|
alpar@461
|
54 |
/// \e
|
|
alpar@462
|
55 |
ClpLp(const ClpLp&);
|
|
alpar@461
|
56 |
/// \e
|
|
alpar@462
|
57 |
~ClpLp();
|
|
alpar@461
|
58 |
|
|
alpar@540
|
59 |
/// \e
|
|
alpar@540
|
60 |
virtual ClpLp* newSolver() const;
|
|
alpar@540
|
61 |
/// \e
|
|
alpar@540
|
62 |
virtual ClpLp* cloneSolver() const;
|
|
alpar@540
|
63 |
|
|
alpar@461
|
64 |
protected:
|
|
alpar@461
|
65 |
|
|
alpar@461
|
66 |
mutable double* _primal_ray;
|
|
alpar@461
|
67 |
mutable double* _dual_ray;
|
|
alpar@461
|
68 |
|
|
alpar@461
|
69 |
void _init_temporals();
|
|
alpar@461
|
70 |
void _clear_temporals();
|
|
alpar@461
|
71 |
|
|
alpar@461
|
72 |
protected:
|
|
alpar@461
|
73 |
|
|
alpar@461
|
74 |
virtual const char* _solverName() const;
|
|
alpar@461
|
75 |
|
|
alpar@461
|
76 |
virtual int _addCol();
|
|
alpar@461
|
77 |
virtual int _addRow();
|
|
deba@746
|
78 |
virtual int _addRow(Value l, ExprIterator b, ExprIterator e, Value u);
|
|
alpar@461
|
79 |
|
|
alpar@461
|
80 |
virtual void _eraseCol(int i);
|
|
alpar@461
|
81 |
virtual void _eraseRow(int i);
|
|
alpar@461
|
82 |
|
|
alpar@461
|
83 |
virtual void _eraseColId(int i);
|
|
alpar@461
|
84 |
virtual void _eraseRowId(int i);
|
|
alpar@461
|
85 |
|
|
alpar@461
|
86 |
virtual void _getColName(int col, std::string& name) const;
|
|
alpar@461
|
87 |
virtual void _setColName(int col, const std::string& name);
|
|
alpar@461
|
88 |
virtual int _colByName(const std::string& name) const;
|
|
alpar@461
|
89 |
|
|
alpar@461
|
90 |
virtual void _getRowName(int row, std::string& name) const;
|
|
alpar@461
|
91 |
virtual void _setRowName(int row, const std::string& name);
|
|
alpar@461
|
92 |
virtual int _rowByName(const std::string& name) const;
|
|
alpar@461
|
93 |
|
|
alpar@461
|
94 |
virtual void _setRowCoeffs(int i, ExprIterator b, ExprIterator e);
|
|
alpar@461
|
95 |
virtual void _getRowCoeffs(int i, InsertIterator b) const;
|
|
alpar@461
|
96 |
|
|
alpar@461
|
97 |
virtual void _setColCoeffs(int i, ExprIterator b, ExprIterator e);
|
|
alpar@461
|
98 |
virtual void _getColCoeffs(int i, InsertIterator b) const;
|
|
alpar@461
|
99 |
|
|
alpar@461
|
100 |
virtual void _setCoeff(int row, int col, Value value);
|
|
alpar@461
|
101 |
virtual Value _getCoeff(int row, int col) const;
|
|
alpar@461
|
102 |
|
|
alpar@461
|
103 |
virtual void _setColLowerBound(int i, Value value);
|
|
alpar@461
|
104 |
virtual Value _getColLowerBound(int i) const;
|
|
alpar@461
|
105 |
virtual void _setColUpperBound(int i, Value value);
|
|
alpar@461
|
106 |
virtual Value _getColUpperBound(int i) const;
|
|
alpar@461
|
107 |
|
|
alpar@461
|
108 |
virtual void _setRowLowerBound(int i, Value value);
|
|
alpar@461
|
109 |
virtual Value _getRowLowerBound(int i) const;
|
|
alpar@461
|
110 |
virtual void _setRowUpperBound(int i, Value value);
|
|
alpar@461
|
111 |
virtual Value _getRowUpperBound(int i) const;
|
|
alpar@461
|
112 |
|
|
alpar@461
|
113 |
virtual void _setObjCoeffs(ExprIterator, ExprIterator);
|
|
alpar@461
|
114 |
virtual void _getObjCoeffs(InsertIterator) const;
|
|
alpar@461
|
115 |
|
|
alpar@461
|
116 |
virtual void _setObjCoeff(int i, Value obj_coef);
|
|
alpar@461
|
117 |
virtual Value _getObjCoeff(int i) const;
|
|
alpar@461
|
118 |
|
|
alpar@461
|
119 |
virtual void _setSense(Sense sense);
|
|
alpar@461
|
120 |
virtual Sense _getSense() const;
|
|
alpar@461
|
121 |
|
|
alpar@461
|
122 |
virtual SolveExitStatus _solve();
|
|
alpar@461
|
123 |
|
|
alpar@461
|
124 |
virtual Value _getPrimal(int i) const;
|
|
alpar@461
|
125 |
virtual Value _getDual(int i) const;
|
|
alpar@461
|
126 |
|
|
alpar@461
|
127 |
virtual Value _getPrimalValue() const;
|
|
alpar@461
|
128 |
|
|
alpar@461
|
129 |
virtual Value _getPrimalRay(int i) const;
|
|
alpar@461
|
130 |
virtual Value _getDualRay(int i) const;
|
|
alpar@461
|
131 |
|
|
alpar@461
|
132 |
virtual VarStatus _getColStatus(int i) const;
|
|
alpar@461
|
133 |
virtual VarStatus _getRowStatus(int i) const;
|
|
alpar@461
|
134 |
|
|
alpar@461
|
135 |
virtual ProblemType _getPrimalType() const;
|
|
alpar@461
|
136 |
virtual ProblemType _getDualType() const;
|
|
alpar@461
|
137 |
|
|
alpar@461
|
138 |
virtual void _clear();
|
|
alpar@461
|
139 |
|
|
deba@576
|
140 |
virtual void _messageLevel(MessageLevel);
|
|
alpar@877
|
141 |
|
|
alpar@461
|
142 |
public:
|
|
alpar@461
|
143 |
|
|
alpar@461
|
144 |
///Solves LP with primal simplex method.
|
|
alpar@461
|
145 |
SolveExitStatus solvePrimal();
|
|
alpar@461
|
146 |
|
|
alpar@461
|
147 |
///Solves LP with dual simplex method.
|
|
alpar@461
|
148 |
SolveExitStatus solveDual();
|
|
alpar@461
|
149 |
|
|
alpar@461
|
150 |
///Solves LP with barrier method.
|
|
alpar@461
|
151 |
SolveExitStatus solveBarrier();
|
|
alpar@461
|
152 |
|
|
alpar@461
|
153 |
///Returns the constraint identifier understood by CLP.
|
|
ggab90@1130
|
154 |
int clpRow(Row r) const { return _rows(id(r)); }
|
|
alpar@461
|
155 |
|
|
alpar@461
|
156 |
///Returns the variable identifier understood by CLP.
|
|
ggab90@1130
|
157 |
int clpCol(Col c) const { return _cols(id(c)); }
|
|
alpar@461
|
158 |
|
|
alpar@461
|
159 |
};
|
|
alpar@461
|
160 |
|
|
alpar@461
|
161 |
} //END OF NAMESPACE LEMON
|
|
alpar@461
|
162 |
|
|
alpar@461
|
163 |
#endif //LEMON_CLP_H
|
|
alpar@461
|
164 |
|