0
4
1
1 |
/* -*- mode: C++; indent-tabs-mode: nil; -*- |
|
2 |
* |
|
3 |
* This file is a part of LEMON, a generic C++ optimization library. |
|
4 |
* |
|
5 |
* Copyright (C) 2003-2009 |
|
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 <iostream> |
|
20 |
#include <sstream> |
|
21 |
|
|
22 |
#include <lemon/smart_graph.h> |
|
23 |
#include <lemon/lgf_reader.h> |
|
24 |
#include <lemon/min_mean_cycle.h> |
|
25 |
#include <lemon/path.h> |
|
26 |
#include <lemon/concepts/digraph.h> |
|
27 |
#include <lemon/concept_check.h> |
|
28 |
|
|
29 |
#include "test_tools.h" |
|
30 |
|
|
31 |
using namespace lemon; |
|
32 |
|
|
33 |
char test_lgf[] = |
|
34 |
"@nodes\n" |
|
35 |
"label\n" |
|
36 |
"1\n" |
|
37 |
"2\n" |
|
38 |
"3\n" |
|
39 |
"4\n" |
|
40 |
"5\n" |
|
41 |
"6\n" |
|
42 |
"7\n" |
|
43 |
"@arcs\n" |
|
44 |
" len1 len2 len3 len4 c1 c2 c3 c4\n" |
|
45 |
"1 2 1 1 1 1 0 0 0 0\n" |
|
46 |
"2 4 5 5 5 5 1 0 0 0\n" |
|
47 |
"2 3 8 8 8 8 0 0 0 0\n" |
|
48 |
"3 2 -2 0 0 0 1 0 0 0\n" |
|
49 |
"3 4 4 4 4 4 0 0 0 0\n" |
|
50 |
"3 7 -4 -4 -4 -4 0 0 0 0\n" |
|
51 |
"4 1 2 2 2 2 0 0 0 0\n" |
|
52 |
"4 3 3 3 3 3 1 0 0 0\n" |
|
53 |
"4 4 3 3 0 0 0 0 1 0\n" |
|
54 |
"5 2 4 4 4 4 0 0 0 0\n" |
|
55 |
"5 6 3 3 3 3 0 1 0 0\n" |
|
56 |
"6 5 2 2 2 2 0 1 0 0\n" |
|
57 |
"6 4 -1 -1 -1 -1 0 0 0 0\n" |
|
58 |
"6 7 1 1 1 1 0 0 0 0\n" |
|
59 |
"7 7 4 4 4 -1 0 0 0 1\n"; |
|
60 |
|
|
61 |
|
|
62 |
// Check the interface of an MMC algorithm |
|
63 |
template <typename GR, typename Value> |
|
64 |
struct MmcClassConcept |
|
65 |
{ |
|
66 |
template <typename MMC> |
|
67 |
struct Constraints { |
|
68 |
void constraints() { |
|
69 |
const Constraints& me = *this; |
|
70 |
|
|
71 |
typedef typename MMC |
|
72 |
::template SetPath<ListPath<GR> > |
|
73 |
::template SetLargeValue<Value> |
|
74 |
::Create MmcAlg; |
|
75 |
MmcAlg mmc(me.g, me.length); |
|
76 |
const MmcAlg& const_mmc = mmc; |
|
77 |
|
|
78 |
b = mmc.cycle(p).run(); |
|
79 |
b = mmc.findMinMean(); |
|
80 |
b = mmc.findCycle(); |
|
81 |
|
|
82 |
v = const_mmc.cycleLength(); |
|
83 |
i = const_mmc.cycleArcNum(); |
|
84 |
d = const_mmc.cycleMean(); |
|
85 |
p = const_mmc.cycle(); |
|
86 |
} |
|
87 |
|
|
88 |
typedef concepts::ReadMap<typename GR::Arc, Value> LM; |
|
89 |
|
|
90 |
GR g; |
|
91 |
LM length; |
|
92 |
ListPath<GR> p; |
|
93 |
Value v; |
|
94 |
int i; |
|
95 |
double d; |
|
96 |
bool b; |
|
97 |
}; |
|
98 |
}; |
|
99 |
|
|
100 |
// Perform a test with the given parameters |
|
101 |
template <typename MMC> |
|
102 |
void checkMmcAlg(const SmartDigraph& gr, |
|
103 |
const SmartDigraph::ArcMap<int>& lm, |
|
104 |
const SmartDigraph::ArcMap<int>& cm, |
|
105 |
int length, int size) { |
|
106 |
MMC alg(gr, lm); |
|
107 |
alg.findMinMean(); |
|
108 |
check(alg.cycleMean() == static_cast<double>(length) / size, |
|
109 |
"Wrong cycle mean"); |
|
110 |
alg.findCycle(); |
|
111 |
check(alg.cycleLength() == length && alg.cycleArcNum() == size, |
|
112 |
"Wrong path"); |
|
113 |
SmartDigraph::ArcMap<int> cycle(gr, 0); |
|
114 |
for (typename MMC::Path::ArcIt a(alg.cycle()); a != INVALID; ++a) { |
|
115 |
++cycle[a]; |
|
116 |
} |
|
117 |
for (SmartDigraph::ArcIt a(gr); a != INVALID; ++a) { |
|
118 |
check(cm[a] == cycle[a], "Wrong path"); |
|
119 |
} |
|
120 |
} |
|
121 |
|
|
122 |
// Class for comparing types |
|
123 |
template <typename T1, typename T2> |
|
124 |
struct IsSameType { |
|
125 |
static const int result = 0; |
|
126 |
}; |
|
127 |
|
|
128 |
template <typename T> |
|
129 |
struct IsSameType<T,T> { |
|
130 |
static const int result = 1; |
|
131 |
}; |
|
132 |
|
|
133 |
|
|
134 |
int main() { |
|
135 |
#ifdef LEMON_HAVE_LONG_LONG |
|
136 |
typedef long long long_int; |
|
137 |
#else |
|
138 |
typedef long long_int; |
|
139 |
#endif |
|
140 |
|
|
141 |
// Check the interface |
|
142 |
{ |
|
143 |
typedef concepts::Digraph GR; |
|
144 |
typedef MinMeanCycle<GR, concepts::ReadMap<GR::Arc, int> > IntMmcAlg; |
|
145 |
typedef MinMeanCycle<GR, concepts::ReadMap<GR::Arc, float> > FloatMmcAlg; |
|
146 |
|
|
147 |
checkConcept<MmcClassConcept<GR, int>, IntMmcAlg>(); |
|
148 |
checkConcept<MmcClassConcept<GR, float>, FloatMmcAlg>(); |
|
149 |
|
|
150 |
if (IsSameType<IntMmcAlg::LargeValue, long_int>::result == 0) |
|
151 |
check(false, "Wrong LargeValue type"); |
|
152 |
if (IsSameType<FloatMmcAlg::LargeValue, double>::result == 0) |
|
153 |
check(false, "Wrong LargeValue type"); |
|
154 |
} |
|
155 |
|
|
156 |
// Run various tests |
|
157 |
{ |
|
158 |
typedef SmartDigraph GR; |
|
159 |
DIGRAPH_TYPEDEFS(GR); |
|
160 |
|
|
161 |
GR gr; |
|
162 |
IntArcMap l1(gr), l2(gr), l3(gr), l4(gr); |
|
163 |
IntArcMap c1(gr), c2(gr), c3(gr), c4(gr); |
|
164 |
|
|
165 |
std::istringstream input(test_lgf); |
|
166 |
digraphReader(gr, input). |
|
167 |
arcMap("len1", l1). |
|
168 |
arcMap("len2", l2). |
|
169 |
arcMap("len3", l3). |
|
170 |
arcMap("len4", l4). |
|
171 |
arcMap("c1", c1). |
|
172 |
arcMap("c2", c2). |
|
173 |
arcMap("c3", c3). |
|
174 |
arcMap("c4", c4). |
|
175 |
run(); |
|
176 |
|
|
177 |
checkMmcAlg<MinMeanCycle<GR, IntArcMap> >(gr, l1, c1, 6, 3); |
|
178 |
checkMmcAlg<MinMeanCycle<GR, IntArcMap> >(gr, l2, c2, 5, 2); |
|
179 |
checkMmcAlg<MinMeanCycle<GR, IntArcMap> >(gr, l3, c3, 0, 1); |
|
180 |
checkMmcAlg<MinMeanCycle<GR, IntArcMap> >(gr, l4, c4, -1, 1); |
|
181 |
} |
|
182 |
|
|
183 |
return 0; |
|
184 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 | 5 |
* Copyright (C) 2003-2008 |
6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
8 | 8 |
* |
9 | 9 |
* Permission to use, modify and distribute this software is granted |
10 | 10 |
* provided that this copyright notice appears in all copies. For |
11 | 11 |
* precise terms see the accompanying LICENSE file. |
12 | 12 |
* |
13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
14 | 14 |
* express or implied, and with no claim as to its suitability for any |
15 | 15 |
* purpose. |
16 | 16 |
* |
17 | 17 |
*/ |
18 | 18 |
|
19 | 19 |
#ifndef LEMON_MIN_MEAN_CYCLE_H |
20 | 20 |
#define LEMON_MIN_MEAN_CYCLE_H |
21 | 21 |
|
22 | 22 |
/// \ingroup shortest_path |
23 | 23 |
/// |
24 | 24 |
/// \file |
25 | 25 |
/// \brief Howard's algorithm for finding a minimum mean cycle. |
26 | 26 |
|
27 | 27 |
#include <vector> |
28 |
#include <limits> |
|
28 | 29 |
#include <lemon/core.h> |
29 | 30 |
#include <lemon/path.h> |
30 | 31 |
#include <lemon/tolerance.h> |
31 | 32 |
#include <lemon/connectivity.h> |
32 | 33 |
|
33 | 34 |
namespace lemon { |
34 | 35 |
|
35 | 36 |
/// \brief Default traits class of MinMeanCycle class. |
36 | 37 |
/// |
37 | 38 |
/// Default traits class of MinMeanCycle class. |
38 | 39 |
/// \tparam GR The type of the digraph. |
39 | 40 |
/// \tparam LEN The type of the length map. |
40 | 41 |
/// It must conform to the \ref concepts::ReadMap "ReadMap" concept. |
41 | 42 |
#ifdef DOXYGEN |
42 | 43 |
template <typename GR, typename LEN> |
43 | 44 |
#else |
44 | 45 |
template <typename GR, typename LEN, |
45 | 46 |
bool integer = std::numeric_limits<typename LEN::Value>::is_integer> |
46 | 47 |
#endif |
47 | 48 |
struct MinMeanCycleDefaultTraits |
48 | 49 |
{ |
49 | 50 |
/// The type of the digraph |
50 | 51 |
typedef GR Digraph; |
51 | 52 |
/// The type of the length map |
52 | 53 |
typedef LEN LengthMap; |
53 | 54 |
/// The type of the arc lengths |
54 | 55 |
typedef typename LengthMap::Value Value; |
55 | 56 |
|
56 | 57 |
/// \brief The large value type used for internal computations |
57 | 58 |
/// |
58 | 59 |
/// The large value type used for internal computations. |
59 | 60 |
/// It is \c long \c long if the \c Value type is integer, |
60 | 61 |
/// otherwise it is \c double. |
61 | 62 |
/// \c Value must be convertible to \c LargeValue. |
62 | 63 |
typedef double LargeValue; |
63 | 64 |
|
64 | 65 |
/// The tolerance type used for internal computations |
65 | 66 |
typedef lemon::Tolerance<LargeValue> Tolerance; |
66 | 67 |
|
67 | 68 |
/// \brief The path type of the found cycles |
68 | 69 |
/// |
69 | 70 |
/// The path type of the found cycles. |
70 | 71 |
/// It must conform to the \ref lemon::concepts::Path "Path" concept |
71 | 72 |
/// and it must have an \c addBack() function. |
72 | 73 |
typedef lemon::Path<Digraph> Path; |
73 | 74 |
}; |
74 | 75 |
|
75 | 76 |
// Default traits class for integer value types |
76 | 77 |
template <typename GR, typename LEN> |
77 | 78 |
struct MinMeanCycleDefaultTraits<GR, LEN, true> |
78 | 79 |
{ |
79 | 80 |
typedef GR Digraph; |
80 | 81 |
typedef LEN LengthMap; |
81 | 82 |
typedef typename LengthMap::Value Value; |
82 | 83 |
#ifdef LEMON_HAVE_LONG_LONG |
83 | 84 |
typedef long long LargeValue; |
84 | 85 |
#else |
85 | 86 |
typedef long LargeValue; |
86 | 87 |
#endif |
87 | 88 |
typedef lemon::Tolerance<LargeValue> Tolerance; |
88 | 89 |
typedef lemon::Path<Digraph> Path; |
89 | 90 |
}; |
90 | 91 |
|
91 | 92 |
|
92 | 93 |
/// \addtogroup shortest_path |
93 | 94 |
/// @{ |
94 | 95 |
|
95 | 96 |
/// \brief Implementation of Howard's algorithm for finding a minimum |
96 | 97 |
/// mean cycle. |
97 | 98 |
/// |
98 | 99 |
/// \ref MinMeanCycle implements Howard's algorithm for finding a |
99 | 100 |
/// directed cycle of minimum mean length (cost) in a digraph. |
100 | 101 |
/// |
101 | 102 |
/// \tparam GR The type of the digraph the algorithm runs on. |
102 | 103 |
/// \tparam LEN The type of the length map. The default |
103 | 104 |
/// map type is \ref concepts::Digraph::ArcMap "GR::ArcMap<int>". |
104 | 105 |
#ifdef DOXYGEN |
105 | 106 |
template <typename GR, typename LEN, typename TR> |
106 | 107 |
#else |
107 | 108 |
template < typename GR, |
108 | 109 |
typename LEN = typename GR::template ArcMap<int>, |
109 | 110 |
typename TR = MinMeanCycleDefaultTraits<GR, LEN> > |
110 | 111 |
#endif |
111 | 112 |
class MinMeanCycle |
112 | 113 |
{ |
113 | 114 |
public: |
114 | 115 |
|
115 | 116 |
/// The type of the digraph |
116 | 117 |
typedef typename TR::Digraph Digraph; |
117 | 118 |
/// The type of the length map |
118 | 119 |
typedef typename TR::LengthMap LengthMap; |
119 | 120 |
/// The type of the arc lengths |
120 | 121 |
typedef typename TR::Value Value; |
121 | 122 |
|
122 | 123 |
/// \brief The large value type |
123 | 124 |
/// |
124 | 125 |
/// The large value type used for internal computations. |
125 | 126 |
/// Using the \ref MinMeanCycleDefaultTraits "default traits class", |
126 | 127 |
/// it is \c long \c long if the \c Value type is integer, |
127 | 128 |
/// otherwise it is \c double. |
128 | 129 |
typedef typename TR::LargeValue LargeValue; |
129 | 130 |
|
130 | 131 |
/// The tolerance type |
131 | 132 |
typedef typename TR::Tolerance Tolerance; |
132 | 133 |
|
133 | 134 |
/// \brief The path type of the found cycles |
134 | 135 |
/// |
135 | 136 |
/// The path type of the found cycles. |
136 | 137 |
/// Using the \ref MinMeanCycleDefaultTraits "default traits class", |
137 | 138 |
/// it is \ref lemon::Path "Path<Digraph>". |
138 | 139 |
typedef typename TR::Path Path; |
139 | 140 |
|
140 | 141 |
/// The \ref MinMeanCycleDefaultTraits "traits class" of the algorithm |
141 | 142 |
typedef TR Traits; |
142 | 143 |
|
143 | 144 |
private: |
144 | 145 |
|
145 | 146 |
TEMPLATE_DIGRAPH_TYPEDEFS(Digraph); |
146 | 147 |
|
147 | 148 |
// The digraph the algorithm runs on |
148 | 149 |
const Digraph &_gr; |
149 | 150 |
// The length of the arcs |
150 | 151 |
const LengthMap &_length; |
151 | 152 |
|
152 | 153 |
// Data for the found cycles |
153 | 154 |
bool _curr_found, _best_found; |
154 | 155 |
LargeValue _curr_length, _best_length; |
155 | 156 |
int _curr_size, _best_size; |
156 | 157 |
Node _curr_node, _best_node; |
157 | 158 |
|
158 | 159 |
Path *_cycle_path; |
159 | 160 |
bool _local_path; |
160 | 161 |
|
161 | 162 |
// Internal data used by the algorithm |
162 | 163 |
typename Digraph::template NodeMap<Arc> _policy; |
163 | 164 |
typename Digraph::template NodeMap<bool> _reached; |
164 | 165 |
typename Digraph::template NodeMap<int> _level; |
165 | 166 |
typename Digraph::template NodeMap<LargeValue> _dist; |
166 | 167 |
|
167 | 168 |
// Data for storing the strongly connected components |
168 | 169 |
int _comp_num; |
169 | 170 |
typename Digraph::template NodeMap<int> _comp; |
170 | 171 |
std::vector<std::vector<Node> > _comp_nodes; |
171 | 172 |
std::vector<Node>* _nodes; |
172 | 173 |
typename Digraph::template NodeMap<std::vector<Arc> > _in_arcs; |
173 | 174 |
|
174 | 175 |
// Queue used for BFS search |
175 | 176 |
std::vector<Node> _queue; |
176 | 177 |
int _qfront, _qback; |
177 | 178 |
|
178 | 179 |
Tolerance _tolerance; |
179 | 180 |
|
180 | 181 |
public: |
181 | 182 |
|
182 | 183 |
/// \name Named Template Parameters |
183 | 184 |
/// @{ |
184 | 185 |
|
185 | 186 |
template <typename T> |
186 | 187 |
struct SetLargeValueTraits : public Traits { |
187 | 188 |
typedef T LargeValue; |
188 | 189 |
typedef lemon::Tolerance<T> Tolerance; |
189 | 190 |
}; |
190 | 191 |
|
191 | 192 |
/// \brief \ref named-templ-param "Named parameter" for setting |
192 | 193 |
/// \c LargeValue type. |
193 | 194 |
/// |
194 | 195 |
/// \ref named-templ-param "Named parameter" for setting \c LargeValue |
195 | 196 |
/// type. It is used for internal computations in the algorithm. |
196 | 197 |
template <typename T> |
197 | 198 |
struct SetLargeValue |
198 | 199 |
: public MinMeanCycle<GR, LEN, SetLargeValueTraits<T> > { |
199 | 200 |
typedef MinMeanCycle<GR, LEN, SetLargeValueTraits<T> > Create; |
200 | 201 |
}; |
201 | 202 |
|
202 | 203 |
template <typename T> |
203 | 204 |
struct SetPathTraits : public Traits { |
204 | 205 |
typedef T Path; |
205 | 206 |
}; |
206 | 207 |
|
207 | 208 |
/// \brief \ref named-templ-param "Named parameter" for setting |
208 | 209 |
/// \c %Path type. |
209 | 210 |
/// |
210 | 211 |
/// \ref named-templ-param "Named parameter" for setting the \c %Path |
211 | 212 |
/// type of the found cycles. |
212 | 213 |
/// It must conform to the \ref lemon::concepts::Path "Path" concept |
213 | 214 |
/// and it must have an \c addBack() function. |
214 | 215 |
template <typename T> |
215 | 216 |
struct SetPath |
216 | 217 |
: public MinMeanCycle<GR, LEN, SetPathTraits<T> > { |
217 | 218 |
typedef MinMeanCycle<GR, LEN, SetPathTraits<T> > Create; |
218 | 219 |
}; |
219 | 220 |
|
220 | 221 |
/// @} |
221 | 222 |
|
222 | 223 |
public: |
223 | 224 |
|
224 | 225 |
/// \brief Constructor. |
225 | 226 |
/// |
226 | 227 |
/// The constructor of the class. |
227 | 228 |
/// |
228 | 229 |
/// \param digraph The digraph the algorithm runs on. |
229 | 230 |
/// \param length The lengths (costs) of the arcs. |
230 | 231 |
MinMeanCycle( const Digraph &digraph, |
231 | 232 |
const LengthMap &length ) : |
232 | 233 |
_gr(digraph), _length(length), _cycle_path(NULL), _local_path(false), |
233 | 234 |
_policy(digraph), _reached(digraph), _level(digraph), _dist(digraph), |
234 | 235 |
_comp(digraph), _in_arcs(digraph) |
235 | 236 |
{} |
236 | 237 |
|
237 | 238 |
/// Destructor. |
238 | 239 |
~MinMeanCycle() { |
239 | 240 |
if (_local_path) delete _cycle_path; |
240 | 241 |
} |
241 | 242 |
|
242 | 243 |
/// \brief Set the path structure for storing the found cycle. |
243 | 244 |
/// |
244 | 245 |
/// This function sets an external path structure for storing the |
245 | 246 |
/// found cycle. |
246 | 247 |
/// |
247 | 248 |
/// If you don't call this function before calling \ref run() or |
248 | 249 |
/// \ref findMinMean(), it will allocate a local \ref Path "path" |
249 | 250 |
/// structure. The destuctor deallocates this automatically |
250 | 251 |
/// allocated object, of course. |
251 | 252 |
/// |
252 | 253 |
/// \note The algorithm calls only the \ref lemon::Path::addBack() |
253 | 254 |
/// "addBack()" function of the given path structure. |
254 | 255 |
/// |
255 | 256 |
/// \return <tt>(*this)</tt> |
256 | 257 |
MinMeanCycle& cycle(Path &path) { |
257 | 258 |
if (_local_path) { |
258 | 259 |
delete _cycle_path; |
259 | 260 |
_local_path = false; |
260 | 261 |
} |
261 | 262 |
_cycle_path = &path; |
262 | 263 |
return *this; |
263 | 264 |
} |
264 | 265 |
|
265 | 266 |
/// \name Execution control |
266 | 267 |
/// The simplest way to execute the algorithm is to call the \ref run() |
267 | 268 |
/// function.\n |
268 | 269 |
/// If you only need the minimum mean length, you may call |
269 | 270 |
/// \ref findMinMean(). |
270 | 271 |
|
271 | 272 |
/// @{ |
272 | 273 |
|
273 | 274 |
/// \brief Run the algorithm. |
274 | 275 |
/// |
275 | 276 |
/// This function runs the algorithm. |
276 | 277 |
/// It can be called more than once (e.g. if the underlying digraph |
277 | 278 |
/// and/or the arc lengths have been modified). |
278 | 279 |
/// |
279 | 280 |
/// \return \c true if a directed cycle exists in the digraph. |
280 | 281 |
/// |
281 | 282 |
/// \note <tt>mmc.run()</tt> is just a shortcut of the following code. |
282 | 283 |
/// \code |
283 | 284 |
/// return mmc.findMinMean() && mmc.findCycle(); |
284 | 285 |
/// \endcode |
285 | 286 |
bool run() { |
286 | 287 |
return findMinMean() && findCycle(); |
287 | 288 |
} |
288 | 289 |
|
289 | 290 |
/// \brief Find the minimum cycle mean. |
290 | 291 |
/// |
291 | 292 |
/// This function finds the minimum mean length of the directed |
292 | 293 |
/// cycles in the digraph. |
293 | 294 |
/// |
294 | 295 |
/// \return \c true if a directed cycle exists in the digraph. |
295 | 296 |
bool findMinMean() { |
296 | 297 |
// Initialize and find strongly connected components |
297 | 298 |
init(); |
298 | 299 |
findComponents(); |
299 | 300 |
|
300 | 301 |
// Find the minimum cycle mean in the components |
301 | 302 |
for (int comp = 0; comp < _comp_num; ++comp) { |
302 | 303 |
// Find the minimum mean cycle in the current component |
303 | 304 |
if (!buildPolicyGraph(comp)) continue; |
304 | 305 |
while (true) { |
305 | 306 |
findPolicyCycle(); |
306 | 307 |
if (!computeNodeDistances()) break; |
307 | 308 |
} |
308 | 309 |
// Update the best cycle (global minimum mean cycle) |
309 | 310 |
if ( !_best_found || (_curr_found && |
310 | 311 |
_curr_length * _best_size < _best_length * _curr_size) ) { |
311 | 312 |
_best_found = true; |
312 | 313 |
_best_length = _curr_length; |
313 | 314 |
_best_size = _curr_size; |
314 | 315 |
_best_node = _curr_node; |
315 | 316 |
} |
316 | 317 |
} |
317 | 318 |
return _best_found; |
318 | 319 |
} |
319 | 320 |
|
320 | 321 |
/// \brief Find a minimum mean directed cycle. |
321 | 322 |
/// |
322 | 323 |
/// This function finds a directed cycle of minimum mean length |
323 | 324 |
/// in the digraph using the data computed by findMinMean(). |
324 | 325 |
/// |
325 | 326 |
/// \return \c true if a directed cycle exists in the digraph. |
326 | 327 |
/// |
327 | 328 |
/// \pre \ref findMinMean() must be called before using this function. |
328 | 329 |
bool findCycle() { |
329 | 330 |
if (!_best_found) return false; |
330 | 331 |
_cycle_path->addBack(_policy[_best_node]); |
331 | 332 |
for ( Node v = _best_node; |
332 | 333 |
(v = _gr.target(_policy[v])) != _best_node; ) { |
333 | 334 |
_cycle_path->addBack(_policy[v]); |
334 | 335 |
} |
335 | 336 |
return true; |
336 | 337 |
} |
337 | 338 |
|
338 | 339 |
/// @} |
339 | 340 |
|
340 | 341 |
/// \name Query Functions |
341 | 342 |
/// The results of the algorithm can be obtained using these |
342 | 343 |
/// functions.\n |
343 | 344 |
/// The algorithm should be executed before using them. |
344 | 345 |
|
345 | 346 |
/// @{ |
346 | 347 |
|
347 | 348 |
/// \brief Return the total length of the found cycle. |
348 | 349 |
/// |
349 | 350 |
/// This function returns the total length of the found cycle. |
350 | 351 |
/// |
351 | 352 |
/// \pre \ref run() or \ref findMinMean() must be called before |
352 | 353 |
/// using this function. |
353 | 354 |
LargeValue cycleLength() const { |
354 | 355 |
return _best_length; |
355 | 356 |
} |
356 | 357 |
|
357 | 358 |
/// \brief Return the number of arcs on the found cycle. |
358 | 359 |
/// |
359 | 360 |
/// This function returns the number of arcs on the found cycle. |
360 | 361 |
/// |
361 | 362 |
/// \pre \ref run() or \ref findMinMean() must be called before |
362 | 363 |
/// using this function. |
363 | 364 |
int cycleArcNum() const { |
364 | 365 |
return _best_size; |
365 | 366 |
} |
366 | 367 |
|
367 | 368 |
/// \brief Return the mean length of the found cycle. |
368 | 369 |
/// |
369 | 370 |
/// This function returns the mean length of the found cycle. |
370 | 371 |
/// |
371 | 372 |
/// \note <tt>alg.cycleMean()</tt> is just a shortcut of the |
372 | 373 |
/// following code. |
373 | 374 |
/// \code |
374 | 375 |
/// return static_cast<double>(alg.cycleLength()) / alg.cycleArcNum(); |
375 | 376 |
/// \endcode |
376 | 377 |
/// |
377 | 378 |
/// \pre \ref run() or \ref findMinMean() must be called before |
378 | 379 |
/// using this function. |
379 | 380 |
double cycleMean() const { |
380 | 381 |
return static_cast<double>(_best_length) / _best_size; |
381 | 382 |
} |
382 | 383 |
|
383 | 384 |
/// \brief Return the found cycle. |
384 | 385 |
/// |
385 | 386 |
/// This function returns a const reference to the path structure |
386 | 387 |
/// storing the found cycle. |
387 | 388 |
/// |
388 | 389 |
/// \pre \ref run() or \ref findCycle() must be called before using |
389 | 390 |
/// this function. |
390 | 391 |
const Path& cycle() const { |
391 | 392 |
return *_cycle_path; |
392 | 393 |
} |
393 | 394 |
|
394 | 395 |
///@} |
395 | 396 |
|
396 | 397 |
private: |
397 | 398 |
|
398 | 399 |
// Initialize |
399 | 400 |
void init() { |
400 | 401 |
if (!_cycle_path) { |
401 | 402 |
_local_path = true; |
402 | 403 |
_cycle_path = new Path; |
403 | 404 |
} |
404 | 405 |
_queue.resize(countNodes(_gr)); |
405 | 406 |
_best_found = false; |
406 | 407 |
_best_length = 0; |
407 | 408 |
_best_size = 1; |
408 | 409 |
_cycle_path->clear(); |
409 | 410 |
} |
410 | 411 |
|
411 | 412 |
// Find strongly connected components and initialize _comp_nodes |
412 | 413 |
// and _in_arcs |
413 | 414 |
void findComponents() { |
414 | 415 |
_comp_num = stronglyConnectedComponents(_gr, _comp); |
415 | 416 |
_comp_nodes.resize(_comp_num); |
416 | 417 |
if (_comp_num == 1) { |
417 | 418 |
_comp_nodes[0].clear(); |
418 | 419 |
for (NodeIt n(_gr); n != INVALID; ++n) { |
419 | 420 |
_comp_nodes[0].push_back(n); |
420 | 421 |
_in_arcs[n].clear(); |
421 | 422 |
for (InArcIt a(_gr, n); a != INVALID; ++a) { |
422 | 423 |
_in_arcs[n].push_back(a); |
423 | 424 |
} |
424 | 425 |
} |
425 | 426 |
} else { |
426 | 427 |
for (int i = 0; i < _comp_num; ++i) |
427 | 428 |
_comp_nodes[i].clear(); |
428 | 429 |
for (NodeIt n(_gr); n != INVALID; ++n) { |
429 | 430 |
int k = _comp[n]; |
430 | 431 |
_comp_nodes[k].push_back(n); |
431 | 432 |
_in_arcs[n].clear(); |
432 | 433 |
for (InArcIt a(_gr, n); a != INVALID; ++a) { |
433 | 434 |
if (_comp[_gr.source(a)] == k) _in_arcs[n].push_back(a); |
434 | 435 |
} |
435 | 436 |
} |
436 | 437 |
} |
437 | 438 |
} |
438 | 439 |
|
439 | 440 |
// Build the policy graph in the given strongly connected component |
440 | 441 |
// (the out-degree of every node is 1) |
441 | 442 |
bool buildPolicyGraph(int comp) { |
442 | 443 |
_nodes = &(_comp_nodes[comp]); |
443 | 444 |
if (_nodes->size() < 1 || |
444 | 445 |
(_nodes->size() == 1 && _in_arcs[(*_nodes)[0]].size() == 0)) { |
445 | 446 |
return false; |
446 | 447 |
} |
447 | 448 |
for (int i = 0; i < int(_nodes->size()); ++i) { |
448 | 449 |
_dist[(*_nodes)[i]] = std::numeric_limits<LargeValue>::max(); |
449 | 450 |
} |
450 | 451 |
Node u, v; |
451 | 452 |
Arc e; |
452 | 453 |
for (int i = 0; i < int(_nodes->size()); ++i) { |
453 | 454 |
v = (*_nodes)[i]; |
454 | 455 |
for (int j = 0; j < int(_in_arcs[v].size()); ++j) { |
455 | 456 |
e = _in_arcs[v][j]; |
456 | 457 |
u = _gr.source(e); |
457 | 458 |
if (_length[e] < _dist[u]) { |
458 | 459 |
_dist[u] = _length[e]; |
459 | 460 |
_policy[u] = e; |
460 | 461 |
} |
461 | 462 |
} |
462 | 463 |
} |
463 | 464 |
return true; |
464 | 465 |
} |
465 | 466 |
|
466 | 467 |
// Find the minimum mean cycle in the policy graph |
467 | 468 |
void findPolicyCycle() { |
468 | 469 |
for (int i = 0; i < int(_nodes->size()); ++i) { |
469 | 470 |
_level[(*_nodes)[i]] = -1; |
470 | 471 |
} |
471 | 472 |
LargeValue clength; |
472 | 473 |
int csize; |
473 | 474 |
Node u, v; |
474 | 475 |
_curr_found = false; |
475 | 476 |
for (int i = 0; i < int(_nodes->size()); ++i) { |
476 | 477 |
u = (*_nodes)[i]; |
477 | 478 |
if (_level[u] >= 0) continue; |
478 | 479 |
for (; _level[u] < 0; u = _gr.target(_policy[u])) { |
479 | 480 |
_level[u] = i; |
480 | 481 |
} |
481 | 482 |
if (_level[u] == i) { |
482 | 483 |
// A cycle is found |
483 | 484 |
clength = _length[_policy[u]]; |
484 | 485 |
csize = 1; |
485 | 486 |
for (v = u; (v = _gr.target(_policy[v])) != u; ) { |
486 | 487 |
clength += _length[_policy[v]]; |
487 | 488 |
++csize; |
488 | 489 |
} |
489 | 490 |
if ( !_curr_found || |
490 | 491 |
(clength * _curr_size < _curr_length * csize) ) { |
491 | 492 |
_curr_found = true; |
492 | 493 |
_curr_length = clength; |
493 | 494 |
_curr_size = csize; |
494 | 495 |
_curr_node = u; |
495 | 496 |
} |
496 | 497 |
} |
497 | 498 |
} |
498 | 499 |
} |
499 | 500 |
|
500 | 501 |
// Contract the policy graph and compute node distances |
501 | 502 |
bool computeNodeDistances() { |
502 | 503 |
// Find the component of the main cycle and compute node distances |
503 | 504 |
// using reverse BFS |
504 | 505 |
for (int i = 0; i < int(_nodes->size()); ++i) { |
505 | 506 |
_reached[(*_nodes)[i]] = false; |
506 | 507 |
} |
507 | 508 |
_qfront = _qback = 0; |
508 | 509 |
_queue[0] = _curr_node; |
509 | 510 |
_reached[_curr_node] = true; |
510 | 511 |
_dist[_curr_node] = 0; |
511 | 512 |
Node u, v; |
512 | 513 |
Arc e; |
513 | 514 |
while (_qfront <= _qback) { |
514 | 515 |
v = _queue[_qfront++]; |
515 | 516 |
for (int j = 0; j < int(_in_arcs[v].size()); ++j) { |
516 | 517 |
e = _in_arcs[v][j]; |
517 | 518 |
u = _gr.source(e); |
518 | 519 |
if (_policy[u] == e && !_reached[u]) { |
519 | 520 |
_reached[u] = true; |
520 | 521 |
_dist[u] = _dist[v] + _length[e] * _curr_size - _curr_length; |
521 | 522 |
_queue[++_qback] = u; |
522 | 523 |
} |
523 | 524 |
} |
524 | 525 |
} |
525 | 526 |
|
526 | 527 |
// Connect all other nodes to this component and compute node |
527 | 528 |
// distances using reverse BFS |
528 | 529 |
_qfront = 0; |
529 | 530 |
while (_qback < int(_nodes->size())-1) { |
530 | 531 |
v = _queue[_qfront++]; |
531 | 532 |
for (int j = 0; j < int(_in_arcs[v].size()); ++j) { |
532 | 533 |
e = _in_arcs[v][j]; |
533 | 534 |
u = _gr.source(e); |
534 | 535 |
if (!_reached[u]) { |
535 | 536 |
_reached[u] = true; |
536 | 537 |
_policy[u] = e; |
537 | 538 |
_dist[u] = _dist[v] + _length[e] * _curr_size - _curr_length; |
538 | 539 |
_queue[++_qback] = u; |
539 | 540 |
} |
1 | 1 |
INCLUDE_DIRECTORIES( |
2 | 2 |
${PROJECT_SOURCE_DIR} |
3 | 3 |
${PROJECT_BINARY_DIR} |
4 | 4 |
) |
5 | 5 |
|
6 | 6 |
LINK_DIRECTORIES( |
7 | 7 |
${PROJECT_BINARY_DIR}/lemon |
8 | 8 |
) |
9 | 9 |
|
10 | 10 |
SET(TESTS |
11 | 11 |
adaptors_test |
12 | 12 |
bfs_test |
13 | 13 |
circulation_test |
14 | 14 |
connectivity_test |
15 | 15 |
counter_test |
16 | 16 |
dfs_test |
17 | 17 |
digraph_test |
18 | 18 |
dijkstra_test |
19 | 19 |
dim_test |
20 | 20 |
edge_set_test |
21 | 21 |
error_test |
22 | 22 |
euler_test |
23 | 23 |
gomory_hu_test |
24 | 24 |
graph_copy_test |
25 | 25 |
graph_test |
26 | 26 |
graph_utils_test |
27 | 27 |
hao_orlin_test |
28 | 28 |
heap_test |
29 | 29 |
kruskal_test |
30 | 30 |
maps_test |
31 | 31 |
matching_test |
32 | 32 |
min_cost_arborescence_test |
33 | 33 |
min_cost_flow_test |
34 |
min_mean_cycle_test |
|
34 | 35 |
path_test |
35 | 36 |
preflow_test |
36 | 37 |
radix_sort_test |
37 | 38 |
random_test |
38 | 39 |
suurballe_test |
39 | 40 |
time_measure_test |
40 | 41 |
unionfind_test |
41 | 42 |
) |
42 | 43 |
|
43 | 44 |
IF(LEMON_HAVE_LP) |
44 | 45 |
ADD_EXECUTABLE(lp_test lp_test.cc) |
45 | 46 |
SET(LP_TEST_LIBS lemon) |
46 | 47 |
|
47 | 48 |
IF(LEMON_HAVE_GLPK) |
48 | 49 |
SET(LP_TEST_LIBS ${LP_TEST_LIBS} ${GLPK_LIBRARIES}) |
49 | 50 |
ENDIF() |
50 | 51 |
IF(LEMON_HAVE_CPLEX) |
51 | 52 |
SET(LP_TEST_LIBS ${LP_TEST_LIBS} ${CPLEX_LIBRARIES}) |
52 | 53 |
ENDIF() |
53 | 54 |
IF(LEMON_HAVE_CLP) |
54 | 55 |
SET(LP_TEST_LIBS ${LP_TEST_LIBS} ${COIN_CLP_LIBRARIES}) |
55 | 56 |
ENDIF() |
56 | 57 |
|
57 | 58 |
TARGET_LINK_LIBRARIES(lp_test ${LP_TEST_LIBS}) |
58 | 59 |
ADD_TEST(lp_test lp_test) |
59 | 60 |
|
60 | 61 |
IF(WIN32 AND LEMON_HAVE_GLPK) |
61 | 62 |
GET_TARGET_PROPERTY(TARGET_LOC lp_test LOCATION) |
62 | 63 |
GET_FILENAME_COMPONENT(TARGET_PATH ${TARGET_LOC} PATH) |
63 | 64 |
ADD_CUSTOM_COMMAND(TARGET lp_test POST_BUILD |
64 | 65 |
COMMAND ${CMAKE_COMMAND} -E copy ${GLPK_BIN_DIR}/glpk.dll ${TARGET_PATH} |
65 | 66 |
COMMAND ${CMAKE_COMMAND} -E copy ${GLPK_BIN_DIR}/libltdl3.dll ${TARGET_PATH} |
66 | 67 |
COMMAND ${CMAKE_COMMAND} -E copy ${GLPK_BIN_DIR}/zlib1.dll ${TARGET_PATH} |
67 | 68 |
) |
68 | 69 |
ENDIF() |
69 | 70 |
|
70 | 71 |
IF(WIN32 AND LEMON_HAVE_CPLEX) |
71 | 72 |
GET_TARGET_PROPERTY(TARGET_LOC lp_test LOCATION) |
72 | 73 |
GET_FILENAME_COMPONENT(TARGET_PATH ${TARGET_LOC} PATH) |
73 | 74 |
ADD_CUSTOM_COMMAND(TARGET lp_test POST_BUILD |
74 | 75 |
COMMAND ${CMAKE_COMMAND} -E copy ${CPLEX_BIN_DIR}/cplex91.dll ${TARGET_PATH} |
75 | 76 |
) |
76 | 77 |
ENDIF() |
77 | 78 |
ENDIF() |
78 | 79 |
|
79 | 80 |
IF(LEMON_HAVE_MIP) |
80 | 81 |
ADD_EXECUTABLE(mip_test mip_test.cc) |
81 | 82 |
SET(MIP_TEST_LIBS lemon) |
82 | 83 |
|
83 | 84 |
IF(LEMON_HAVE_GLPK) |
84 | 85 |
SET(MIP_TEST_LIBS ${MIP_TEST_LIBS} ${GLPK_LIBRARIES}) |
85 | 86 |
ENDIF() |
86 | 87 |
IF(LEMON_HAVE_CPLEX) |
87 | 88 |
SET(MIP_TEST_LIBS ${MIP_TEST_LIBS} ${CPLEX_LIBRARIES}) |
88 | 89 |
ENDIF() |
89 | 90 |
IF(LEMON_HAVE_CBC) |
90 | 91 |
SET(MIP_TEST_LIBS ${MIP_TEST_LIBS} ${COIN_CBC_LIBRARIES}) |
91 | 92 |
ENDIF() |
92 | 93 |
|
93 | 94 |
TARGET_LINK_LIBRARIES(mip_test ${MIP_TEST_LIBS}) |
94 | 95 |
ADD_TEST(mip_test mip_test) |
95 | 96 |
|
96 | 97 |
IF(WIN32 AND LEMON_HAVE_GLPK) |
97 | 98 |
GET_TARGET_PROPERTY(TARGET_LOC mip_test LOCATION) |
98 | 99 |
GET_FILENAME_COMPONENT(TARGET_PATH ${TARGET_LOC} PATH) |
99 | 100 |
ADD_CUSTOM_COMMAND(TARGET mip_test POST_BUILD |
100 | 101 |
COMMAND ${CMAKE_COMMAND} -E copy ${GLPK_BIN_DIR}/glpk.dll ${TARGET_PATH} |
101 | 102 |
COMMAND ${CMAKE_COMMAND} -E copy ${GLPK_BIN_DIR}/libltdl3.dll ${TARGET_PATH} |
102 | 103 |
COMMAND ${CMAKE_COMMAND} -E copy ${GLPK_BIN_DIR}/zlib1.dll ${TARGET_PATH} |
103 | 104 |
) |
104 | 105 |
ENDIF() |
105 | 106 |
|
106 | 107 |
IF(WIN32 AND LEMON_HAVE_CPLEX) |
107 | 108 |
GET_TARGET_PROPERTY(TARGET_LOC mip_test LOCATION) |
108 | 109 |
GET_FILENAME_COMPONENT(TARGET_PATH ${TARGET_LOC} PATH) |
109 | 110 |
ADD_CUSTOM_COMMAND(TARGET mip_test POST_BUILD |
110 | 111 |
COMMAND ${CMAKE_COMMAND} -E copy ${CPLEX_BIN_DIR}/cplex91.dll ${TARGET_PATH} |
111 | 112 |
) |
112 | 113 |
ENDIF() |
113 | 114 |
ENDIF() |
114 | 115 |
|
115 | 116 |
FOREACH(TEST_NAME ${TESTS}) |
116 | 117 |
ADD_EXECUTABLE(${TEST_NAME} ${TEST_NAME}.cc) |
117 | 118 |
TARGET_LINK_LIBRARIES(${TEST_NAME} lemon) |
118 | 119 |
ADD_TEST(${TEST_NAME} ${TEST_NAME}) |
119 | 120 |
ENDFOREACH() |
1 | 1 |
EXTRA_DIST += \ |
2 | 2 |
test/CMakeLists.txt |
3 | 3 |
|
4 | 4 |
noinst_HEADERS += \ |
5 | 5 |
test/graph_test.h \ |
6 | 6 |
test/test_tools.h |
7 | 7 |
|
8 | 8 |
check_PROGRAMS += \ |
9 | 9 |
test/adaptors_test \ |
10 | 10 |
test/bfs_test \ |
11 | 11 |
test/circulation_test \ |
12 | 12 |
test/connectivity_test \ |
13 | 13 |
test/counter_test \ |
14 | 14 |
test/dfs_test \ |
15 | 15 |
test/digraph_test \ |
16 | 16 |
test/dijkstra_test \ |
17 | 17 |
test/dim_test \ |
18 | 18 |
test/edge_set_test \ |
19 | 19 |
test/error_test \ |
20 | 20 |
test/euler_test \ |
21 | 21 |
test/gomory_hu_test \ |
22 | 22 |
test/graph_copy_test \ |
23 | 23 |
test/graph_test \ |
24 | 24 |
test/graph_utils_test \ |
25 | 25 |
test/hao_orlin_test \ |
26 | 26 |
test/heap_test \ |
27 | 27 |
test/kruskal_test \ |
28 | 28 |
test/maps_test \ |
29 | 29 |
test/matching_test \ |
30 | 30 |
test/min_cost_arborescence_test \ |
31 | 31 |
test/min_cost_flow_test \ |
32 |
test/min_mean_cycle_test \ |
|
32 | 33 |
test/path_test \ |
33 | 34 |
test/preflow_test \ |
34 | 35 |
test/radix_sort_test \ |
35 | 36 |
test/random_test \ |
36 | 37 |
test/suurballe_test \ |
37 | 38 |
test/test_tools_fail \ |
38 | 39 |
test/test_tools_pass \ |
39 | 40 |
test/time_measure_test \ |
40 | 41 |
test/unionfind_test |
41 | 42 |
|
42 | 43 |
test_test_tools_pass_DEPENDENCIES = demo |
43 | 44 |
|
44 | 45 |
if HAVE_LP |
45 | 46 |
check_PROGRAMS += test/lp_test |
46 | 47 |
endif HAVE_LP |
47 | 48 |
if HAVE_MIP |
48 | 49 |
check_PROGRAMS += test/mip_test |
49 | 50 |
endif HAVE_MIP |
50 | 51 |
|
51 | 52 |
TESTS += $(check_PROGRAMS) |
52 | 53 |
XFAIL_TESTS += test/test_tools_fail$(EXEEXT) |
53 | 54 |
|
54 | 55 |
test_adaptors_test_SOURCES = test/adaptors_test.cc |
55 | 56 |
test_bfs_test_SOURCES = test/bfs_test.cc |
56 | 57 |
test_circulation_test_SOURCES = test/circulation_test.cc |
57 | 58 |
test_counter_test_SOURCES = test/counter_test.cc |
58 | 59 |
test_connectivity_test_SOURCES = test/connectivity_test.cc |
59 | 60 |
test_dfs_test_SOURCES = test/dfs_test.cc |
60 | 61 |
test_digraph_test_SOURCES = test/digraph_test.cc |
61 | 62 |
test_dijkstra_test_SOURCES = test/dijkstra_test.cc |
62 | 63 |
test_dim_test_SOURCES = test/dim_test.cc |
63 | 64 |
test_edge_set_test_SOURCES = test/edge_set_test.cc |
64 | 65 |
test_error_test_SOURCES = test/error_test.cc |
65 | 66 |
test_euler_test_SOURCES = test/euler_test.cc |
66 | 67 |
test_gomory_hu_test_SOURCES = test/gomory_hu_test.cc |
67 | 68 |
test_graph_copy_test_SOURCES = test/graph_copy_test.cc |
68 | 69 |
test_graph_test_SOURCES = test/graph_test.cc |
69 | 70 |
test_graph_utils_test_SOURCES = test/graph_utils_test.cc |
70 | 71 |
test_heap_test_SOURCES = test/heap_test.cc |
71 | 72 |
test_kruskal_test_SOURCES = test/kruskal_test.cc |
72 | 73 |
test_hao_orlin_test_SOURCES = test/hao_orlin_test.cc |
73 | 74 |
test_lp_test_SOURCES = test/lp_test.cc |
74 | 75 |
test_maps_test_SOURCES = test/maps_test.cc |
75 | 76 |
test_mip_test_SOURCES = test/mip_test.cc |
76 | 77 |
test_matching_test_SOURCES = test/matching_test.cc |
77 | 78 |
test_min_cost_arborescence_test_SOURCES = test/min_cost_arborescence_test.cc |
78 | 79 |
test_min_cost_flow_test_SOURCES = test/min_cost_flow_test.cc |
80 |
test_min_mean_cycle_test_SOURCES = test/min_mean_cycle_test.cc |
|
79 | 81 |
test_path_test_SOURCES = test/path_test.cc |
80 | 82 |
test_preflow_test_SOURCES = test/preflow_test.cc |
81 | 83 |
test_radix_sort_test_SOURCES = test/radix_sort_test.cc |
82 | 84 |
test_suurballe_test_SOURCES = test/suurballe_test.cc |
83 | 85 |
test_random_test_SOURCES = test/random_test.cc |
84 | 86 |
test_test_tools_fail_SOURCES = test/test_tools_fail.cc |
85 | 87 |
test_test_tools_pass_SOURCES = test/test_tools_pass.cc |
86 | 88 |
test_time_measure_test_SOURCES = test/time_measure_test.cc |
87 | 89 |
test_unionfind_test_SOURCES = test/unionfind_test.cc |
1 | 1 |
/* -*- mode: C++; indent-tabs-mode: nil; -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library. |
4 | 4 |
* |
5 | 5 |
* Copyright (C) 2003-2009 |
6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
8 | 8 |
* |
9 | 9 |
* Permission to use, modify and distribute this software is granted |
10 | 10 |
* provided that this copyright notice appears in all copies. For |
11 | 11 |
* precise terms see the accompanying LICENSE file. |
12 | 12 |
* |
13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
14 | 14 |
* express or implied, and with no claim as to its suitability for any |
15 | 15 |
* purpose. |
16 | 16 |
* |
17 | 17 |
*/ |
18 | 18 |
|
19 | 19 |
#ifndef LEMON_TEST_TEST_TOOLS_H |
20 | 20 |
#define LEMON_TEST_TEST_TOOLS_H |
21 | 21 |
|
22 | 22 |
///\ingroup misc |
23 | 23 |
///\file |
24 | 24 |
///\brief Some utilities to write test programs. |
25 | 25 |
|
26 | 26 |
#include <iostream> |
27 | 27 |
#include <stdlib.h> |
28 | 28 |
|
29 | 29 |
///If \c rc is fail, writes an error message and exits. |
30 | 30 |
|
31 | 31 |
///If \c rc is fail, writes an error message and exits. |
32 | 32 |
///The error message contains the file name and the line number of the |
33 | 33 |
///source code in a standard from, which makes it possible to go there |
34 | 34 |
///using good source browsers like e.g. \c emacs. |
35 | 35 |
/// |
36 | 36 |
///For example |
37 | 37 |
///\code check(0==1,"This is obviously false.");\endcode will |
38 | 38 |
///print something like this (and then exits). |
39 | 39 |
///\verbatim file_name.cc:123: error: This is obviously false. \endverbatim |
40 |
#define check(rc, msg) \ |
|
41 |
if(!(rc)) { \ |
|
42 |
std::cerr << __FILE__ ":" << __LINE__ << ": error: " << msg << std::endl; \ |
|
43 |
abort(); \ |
|
44 |
|
|
40 |
#define check(rc, msg) \ |
|
41 |
{ \ |
|
42 |
if(!(rc)) { \ |
|
43 |
std::cerr << __FILE__ ":" << __LINE__ << ": error: " \ |
|
44 |
<< msg << std::endl; \ |
|
45 |
abort(); \ |
|
46 |
} else { } \ |
|
47 |
} \ |
|
48 |
|
|
45 | 49 |
|
46 | 50 |
#endif |
0 comments (0 inline)