0
5
0
... | ... |
@@ -1233,64 +1233,69 @@ |
1233 | 1233 |
|
1234 | 1234 |
/// \brief The type of the digraph the algorithm runs on. |
1235 | 1235 |
typedef _Digraph Digraph; |
1236 | 1236 |
|
1237 | 1237 |
/// \brief The type of the map that indicates which nodes are reached. |
1238 | 1238 |
/// |
1239 | 1239 |
/// The type of the map that indicates which nodes are reached. |
1240 | 1240 |
/// It must meet the \ref concepts::ReadWriteMap "ReadWriteMap" concept. |
1241 | 1241 |
typedef typename Digraph::template NodeMap<bool> ReachedMap; |
1242 | 1242 |
|
1243 | 1243 |
/// \brief Instantiates a \ref ReachedMap. |
1244 | 1244 |
/// |
1245 | 1245 |
/// This function instantiates a \ref ReachedMap. |
1246 | 1246 |
/// \param digraph is the digraph, to which |
1247 | 1247 |
/// we would like to define the \ref ReachedMap. |
1248 | 1248 |
static ReachedMap *createReachedMap(const Digraph &digraph) { |
1249 | 1249 |
return new ReachedMap(digraph); |
1250 | 1250 |
} |
1251 | 1251 |
|
1252 | 1252 |
}; |
1253 | 1253 |
|
1254 | 1254 |
/// \ingroup search |
1255 | 1255 |
/// |
1256 | 1256 |
/// \brief %BFS algorithm class with visitor interface. |
1257 | 1257 |
/// |
1258 | 1258 |
/// This class provides an efficient implementation of the %BFS algorithm |
1259 | 1259 |
/// with visitor interface. |
1260 | 1260 |
/// |
1261 | 1261 |
/// The %BfsVisit class provides an alternative interface to the Bfs |
1262 | 1262 |
/// class. It works with callback mechanism, the BfsVisit object calls |
1263 | 1263 |
/// the member functions of the \c Visitor class on every BFS event. |
1264 | 1264 |
/// |
1265 |
/// This interface of the BFS algorithm should be used in special cases |
|
1266 |
/// when extra actions have to be performed in connection with certain |
|
1267 |
/// events of the BFS algorithm. Otherwise consider to use Bfs or bfs() |
|
1268 |
/// instead. |
|
1269 |
/// |
|
1265 | 1270 |
/// \tparam _Digraph The type of the digraph the algorithm runs on. |
1266 | 1271 |
/// The default value is |
1267 | 1272 |
/// \ref ListDigraph. The value of _Digraph is not used directly by |
1268 | 1273 |
/// \ref BfsVisit, it is only passed to \ref BfsVisitDefaultTraits. |
1269 | 1274 |
/// \tparam _Visitor The Visitor type that is used by the algorithm. |
1270 | 1275 |
/// \ref BfsVisitor "BfsVisitor<_Digraph>" is an empty visitor, which |
1271 | 1276 |
/// does not observe the BFS events. If you want to observe the BFS |
1272 | 1277 |
/// events, you should implement your own visitor class. |
1273 | 1278 |
/// \tparam _Traits Traits class to set various data types used by the |
1274 | 1279 |
/// algorithm. The default traits class is |
1275 | 1280 |
/// \ref BfsVisitDefaultTraits "BfsVisitDefaultTraits<_Digraph>". |
1276 | 1281 |
/// See \ref BfsVisitDefaultTraits for the documentation of |
1277 | 1282 |
/// a BFS visit traits class. |
1278 | 1283 |
#ifdef DOXYGEN |
1279 | 1284 |
template <typename _Digraph, typename _Visitor, typename _Traits> |
1280 | 1285 |
#else |
1281 | 1286 |
template <typename _Digraph = ListDigraph, |
1282 | 1287 |
typename _Visitor = BfsVisitor<_Digraph>, |
1283 | 1288 |
typename _Traits = BfsDefaultTraits<_Digraph> > |
1284 | 1289 |
#endif |
1285 | 1290 |
class BfsVisit { |
1286 | 1291 |
public: |
1287 | 1292 |
|
1288 | 1293 |
/// \brief \ref Exception for uninitialized parameters. |
1289 | 1294 |
/// |
1290 | 1295 |
/// This error represents problems in the initialization |
1291 | 1296 |
/// of the parameters of the algorithm. |
1292 | 1297 |
class UninitializedParameter : public lemon::UninitializedParameter { |
1293 | 1298 |
public: |
1294 | 1299 |
virtual const char* what() const throw() |
1295 | 1300 |
{ |
1296 | 1301 |
return "lemon::BfsVisit::UninitializedParameter"; |
... | ... |
@@ -1180,64 +1180,69 @@ |
1180 | 1180 |
|
1181 | 1181 |
/// \brief The type of the digraph the algorithm runs on. |
1182 | 1182 |
typedef _Digraph Digraph; |
1183 | 1183 |
|
1184 | 1184 |
/// \brief The type of the map that indicates which nodes are reached. |
1185 | 1185 |
/// |
1186 | 1186 |
/// The type of the map that indicates which nodes are reached. |
1187 | 1187 |
/// It must meet the \ref concepts::ReadWriteMap "ReadWriteMap" concept. |
1188 | 1188 |
typedef typename Digraph::template NodeMap<bool> ReachedMap; |
1189 | 1189 |
|
1190 | 1190 |
/// \brief Instantiates a \ref ReachedMap. |
1191 | 1191 |
/// |
1192 | 1192 |
/// This function instantiates a \ref ReachedMap. |
1193 | 1193 |
/// \param digraph is the digraph, to which |
1194 | 1194 |
/// we would like to define the \ref ReachedMap. |
1195 | 1195 |
static ReachedMap *createReachedMap(const Digraph &digraph) { |
1196 | 1196 |
return new ReachedMap(digraph); |
1197 | 1197 |
} |
1198 | 1198 |
|
1199 | 1199 |
}; |
1200 | 1200 |
|
1201 | 1201 |
/// \ingroup search |
1202 | 1202 |
/// |
1203 | 1203 |
/// \brief %DFS algorithm class with visitor interface. |
1204 | 1204 |
/// |
1205 | 1205 |
/// This class provides an efficient implementation of the %DFS algorithm |
1206 | 1206 |
/// with visitor interface. |
1207 | 1207 |
/// |
1208 | 1208 |
/// The %DfsVisit class provides an alternative interface to the Dfs |
1209 | 1209 |
/// class. It works with callback mechanism, the DfsVisit object calls |
1210 | 1210 |
/// the member functions of the \c Visitor class on every DFS event. |
1211 | 1211 |
/// |
1212 |
/// This interface of the DFS algorithm should be used in special cases |
|
1213 |
/// when extra actions have to be performed in connection with certain |
|
1214 |
/// events of the DFS algorithm. Otherwise consider to use Dfs or dfs() |
|
1215 |
/// instead. |
|
1216 |
/// |
|
1212 | 1217 |
/// \tparam _Digraph The type of the digraph the algorithm runs on. |
1213 | 1218 |
/// The default value is |
1214 | 1219 |
/// \ref ListDigraph. The value of _Digraph is not used directly by |
1215 | 1220 |
/// \ref DfsVisit, it is only passed to \ref DfsVisitDefaultTraits. |
1216 | 1221 |
/// \tparam _Visitor The Visitor type that is used by the algorithm. |
1217 | 1222 |
/// \ref DfsVisitor "DfsVisitor<_Digraph>" is an empty visitor, which |
1218 | 1223 |
/// does not observe the DFS events. If you want to observe the DFS |
1219 | 1224 |
/// events, you should implement your own visitor class. |
1220 | 1225 |
/// \tparam _Traits Traits class to set various data types used by the |
1221 | 1226 |
/// algorithm. The default traits class is |
1222 | 1227 |
/// \ref DfsVisitDefaultTraits "DfsVisitDefaultTraits<_Digraph>". |
1223 | 1228 |
/// See \ref DfsVisitDefaultTraits for the documentation of |
1224 | 1229 |
/// a DFS visit traits class. |
1225 | 1230 |
#ifdef DOXYGEN |
1226 | 1231 |
template <typename _Digraph, typename _Visitor, typename _Traits> |
1227 | 1232 |
#else |
1228 | 1233 |
template <typename _Digraph = ListDigraph, |
1229 | 1234 |
typename _Visitor = DfsVisitor<_Digraph>, |
1230 | 1235 |
typename _Traits = DfsDefaultTraits<_Digraph> > |
1231 | 1236 |
#endif |
1232 | 1237 |
class DfsVisit { |
1233 | 1238 |
public: |
1234 | 1239 |
|
1235 | 1240 |
/// \brief \ref Exception for uninitialized parameters. |
1236 | 1241 |
/// |
1237 | 1242 |
/// This error represents problems in the initialization |
1238 | 1243 |
/// of the parameters of the algorithm. |
1239 | 1244 |
class UninitializedParameter : public lemon::UninitializedParameter { |
1240 | 1245 |
public: |
1241 | 1246 |
virtual const char* what() const throw() |
1242 | 1247 |
{ |
1243 | 1248 |
return "lemon::DfsVisit::UninitializedParameter"; |
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-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_DIM2_H |
20 | 20 |
#define LEMON_DIM2_H |
21 | 21 |
|
22 | 22 |
#include <iostream> |
23 | 23 |
|
24 | 24 |
///\ingroup misc |
25 | 25 |
///\file |
26 | 26 |
///\brief A simple two dimensional vector and a bounding box implementation |
27 | 27 |
/// |
28 | 28 |
/// The class \ref lemon::dim2::Point "dim2::Point" implements |
29 | 29 |
/// a two dimensional vector with the usual operations. |
30 | 30 |
/// |
31 |
/// The class \ref lemon::dim2::BoundingBox "dim2::BoundingBox" |
|
32 |
/// can be used to determine |
|
31 |
/// The class \ref lemon::dim2::Box "dim2::Box" can be used to determine |
|
33 | 32 |
/// the rectangular bounding box of a set of |
34 | 33 |
/// \ref lemon::dim2::Point "dim2::Point"'s. |
35 | 34 |
|
36 | 35 |
namespace lemon { |
37 | 36 |
|
38 | 37 |
///Tools for handling two dimensional coordinates |
39 | 38 |
|
40 | 39 |
///This namespace is a storage of several |
41 | 40 |
///tools for handling two dimensional coordinates |
42 | 41 |
namespace dim2 { |
43 | 42 |
|
44 | 43 |
/// \addtogroup misc |
45 | 44 |
/// @{ |
46 | 45 |
|
47 |
/// |
|
46 |
/// Two dimensional vector (plain vector) |
|
48 | 47 |
|
49 | 48 |
/// A simple two dimensional vector (plain vector) implementation |
50 | 49 |
/// with the usual vector operations. |
51 | 50 |
template<typename T> |
52 | 51 |
class Point { |
53 | 52 |
|
54 | 53 |
public: |
55 | 54 |
|
56 | 55 |
typedef T Value; |
57 | 56 |
|
58 | 57 |
///First coordinate |
59 | 58 |
T x; |
60 | 59 |
///Second coordinate |
61 | 60 |
T y; |
62 | 61 |
|
63 | 62 |
///Default constructor |
64 | 63 |
Point() {} |
65 | 64 |
|
66 | 65 |
///Construct an instance from coordinates |
67 | 66 |
Point(T a, T b) : x(a), y(b) { } |
68 | 67 |
|
69 | 68 |
///Returns the dimension of the vector (i.e. returns 2). |
70 | 69 |
|
71 | 70 |
///The dimension of the vector. |
72 | 71 |
///This function always returns 2. |
73 | 72 |
int size() const { return 2; } |
74 | 73 |
|
75 | 74 |
///Subscripting operator |
76 | 75 |
|
77 | 76 |
///\c p[0] is \c p.x and \c p[1] is \c p.y |
78 | 77 |
/// |
79 | 78 |
T& operator[](int idx) { return idx == 0 ? x : y; } |
... | ... |
@@ -192,376 +191,417 @@ |
192 | 191 |
/// |
193 | 192 |
template<typename T> |
194 | 193 |
inline std::istream& operator>>(std::istream &is, Point<T> &z) { |
195 | 194 |
char c; |
196 | 195 |
if (is >> c) { |
197 | 196 |
if (c != '(') is.putback(c); |
198 | 197 |
} else { |
199 | 198 |
is.clear(); |
200 | 199 |
} |
201 | 200 |
if (!(is >> z.x)) return is; |
202 | 201 |
if (is >> c) { |
203 | 202 |
if (c != ',') is.putback(c); |
204 | 203 |
} else { |
205 | 204 |
is.clear(); |
206 | 205 |
} |
207 | 206 |
if (!(is >> z.y)) return is; |
208 | 207 |
if (is >> c) { |
209 | 208 |
if (c != ')') is.putback(c); |
210 | 209 |
} else { |
211 | 210 |
is.clear(); |
212 | 211 |
} |
213 | 212 |
return is; |
214 | 213 |
} |
215 | 214 |
|
216 | 215 |
///Write a plain vector to a stream |
217 | 216 |
|
218 | 217 |
///Write a plain vector to a stream. |
219 | 218 |
///\relates Point |
220 | 219 |
/// |
221 | 220 |
template<typename T> |
222 | 221 |
inline std::ostream& operator<<(std::ostream &os, const Point<T>& z) |
223 | 222 |
{ |
224 |
os << "(" << z.x << ", |
|
223 |
os << "(" << z.x << "," << z.y << ")"; |
|
225 | 224 |
return os; |
226 | 225 |
} |
227 | 226 |
|
228 | 227 |
///Rotate by 90 degrees |
229 | 228 |
|
230 | 229 |
///Returns the parameter rotated by 90 degrees in positive direction. |
231 | 230 |
///\relates Point |
232 | 231 |
/// |
233 | 232 |
template<typename T> |
234 | 233 |
inline Point<T> rot90(const Point<T> &z) |
235 | 234 |
{ |
236 | 235 |
return Point<T>(-z.y,z.x); |
237 | 236 |
} |
238 | 237 |
|
239 | 238 |
///Rotate by 180 degrees |
240 | 239 |
|
241 | 240 |
///Returns the parameter rotated by 180 degrees. |
242 | 241 |
///\relates Point |
243 | 242 |
/// |
244 | 243 |
template<typename T> |
245 | 244 |
inline Point<T> rot180(const Point<T> &z) |
246 | 245 |
{ |
247 | 246 |
return Point<T>(-z.x,-z.y); |
248 | 247 |
} |
249 | 248 |
|
250 | 249 |
///Rotate by 270 degrees |
251 | 250 |
|
252 | 251 |
///Returns the parameter rotated by 90 degrees in negative direction. |
253 | 252 |
///\relates Point |
254 | 253 |
/// |
255 | 254 |
template<typename T> |
256 | 255 |
inline Point<T> rot270(const Point<T> &z) |
257 | 256 |
{ |
258 | 257 |
return Point<T>(z.y,-z.x); |
259 | 258 |
} |
260 | 259 |
|
261 | 260 |
|
262 | 261 |
|
263 |
|
|
262 |
/// Bounding box of plain vectors (\ref Point points). |
|
264 | 263 |
|
265 |
/// A class to calculate or store the bounding box of plain vectors. |
|
266 |
/// |
|
267 |
template<typename T> |
|
268 |
class BoundingBox { |
|
264 |
/// A class to calculate or store the bounding box of plain vectors |
|
265 |
/// (\ref Point points). |
|
266 |
template<typename T> |
|
267 |
class Box { |
|
269 | 268 |
Point<T> _bottom_left, _top_right; |
270 | 269 |
bool _empty; |
271 | 270 |
public: |
272 | 271 |
|
273 |
///Default constructor: creates an empty bounding box |
|
274 |
BoundingBox() { _empty = true; } |
|
272 |
///Default constructor: creates an empty box |
|
273 |
Box() { _empty = true; } |
|
275 | 274 |
|
276 |
///Construct an instance from one point |
|
277 |
BoundingBox(Point<T> a) { |
|
275 |
///Construct a box from one point |
|
276 |
Box(Point<T> a) { |
|
278 | 277 |
_bottom_left = _top_right = a; |
279 | 278 |
_empty = false; |
280 | 279 |
} |
281 | 280 |
|
282 |
///Construct |
|
281 |
///Construct a box from two points |
|
283 | 282 |
|
284 |
///Construct |
|
283 |
///Construct a box from two points. |
|
285 | 284 |
///\param a The bottom left corner. |
286 | 285 |
///\param b The top right corner. |
287 | 286 |
///\warning The coordinates of the bottom left corner must be no more |
288 | 287 |
///than those of the top right one. |
289 |
|
|
288 |
Box(Point<T> a,Point<T> b) |
|
290 | 289 |
{ |
291 | 290 |
_bottom_left = a; |
292 | 291 |
_top_right = b; |
293 | 292 |
_empty = false; |
294 | 293 |
} |
295 | 294 |
|
296 |
///Construct |
|
295 |
///Construct a box from four numbers |
|
297 | 296 |
|
298 |
///Construct |
|
297 |
///Construct a box from four numbers. |
|
299 | 298 |
///\param l The left side of the box. |
300 | 299 |
///\param b The bottom of the box. |
301 | 300 |
///\param r The right side of the box. |
302 | 301 |
///\param t The top of the box. |
303 | 302 |
///\warning The left side must be no more than the right side and |
304 | 303 |
///bottom must be no more than the top. |
305 |
|
|
304 |
Box(T l,T b,T r,T t) |
|
306 | 305 |
{ |
307 | 306 |
_bottom_left=Point<T>(l,b); |
308 | 307 |
_top_right=Point<T>(r,t); |
309 | 308 |
_empty = false; |
310 | 309 |
} |
311 | 310 |
|
312 |
///Return \c true if the |
|
311 |
///Return \c true if the box is empty. |
|
313 | 312 |
|
314 |
///Return \c true if the |
|
313 |
///Return \c true if the box is empty (i.e. return \c false |
|
315 | 314 |
///if at least one point was added to the box or the coordinates of |
316 | 315 |
///the box were set). |
317 | 316 |
/// |
318 |
///The coordinates of an empty |
|
317 |
///The coordinates of an empty box are not defined. |
|
319 | 318 |
bool empty() const { |
320 | 319 |
return _empty; |
321 | 320 |
} |
322 | 321 |
|
323 |
///Make the |
|
322 |
///Make the box empty |
|
324 | 323 |
void clear() { |
325 | 324 |
_empty = true; |
326 | 325 |
} |
327 | 326 |
|
328 | 327 |
///Give back the bottom left corner of the box |
329 | 328 |
|
330 | 329 |
///Give back the bottom left corner of the box. |
331 |
///If the |
|
330 |
///If the box is empty, then the return value is not defined. |
|
332 | 331 |
Point<T> bottomLeft() const { |
333 | 332 |
return _bottom_left; |
334 | 333 |
} |
335 | 334 |
|
336 | 335 |
///Set the bottom left corner of the box |
337 | 336 |
|
338 | 337 |
///Set the bottom left corner of the box. |
339 | 338 |
///\pre The box must not be empty. |
340 | 339 |
void bottomLeft(Point<T> p) { |
341 | 340 |
_bottom_left = p; |
342 | 341 |
} |
343 | 342 |
|
344 | 343 |
///Give back the top right corner of the box |
345 | 344 |
|
346 | 345 |
///Give back the top right corner of the box. |
347 |
///If the |
|
346 |
///If the box is empty, then the return value is not defined. |
|
348 | 347 |
Point<T> topRight() const { |
349 | 348 |
return _top_right; |
350 | 349 |
} |
351 | 350 |
|
352 | 351 |
///Set the top right corner of the box |
353 | 352 |
|
354 | 353 |
///Set the top right corner of the box. |
355 | 354 |
///\pre The box must not be empty. |
356 | 355 |
void topRight(Point<T> p) { |
357 | 356 |
_top_right = p; |
358 | 357 |
} |
359 | 358 |
|
360 | 359 |
///Give back the bottom right corner of the box |
361 | 360 |
|
362 | 361 |
///Give back the bottom right corner of the box. |
363 |
///If the |
|
362 |
///If the box is empty, then the return value is not defined. |
|
364 | 363 |
Point<T> bottomRight() const { |
365 | 364 |
return Point<T>(_top_right.x,_bottom_left.y); |
366 | 365 |
} |
367 | 366 |
|
368 | 367 |
///Set the bottom right corner of the box |
369 | 368 |
|
370 | 369 |
///Set the bottom right corner of the box. |
371 | 370 |
///\pre The box must not be empty. |
372 | 371 |
void bottomRight(Point<T> p) { |
373 | 372 |
_top_right.x = p.x; |
374 | 373 |
_bottom_left.y = p.y; |
375 | 374 |
} |
376 | 375 |
|
377 | 376 |
///Give back the top left corner of the box |
378 | 377 |
|
379 | 378 |
///Give back the top left corner of the box. |
380 |
///If the |
|
379 |
///If the box is empty, then the return value is not defined. |
|
381 | 380 |
Point<T> topLeft() const { |
382 | 381 |
return Point<T>(_bottom_left.x,_top_right.y); |
383 | 382 |
} |
384 | 383 |
|
385 | 384 |
///Set the top left corner of the box |
386 | 385 |
|
387 | 386 |
///Set the top left corner of the box. |
388 | 387 |
///\pre The box must not be empty. |
389 | 388 |
void topLeft(Point<T> p) { |
390 | 389 |
_top_right.y = p.y; |
391 | 390 |
_bottom_left.x = p.x; |
392 | 391 |
} |
393 | 392 |
|
394 | 393 |
///Give back the bottom of the box |
395 | 394 |
|
396 | 395 |
///Give back the bottom of the box. |
397 |
///If the |
|
396 |
///If the box is empty, then the return value is not defined. |
|
398 | 397 |
T bottom() const { |
399 | 398 |
return _bottom_left.y; |
400 | 399 |
} |
401 | 400 |
|
402 | 401 |
///Set the bottom of the box |
403 | 402 |
|
404 | 403 |
///Set the bottom of the box. |
405 | 404 |
///\pre The box must not be empty. |
406 | 405 |
void bottom(T t) { |
407 | 406 |
_bottom_left.y = t; |
408 | 407 |
} |
409 | 408 |
|
410 | 409 |
///Give back the top of the box |
411 | 410 |
|
412 | 411 |
///Give back the top of the box. |
413 |
///If the |
|
412 |
///If the box is empty, then the return value is not defined. |
|
414 | 413 |
T top() const { |
415 | 414 |
return _top_right.y; |
416 | 415 |
} |
417 | 416 |
|
418 | 417 |
///Set the top of the box |
419 | 418 |
|
420 | 419 |
///Set the top of the box. |
421 | 420 |
///\pre The box must not be empty. |
422 | 421 |
void top(T t) { |
423 | 422 |
_top_right.y = t; |
424 | 423 |
} |
425 | 424 |
|
426 | 425 |
///Give back the left side of the box |
427 | 426 |
|
428 | 427 |
///Give back the left side of the box. |
429 |
///If the |
|
428 |
///If the box is empty, then the return value is not defined. |
|
430 | 429 |
T left() const { |
431 | 430 |
return _bottom_left.x; |
432 | 431 |
} |
433 | 432 |
|
434 | 433 |
///Set the left side of the box |
435 | 434 |
|
436 | 435 |
///Set the left side of the box. |
437 | 436 |
///\pre The box must not be empty. |
438 | 437 |
void left(T t) { |
439 | 438 |
_bottom_left.x = t; |
440 | 439 |
} |
441 | 440 |
|
442 | 441 |
/// Give back the right side of the box |
443 | 442 |
|
444 | 443 |
/// Give back the right side of the box. |
445 |
///If the |
|
444 |
///If the box is empty, then the return value is not defined. |
|
446 | 445 |
T right() const { |
447 | 446 |
return _top_right.x; |
448 | 447 |
} |
449 | 448 |
|
450 | 449 |
///Set the right side of the box |
451 | 450 |
|
452 | 451 |
///Set the right side of the box. |
453 | 452 |
///\pre The box must not be empty. |
454 | 453 |
void right(T t) { |
455 | 454 |
_top_right.x = t; |
456 | 455 |
} |
457 | 456 |
|
458 | 457 |
///Give back the height of the box |
459 | 458 |
|
460 | 459 |
///Give back the height of the box. |
461 |
///If the |
|
460 |
///If the box is empty, then the return value is not defined. |
|
462 | 461 |
T height() const { |
463 | 462 |
return _top_right.y-_bottom_left.y; |
464 | 463 |
} |
465 | 464 |
|
466 | 465 |
///Give back the width of the box |
467 | 466 |
|
468 | 467 |
///Give back the width of the box. |
469 |
///If the |
|
468 |
///If the box is empty, then the return value is not defined. |
|
470 | 469 |
T width() const { |
471 | 470 |
return _top_right.x-_bottom_left.x; |
472 | 471 |
} |
473 | 472 |
|
474 |
///Checks whether a point is inside |
|
473 |
///Checks whether a point is inside the box |
|
475 | 474 |
bool inside(const Point<T>& u) const { |
476 | 475 |
if (_empty) |
477 | 476 |
return false; |
478 | 477 |
else { |
479 | 478 |
return ( (u.x-_bottom_left.x)*(_top_right.x-u.x) >= 0 && |
480 | 479 |
(u.y-_bottom_left.y)*(_top_right.y-u.y) >= 0 ); |
481 | 480 |
} |
482 | 481 |
} |
483 | 482 |
|
484 |
///Increments |
|
483 |
///Increments the box with a point |
|
485 | 484 |
|
486 |
///Increments |
|
485 |
///Increments the box with a point. |
|
487 | 486 |
/// |
488 |
|
|
487 |
Box& add(const Point<T>& u){ |
|
489 | 488 |
if (_empty) { |
490 | 489 |
_bottom_left = _top_right = u; |
491 | 490 |
_empty = false; |
492 | 491 |
} |
493 | 492 |
else { |
494 | 493 |
if (_bottom_left.x > u.x) _bottom_left.x = u.x; |
495 | 494 |
if (_bottom_left.y > u.y) _bottom_left.y = u.y; |
496 | 495 |
if (_top_right.x < u.x) _top_right.x = u.x; |
497 | 496 |
if (_top_right.y < u.y) _top_right.y = u.y; |
498 | 497 |
} |
499 | 498 |
return *this; |
500 | 499 |
} |
501 | 500 |
|
502 |
///Increments |
|
501 |
///Increments the box to contain another box |
|
503 | 502 |
|
504 |
///Increments |
|
503 |
///Increments the box to contain another box. |
|
505 | 504 |
/// |
506 |
|
|
505 |
Box& add(const Box &u){ |
|
507 | 506 |
if ( !u.empty() ){ |
508 | 507 |
add(u._bottom_left); |
509 | 508 |
add(u._top_right); |
510 | 509 |
} |
511 | 510 |
return *this; |
512 | 511 |
} |
513 | 512 |
|
514 |
///Intersection of two |
|
513 |
///Intersection of two boxes |
|
515 | 514 |
|
516 |
///Intersection of two |
|
515 |
///Intersection of two boxes. |
|
517 | 516 |
/// |
518 |
BoundingBox operator&(const BoundingBox& u) const { |
|
519 |
BoundingBox b; |
|
517 |
Box operator&(const Box& u) const { |
|
518 |
Box b; |
|
520 | 519 |
if (_empty || u._empty) { |
521 | 520 |
b._empty = true; |
522 | 521 |
} else { |
523 | 522 |
b._bottom_left.x = std::max(_bottom_left.x, u._bottom_left.x); |
524 | 523 |
b._bottom_left.y = std::max(_bottom_left.y, u._bottom_left.y); |
525 | 524 |
b._top_right.x = std::min(_top_right.x, u._top_right.x); |
526 | 525 |
b._top_right.y = std::min(_top_right.y, u._top_right.y); |
527 | 526 |
b._empty = b._bottom_left.x > b._top_right.x || |
528 | 527 |
b._bottom_left.y > b._top_right.y; |
529 | 528 |
} |
530 | 529 |
return b; |
531 | 530 |
} |
532 | 531 |
|
533 |
|
|
532 |
};//class Box |
|
534 | 533 |
|
535 | 534 |
|
535 |
///Read a box from a stream |
|
536 |
|
|
537 |
///Read a box from a stream. |
|
538 |
///\relates Box |
|
539 |
template<typename T> |
|
540 |
inline std::istream& operator>>(std::istream &is, Box<T>& b) { |
|
541 |
char c; |
|
542 |
Point<T> p; |
|
543 |
if (is >> c) { |
|
544 |
if (c != '(') is.putback(c); |
|
545 |
} else { |
|
546 |
is.clear(); |
|
547 |
} |
|
548 |
if (!(is >> p)) return is; |
|
549 |
b.bottomLeft(p); |
|
550 |
if (is >> c) { |
|
551 |
if (c != ',') is.putback(c); |
|
552 |
} else { |
|
553 |
is.clear(); |
|
554 |
} |
|
555 |
if (!(is >> p)) return is; |
|
556 |
b.topRight(p); |
|
557 |
if (is >> c) { |
|
558 |
if (c != ')') is.putback(c); |
|
559 |
} else { |
|
560 |
is.clear(); |
|
561 |
} |
|
562 |
return is; |
|
563 |
} |
|
564 |
|
|
565 |
///Write a box to a stream |
|
566 |
|
|
567 |
///Write a box to a stream. |
|
568 |
///\relates Box |
|
569 |
template<typename T> |
|
570 |
inline std::ostream& operator<<(std::ostream &os, const Box<T>& b) |
|
571 |
{ |
|
572 |
os << "(" << b.bottomLeft() << "," << b.topRight() << ")"; |
|
573 |
return os; |
|
574 |
} |
|
575 |
|
|
536 | 576 |
///Map of x-coordinates of a \ref Point "Point"-map |
537 | 577 |
|
538 | 578 |
///\ingroup maps |
539 | 579 |
///Map of x-coordinates of a \ref Point "Point"-map. |
540 | 580 |
/// |
541 | 581 |
template<class M> |
542 | 582 |
class XMap |
543 | 583 |
{ |
544 | 584 |
M& _map; |
545 | 585 |
public: |
546 | 586 |
|
547 | 587 |
typedef typename M::Value::Value Value; |
548 | 588 |
typedef typename M::Key Key; |
549 | 589 |
///\e |
550 | 590 |
XMap(M& map) : _map(map) {} |
551 | 591 |
Value operator[](Key k) const {return _map[k].x;} |
552 | 592 |
void set(Key k,Value v) {_map.set(k,typename M::Value(v,_map[k].y));} |
553 | 593 |
}; |
554 | 594 |
|
555 | 595 |
///Returns an \ref XMap class |
556 | 596 |
|
557 | 597 |
///This function just returns an \ref XMap class. |
558 | 598 |
/// |
559 | 599 |
///\ingroup maps |
560 | 600 |
///\relates XMap |
561 | 601 |
template<class M> |
562 | 602 |
inline XMap<M> xMap(M &m) |
563 | 603 |
{ |
564 | 604 |
return XMap<M>(m); |
565 | 605 |
} |
566 | 606 |
|
567 | 607 |
template<class M> |
... | ... |
@@ -696,98 +696,98 @@ |
696 | 696 |
GetTimeFormat(LOCALE_USER_DEFAULT, 0, &time, |
697 | 697 |
"HH':'mm':'ss", buf2, 9) && |
698 | 698 |
GetDateFormat(LOCALE_USER_DEFAULT, 0, &time, |
699 | 699 |
"yyyy", buf3, 5)) { |
700 | 700 |
os << "%%CreationDate: " << buf1 << ' ' |
701 | 701 |
<< buf2 << ' ' << buf3 << std::endl; |
702 | 702 |
} |
703 | 703 |
#endif |
704 | 704 |
} |
705 | 705 |
|
706 | 706 |
if (_autoArcWidthScale) { |
707 | 707 |
double max_w=0; |
708 | 708 |
for(ArcIt e(g);e!=INVALID;++e) |
709 | 709 |
max_w=std::max(double(_arcWidths[e]),max_w); |
710 | 710 |
//\todo better 'epsilon' would be nice here. |
711 | 711 |
if(max_w>EPSILON) { |
712 | 712 |
_arcWidthScale/=max_w; |
713 | 713 |
} |
714 | 714 |
} |
715 | 715 |
|
716 | 716 |
if (_autoNodeScale) { |
717 | 717 |
double max_s=0; |
718 | 718 |
for(NodeIt n(g);n!=INVALID;++n) |
719 | 719 |
max_s=std::max(double(_nodeSizes[n]),max_s); |
720 | 720 |
//\todo better 'epsilon' would be nice here. |
721 | 721 |
if(max_s>EPSILON) { |
722 | 722 |
_nodeScale/=max_s; |
723 | 723 |
} |
724 | 724 |
} |
725 | 725 |
|
726 | 726 |
double diag_len = 1; |
727 | 727 |
if(!(_absoluteNodeSizes&&_absoluteArcWidths)) { |
728 |
dim2:: |
|
728 |
dim2::Box<double> bb; |
|
729 | 729 |
for(NodeIt n(g);n!=INVALID;++n) bb.add(mycoords[n]); |
730 | 730 |
if (bb.empty()) { |
731 |
bb = dim2:: |
|
731 |
bb = dim2::Box<double>(dim2::Point<double>(0,0)); |
|
732 | 732 |
} |
733 | 733 |
diag_len = std::sqrt((bb.bottomLeft()-bb.topRight()).normSquare()); |
734 | 734 |
if(diag_len<EPSILON) diag_len = 1; |
735 | 735 |
if(!_absoluteNodeSizes) _nodeScale*=diag_len; |
736 | 736 |
if(!_absoluteArcWidths) _arcWidthScale*=diag_len; |
737 | 737 |
} |
738 | 738 |
|
739 |
dim2:: |
|
739 |
dim2::Box<double> bb; |
|
740 | 740 |
for(NodeIt n(g);n!=INVALID;++n) { |
741 | 741 |
double ns=_nodeSizes[n]*_nodeScale; |
742 | 742 |
dim2::Point<double> p(ns,ns); |
743 | 743 |
switch(_nodeShapes[n]) { |
744 | 744 |
case CIRCLE: |
745 | 745 |
case SQUARE: |
746 | 746 |
case DIAMOND: |
747 | 747 |
bb.add(p+mycoords[n]); |
748 | 748 |
bb.add(-p+mycoords[n]); |
749 | 749 |
break; |
750 | 750 |
case MALE: |
751 | 751 |
bb.add(-p+mycoords[n]); |
752 | 752 |
bb.add(dim2::Point<double>(1.5*ns,1.5*std::sqrt(3.0)*ns)+mycoords[n]); |
753 | 753 |
break; |
754 | 754 |
case FEMALE: |
755 | 755 |
bb.add(p+mycoords[n]); |
756 | 756 |
bb.add(dim2::Point<double>(-ns,-3.01*ns)+mycoords[n]); |
757 | 757 |
break; |
758 | 758 |
} |
759 | 759 |
} |
760 | 760 |
if (bb.empty()) { |
761 |
bb = dim2:: |
|
761 |
bb = dim2::Box<double>(dim2::Point<double>(0,0)); |
|
762 | 762 |
} |
763 | 763 |
|
764 | 764 |
if(_scaleToA4) |
765 | 765 |
os <<"%%BoundingBox: 0 0 596 842\n%%DocumentPaperSizes: a4\n"; |
766 | 766 |
else { |
767 | 767 |
if(_preScale) { |
768 | 768 |
//Rescale so that BoundingBox won't be neither to big nor too small. |
769 | 769 |
while(bb.height()*_scale>1000||bb.width()*_scale>1000) _scale/=10; |
770 | 770 |
while(bb.height()*_scale<100||bb.width()*_scale<100) _scale*=10; |
771 | 771 |
} |
772 | 772 |
|
773 | 773 |
os << "%%BoundingBox: " |
774 | 774 |
<< int(floor(bb.left() * _scale - _xBorder)) << ' ' |
775 | 775 |
<< int(floor(bb.bottom() * _scale - _yBorder)) << ' ' |
776 | 776 |
<< int(ceil(bb.right() * _scale + _xBorder)) << ' ' |
777 | 777 |
<< int(ceil(bb.top() * _scale + _yBorder)) << '\n'; |
778 | 778 |
} |
779 | 779 |
|
780 | 780 |
os << "%%EndComments\n"; |
781 | 781 |
|
782 | 782 |
//x1 y1 x2 y2 x3 y3 cr cg cb w |
783 | 783 |
os << "/lb { setlinewidth setrgbcolor newpath moveto\n" |
784 | 784 |
<< " 4 2 roll 1 index 1 index curveto stroke } bind def\n"; |
785 | 785 |
os << "/l { setlinewidth setrgbcolor newpath moveto lineto stroke }" |
786 | 786 |
<< " bind def\n"; |
787 | 787 |
//x y r |
788 | 788 |
os << "/c { newpath dup 3 index add 2 index moveto 0 360 arc closepath }" |
789 | 789 |
<< " bind def\n"; |
790 | 790 |
//x y r |
791 | 791 |
os << "/sq { newpath 2 index 1 index add 2 index 2 index add moveto\n" |
792 | 792 |
<< " 2 index 1 index sub 2 index 2 index add lineto\n" |
793 | 793 |
<< " 2 index 1 index sub 2 index 2 index sub lineto\n" |
... | ... |
@@ -21,67 +21,67 @@ |
21 | 21 |
#include "test_tools.h" |
22 | 22 |
|
23 | 23 |
using namespace std; |
24 | 24 |
using namespace lemon; |
25 | 25 |
|
26 | 26 |
int main() |
27 | 27 |
{ |
28 | 28 |
typedef dim2::Point<int> Point; |
29 | 29 |
|
30 | 30 |
Point p; |
31 | 31 |
check(p.size()==2, "Wrong dim2::Point initialization."); |
32 | 32 |
|
33 | 33 |
Point a(1,2); |
34 | 34 |
Point b(3,4); |
35 | 35 |
check(a[0]==1 && a[1]==2, "Wrong dim2::Point initialization."); |
36 | 36 |
|
37 | 37 |
p = a+b; |
38 | 38 |
check(p.x==4 && p.y==6, "Wrong dim2::Point addition."); |
39 | 39 |
|
40 | 40 |
p = a-b; |
41 | 41 |
check(p.x==-2 && p.y==-2, "Wrong dim2::Point subtraction."); |
42 | 42 |
|
43 | 43 |
check(a.normSquare()==5,"Wrong dim2::Point norm calculation."); |
44 | 44 |
check(a*b==11, "Wrong dim2::Point scalar product."); |
45 | 45 |
|
46 | 46 |
int l=2; |
47 | 47 |
p = a*l; |
48 | 48 |
check(p.x==2 && p.y==4, "Wrong dim2::Point multiplication by a scalar."); |
49 | 49 |
|
50 | 50 |
p = b/l; |
51 | 51 |
check(p.x==1 && p.y==2, "Wrong dim2::Point division by a scalar."); |
52 | 52 |
|
53 |
typedef dim2::BoundingBox<int> BB; |
|
54 |
BB box1; |
|
55 |
|
|
53 |
typedef dim2::Box<int> Box; |
|
54 |
Box box1; |
|
55 |
check(box1.empty(), "Wrong empty() in dim2::Box."); |
|
56 | 56 |
|
57 | 57 |
box1.add(a); |
58 |
check(!box1.empty(), "Wrong empty() in dim2:: |
|
58 |
check(!box1.empty(), "Wrong empty() in dim2::Box."); |
|
59 | 59 |
box1.add(b); |
60 | 60 |
|
61 | 61 |
check(box1.left()==1 && box1.bottom()==2 && |
62 | 62 |
box1.right()==3 && box1.top()==4, |
63 |
"Wrong addition of points to dim2:: |
|
63 |
"Wrong addition of points to dim2::Box."); |
|
64 | 64 |
|
65 |
check(box1.inside(Point(2,3)), "Wrong inside() in dim2::BoundingBox."); |
|
66 |
check(box1.inside(Point(1,3)), "Wrong inside() in dim2::BoundingBox."); |
|
67 |
check( |
|
65 |
check(box1.inside(Point(2,3)), "Wrong inside() in dim2::Box."); |
|
66 |
check(box1.inside(Point(1,3)), "Wrong inside() in dim2::Box."); |
|
67 |
check(!box1.inside(Point(0,3)), "Wrong inside() in dim2::Box."); |
|
68 | 68 |
|
69 |
BB box2(Point(2,2)); |
|
70 |
check(!box2.empty(), "Wrong empty() in dim2::BoundingBox."); |
|
71 |
|
|
69 |
Box box2(Point(2,2)); |
|
70 |
check(!box2.empty(), "Wrong empty() in dim2::Box."); |
|
71 |
|
|
72 | 72 |
box2.bottomLeft(Point(2,0)); |
73 | 73 |
box2.topRight(Point(5,3)); |
74 |
|
|
74 |
Box box3 = box1 & box2; |
|
75 | 75 |
check(!box3.empty() && |
76 |
box3.left()==2 && box3.bottom()==2 && |
|
76 |
box3.left()==2 && box3.bottom()==2 && |
|
77 | 77 |
box3.right()==3 && box3.top()==3, |
78 |
"Wrong intersection of two dim2::BoundingBox objects."); |
|
79 |
|
|
78 |
"Wrong intersection of two dim2::Box objects."); |
|
79 |
|
|
80 | 80 |
box1.add(box2); |
81 | 81 |
check(!box1.empty() && |
82 | 82 |
box1.left()==1 && box1.bottom()==0 && |
83 | 83 |
box1.right()==5 && box1.top()==4, |
84 |
"Wrong addition of two dim2:: |
|
84 |
"Wrong addition of two dim2::Box objects."); |
|
85 | 85 |
|
86 | 86 |
return 0; |
87 | 87 |
} |
0 comments (0 inline)