0
17
0
1
1
1
1
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
///\file |
20 | 20 |
///\brief Some basic non inline function and static global data. |
21 | 21 |
|
22 | 22 |
#include<lemon/tolerance.h> |
23 | 23 |
#include<lemon/bits/invalid.h> |
24 | 24 |
namespace lemon { |
25 | 25 |
|
26 | 26 |
float Tolerance<float>::def_epsilon = 1e-4; |
27 | 27 |
double Tolerance<double>::def_epsilon = 1e-10; |
28 | 28 |
long double Tolerance<long double>::def_epsilon = 1e-14; |
29 | 29 |
|
30 | 30 |
#ifndef LEMON_ONLY_TEMPLATES |
31 | 31 |
const Invalid INVALID = Invalid(); |
32 | 32 |
#endif |
33 | 33 |
|
34 | 34 |
} //namespace lemon |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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_BITS_INVALID_H |
20 | 20 |
#define LEMON_BITS_INVALID_H |
21 | 21 |
|
22 | 22 |
///\file |
23 | 23 |
///\brief Definition of INVALID. |
24 | 24 |
|
25 | 25 |
namespace lemon { |
26 | 26 |
|
27 | 27 |
/// \brief Dummy type to make it easier to create invalid iterators. |
28 | 28 |
/// |
29 | 29 |
/// See \ref INVALID for the usage. |
30 | 30 |
struct Invalid { |
31 | 31 |
public: |
32 | 32 |
bool operator==(Invalid) { return true; } |
33 | 33 |
bool operator!=(Invalid) { return false; } |
34 | 34 |
bool operator< (Invalid) { return false; } |
35 | 35 |
}; |
36 | 36 |
|
37 | 37 |
/// \brief Invalid iterators. |
38 | 38 |
/// |
39 | 39 |
/// \ref Invalid is a global type that converts to each iterator |
40 | 40 |
/// in such a way that the value of the target iterator will be invalid. |
41 | 41 |
|
42 | 42 |
//Some people didn't like this: |
43 | 43 |
//const Invalid &INVALID = *(Invalid *)0; |
44 | 44 |
|
45 | 45 |
#ifdef LEMON_ONLY_TEMPLATES |
46 | 46 |
const Invalid INVALID = Invalid(); |
47 | 47 |
#else |
48 | 48 |
extern const Invalid INVALID; |
49 | 49 |
#endif |
50 | 50 |
|
51 | 51 |
} //namespace lemon |
52 | 52 |
|
53 | 53 |
#endif |
54 | 54 |
|
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
// This file contains a modified version of the enable_if library from BOOST. |
20 | 20 |
// See the appropriate copyright notice below. |
21 | 21 |
|
22 | 22 |
// Boost enable_if library |
23 | 23 |
|
24 | 24 |
// Copyright 2003 � The Trustees of Indiana University. |
25 | 25 |
|
26 | 26 |
// Use, modification, and distribution is subject to the Boost Software |
27 | 27 |
// License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at |
28 | 28 |
// http://www.boost.org/LICENSE_1_0.txt) |
29 | 29 |
|
30 | 30 |
// Authors: Jaakko J�rvi (jajarvi at osl.iu.edu) |
31 | 31 |
// Jeremiah Willcock (jewillco at osl.iu.edu) |
32 | 32 |
// Andrew Lumsdaine (lums at osl.iu.edu) |
33 | 33 |
|
34 | 34 |
|
35 | 35 |
#ifndef LEMON_BITS_UTILITY_H |
36 | 36 |
#define LEMON_BITS_UTILITY_H |
37 | 37 |
|
38 | 38 |
///\file |
39 | 39 |
///\brief Miscellaneous basic utilities |
40 | 40 |
/// |
41 | 41 |
///\todo Please rethink the organisation of the basic files like this. |
42 | 42 |
///E.g. this file might be merged with invalid.h. |
43 | 43 |
|
44 | 44 |
|
45 | 45 |
namespace lemon |
46 | 46 |
{ |
47 | 47 |
|
48 | 48 |
/// Basic type for defining "tags". A "YES" condition for \c enable_if. |
49 | 49 |
|
50 | 50 |
/// Basic type for defining "tags". A "YES" condition for \c enable_if. |
51 | 51 |
/// |
52 | 52 |
///\sa False |
53 | 53 |
/// |
54 | 54 |
/// \todo This should go to a separate "basic_types.h" (or something) |
55 | 55 |
/// file. |
56 | 56 |
struct True { |
57 | 57 |
///\e |
58 | 58 |
static const bool value = true; |
59 | 59 |
}; |
60 | 60 |
|
61 | 61 |
/// Basic type for defining "tags". A "NO" condition for \c enable_if. |
62 | 62 |
|
63 | 63 |
/// Basic type for defining "tags". A "NO" condition for \c enable_if. |
64 | 64 |
/// |
65 | 65 |
///\sa True |
66 | 66 |
struct False { |
67 | 67 |
///\e |
68 | 68 |
static const bool value = false; |
69 | 69 |
}; |
70 | 70 |
|
71 | 71 |
|
72 | 72 |
struct InvalidType { |
73 | 73 |
}; |
74 | 74 |
|
75 | 75 |
template <typename T> |
76 | 76 |
struct Wrap { |
77 | 77 |
const T &value; |
78 | 78 |
Wrap(const T &t) : value(t) {} |
79 | 79 |
}; |
80 | 80 |
|
81 | 81 |
/**************** dummy class to avoid ambiguity ****************/ |
82 | 82 |
|
83 | 83 |
template<int T> struct dummy { dummy(int) {} }; |
84 | 84 |
|
85 | 85 |
/**************** enable_if from BOOST ****************/ |
86 | 86 |
|
87 | 87 |
template <typename Type, typename T = void> |
88 | 88 |
struct exists { |
89 | 89 |
typedef T type; |
90 | 90 |
}; |
91 | 91 |
|
92 | 92 |
|
93 | 93 |
template <bool B, class T = void> |
94 | 94 |
struct enable_if_c { |
95 | 95 |
typedef T type; |
96 | 96 |
}; |
97 | 97 |
|
98 | 98 |
template <class T> |
99 | 99 |
struct enable_if_c<false, T> {}; |
100 | 100 |
|
101 | 101 |
template <class Cond, class T = void> |
102 | 102 |
struct enable_if : public enable_if_c<Cond::value, T> {}; |
103 | 103 |
|
104 | 104 |
template <bool B, class T> |
105 | 105 |
struct lazy_enable_if_c { |
106 | 106 |
typedef typename T::type type; |
107 | 107 |
}; |
108 | 108 |
|
109 | 109 |
template <class T> |
110 | 110 |
struct lazy_enable_if_c<false, T> {}; |
111 | 111 |
|
112 | 112 |
template <class Cond, class T> |
113 | 113 |
struct lazy_enable_if : public lazy_enable_if_c<Cond::value, T> {}; |
114 | 114 |
|
115 | 115 |
|
116 | 116 |
template <bool B, class T = void> |
117 | 117 |
struct disable_if_c { |
118 | 118 |
typedef T type; |
119 | 119 |
}; |
120 | 120 |
|
121 | 121 |
template <class T> |
122 | 122 |
struct disable_if_c<true, T> {}; |
123 | 123 |
|
124 | 124 |
template <class Cond, class T = void> |
125 | 125 |
struct disable_if : public disable_if_c<Cond::value, T> {}; |
126 | 126 |
|
127 | 127 |
template <bool B, class T> |
128 | 128 |
struct lazy_disable_if_c { |
129 | 129 |
typedef typename T::type type; |
130 | 130 |
}; |
131 | 131 |
|
132 | 132 |
template <class T> |
133 | 133 |
struct lazy_disable_if_c<true, T> {}; |
134 | 134 |
|
135 | 135 |
template <class Cond, class T> |
136 | 136 |
struct lazy_disable_if : public lazy_disable_if_c<Cond::value, T> {}; |
137 | 137 |
|
138 | 138 |
} // namespace lemon |
139 | 139 |
|
140 | 140 |
#endif |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
// This file contains a modified version of the concept checking |
20 | 20 |
// utility from BOOST. |
21 | 21 |
// See the appropriate copyright notice below. |
22 | 22 |
|
23 | 23 |
// (C) Copyright Jeremy Siek 2000. |
24 | 24 |
// Distributed under the Boost Software License, Version 1.0. (See |
25 | 25 |
// accompanying file LICENSE_1_0.txt or copy at |
26 | 26 |
// http://www.boost.org/LICENSE_1_0.txt) |
27 | 27 |
// |
28 | 28 |
// Revision History: |
29 | 29 |
// 05 May 2001: Workarounds for HP aCC from Thomas Matelich. (Jeremy Siek) |
30 | 30 |
// 02 April 2001: Removed limits header altogether. (Jeremy Siek) |
31 | 31 |
// 01 April 2001: Modified to use new <boost/limits.hpp> header. (JMaddock) |
32 | 32 |
// |
33 | 33 |
|
34 | 34 |
// See http://www.boost.org/libs/concept_check for documentation. |
35 | 35 |
|
36 | 36 |
#ifndef LEMON_CONCEPT_CHECKS_H |
37 | 37 |
#define LEMON_CONCEPT_CHECKS_H |
38 | 38 |
|
39 | 39 |
namespace lemon { |
40 | 40 |
|
41 | 41 |
/* |
42 | 42 |
"inline" is used for ignore_unused_variable_warning() |
43 | 43 |
and function_requires() to make sure there is no |
44 | 44 |
overtarget with g++. |
45 | 45 |
*/ |
46 | 46 |
|
47 | 47 |
template <class T> inline void ignore_unused_variable_warning(const T&) { } |
48 | 48 |
|
49 | 49 |
template <class Concept> |
50 | 50 |
inline void function_requires() |
51 | 51 |
{ |
52 | 52 |
#if !defined(NDEBUG) |
53 | 53 |
void (Concept::*x)() = & Concept::constraints; |
54 | 54 |
ignore_unused_variable_warning(x); |
55 | 55 |
#endif |
56 | 56 |
} |
57 | 57 |
|
58 | 58 |
template <typename Concept, typename Type> |
59 | 59 |
inline void checkConcept() { |
60 | 60 |
#if !defined(NDEBUG) |
61 | 61 |
typedef typename Concept::template Constraints<Type> ConceptCheck; |
62 | 62 |
void (ConceptCheck::*x)() = & ConceptCheck::constraints; |
63 | 63 |
ignore_unused_variable_warning(x); |
64 | 64 |
#endif |
65 | 65 |
} |
66 | 66 |
|
67 | 67 |
#define BOOST_CLASS_REQUIRE(type_var, ns, concept) \ |
68 | 68 |
typedef void (ns::concept <type_var>::* func##type_var##concept)(); \ |
69 | 69 |
template <func##type_var##concept Tp1_> \ |
70 | 70 |
struct concept_checking_##type_var##concept { }; \ |
71 | 71 |
typedef concept_checking_##type_var##concept< \ |
72 | 72 |
BOOST_FPTR ns::concept<type_var>::constraints> \ |
73 | 73 |
concept_checking_typedef_##type_var##concept |
74 | 74 |
|
75 | 75 |
#define BOOST_CLASS_REQUIRE2(type_var1, type_var2, ns, concept) \ |
76 | 76 |
typedef void (ns::concept <type_var1,type_var2>::* \ |
77 | 77 |
func##type_var1##type_var2##concept)(); \ |
78 | 78 |
template <func##type_var1##type_var2##concept Tp1_> \ |
79 | 79 |
struct concept_checking_##type_var1##type_var2##concept { }; \ |
80 | 80 |
typedef concept_checking_##type_var1##type_var2##concept< \ |
81 | 81 |
BOOST_FPTR ns::concept<type_var1,type_var2>::constraints> \ |
82 | 82 |
concept_checking_typedef_##type_var1##type_var2##concept |
83 | 83 |
|
84 | 84 |
#define BOOST_CLASS_REQUIRE3(tv1, tv2, tv3, ns, concept) \ |
85 | 85 |
typedef void (ns::concept <tv1,tv2,tv3>::* \ |
86 | 86 |
func##tv1##tv2##tv3##concept)(); \ |
87 | 87 |
template <func##tv1##tv2##tv3##concept Tp1_> \ |
88 | 88 |
struct concept_checking_##tv1##tv2##tv3##concept { }; \ |
89 | 89 |
typedef concept_checking_##tv1##tv2##tv3##concept< \ |
90 | 90 |
BOOST_FPTR ns::concept<tv1,tv2,tv3>::constraints> \ |
91 | 91 |
concept_checking_typedef_##tv1##tv2##tv3##concept |
92 | 92 |
|
93 | 93 |
#define BOOST_CLASS_REQUIRE4(tv1, tv2, tv3, tv4, ns, concept) \ |
94 | 94 |
typedef void (ns::concept <tv1,tv2,tv3,tv4>::* \ |
95 | 95 |
func##tv1##tv2##tv3##tv4##concept)(); \ |
96 | 96 |
template <func##tv1##tv2##tv3##tv4##concept Tp1_> \ |
97 | 97 |
struct concept_checking_##tv1##tv2##tv3##tv4##concept { }; \ |
98 | 98 |
typedef concept_checking_##tv1##tv2##tv3##tv4##concept< \ |
99 | 99 |
BOOST_FPTR ns::concept<tv1,tv2,tv3,tv4>::constraints> \ |
100 | 100 |
concept_checking_typedef_##tv1##tv2##tv3##tv4##concept |
101 | 101 |
|
102 | 102 |
|
103 | 103 |
} // namespace lemon |
104 | 104 |
|
105 | 105 |
#endif // LEMON_CONCEPT_CHECKS_H |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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_CONCEPT_MAPS_H |
20 | 20 |
#define LEMON_CONCEPT_MAPS_H |
21 | 21 |
|
22 | 22 |
#include <lemon/bits/utility.h> |
23 | 23 |
#include <lemon/concept_check.h> |
24 | 24 |
|
25 | 25 |
///\ingroup concept |
26 | 26 |
///\file |
27 | 27 |
///\brief Map concepts checking classes for testing and documenting. |
28 | 28 |
|
29 | 29 |
namespace lemon { |
30 | 30 |
|
31 | 31 |
namespace concepts { |
32 | 32 |
|
33 | 33 |
/// \addtogroup concept |
34 | 34 |
/// @{ |
35 | 35 |
|
36 | 36 |
/// Readable map concept |
37 | 37 |
|
38 | 38 |
/// Readable map concept. |
39 | 39 |
/// |
40 | 40 |
template<typename K, typename T> |
41 | 41 |
class ReadMap |
42 | 42 |
{ |
43 | 43 |
public: |
44 | 44 |
/// The key type of the map. |
45 | 45 |
typedef K Key; |
46 | 46 |
/// The value type of the map. (The type of objects associated with the keys). |
47 | 47 |
typedef T Value; |
48 | 48 |
|
49 | 49 |
/// Returns the value associated with a key. |
50 | 50 |
|
51 | 51 |
/// \bug Value shouldn't need to be default constructible. |
52 | 52 |
/// |
53 | 53 |
Value operator[](const Key &) const {return Value();} |
54 | 54 |
|
55 | 55 |
template<typename _ReadMap> |
56 | 56 |
struct Constraints { |
57 | 57 |
|
58 | 58 |
void constraints() { |
59 | 59 |
Value val = m[key]; |
60 | 60 |
val = m[key]; |
61 | 61 |
typename _ReadMap::Value own_val = m[own_key]; |
62 | 62 |
own_val = m[own_key]; |
63 | 63 |
|
64 | 64 |
ignore_unused_variable_warning(val); |
65 | 65 |
ignore_unused_variable_warning(own_val); |
66 | 66 |
ignore_unused_variable_warning(key); |
67 | 67 |
} |
68 | 68 |
Key& key; |
69 | 69 |
typename _ReadMap::Key& own_key; |
70 | 70 |
_ReadMap& m; |
71 | 71 |
}; |
72 | 72 |
|
73 | 73 |
}; |
74 | 74 |
|
75 | 75 |
|
76 | 76 |
/// Writable map concept |
77 | 77 |
|
78 | 78 |
/// Writable map concept. |
79 | 79 |
/// |
80 | 80 |
template<typename K, typename T> |
81 | 81 |
class WriteMap |
82 | 82 |
{ |
83 | 83 |
public: |
84 | 84 |
/// The key type of the map. |
85 | 85 |
typedef K Key; |
86 | 86 |
/// The value type of the map. (The type of objects associated with the keys). |
87 | 87 |
typedef T Value; |
88 | 88 |
|
89 | 89 |
/// Sets the value associated with a key. |
90 | 90 |
void set(const Key &,const Value &) {} |
91 | 91 |
|
92 | 92 |
///Default constructor |
93 | 93 |
WriteMap() {} |
94 | 94 |
|
95 | 95 |
template <typename _WriteMap> |
96 | 96 |
struct Constraints { |
97 | 97 |
void constraints() { |
98 | 98 |
// No constraints for constructor. |
99 | 99 |
m.set(key, val); |
100 | 100 |
m.set(own_key, own_val); |
101 | 101 |
ignore_unused_variable_warning(key); |
102 | 102 |
ignore_unused_variable_warning(val); |
103 | 103 |
ignore_unused_variable_warning(own_key); |
104 | 104 |
ignore_unused_variable_warning(own_val); |
105 | 105 |
} |
106 | 106 |
|
107 | 107 |
Value& val; |
108 | 108 |
typename _WriteMap::Value own_val; |
109 | 109 |
Key& key; |
110 | 110 |
typename _WriteMap::Key& own_key; |
111 | 111 |
_WriteMap& m; |
112 | 112 |
|
113 | 113 |
}; |
114 | 114 |
}; |
115 | 115 |
|
116 | 116 |
/// Read/Writable map concept |
117 | 117 |
|
118 | 118 |
/// Read/writable map concept. |
119 | 119 |
/// |
120 | 120 |
template<typename K, typename T> |
121 | 121 |
class ReadWriteMap : public ReadMap<K,T>, |
122 | 122 |
public WriteMap<K,T> |
123 | 123 |
{ |
124 | 124 |
public: |
125 | 125 |
/// The key type of the map. |
126 | 126 |
typedef K Key; |
127 | 127 |
/// The value type of the map. (The type of objects associated with the keys). |
128 | 128 |
typedef T Value; |
129 | 129 |
|
130 | 130 |
/// Returns the value associated with a key. |
131 | 131 |
Value operator[](const Key &) const {return Value();} |
132 | 132 |
/// Sets the value associated with a key. |
133 | 133 |
void set(const Key & ,const Value &) {} |
134 | 134 |
|
135 | 135 |
template<typename _ReadWriteMap> |
136 | 136 |
struct Constraints { |
137 | 137 |
void constraints() { |
138 | 138 |
checkConcept<ReadMap<K, T>, _ReadWriteMap >(); |
139 | 139 |
checkConcept<WriteMap<K, T>, _ReadWriteMap >(); |
140 | 140 |
} |
141 | 141 |
}; |
142 | 142 |
}; |
143 | 143 |
|
144 | 144 |
|
145 | 145 |
/// Dereferable map concept |
146 | 146 |
|
147 | 147 |
/// Dereferable map concept. |
148 | 148 |
/// |
149 | 149 |
template<typename K, typename T, typename R, typename CR> |
150 | 150 |
class ReferenceMap : public ReadWriteMap<K,T> |
151 | 151 |
{ |
152 | 152 |
public: |
153 | 153 |
/// Tag for reference maps. |
154 | 154 |
typedef True ReferenceMapTag; |
155 | 155 |
/// The key type of the map. |
156 | 156 |
typedef K Key; |
157 | 157 |
/// The value type of the map. (The type of objects associated with the keys). |
158 | 158 |
typedef T Value; |
159 | 159 |
/// The reference type of the map. |
160 | 160 |
typedef R Reference; |
161 | 161 |
/// The const reference type of the map. |
162 | 162 |
typedef CR ConstReference; |
163 | 163 |
|
164 | 164 |
protected: |
165 | 165 |
Value tmp; |
166 | 166 |
public: |
167 | 167 |
|
168 | 168 |
///Returns a reference to the value associated to a key. |
169 | 169 |
Reference operator[](const Key &) { return tmp; } |
170 | 170 |
///Returns a const reference to the value associated to a key. |
171 | 171 |
ConstReference operator[](const Key &) const { return tmp; } |
172 | 172 |
/// Sets the value associated with a key. |
173 | 173 |
void set(const Key &k,const Value &t) { operator[](k)=t; } |
174 | 174 |
|
175 | 175 |
/// \todo Rethink this concept. |
176 | 176 |
template<typename _ReferenceMap> |
177 | 177 |
struct ReferenceMapConcept { |
178 | 178 |
|
179 | 179 |
void constraints() { |
180 | 180 |
checkConcept<ReadWriteMap, _ReferenceMap >(); |
181 | 181 |
m[key] = val; |
182 | 182 |
val = m[key]; |
183 | 183 |
m[key] = ref; |
184 | 184 |
ref = m[key]; |
185 | 185 |
m[own_key] = own_val; |
186 | 186 |
own_val = m[own_key]; |
187 | 187 |
m[own_key] = own_ref; |
188 | 188 |
own_ref = m[own_key]; |
189 | 189 |
} |
190 | 190 |
|
191 | 191 |
typename _ReferenceMap::Key& own_key; |
192 | 192 |
typename _ReferenceMap::Value& own_val; |
193 | 193 |
typename _ReferenceMap::Reference& own_ref; |
194 | 194 |
Key& key; |
195 | 195 |
Value& val; |
196 | 196 |
Reference& ref; |
197 | 197 |
_ReferenceMap& m; |
198 | 198 |
}; |
199 | 199 |
}; |
200 | 200 |
|
201 | 201 |
// @} |
202 | 202 |
|
203 | 203 |
} //namespace concepts |
204 | 204 |
|
205 | 205 |
} //namespace lemon |
206 | 206 |
|
207 | 207 |
#endif // LEMON_CONCEPT_MAPS_H |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
#include <lemon/bits/utility.h> |
24 | 24 |
|
25 | 25 |
///\ingroup misc |
26 | 26 |
///\file |
27 | 27 |
///\brief A simple two dimensional vector and a bounding box implementation |
28 | 28 |
/// |
29 | 29 |
/// The class \ref lemon::dim2::Point "dim2::Point" implements |
30 | 30 |
///a two dimensional vector with the usual |
31 | 31 |
/// operations. |
32 | 32 |
/// |
33 | 33 |
/// The class \ref lemon::dim2::BoundingBox "dim2::BoundingBox" |
34 | 34 |
/// can be used to determine |
35 | 35 |
/// the rectangular bounding box of a set of |
36 | 36 |
/// \ref lemon::dim2::Point "dim2::Point"'s. |
37 | 37 |
|
38 | 38 |
namespace lemon { |
39 | 39 |
|
40 | 40 |
///Tools for handling two dimensional coordinates |
41 | 41 |
|
42 | 42 |
///This namespace is a storage of several |
43 | 43 |
///tools for handling two dimensional coordinates |
44 | 44 |
namespace dim2 { |
45 | 45 |
|
46 | 46 |
/// \addtogroup misc |
47 | 47 |
/// @{ |
48 | 48 |
|
49 | 49 |
/// A simple two dimensional vector (plainvector) implementation |
50 | 50 |
|
51 | 51 |
/// A simple two dimensional vector (plainvector) implementation |
52 | 52 |
///with the usual vector |
53 | 53 |
/// operators. |
54 | 54 |
/// |
55 | 55 |
template<typename T> |
56 | 56 |
class Point { |
57 | 57 |
|
58 | 58 |
public: |
59 | 59 |
|
60 | 60 |
typedef T Value; |
61 | 61 |
|
62 | 62 |
///First coordinate |
63 | 63 |
T x; |
64 | 64 |
///Second coordinate |
65 | 65 |
T y; |
66 | 66 |
|
67 | 67 |
///Default constructor |
68 | 68 |
Point() {} |
69 | 69 |
|
70 | 70 |
///Construct an instance from coordinates |
71 | 71 |
Point(T a, T b) : x(a), y(b) { } |
72 | 72 |
|
73 | 73 |
///The dimension of the vector. |
74 | 74 |
|
75 | 75 |
///The dimension of the vector. |
76 | 76 |
///This function always returns 2. |
77 | 77 |
int size() const { return 2; } |
78 | 78 |
|
79 | 79 |
///Subscripting operator |
80 | 80 |
|
81 | 81 |
///\c p[0] is \c p.x and \c p[1] is \c p.y |
82 | 82 |
/// |
83 | 83 |
T& operator[](int idx) { return idx == 0 ? x : y; } |
84 | 84 |
|
85 | 85 |
///Const subscripting operator |
86 | 86 |
|
87 | 87 |
///\c p[0] is \c p.x and \c p[1] is \c p.y |
88 | 88 |
/// |
89 | 89 |
const T& operator[](int idx) const { return idx == 0 ? x : y; } |
90 | 90 |
|
91 | 91 |
///Conversion constructor |
92 | 92 |
template<class TT> Point(const Point<TT> &p) : x(p.x), y(p.y) {} |
93 | 93 |
|
94 | 94 |
///Give back the square of the norm of the vector |
95 | 95 |
T normSquare() const { |
96 | 96 |
return x*x+y*y; |
97 | 97 |
} |
98 | 98 |
|
99 | 99 |
///Increment the left hand side by u |
100 | 100 |
Point<T>& operator +=(const Point<T>& u) { |
101 | 101 |
x += u.x; |
102 | 102 |
y += u.y; |
103 | 103 |
return *this; |
104 | 104 |
} |
105 | 105 |
|
106 | 106 |
///Decrement the left hand side by u |
107 | 107 |
Point<T>& operator -=(const Point<T>& u) { |
108 | 108 |
x -= u.x; |
109 | 109 |
y -= u.y; |
110 | 110 |
return *this; |
111 | 111 |
} |
112 | 112 |
|
113 | 113 |
///Multiply the left hand side with a scalar |
114 | 114 |
Point<T>& operator *=(const T &u) { |
115 | 115 |
x *= u; |
116 | 116 |
y *= u; |
117 | 117 |
return *this; |
118 | 118 |
} |
119 | 119 |
|
120 | 120 |
///Divide the left hand side by a scalar |
121 | 121 |
Point<T>& operator /=(const T &u) { |
122 | 122 |
x /= u; |
123 | 123 |
y /= u; |
124 | 124 |
return *this; |
125 | 125 |
} |
126 | 126 |
|
127 | 127 |
///Return the scalar product of two vectors |
128 | 128 |
T operator *(const Point<T>& u) const { |
129 | 129 |
return x*u.x+y*u.y; |
130 | 130 |
} |
131 | 131 |
|
132 | 132 |
///Return the sum of two vectors |
133 | 133 |
Point<T> operator+(const Point<T> &u) const { |
134 | 134 |
Point<T> b=*this; |
135 | 135 |
return b+=u; |
136 | 136 |
} |
137 | 137 |
|
138 | 138 |
///Return the negative of the vector |
139 | 139 |
Point<T> operator-() const { |
140 | 140 |
Point<T> b=*this; |
141 | 141 |
b.x=-b.x; b.y=-b.y; |
142 | 142 |
return b; |
143 | 143 |
} |
144 | 144 |
|
145 | 145 |
///Return the difference of two vectors |
146 | 146 |
Point<T> operator-(const Point<T> &u) const { |
147 | 147 |
Point<T> b=*this; |
148 | 148 |
return b-=u; |
149 | 149 |
} |
150 | 150 |
|
151 | 151 |
///Return a vector multiplied by a scalar |
152 | 152 |
Point<T> operator*(const T &u) const { |
153 | 153 |
Point<T> b=*this; |
154 | 154 |
return b*=u; |
155 | 155 |
} |
156 | 156 |
|
157 | 157 |
///Return a vector divided by a scalar |
158 | 158 |
Point<T> operator/(const T &u) const { |
159 | 159 |
Point<T> b=*this; |
160 | 160 |
return b/=u; |
161 | 161 |
} |
162 | 162 |
|
163 | 163 |
///Test equality |
164 | 164 |
bool operator==(const Point<T> &u) const { |
165 | 165 |
return (x==u.x) && (y==u.y); |
166 | 166 |
} |
167 | 167 |
|
168 | 168 |
///Test inequality |
169 | 169 |
bool operator!=(Point u) const { |
170 | 170 |
return (x!=u.x) || (y!=u.y); |
171 | 171 |
} |
172 | 172 |
|
173 | 173 |
}; |
174 | 174 |
|
175 | 175 |
///Return a Point |
176 | 176 |
|
177 | 177 |
///Return a Point. |
178 | 178 |
///\relates Point |
179 | 179 |
template <typename T> |
180 | 180 |
inline Point<T> makePoint(const T& x, const T& y) { |
181 | 181 |
return Point<T>(x, y); |
182 | 182 |
} |
183 | 183 |
|
184 | 184 |
///Return a vector multiplied by a scalar |
185 | 185 |
|
186 | 186 |
///Return a vector multiplied by a scalar. |
187 | 187 |
///\relates Point |
188 | 188 |
template<typename T> Point<T> operator*(const T &u,const Point<T> &x) { |
189 | 189 |
return x*u; |
190 | 190 |
} |
191 | 191 |
|
192 | 192 |
///Read a plainvector from a stream |
193 | 193 |
|
194 | 194 |
///Read a plainvector from a stream. |
195 | 195 |
///\relates Point |
196 | 196 |
/// |
197 | 197 |
template<typename T> |
198 | 198 |
inline std::istream& operator>>(std::istream &is, Point<T> &z) { |
199 | 199 |
char c; |
200 | 200 |
if (is >> c) { |
201 | 201 |
if (c != '(') is.putback(c); |
202 | 202 |
} else { |
203 | 203 |
is.clear(); |
204 | 204 |
} |
205 | 205 |
if (!(is >> z.x)) return is; |
206 | 206 |
if (is >> c) { |
207 | 207 |
if (c != ',') is.putback(c); |
208 | 208 |
} else { |
209 | 209 |
is.clear(); |
210 | 210 |
} |
211 | 211 |
if (!(is >> z.y)) return is; |
212 | 212 |
if (is >> c) { |
213 | 213 |
if (c != ')') is.putback(c); |
214 | 214 |
} else { |
215 | 215 |
is.clear(); |
216 | 216 |
} |
217 | 217 |
return is; |
218 | 218 |
} |
219 | 219 |
|
220 | 220 |
///Write a plainvector to a stream |
221 | 221 |
|
222 | 222 |
///Write a plainvector to a stream. |
223 | 223 |
///\relates Point |
224 | 224 |
/// |
225 | 225 |
template<typename T> |
226 | 226 |
inline std::ostream& operator<<(std::ostream &os, const Point<T>& z) |
227 | 227 |
{ |
228 | 228 |
os << "(" << z.x << ", " << z.y << ")"; |
229 | 229 |
return os; |
230 | 230 |
} |
231 | 231 |
|
232 | 232 |
///Rotate by 90 degrees |
233 | 233 |
|
234 | 234 |
///Returns the parameter rotated by 90 degrees in positive direction. |
235 | 235 |
///\relates Point |
236 | 236 |
/// |
237 | 237 |
template<typename T> |
238 | 238 |
inline Point<T> rot90(const Point<T> &z) |
239 | 239 |
{ |
240 | 240 |
return Point<T>(-z.y,z.x); |
241 | 241 |
} |
242 | 242 |
|
243 | 243 |
///Rotate by 180 degrees |
244 | 244 |
|
245 | 245 |
///Returns the parameter rotated by 180 degrees. |
246 | 246 |
///\relates Point |
247 | 247 |
/// |
248 | 248 |
template<typename T> |
249 | 249 |
inline Point<T> rot180(const Point<T> &z) |
250 | 250 |
{ |
251 | 251 |
return Point<T>(-z.x,-z.y); |
252 | 252 |
} |
253 | 253 |
|
254 | 254 |
///Rotate by 270 degrees |
255 | 255 |
|
256 | 256 |
///Returns the parameter rotated by 90 degrees in negative direction. |
257 | 257 |
///\relates Point |
258 | 258 |
/// |
259 | 259 |
template<typename T> |
260 | 260 |
inline Point<T> rot270(const Point<T> &z) |
261 | 261 |
{ |
262 | 262 |
return Point<T>(z.y,-z.x); |
263 | 263 |
} |
264 | 264 |
|
265 | 265 |
|
266 | 266 |
|
267 | 267 |
/// A class to calculate or store the bounding box of plainvectors. |
268 | 268 |
|
269 | 269 |
/// A class to calculate or store the bounding box of plainvectors. |
270 | 270 |
/// |
271 | 271 |
template<typename T> |
272 | 272 |
class BoundingBox { |
273 | 273 |
Point<T> bottom_left, top_right; |
274 | 274 |
bool _empty; |
275 | 275 |
public: |
276 | 276 |
|
277 | 277 |
///Default constructor: creates an empty bounding box |
278 | 278 |
BoundingBox() { _empty = true; } |
279 | 279 |
|
280 | 280 |
///Construct an instance from one point |
281 | 281 |
BoundingBox(Point<T> a) { bottom_left=top_right=a; _empty = false; } |
282 | 282 |
|
283 | 283 |
///Construct an instance from two points |
284 | 284 |
|
285 | 285 |
///Construct an instance from two points. |
286 | 286 |
///\param a The bottom left corner. |
287 | 287 |
///\param b The top right corner. |
288 | 288 |
///\warning The coordinates of the bottom left corner must be no more |
289 | 289 |
///than those of the top right one. |
290 | 290 |
BoundingBox(Point<T> a,Point<T> b) |
291 | 291 |
{ |
292 | 292 |
bottom_left=a; |
293 | 293 |
top_right=b; |
294 | 294 |
_empty = false; |
295 | 295 |
} |
296 | 296 |
|
297 | 297 |
///Construct an instance from four numbers |
298 | 298 |
|
299 | 299 |
///Construct an instance from four numbers. |
300 | 300 |
///\param l The left side of the box. |
301 | 301 |
///\param b The bottom of the box. |
302 | 302 |
///\param r The right side of the box. |
303 | 303 |
///\param t The top of the box. |
304 | 304 |
///\warning The left side must be no more than the right side and |
305 | 305 |
///bottom must be no more than the top. |
306 | 306 |
BoundingBox(T l,T b,T r,T t) |
307 | 307 |
{ |
308 | 308 |
bottom_left=Point<T>(l,b); |
309 | 309 |
top_right=Point<T>(r,t); |
310 | 310 |
_empty = false; |
311 | 311 |
} |
312 | 312 |
|
313 | 313 |
///Return \c true if the bounding box is empty. |
314 | 314 |
|
315 | 315 |
///Return \c true if the bounding box is empty (i.e. return \c false |
316 | 316 |
///if at least one point was added to the box or the coordinates of |
317 | 317 |
///the box were set). |
318 | 318 |
///The coordinates of an empty bounding box are not defined. |
319 | 319 |
bool empty() const { |
320 | 320 |
return _empty; |
321 | 321 |
} |
322 | 322 |
|
323 | 323 |
///Make the BoundingBox empty |
324 | 324 |
void clear() { |
325 | 325 |
_empty=1; |
326 | 326 |
} |
327 | 327 |
|
328 | 328 |
///Give back the bottom left corner |
329 | 329 |
|
330 | 330 |
///Give back the bottom left corner. |
331 | 331 |
///If the bounding box is empty, then the return value is not defined. |
332 | 332 |
Point<T> bottomLeft() const { |
333 | 333 |
return bottom_left; |
334 | 334 |
} |
335 | 335 |
|
336 | 336 |
///Set the bottom left corner |
337 | 337 |
|
338 | 338 |
///Set the bottom left corner. |
339 | 339 |
///It should only be used for non-empty box. |
340 | 340 |
void bottomLeft(Point<T> p) { |
341 | 341 |
bottom_left = p; |
342 | 342 |
} |
343 | 343 |
|
344 | 344 |
///Give back the top right corner |
345 | 345 |
|
346 | 346 |
///Give back the top right corner. |
347 | 347 |
///If the bounding box is empty, then the return value is not defined. |
348 | 348 |
Point<T> topRight() const { |
349 | 349 |
return top_right; |
350 | 350 |
} |
351 | 351 |
|
352 | 352 |
///Set the top right corner |
353 | 353 |
|
354 | 354 |
///Set the top right corner. |
355 | 355 |
///It should only be used for non-empty box. |
356 | 356 |
void topRight(Point<T> p) { |
357 | 357 |
top_right = p; |
358 | 358 |
} |
359 | 359 |
|
360 | 360 |
///Give back the bottom right corner |
361 | 361 |
|
362 | 362 |
///Give back the bottom right corner. |
363 | 363 |
///If the bounding box is empty, then the return value is not defined. |
364 | 364 |
Point<T> bottomRight() const { |
365 | 365 |
return Point<T>(top_right.x,bottom_left.y); |
366 | 366 |
} |
367 | 367 |
|
368 | 368 |
///Set the bottom right corner |
369 | 369 |
|
370 | 370 |
///Set the bottom right corner. |
371 | 371 |
///It should only be used for non-empty box. |
372 | 372 |
void bottomRight(Point<T> p) { |
373 | 373 |
top_right.x = p.x; |
374 | 374 |
bottom_left.y = p.y; |
375 | 375 |
} |
376 | 376 |
|
377 | 377 |
///Give back the top left corner |
378 | 378 |
|
379 | 379 |
///Give back the top left corner. |
380 | 380 |
///If the bounding box is empty, then the return value is not defined. |
381 | 381 |
Point<T> topLeft() const { |
382 | 382 |
return Point<T>(bottom_left.x,top_right.y); |
383 | 383 |
} |
384 | 384 |
|
385 | 385 |
///Set the top left corner |
386 | 386 |
|
387 | 387 |
///Set the top left corner. |
388 | 388 |
///It should only be used for non-empty box. |
389 | 389 |
void topLeft(Point<T> p) { |
1 |
/* -*- C++ -*- |
|
2 |
* |
|
3 |
* This file is a part of LEMON, a generic C++ optimization library |
|
4 |
* |
|
5 |
* Copyright (C) 2003-2008 |
|
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 |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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_MAPS_H |
20 | 20 |
#define LEMON_MAPS_H |
21 | 21 |
|
22 | 22 |
#include <iterator> |
23 | 23 |
#include <functional> |
24 | 24 |
#include <vector> |
25 | 25 |
|
26 | 26 |
#include <lemon/bits/utility.h> |
27 | 27 |
// #include <lemon/bits/traits.h> |
28 | 28 |
|
29 | 29 |
///\file |
30 | 30 |
///\ingroup maps |
31 | 31 |
///\brief Miscellaneous property maps |
32 | 32 |
/// |
33 | 33 |
#include <map> |
34 | 34 |
|
35 | 35 |
namespace lemon { |
36 | 36 |
|
37 | 37 |
/// \addtogroup maps |
38 | 38 |
/// @{ |
39 | 39 |
|
40 | 40 |
/// Base class of maps. |
41 | 41 |
|
42 | 42 |
/// Base class of maps. |
43 | 43 |
/// It provides the necessary <tt>typedef</tt>s required by the map concept. |
44 | 44 |
template<typename K, typename T> |
45 | 45 |
class MapBase { |
46 | 46 |
public: |
47 | 47 |
/// The key type of the map. |
48 | 48 |
typedef K Key; |
49 | 49 |
/// The value type of the map. (The type of objects associated with the keys). |
50 | 50 |
typedef T Value; |
51 | 51 |
}; |
52 | 52 |
|
53 | 53 |
/// Null map. (a.k.a. DoNothingMap) |
54 | 54 |
|
55 | 55 |
/// This map can be used if you have to provide a map only for |
56 | 56 |
/// its type definitions, or if you have to provide a writable map, |
57 | 57 |
/// but data written to it is not required (i.e. it will be sent to |
58 | 58 |
/// <tt>/dev/null</tt>). |
59 | 59 |
template<typename K, typename T> |
60 | 60 |
class NullMap : public MapBase<K, T> { |
61 | 61 |
public: |
62 | 62 |
typedef MapBase<K, T> Parent; |
63 | 63 |
typedef typename Parent::Key Key; |
64 | 64 |
typedef typename Parent::Value Value; |
65 | 65 |
|
66 | 66 |
/// Gives back a default constructed element. |
67 | 67 |
T operator[](const K&) const { return T(); } |
68 | 68 |
/// Absorbs the value. |
69 | 69 |
void set(const K&, const T&) {} |
70 | 70 |
}; |
71 | 71 |
|
72 | 72 |
///Returns a \c NullMap class |
73 | 73 |
|
74 | 74 |
///This function just returns a \c NullMap class. |
75 | 75 |
///\relates NullMap |
76 | 76 |
template <typename K, typename V> |
77 | 77 |
NullMap<K, V> nullMap() { |
78 | 78 |
return NullMap<K, V>(); |
79 | 79 |
} |
80 | 80 |
|
81 | 81 |
|
82 | 82 |
/// Constant map. |
83 | 83 |
|
84 | 84 |
/// This is a readable map which assigns a specified value to each key. |
85 | 85 |
/// In other aspects it is equivalent to the \c NullMap. |
86 | 86 |
template<typename K, typename T> |
87 | 87 |
class ConstMap : public MapBase<K, T> { |
88 | 88 |
private: |
89 | 89 |
T v; |
90 | 90 |
public: |
91 | 91 |
|
92 | 92 |
typedef MapBase<K, T> Parent; |
93 | 93 |
typedef typename Parent::Key Key; |
94 | 94 |
typedef typename Parent::Value Value; |
95 | 95 |
|
96 | 96 |
/// Default constructor |
97 | 97 |
|
98 | 98 |
/// Default constructor. |
99 | 99 |
/// The value of the map will be uninitialized. |
100 | 100 |
/// (More exactly it will be default constructed.) |
101 | 101 |
ConstMap() {} |
102 | 102 |
|
103 | 103 |
/// Constructor with specified initial value |
104 | 104 |
|
105 | 105 |
/// Constructor with specified initial value. |
106 | 106 |
/// \param _v is the initial value of the map. |
107 | 107 |
ConstMap(const T &_v) : v(_v) {} |
108 | 108 |
|
109 | 109 |
///\e |
110 | 110 |
T operator[](const K&) const { return v; } |
111 | 111 |
|
112 | 112 |
///\e |
113 | 113 |
void setAll(const T &t) { |
114 | 114 |
v = t; |
115 | 115 |
} |
116 | 116 |
|
117 | 117 |
template<typename T1> |
118 | 118 |
struct rebind { |
119 | 119 |
typedef ConstMap<K, T1> other; |
120 | 120 |
}; |
121 | 121 |
|
122 | 122 |
template<typename T1> |
123 | 123 |
ConstMap(const ConstMap<K, T1> &, const T &_v) : v(_v) {} |
124 | 124 |
}; |
125 | 125 |
|
126 | 126 |
///Returns a \c ConstMap class |
127 | 127 |
|
128 | 128 |
///This function just returns a \c ConstMap class. |
129 | 129 |
///\relates ConstMap |
130 | 130 |
template<typename K, typename V> |
131 | 131 |
inline ConstMap<K, V> constMap(const V &v) { |
132 | 132 |
return ConstMap<K, V>(v); |
133 | 133 |
} |
134 | 134 |
|
135 | 135 |
|
136 | 136 |
template<typename T, T v> |
137 | 137 |
struct Const { }; |
138 | 138 |
|
139 | 139 |
/// Constant map with inlined constant value. |
140 | 140 |
|
141 | 141 |
/// This is a readable map which assigns a specified value to each key. |
142 | 142 |
/// In other aspects it is equivalent to the \c NullMap. |
143 | 143 |
template<typename K, typename V, V v> |
144 | 144 |
class ConstMap<K, Const<V, v> > : public MapBase<K, V> { |
145 | 145 |
public: |
146 | 146 |
typedef MapBase<K, V> Parent; |
147 | 147 |
typedef typename Parent::Key Key; |
148 | 148 |
typedef typename Parent::Value Value; |
149 | 149 |
|
150 | 150 |
ConstMap() { } |
151 | 151 |
///\e |
152 | 152 |
V operator[](const K&) const { return v; } |
153 | 153 |
///\e |
154 | 154 |
void set(const K&, const V&) { } |
155 | 155 |
}; |
156 | 156 |
|
157 | 157 |
///Returns a \c ConstMap class |
158 | 158 |
|
159 | 159 |
///This function just returns a \c ConstMap class with inlined value. |
160 | 160 |
///\relates ConstMap |
161 | 161 |
template<typename K, typename V, V v> |
162 | 162 |
inline ConstMap<K, Const<V, v> > constMap() { |
163 | 163 |
return ConstMap<K, Const<V, v> >(); |
164 | 164 |
} |
165 | 165 |
|
166 | 166 |
///Map based on std::map |
167 | 167 |
|
168 | 168 |
///This is essentially a wrapper for \c std::map with addition that |
169 | 169 |
///you can specify a default value different from \c Value(). |
170 | 170 |
template <typename K, typename T, typename Compare = std::less<K> > |
171 | 171 |
class StdMap { |
172 | 172 |
template <typename K1, typename T1, typename C1> |
173 | 173 |
friend class StdMap; |
174 | 174 |
public: |
175 | 175 |
|
176 | 176 |
typedef True ReferenceMapTag; |
177 | 177 |
///\e |
178 | 178 |
typedef K Key; |
179 | 179 |
///\e |
180 | 180 |
typedef T Value; |
181 | 181 |
///\e |
182 | 182 |
typedef T& Reference; |
183 | 183 |
///\e |
184 | 184 |
typedef const T& ConstReference; |
185 | 185 |
|
186 | 186 |
private: |
187 | 187 |
|
188 | 188 |
typedef std::map<K, T, Compare> Map; |
189 | 189 |
Value _value; |
190 | 190 |
Map _map; |
191 | 191 |
|
192 | 192 |
public: |
193 | 193 |
|
194 | 194 |
/// Constructor with specified default value |
195 | 195 |
StdMap(const T& value = T()) : _value(value) {} |
196 | 196 |
/// \brief Constructs the map from an appropriate std::map, and explicitly |
197 | 197 |
/// specifies a default value. |
198 | 198 |
template <typename T1, typename Comp1> |
199 | 199 |
StdMap(const std::map<Key, T1, Comp1> &map, const T& value = T()) |
200 | 200 |
: _map(map.begin(), map.end()), _value(value) {} |
201 | 201 |
|
202 | 202 |
/// \brief Constructs a map from an other StdMap. |
203 | 203 |
template<typename T1, typename Comp1> |
204 | 204 |
StdMap(const StdMap<Key, T1, Comp1> &c) |
205 | 205 |
: _map(c._map.begin(), c._map.end()), _value(c._value) {} |
206 | 206 |
|
207 | 207 |
private: |
208 | 208 |
|
209 | 209 |
StdMap& operator=(const StdMap&); |
210 | 210 |
|
211 | 211 |
public: |
212 | 212 |
|
213 | 213 |
///\e |
214 | 214 |
Reference operator[](const Key &k) { |
215 | 215 |
typename Map::iterator it = _map.lower_bound(k); |
216 | 216 |
if (it != _map.end() && !_map.key_comp()(k, it->first)) |
217 | 217 |
return it->second; |
218 | 218 |
else |
219 | 219 |
return _map.insert(it, std::make_pair(k, _value))->second; |
220 | 220 |
} |
221 | 221 |
|
222 | 222 |
/// \e |
223 | 223 |
ConstReference operator[](const Key &k) const { |
224 | 224 |
typename Map::const_iterator it = _map.find(k); |
225 | 225 |
if (it != _map.end()) |
226 | 226 |
return it->second; |
227 | 227 |
else |
228 | 228 |
return _value; |
229 | 229 |
} |
230 | 230 |
|
231 | 231 |
/// \e |
232 | 232 |
void set(const Key &k, const T &t) { |
233 | 233 |
typename Map::iterator it = _map.lower_bound(k); |
234 | 234 |
if (it != _map.end() && !_map.key_comp()(k, it->first)) |
235 | 235 |
it->second = t; |
236 | 236 |
else |
237 | 237 |
_map.insert(it, std::make_pair(k, t)); |
238 | 238 |
} |
239 | 239 |
|
240 | 240 |
/// \e |
241 | 241 |
void setAll(const T &t) { |
242 | 242 |
_value = t; |
243 | 243 |
_map.clear(); |
244 | 244 |
} |
245 | 245 |
|
246 | 246 |
template <typename T1, typename C1 = std::less<T1> > |
247 | 247 |
struct rebind { |
248 | 248 |
typedef StdMap<Key, T1, C1> other; |
249 | 249 |
}; |
250 | 250 |
}; |
251 | 251 |
|
252 | 252 |
/// \brief Map for storing values for keys from the range <tt>[0..size-1]</tt> |
253 | 253 |
/// |
254 | 254 |
/// The current map has the <tt>[0..size-1]</tt> keyset and the values |
255 | 255 |
/// are stored in a \c std::vector<T> container. It can be used with |
256 | 256 |
/// some data structures, for example \c UnionFind, \c BinHeap, when |
257 | 257 |
/// the used items are small integer numbers. |
258 | 258 |
/// |
259 | 259 |
/// \todo Revise its name |
260 | 260 |
template <typename T> |
261 | 261 |
class IntegerMap { |
262 | 262 |
|
263 | 263 |
template <typename T1> |
264 | 264 |
friend class IntegerMap; |
265 | 265 |
|
266 | 266 |
public: |
267 | 267 |
|
268 | 268 |
typedef True ReferenceMapTag; |
269 | 269 |
///\e |
270 | 270 |
typedef int Key; |
271 | 271 |
///\e |
272 | 272 |
typedef T Value; |
273 | 273 |
///\e |
274 | 274 |
typedef T& Reference; |
275 | 275 |
///\e |
276 | 276 |
typedef const T& ConstReference; |
277 | 277 |
|
278 | 278 |
private: |
279 | 279 |
|
280 | 280 |
typedef std::vector<T> Vector; |
281 | 281 |
Vector _vector; |
282 | 282 |
|
283 | 283 |
public: |
284 | 284 |
|
285 | 285 |
/// Constructor with specified default value |
286 | 286 |
IntegerMap(int size = 0, const T& value = T()) : _vector(size, value) {} |
287 | 287 |
|
288 | 288 |
/// \brief Constructs the map from an appropriate std::vector. |
289 | 289 |
template <typename T1> |
290 | 290 |
IntegerMap(const std::vector<T1>& vector) |
291 | 291 |
: _vector(vector.begin(), vector.end()) {} |
292 | 292 |
|
293 | 293 |
/// \brief Constructs a map from an other IntegerMap. |
294 | 294 |
template <typename T1> |
295 | 295 |
IntegerMap(const IntegerMap<T1> &c) |
296 | 296 |
: _vector(c._vector.begin(), c._vector.end()) {} |
297 | 297 |
|
298 | 298 |
/// \brief Resize the container |
299 | 299 |
void resize(int size, const T& value = T()) { |
300 | 300 |
_vector.resize(size, value); |
301 | 301 |
} |
302 | 302 |
|
303 | 303 |
private: |
304 | 304 |
|
305 | 305 |
IntegerMap& operator=(const IntegerMap&); |
306 | 306 |
|
307 | 307 |
public: |
308 | 308 |
|
309 | 309 |
///\e |
310 | 310 |
Reference operator[](Key k) { |
311 | 311 |
return _vector[k]; |
312 | 312 |
} |
313 | 313 |
|
314 | 314 |
/// \e |
315 | 315 |
ConstReference operator[](Key k) const { |
316 | 316 |
return _vector[k]; |
317 | 317 |
} |
318 | 318 |
|
319 | 319 |
/// \e |
320 | 320 |
void set(const Key &k, const T& t) { |
321 | 321 |
_vector[k] = t; |
322 | 322 |
} |
323 | 323 |
|
324 | 324 |
}; |
325 | 325 |
|
326 | 326 |
/// @} |
327 | 327 |
|
328 | 328 |
/// \addtogroup map_adaptors |
329 | 329 |
/// @{ |
330 | 330 |
|
331 | 331 |
/// \brief Identity map. |
332 | 332 |
/// |
333 | 333 |
/// This map gives back the given key as value without any |
334 | 334 |
/// modification. |
335 | 335 |
template <typename T> |
336 | 336 |
class IdentityMap : public MapBase<T, T> { |
337 | 337 |
public: |
338 | 338 |
typedef MapBase<T, T> Parent; |
339 | 339 |
typedef typename Parent::Key Key; |
340 | 340 |
typedef typename Parent::Value Value; |
341 | 341 |
|
342 | 342 |
/// \e |
343 | 343 |
const T& operator[](const T& t) const { |
344 | 344 |
return t; |
345 | 345 |
} |
346 | 346 |
}; |
347 | 347 |
|
348 | 348 |
///Returns an \c IdentityMap class |
349 | 349 |
|
350 | 350 |
///This function just returns an \c IdentityMap class. |
351 | 351 |
///\relates IdentityMap |
352 | 352 |
template<typename T> |
353 | 353 |
inline IdentityMap<T> identityMap() { |
354 | 354 |
return IdentityMap<T>(); |
355 | 355 |
} |
356 | 356 |
|
357 | 357 |
|
358 | 358 |
///\brief Convert the \c Value of a map to another type using |
359 | 359 |
///the default conversion. |
360 | 360 |
/// |
361 | 361 |
///This \c concepts::ReadMap "read only map" |
362 | 362 |
///converts the \c Value of a map to type \c T. |
363 | 363 |
///Its \c Key is inherited from \c M. |
364 | 364 |
template <typename M, typename T> |
365 | 365 |
class ConvertMap : public MapBase<typename M::Key, T> { |
366 | 366 |
const M& m; |
367 | 367 |
public: |
368 | 368 |
typedef MapBase<typename M::Key, T> Parent; |
369 | 369 |
typedef typename Parent::Key Key; |
370 | 370 |
typedef typename Parent::Value Value; |
371 | 371 |
|
372 | 372 |
///Constructor |
373 | 373 |
|
374 | 374 |
///Constructor. |
375 | 375 |
///\param _m is the underlying map. |
376 | 376 |
ConvertMap(const M &_m) : m(_m) {}; |
377 | 377 |
|
378 | 378 |
/// \brief The subscript operator. |
379 | 379 |
/// |
380 | 380 |
/// The subscript operator. |
381 | 381 |
Value operator[](const Key& k) const {return m[k];} |
382 | 382 |
}; |
383 | 383 |
|
384 | 384 |
///Returns a \c ConvertMap class |
385 | 385 |
|
386 | 386 |
///This function just returns a \c ConvertMap class. |
387 | 387 |
///\relates ConvertMap |
388 | 388 |
template<typename T, typename M> |
389 | 389 |
inline ConvertMap<M, T> convertMap(const M &m) { |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
///\file |
20 | 20 |
///\brief Instantiation of the Random class. |
21 | 21 |
|
22 | 22 |
#include <lemon/random.h> |
23 | 23 |
|
24 | 24 |
namespace lemon { |
25 | 25 |
/// \brief Global random number generator instance |
26 | 26 |
/// |
27 | 27 |
/// A global Mersenne Twister random number generator instance. |
28 | 28 |
Random rnd; |
29 | 29 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
/* |
20 | 20 |
* This file contains the reimplemented version of the Mersenne Twister |
21 | 21 |
* Generator of Matsumoto and Nishimura. |
22 | 22 |
* |
23 | 23 |
* See the appropriate copyright notice below. |
24 | 24 |
* |
25 | 25 |
* Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, |
26 | 26 |
* All rights reserved. |
27 | 27 |
* |
28 | 28 |
* Redistribution and use in source and binary forms, with or without |
29 | 29 |
* modification, are permitted provided that the following conditions |
30 | 30 |
* are met: |
31 | 31 |
* |
32 | 32 |
* 1. Redistributions of source code must retain the above copyright |
33 | 33 |
* notice, this list of conditions and the following disclaimer. |
34 | 34 |
* |
35 | 35 |
* 2. Redistributions in binary form must reproduce the above copyright |
36 | 36 |
* notice, this list of conditions and the following disclaimer in the |
37 | 37 |
* documentation and/or other materials provided with the distribution. |
38 | 38 |
* |
39 | 39 |
* 3. The names of its contributors may not be used to endorse or promote |
40 | 40 |
* products derived from this software without specific prior written |
41 | 41 |
* permission. |
42 | 42 |
* |
43 | 43 |
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
44 | 44 |
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
45 | 45 |
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS |
46 | 46 |
* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE |
47 | 47 |
* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, |
48 | 48 |
* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES |
49 | 49 |
* (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR |
50 | 50 |
* SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) |
51 | 51 |
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, |
52 | 52 |
* STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) |
53 | 53 |
* ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED |
54 | 54 |
* OF THE POSSIBILITY OF SUCH DAMAGE. |
55 | 55 |
* |
56 | 56 |
* |
57 | 57 |
* Any feedback is very welcome. |
58 | 58 |
* http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html |
59 | 59 |
* email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space) |
60 | 60 |
*/ |
61 | 61 |
|
62 | 62 |
#ifndef LEMON_RANDOM_H |
63 | 63 |
#define LEMON_RANDOM_H |
64 | 64 |
|
65 | 65 |
#include <algorithm> |
66 | 66 |
#include <iterator> |
67 | 67 |
#include <vector> |
68 | 68 |
|
69 | 69 |
#include <ctime> |
70 | 70 |
#include <cmath> |
71 | 71 |
|
72 | 72 |
#include <lemon/dim2.h> |
73 | 73 |
///\ingroup misc |
74 | 74 |
///\file |
75 | 75 |
///\brief Mersenne Twister random number generator |
76 | 76 |
|
77 | 77 |
namespace lemon { |
78 | 78 |
|
79 | 79 |
namespace _random_bits { |
80 | 80 |
|
81 | 81 |
template <typename _Word, int _bits = std::numeric_limits<_Word>::digits> |
82 | 82 |
struct RandomTraits {}; |
83 | 83 |
|
84 | 84 |
template <typename _Word> |
85 | 85 |
struct RandomTraits<_Word, 32> { |
86 | 86 |
|
87 | 87 |
typedef _Word Word; |
88 | 88 |
static const int bits = 32; |
89 | 89 |
|
90 | 90 |
static const int length = 624; |
91 | 91 |
static const int shift = 397; |
92 | 92 |
|
93 | 93 |
static const Word mul = 0x6c078965u; |
94 | 94 |
static const Word arrayInit = 0x012BD6AAu; |
95 | 95 |
static const Word arrayMul1 = 0x0019660Du; |
96 | 96 |
static const Word arrayMul2 = 0x5D588B65u; |
97 | 97 |
|
98 | 98 |
static const Word mask = 0x9908B0DFu; |
99 | 99 |
static const Word loMask = (1u << 31) - 1; |
100 | 100 |
static const Word hiMask = ~loMask; |
101 | 101 |
|
102 | 102 |
|
103 | 103 |
static Word tempering(Word rnd) { |
104 | 104 |
rnd ^= (rnd >> 11); |
105 | 105 |
rnd ^= (rnd << 7) & 0x9D2C5680u; |
106 | 106 |
rnd ^= (rnd << 15) & 0xEFC60000u; |
107 | 107 |
rnd ^= (rnd >> 18); |
108 | 108 |
return rnd; |
109 | 109 |
} |
110 | 110 |
|
111 | 111 |
}; |
112 | 112 |
|
113 | 113 |
template <typename _Word> |
114 | 114 |
struct RandomTraits<_Word, 64> { |
115 | 115 |
|
116 | 116 |
typedef _Word Word; |
117 | 117 |
static const int bits = 64; |
118 | 118 |
|
119 | 119 |
static const int length = 312; |
120 | 120 |
static const int shift = 156; |
121 | 121 |
|
122 | 122 |
static const Word mul = Word(0x5851F42Du) << 32 | Word(0x4C957F2Du); |
123 | 123 |
static const Word arrayInit = Word(0x00000000u) << 32 |Word(0x012BD6AAu); |
124 | 124 |
static const Word arrayMul1 = Word(0x369DEA0Fu) << 32 |Word(0x31A53F85u); |
125 | 125 |
static const Word arrayMul2 = Word(0x27BB2EE6u) << 32 |Word(0x87B0B0FDu); |
126 | 126 |
|
127 | 127 |
static const Word mask = Word(0xB5026F5Au) << 32 | Word(0xA96619E9u); |
128 | 128 |
static const Word loMask = (Word(1u) << 31) - 1; |
129 | 129 |
static const Word hiMask = ~loMask; |
130 | 130 |
|
131 | 131 |
static Word tempering(Word rnd) { |
132 | 132 |
rnd ^= (rnd >> 29) & (Word(0x55555555u) << 32 | Word(0x55555555u)); |
133 | 133 |
rnd ^= (rnd << 17) & (Word(0x71D67FFFu) << 32 | Word(0xEDA60000u)); |
134 | 134 |
rnd ^= (rnd << 37) & (Word(0xFFF7EEE0u) << 32 | Word(0x00000000u)); |
135 | 135 |
rnd ^= (rnd >> 43); |
136 | 136 |
return rnd; |
137 | 137 |
} |
138 | 138 |
|
139 | 139 |
}; |
140 | 140 |
|
141 | 141 |
template <typename _Word> |
142 | 142 |
class RandomCore { |
143 | 143 |
public: |
144 | 144 |
|
145 | 145 |
typedef _Word Word; |
146 | 146 |
|
147 | 147 |
private: |
148 | 148 |
|
149 | 149 |
static const int bits = RandomTraits<Word>::bits; |
150 | 150 |
|
151 | 151 |
static const int length = RandomTraits<Word>::length; |
152 | 152 |
static const int shift = RandomTraits<Word>::shift; |
153 | 153 |
|
154 | 154 |
public: |
155 | 155 |
|
156 | 156 |
void initState() { |
157 | 157 |
static const Word seedArray[4] = { |
158 | 158 |
0x12345u, 0x23456u, 0x34567u, 0x45678u |
159 | 159 |
}; |
160 | 160 |
|
161 | 161 |
initState(seedArray, seedArray + 4); |
162 | 162 |
} |
163 | 163 |
|
164 | 164 |
void initState(Word seed) { |
165 | 165 |
|
166 | 166 |
static const Word mul = RandomTraits<Word>::mul; |
167 | 167 |
|
168 | 168 |
current = state; |
169 | 169 |
|
170 | 170 |
Word *curr = state + length - 1; |
171 | 171 |
curr[0] = seed; --curr; |
172 | 172 |
for (int i = 1; i < length; ++i) { |
173 | 173 |
curr[0] = (mul * ( curr[1] ^ (curr[1] >> (bits - 2)) ) + i); |
174 | 174 |
--curr; |
175 | 175 |
} |
176 | 176 |
} |
177 | 177 |
|
178 | 178 |
template <typename Iterator> |
179 | 179 |
void initState(Iterator begin, Iterator end) { |
180 | 180 |
|
181 | 181 |
static const Word init = RandomTraits<Word>::arrayInit; |
182 | 182 |
static const Word mul1 = RandomTraits<Word>::arrayMul1; |
183 | 183 |
static const Word mul2 = RandomTraits<Word>::arrayMul2; |
184 | 184 |
|
185 | 185 |
|
186 | 186 |
Word *curr = state + length - 1; --curr; |
187 | 187 |
Iterator it = begin; int cnt = 0; |
188 | 188 |
int num; |
189 | 189 |
|
190 | 190 |
initState(init); |
191 | 191 |
|
192 | 192 |
num = length > end - begin ? length : end - begin; |
193 | 193 |
while (num--) { |
194 | 194 |
curr[0] = (curr[0] ^ ((curr[1] ^ (curr[1] >> (bits - 2))) * mul1)) |
195 | 195 |
+ *it + cnt; |
196 | 196 |
++it; ++cnt; |
197 | 197 |
if (it == end) { |
198 | 198 |
it = begin; cnt = 0; |
199 | 199 |
} |
200 | 200 |
if (curr == state) { |
201 | 201 |
curr = state + length - 1; curr[0] = state[0]; |
202 | 202 |
} |
203 | 203 |
--curr; |
204 | 204 |
} |
205 | 205 |
|
206 | 206 |
num = length - 1; cnt = length - (curr - state) - 1; |
207 | 207 |
while (num--) { |
208 | 208 |
curr[0] = (curr[0] ^ ((curr[1] ^ (curr[1] >> (bits - 2))) * mul2)) |
209 | 209 |
- cnt; |
210 | 210 |
--curr; ++cnt; |
211 | 211 |
if (curr == state) { |
212 | 212 |
curr = state + length - 1; curr[0] = state[0]; --curr; |
213 | 213 |
cnt = 1; |
214 | 214 |
} |
215 | 215 |
} |
216 | 216 |
|
217 | 217 |
state[length - 1] = Word(1) << (bits - 1); |
218 | 218 |
} |
219 | 219 |
|
220 | 220 |
void copyState(const RandomCore& other) { |
221 | 221 |
std::copy(other.state, other.state + length, state); |
222 | 222 |
current = state + (other.current - other.state); |
223 | 223 |
} |
224 | 224 |
|
225 | 225 |
Word operator()() { |
226 | 226 |
if (current == state) fillState(); |
227 | 227 |
--current; |
228 | 228 |
Word rnd = *current; |
229 | 229 |
return RandomTraits<Word>::tempering(rnd); |
230 | 230 |
} |
231 | 231 |
|
232 | 232 |
private: |
233 | 233 |
|
234 | 234 |
|
235 | 235 |
void fillState() { |
236 | 236 |
static const Word mask[2] = { 0x0ul, RandomTraits<Word>::mask }; |
237 | 237 |
static const Word loMask = RandomTraits<Word>::loMask; |
238 | 238 |
static const Word hiMask = RandomTraits<Word>::hiMask; |
239 | 239 |
|
240 | 240 |
current = state + length; |
241 | 241 |
|
242 | 242 |
register Word *curr = state + length - 1; |
243 | 243 |
register long num; |
244 | 244 |
|
245 | 245 |
num = length - shift; |
246 | 246 |
while (num--) { |
247 | 247 |
curr[0] = (((curr[0] & hiMask) | (curr[-1] & loMask)) >> 1) ^ |
248 | 248 |
curr[- shift] ^ mask[curr[-1] & 1ul]; |
249 | 249 |
--curr; |
250 | 250 |
} |
251 | 251 |
num = shift - 1; |
252 | 252 |
while (num--) { |
253 | 253 |
curr[0] = (((curr[0] & hiMask) | (curr[-1] & loMask)) >> 1) ^ |
254 | 254 |
curr[length - shift] ^ mask[curr[-1] & 1ul]; |
255 | 255 |
--curr; |
256 | 256 |
} |
257 | 257 |
curr[0] = (((curr[0] & hiMask) | (curr[length - 1] & loMask)) >> 1) ^ |
258 | 258 |
curr[length - shift] ^ mask[curr[length - 1] & 1ul]; |
259 | 259 |
|
260 | 260 |
} |
261 | 261 |
|
262 | 262 |
|
263 | 263 |
Word *current; |
264 | 264 |
Word state[length]; |
265 | 265 |
|
266 | 266 |
}; |
267 | 267 |
|
268 | 268 |
|
269 | 269 |
template <typename Result, |
270 | 270 |
int shift = (std::numeric_limits<Result>::digits + 1) / 2> |
271 | 271 |
struct Masker { |
272 | 272 |
static Result mask(const Result& result) { |
273 | 273 |
return Masker<Result, (shift + 1) / 2>:: |
274 | 274 |
mask(static_cast<Result>(result | (result >> shift))); |
275 | 275 |
} |
276 | 276 |
}; |
277 | 277 |
|
278 | 278 |
template <typename Result> |
279 | 279 |
struct Masker<Result, 1> { |
280 | 280 |
static Result mask(const Result& result) { |
281 | 281 |
return static_cast<Result>(result | (result >> 1)); |
282 | 282 |
} |
283 | 283 |
}; |
284 | 284 |
|
285 | 285 |
template <typename Result, typename Word, |
286 | 286 |
int rest = std::numeric_limits<Result>::digits, int shift = 0, |
287 | 287 |
bool last = rest <= std::numeric_limits<Word>::digits> |
288 | 288 |
struct IntConversion { |
289 | 289 |
static const int bits = std::numeric_limits<Word>::digits; |
290 | 290 |
|
291 | 291 |
static Result convert(RandomCore<Word>& rnd) { |
292 | 292 |
return static_cast<Result>(rnd() >> (bits - rest)) << shift; |
293 | 293 |
} |
294 | 294 |
|
295 | 295 |
}; |
296 | 296 |
|
297 | 297 |
template <typename Result, typename Word, int rest, int shift> |
298 | 298 |
struct IntConversion<Result, Word, rest, shift, false> { |
299 | 299 |
static const int bits = std::numeric_limits<Word>::digits; |
300 | 300 |
|
301 | 301 |
static Result convert(RandomCore<Word>& rnd) { |
302 | 302 |
return (static_cast<Result>(rnd()) << shift) | |
303 | 303 |
IntConversion<Result, Word, rest - bits, shift + bits>::convert(rnd); |
304 | 304 |
} |
305 | 305 |
}; |
306 | 306 |
|
307 | 307 |
|
308 | 308 |
template <typename Result, typename Word, |
309 | 309 |
bool one_word = (std::numeric_limits<Word>::digits < |
310 | 310 |
std::numeric_limits<Result>::digits) > |
311 | 311 |
struct Mapping { |
312 | 312 |
static Result map(RandomCore<Word>& rnd, const Result& bound) { |
313 | 313 |
Word max = Word(bound - 1); |
314 | 314 |
Result mask = Masker<Result>::mask(bound - 1); |
315 | 315 |
Result num; |
316 | 316 |
do { |
317 | 317 |
num = IntConversion<Result, Word>::convert(rnd) & mask; |
318 | 318 |
} while (num > max); |
319 | 319 |
return num; |
320 | 320 |
} |
321 | 321 |
}; |
322 | 322 |
|
323 | 323 |
template <typename Result, typename Word> |
324 | 324 |
struct Mapping<Result, Word, false> { |
325 | 325 |
static Result map(RandomCore<Word>& rnd, const Result& bound) { |
326 | 326 |
Word max = Word(bound - 1); |
327 | 327 |
Word mask = Masker<Word, (std::numeric_limits<Result>::digits + 1) / 2> |
328 | 328 |
::mask(max); |
329 | 329 |
Word num; |
330 | 330 |
do { |
331 | 331 |
num = rnd() & mask; |
332 | 332 |
} while (num > max); |
333 | 333 |
return num; |
334 | 334 |
} |
335 | 335 |
}; |
336 | 336 |
|
337 | 337 |
template <typename Result, int exp, bool pos = (exp >= 0)> |
338 | 338 |
struct ShiftMultiplier { |
339 | 339 |
static const Result multiplier() { |
340 | 340 |
Result res = ShiftMultiplier<Result, exp / 2>::multiplier(); |
341 | 341 |
res *= res; |
342 | 342 |
if ((exp & 1) == 1) res *= static_cast<Result>(2.0); |
343 | 343 |
return res; |
344 | 344 |
} |
345 | 345 |
}; |
346 | 346 |
|
347 | 347 |
template <typename Result, int exp> |
348 | 348 |
struct ShiftMultiplier<Result, exp, false> { |
349 | 349 |
static const Result multiplier() { |
350 | 350 |
Result res = ShiftMultiplier<Result, exp / 2>::multiplier(); |
351 | 351 |
res *= res; |
352 | 352 |
if ((exp & 1) == 1) res *= static_cast<Result>(0.5); |
353 | 353 |
return res; |
354 | 354 |
} |
355 | 355 |
}; |
356 | 356 |
|
357 | 357 |
template <typename Result> |
358 | 358 |
struct ShiftMultiplier<Result, 0, true> { |
359 | 359 |
static const Result multiplier() { |
360 | 360 |
return static_cast<Result>(1.0); |
361 | 361 |
} |
362 | 362 |
}; |
363 | 363 |
|
364 | 364 |
template <typename Result> |
365 | 365 |
struct ShiftMultiplier<Result, -20, true> { |
366 | 366 |
static const Result multiplier() { |
367 | 367 |
return static_cast<Result>(1.0/1048576.0); |
368 | 368 |
} |
369 | 369 |
}; |
370 | 370 |
|
371 | 371 |
template <typename Result> |
372 | 372 |
struct ShiftMultiplier<Result, -32, true> { |
373 | 373 |
static const Result multiplier() { |
374 | 374 |
return static_cast<Result>(1.0/424967296.0); |
375 | 375 |
} |
376 | 376 |
}; |
377 | 377 |
|
378 | 378 |
template <typename Result> |
379 | 379 |
struct ShiftMultiplier<Result, -53, true> { |
380 | 380 |
static const Result multiplier() { |
381 | 381 |
return static_cast<Result>(1.0/9007199254740992.0); |
382 | 382 |
} |
383 | 383 |
}; |
384 | 384 |
|
385 | 385 |
template <typename Result> |
386 | 386 |
struct ShiftMultiplier<Result, -64, true> { |
387 | 387 |
static const Result multiplier() { |
388 | 388 |
return static_cast<Result>(1.0/18446744073709551616.0); |
389 | 389 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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_TOLERANCE_H |
20 | 20 |
#define LEMON_TOLERANCE_H |
21 | 21 |
|
22 | 22 |
///\ingroup misc |
23 | 23 |
///\file |
24 | 24 |
///\brief A basic tool to handle the anomalies of calculation with |
25 | 25 |
///floating point numbers. |
26 | 26 |
/// |
27 | 27 |
///\todo It should be in a module like "Basic tools" |
28 | 28 |
|
29 | 29 |
|
30 | 30 |
namespace lemon { |
31 | 31 |
|
32 | 32 |
/// \addtogroup misc |
33 | 33 |
/// @{ |
34 | 34 |
|
35 | 35 |
///\brief A class to provide a basic way to |
36 | 36 |
///handle the comparison of numbers that are obtained |
37 | 37 |
///as a result of a probably inexact computation. |
38 | 38 |
/// |
39 | 39 |
///Tolerance is a class to provide a basic way to |
40 | 40 |
///handle the comparison of numbers that are obtained |
41 | 41 |
///as a result of a probably inexact computation. |
42 | 42 |
/// |
43 | 43 |
///This is an abstract class, it should be specialized for all numerical |
44 | 44 |
///data types. These specialized classes like \ref Tolerance\<double\> |
45 | 45 |
///may offer additional tuning parameters. |
46 | 46 |
/// |
47 | 47 |
///\sa Tolerance<float> |
48 | 48 |
///\sa Tolerance<double> |
49 | 49 |
///\sa Tolerance<long double> |
50 | 50 |
///\sa Tolerance<int> |
51 | 51 |
#if defined __GNUC__ && !defined __STRICT_ANSI__ |
52 | 52 |
///\sa Tolerance<long long int> |
53 | 53 |
#endif |
54 | 54 |
///\sa Tolerance<unsigned int> |
55 | 55 |
#if defined __GNUC__ && !defined __STRICT_ANSI__ |
56 | 56 |
///\sa Tolerance<unsigned long long int> |
57 | 57 |
#endif |
58 | 58 |
|
59 | 59 |
template<class T> |
60 | 60 |
class Tolerance |
61 | 61 |
{ |
62 | 62 |
public: |
63 | 63 |
typedef T Value; |
64 | 64 |
|
65 | 65 |
///\name Comparisons |
66 | 66 |
///The concept is that these bool functions return with \c true only if |
67 | 67 |
///the related comparisons hold even if some numerical error appeared |
68 | 68 |
///during the computations. |
69 | 69 |
|
70 | 70 |
///@{ |
71 | 71 |
|
72 | 72 |
///Returns \c true if \c a is \e surely strictly less than \c b |
73 | 73 |
static bool less(Value a,Value b) {return false;} |
74 | 74 |
///Returns \c true if \c a is \e surely different from \c b |
75 | 75 |
static bool different(Value a,Value b) {return false;} |
76 | 76 |
///Returns \c true if \c a is \e surely positive |
77 | 77 |
static bool positive(Value a) {return false;} |
78 | 78 |
///Returns \c true if \c a is \e surely negative |
79 | 79 |
static bool negative(Value a) {return false;} |
80 | 80 |
///Returns \c true if \c a is \e surely non-zero |
81 | 81 |
static bool nonZero(Value a) {return false;} |
82 | 82 |
|
83 | 83 |
///@} |
84 | 84 |
|
85 | 85 |
///Returns the zero value. |
86 | 86 |
static Value zero() {return T();} |
87 | 87 |
|
88 | 88 |
// static bool finite(Value a) {} |
89 | 89 |
// static Value big() {} |
90 | 90 |
// static Value negativeBig() {} |
91 | 91 |
}; |
92 | 92 |
|
93 | 93 |
|
94 | 94 |
///Float specialization of \ref Tolerance. |
95 | 95 |
|
96 | 96 |
///Float specialization of \ref Tolerance. |
97 | 97 |
///\sa Tolerance |
98 | 98 |
///\relates Tolerance |
99 | 99 |
template<> |
100 | 100 |
class Tolerance<float> |
101 | 101 |
{ |
102 | 102 |
static float def_epsilon; |
103 | 103 |
float _epsilon; |
104 | 104 |
public: |
105 | 105 |
///\e |
106 | 106 |
typedef float Value; |
107 | 107 |
|
108 | 108 |
///Constructor setting the epsilon tolerance to the default value. |
109 | 109 |
Tolerance() : _epsilon(def_epsilon) {} |
110 | 110 |
///Constructor setting the epsilon tolerance. |
111 | 111 |
Tolerance(float e) : _epsilon(e) {} |
112 | 112 |
|
113 | 113 |
///Return the epsilon value. |
114 | 114 |
Value epsilon() const {return _epsilon;} |
115 | 115 |
///Set the epsilon value. |
116 | 116 |
void epsilon(Value e) {_epsilon=e;} |
117 | 117 |
|
118 | 118 |
///Return the default epsilon value. |
119 | 119 |
static Value defaultEpsilon() {return def_epsilon;} |
120 | 120 |
///Set the default epsilon value. |
121 | 121 |
static void defaultEpsilon(Value e) {def_epsilon=e;} |
122 | 122 |
|
123 | 123 |
///\name Comparisons |
124 | 124 |
///See class Tolerance for more details. |
125 | 125 |
|
126 | 126 |
///@{ |
127 | 127 |
|
128 | 128 |
///Returns \c true if \c a is \e surely strictly less than \c b |
129 | 129 |
bool less(Value a,Value b) const {return a+_epsilon<b;} |
130 | 130 |
///Returns \c true if \c a is \e surely different from \c b |
131 | 131 |
bool different(Value a,Value b) const { return less(a,b)||less(b,a); } |
132 | 132 |
///Returns \c true if \c a is \e surely positive |
133 | 133 |
bool positive(Value a) const { return _epsilon<a; } |
134 | 134 |
///Returns \c true if \c a is \e surely negative |
135 | 135 |
bool negative(Value a) const { return -_epsilon>a; } |
136 | 136 |
///Returns \c true if \c a is \e surely non-zero |
137 | 137 |
bool nonZero(Value a) const { return positive(a)||negative(a); } |
138 | 138 |
|
139 | 139 |
///@} |
140 | 140 |
|
141 | 141 |
///Returns zero |
142 | 142 |
static Value zero() {return 0;} |
143 | 143 |
}; |
144 | 144 |
|
145 | 145 |
///Double specialization of \ref Tolerance. |
146 | 146 |
|
147 | 147 |
///Double specialization of \ref Tolerance. |
148 | 148 |
///\sa Tolerance |
149 | 149 |
///\relates Tolerance |
150 | 150 |
template<> |
151 | 151 |
class Tolerance<double> |
152 | 152 |
{ |
153 | 153 |
static double def_epsilon; |
154 | 154 |
double _epsilon; |
155 | 155 |
public: |
156 | 156 |
///\e |
157 | 157 |
typedef double Value; |
158 | 158 |
|
159 | 159 |
///Constructor setting the epsilon tolerance to the default value. |
160 | 160 |
Tolerance() : _epsilon(def_epsilon) {} |
161 | 161 |
///Constructor setting the epsilon tolerance. |
162 | 162 |
Tolerance(double e) : _epsilon(e) {} |
163 | 163 |
|
164 | 164 |
///Return the epsilon value. |
165 | 165 |
Value epsilon() const {return _epsilon;} |
166 | 166 |
///Set the epsilon value. |
167 | 167 |
void epsilon(Value e) {_epsilon=e;} |
168 | 168 |
|
169 | 169 |
///Return the default epsilon value. |
170 | 170 |
static Value defaultEpsilon() {return def_epsilon;} |
171 | 171 |
///Set the default epsilon value. |
172 | 172 |
static void defaultEpsilon(Value e) {def_epsilon=e;} |
173 | 173 |
|
174 | 174 |
///\name Comparisons |
175 | 175 |
///See class Tolerance for more details. |
176 | 176 |
|
177 | 177 |
///@{ |
178 | 178 |
|
179 | 179 |
///Returns \c true if \c a is \e surely strictly less than \c b |
180 | 180 |
bool less(Value a,Value b) const {return a+_epsilon<b;} |
181 | 181 |
///Returns \c true if \c a is \e surely different from \c b |
182 | 182 |
bool different(Value a,Value b) const { return less(a,b)||less(b,a); } |
183 | 183 |
///Returns \c true if \c a is \e surely positive |
184 | 184 |
bool positive(Value a) const { return _epsilon<a; } |
185 | 185 |
///Returns \c true if \c a is \e surely negative |
186 | 186 |
bool negative(Value a) const { return -_epsilon>a; } |
187 | 187 |
///Returns \c true if \c a is \e surely non-zero |
188 | 188 |
bool nonZero(Value a) const { return positive(a)||negative(a); } |
189 | 189 |
|
190 | 190 |
///@} |
191 | 191 |
|
192 | 192 |
///Returns zero |
193 | 193 |
static Value zero() {return 0;} |
194 | 194 |
}; |
195 | 195 |
|
196 | 196 |
///Long double specialization of \ref Tolerance. |
197 | 197 |
|
198 | 198 |
///Long double specialization of \ref Tolerance. |
199 | 199 |
///\sa Tolerance |
200 | 200 |
///\relates Tolerance |
201 | 201 |
template<> |
202 | 202 |
class Tolerance<long double> |
203 | 203 |
{ |
204 | 204 |
static long double def_epsilon; |
205 | 205 |
long double _epsilon; |
206 | 206 |
public: |
207 | 207 |
///\e |
208 | 208 |
typedef long double Value; |
209 | 209 |
|
210 | 210 |
///Constructor setting the epsilon tolerance to the default value. |
211 | 211 |
Tolerance() : _epsilon(def_epsilon) {} |
212 | 212 |
///Constructor setting the epsilon tolerance. |
213 | 213 |
Tolerance(long double e) : _epsilon(e) {} |
214 | 214 |
|
215 | 215 |
///Return the epsilon value. |
216 | 216 |
Value epsilon() const {return _epsilon;} |
217 | 217 |
///Set the epsilon value. |
218 | 218 |
void epsilon(Value e) {_epsilon=e;} |
219 | 219 |
|
220 | 220 |
///Return the default epsilon value. |
221 | 221 |
static Value defaultEpsilon() {return def_epsilon;} |
222 | 222 |
///Set the default epsilon value. |
223 | 223 |
static void defaultEpsilon(Value e) {def_epsilon=e;} |
224 | 224 |
|
225 | 225 |
///\name Comparisons |
226 | 226 |
///See class Tolerance for more details. |
227 | 227 |
|
228 | 228 |
///@{ |
229 | 229 |
|
230 | 230 |
///Returns \c true if \c a is \e surely strictly less than \c b |
231 | 231 |
bool less(Value a,Value b) const {return a+_epsilon<b;} |
232 | 232 |
///Returns \c true if \c a is \e surely different from \c b |
233 | 233 |
bool different(Value a,Value b) const { return less(a,b)||less(b,a); } |
234 | 234 |
///Returns \c true if \c a is \e surely positive |
235 | 235 |
bool positive(Value a) const { return _epsilon<a; } |
236 | 236 |
///Returns \c true if \c a is \e surely negative |
237 | 237 |
bool negative(Value a) const { return -_epsilon>a; } |
238 | 238 |
///Returns \c true if \c a is \e surely non-zero |
239 | 239 |
bool nonZero(Value a) const { return positive(a)||negative(a); } |
240 | 240 |
|
241 | 241 |
///@} |
242 | 242 |
|
243 | 243 |
///Returns zero |
244 | 244 |
static Value zero() {return 0;} |
245 | 245 |
}; |
246 | 246 |
|
247 | 247 |
///Integer specialization of \ref Tolerance. |
248 | 248 |
|
249 | 249 |
///Integer specialization of \ref Tolerance. |
250 | 250 |
///\sa Tolerance |
251 | 251 |
template<> |
252 | 252 |
class Tolerance<int> |
253 | 253 |
{ |
254 | 254 |
public: |
255 | 255 |
///\e |
256 | 256 |
typedef int Value; |
257 | 257 |
|
258 | 258 |
///\name Comparisons |
259 | 259 |
///See \ref Tolerance for more details. |
260 | 260 |
|
261 | 261 |
///@{ |
262 | 262 |
|
263 | 263 |
///Returns \c true if \c a is \e surely strictly less than \c b |
264 | 264 |
static bool less(Value a,Value b) { return a<b;} |
265 | 265 |
///Returns \c true if \c a is \e surely different from \c b |
266 | 266 |
static bool different(Value a,Value b) { return a!=b; } |
267 | 267 |
///Returns \c true if \c a is \e surely positive |
268 | 268 |
static bool positive(Value a) { return 0<a; } |
269 | 269 |
///Returns \c true if \c a is \e surely negative |
270 | 270 |
static bool negative(Value a) { return 0>a; } |
271 | 271 |
///Returns \c true if \c a is \e surely non-zero |
272 | 272 |
static bool nonZero(Value a) { return a!=0; } |
273 | 273 |
|
274 | 274 |
///@} |
275 | 275 |
|
276 | 276 |
///Returns zero |
277 | 277 |
static Value zero() {return 0;} |
278 | 278 |
}; |
279 | 279 |
|
280 | 280 |
///Unsigned integer specialization of \ref Tolerance. |
281 | 281 |
|
282 | 282 |
///Unsigned integer specialization of \ref Tolerance. |
283 | 283 |
///\sa Tolerance |
284 | 284 |
template<> |
285 | 285 |
class Tolerance<unsigned int> |
286 | 286 |
{ |
287 | 287 |
public: |
288 | 288 |
///\e |
289 | 289 |
typedef unsigned int Value; |
290 | 290 |
|
291 | 291 |
///\name Comparisons |
292 | 292 |
///See \ref Tolerance for more details. |
293 | 293 |
|
294 | 294 |
///@{ |
295 | 295 |
|
296 | 296 |
///Returns \c true if \c a is \e surely strictly less than \c b |
297 | 297 |
static bool less(Value a,Value b) { return a<b;} |
298 | 298 |
///Returns \c true if \c a is \e surely different from \c b |
299 | 299 |
static bool different(Value a,Value b) { return a!=b; } |
300 | 300 |
///Returns \c true if \c a is \e surely positive |
301 | 301 |
static bool positive(Value a) { return 0<a; } |
302 | 302 |
///Returns \c true if \c a is \e surely negative |
303 | 303 |
static bool negative(Value) { return false; } |
304 | 304 |
///Returns \c true if \c a is \e surely non-zero |
305 | 305 |
static bool nonZero(Value a) { return a!=0; } |
306 | 306 |
|
307 | 307 |
///@} |
308 | 308 |
|
309 | 309 |
///Returns zero |
310 | 310 |
static Value zero() {return 0;} |
311 | 311 |
}; |
312 | 312 |
|
313 | 313 |
|
314 | 314 |
///Long integer specialization of \ref Tolerance. |
315 | 315 |
|
316 | 316 |
///Long integer specialization of \ref Tolerance. |
317 | 317 |
///\sa Tolerance |
318 | 318 |
template<> |
319 | 319 |
class Tolerance<long int> |
320 | 320 |
{ |
321 | 321 |
public: |
322 | 322 |
///\e |
323 | 323 |
typedef long int Value; |
324 | 324 |
|
325 | 325 |
///\name Comparisons |
326 | 326 |
///See \ref Tolerance for more details. |
327 | 327 |
|
328 | 328 |
///@{ |
329 | 329 |
|
330 | 330 |
///Returns \c true if \c a is \e surely strictly less than \c b |
331 | 331 |
static bool less(Value a,Value b) { return a<b;} |
332 | 332 |
///Returns \c true if \c a is \e surely different from \c b |
333 | 333 |
static bool different(Value a,Value b) { return a!=b; } |
334 | 334 |
///Returns \c true if \c a is \e surely positive |
335 | 335 |
static bool positive(Value a) { return 0<a; } |
336 | 336 |
///Returns \c true if \c a is \e surely negative |
337 | 337 |
static bool negative(Value a) { return 0>a; } |
338 | 338 |
///Returns \c true if \c a is \e surely non-zero |
339 | 339 |
static bool nonZero(Value a) { return a!=0;} |
340 | 340 |
|
341 | 341 |
///@} |
342 | 342 |
|
343 | 343 |
///Returns zero |
344 | 344 |
static Value zero() {return 0;} |
345 | 345 |
}; |
346 | 346 |
|
347 | 347 |
///Unsigned long integer specialization of \ref Tolerance. |
348 | 348 |
|
349 | 349 |
///Unsigned long integer specialization of \ref Tolerance. |
350 | 350 |
///\sa Tolerance |
351 | 351 |
template<> |
352 | 352 |
class Tolerance<unsigned long int> |
353 | 353 |
{ |
354 | 354 |
public: |
355 | 355 |
///\e |
356 | 356 |
typedef unsigned long int Value; |
357 | 357 |
|
358 | 358 |
///\name Comparisons |
359 | 359 |
///See \ref Tolerance for more details. |
360 | 360 |
|
361 | 361 |
///@{ |
362 | 362 |
|
363 | 363 |
///Returns \c true if \c a is \e surely strictly less than \c b |
364 | 364 |
static bool less(Value a,Value b) { return a<b;} |
365 | 365 |
///Returns \c true if \c a is \e surely different from \c b |
366 | 366 |
static bool different(Value a,Value b) { return a!=b; } |
367 | 367 |
///Returns \c true if \c a is \e surely positive |
368 | 368 |
static bool positive(Value a) { return 0<a; } |
369 | 369 |
///Returns \c true if \c a is \e surely negative |
370 | 370 |
static bool negative(Value) { return false; } |
371 | 371 |
///Returns \c true if \c a is \e surely non-zero |
372 | 372 |
static bool nonZero(Value a) { return a!=0;} |
373 | 373 |
|
374 | 374 |
///@} |
375 | 375 |
|
376 | 376 |
///Returns zero |
377 | 377 |
static Value zero() {return 0;} |
378 | 378 |
}; |
379 | 379 |
|
380 | 380 |
#if defined __GNUC__ && !defined __STRICT_ANSI__ |
381 | 381 |
|
382 | 382 |
///Long long integer specialization of \ref Tolerance. |
383 | 383 |
|
384 | 384 |
///Long long integer specialization of \ref Tolerance. |
385 | 385 |
///\warning This class (more exactly, type <tt>long long</tt>) |
386 | 386 |
///is not ansi compatible. |
387 | 387 |
///\sa Tolerance |
388 | 388 |
template<> |
389 | 389 |
class Tolerance<long long int> |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
#include <lemon/dim2.h> |
20 | 20 |
#include <iostream> |
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 |
cout << "Testing classes 'dim2::Point' and 'dim2::BoundingBox'." << endl; |
29 | 29 |
|
30 | 30 |
typedef dim2::Point<int> Point; |
31 | 31 |
|
32 | 32 |
Point p; |
33 | 33 |
check(p.size()==2, "Wrong vector initialization."); |
34 | 34 |
|
35 | 35 |
Point a(1,2); |
36 | 36 |
Point b(3,4); |
37 | 37 |
check(a[0]==1 && a[1]==2, "Wrong vector initialization."); |
38 | 38 |
|
39 | 39 |
p = a+b; |
40 | 40 |
check(p.x==4 && p.y==6, "Wrong vector addition."); |
41 | 41 |
|
42 | 42 |
p = a-b; |
43 | 43 |
check(p.x==-2 && p.y==-2, "Wrong vector subtraction."); |
44 | 44 |
|
45 | 45 |
check(a.normSquare()==5,"Wrong vector norm calculation."); |
46 | 46 |
check(a*b==11, "Wrong vector scalar product."); |
47 | 47 |
|
48 | 48 |
int l=2; |
49 | 49 |
p = a*l; |
50 | 50 |
check(p.x==2 && p.y==4, "Wrong vector multiplication by a scalar."); |
51 | 51 |
|
52 | 52 |
p = b/l; |
53 | 53 |
check(p.x==1 && p.y==2, "Wrong vector division by a scalar."); |
54 | 54 |
|
55 | 55 |
typedef dim2::BoundingBox<int> BB; |
56 | 56 |
BB box1; |
57 | 57 |
check(box1.empty(), "It should be empty."); |
58 | 58 |
|
59 | 59 |
box1.add(a); |
60 | 60 |
check(!box1.empty(), "It should not be empty."); |
61 | 61 |
box1.add(b); |
62 | 62 |
|
63 | 63 |
check(box1.bottomLeft().x==1 && |
64 | 64 |
box1.bottomLeft().y==2 && |
65 | 65 |
box1.topRight().x==3 && |
66 | 66 |
box1.topRight().y==4, |
67 | 67 |
"Wrong addition of points to box."); |
68 | 68 |
|
69 | 69 |
p.x=2; p.y=3; |
70 | 70 |
check(box1.inside(p), "It should be inside."); |
71 | 71 |
|
72 | 72 |
p.x=1; p.y=3; |
73 | 73 |
check(box1.inside(p), "It should be inside."); |
74 | 74 |
|
75 | 75 |
p.x=0; p.y=3; |
76 | 76 |
check(!box1.inside(p), "It should not be inside."); |
77 | 77 |
|
78 | 78 |
BB box2(p); |
79 | 79 |
check(!box2.empty(), |
80 | 80 |
"It should not be empty. Constructed from 1 point."); |
81 | 81 |
|
82 | 82 |
box2.add(box1); |
83 | 83 |
check(box2.inside(p), |
84 | 84 |
"It should be inside. Incremented a box with another one."); |
85 | 85 |
|
86 | 86 |
return 0; |
87 | 87 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
#include <deque> |
20 | 20 |
#include <set> |
21 | 21 |
|
22 | 22 |
#include <lemon/concept_check.h> |
23 | 23 |
#include <lemon/concepts/maps.h> |
24 | 24 |
#include <lemon/maps.h> |
25 | 25 |
|
26 | 26 |
#include "test_tools.h" |
27 | 27 |
|
28 | 28 |
using namespace lemon; |
29 | 29 |
using namespace lemon::concepts; |
30 | 30 |
|
31 | 31 |
struct A {}; |
32 | 32 |
inline bool operator<(A, A) { return true; } |
33 | 33 |
struct B {}; |
34 | 34 |
|
35 | 35 |
class F { |
36 | 36 |
public: |
37 | 37 |
typedef A argument_type; |
38 | 38 |
typedef B result_type; |
39 | 39 |
|
40 | 40 |
B operator()(const A &) const {return B();} |
41 | 41 |
}; |
42 | 42 |
|
43 | 43 |
int func(A) {return 3;} |
44 | 44 |
|
45 | 45 |
int binc(int, B) {return 4;} |
46 | 46 |
|
47 | 47 |
typedef ReadMap<A,double> DoubleMap; |
48 | 48 |
typedef ReadWriteMap<A, double> WriteDoubleMap; |
49 | 49 |
|
50 | 50 |
typedef ReadMap<A,bool> BoolMap; |
51 | 51 |
typedef ReadWriteMap<A, bool> BoolWriteMap; |
52 | 52 |
|
53 | 53 |
int main() |
54 | 54 |
{ // checking graph components |
55 | 55 |
|
56 | 56 |
checkConcept<ReadMap<A,B>, ReadMap<A,B> >(); |
57 | 57 |
checkConcept<WriteMap<A,B>, WriteMap<A,B> >(); |
58 | 58 |
checkConcept<ReadWriteMap<A,B>, ReadWriteMap<A,B> >(); |
59 | 59 |
checkConcept<ReferenceMap<A,B,B&,const B&>, ReferenceMap<A,B,B&,const B&> >(); |
60 | 60 |
|
61 | 61 |
checkConcept<ReadMap<A,double>, AddMap<DoubleMap,DoubleMap> >(); |
62 | 62 |
checkConcept<ReadMap<A,double>, SubMap<DoubleMap,DoubleMap> >(); |
63 | 63 |
checkConcept<ReadMap<A,double>, MulMap<DoubleMap,DoubleMap> >(); |
64 | 64 |
checkConcept<ReadMap<A,double>, DivMap<DoubleMap,DoubleMap> >(); |
65 | 65 |
checkConcept<ReadMap<A,double>, NegMap<DoubleMap> >(); |
66 | 66 |
checkConcept<ReadWriteMap<A,double>, NegWriteMap<WriteDoubleMap> >(); |
67 | 67 |
checkConcept<ReadMap<A,double>, AbsMap<DoubleMap> >(); |
68 | 68 |
checkConcept<ReadMap<A,double>, ShiftMap<DoubleMap> >(); |
69 | 69 |
checkConcept<ReadWriteMap<A,double>, ShiftWriteMap<WriteDoubleMap> >(); |
70 | 70 |
checkConcept<ReadMap<A,double>, ScaleMap<DoubleMap> >(); |
71 | 71 |
checkConcept<ReadWriteMap<A,double>, ScaleWriteMap<WriteDoubleMap> >(); |
72 | 72 |
checkConcept<ReadMap<A,double>, ForkMap<DoubleMap, DoubleMap> >(); |
73 | 73 |
checkConcept<ReadWriteMap<A,double>, |
74 | 74 |
ForkWriteMap<WriteDoubleMap, WriteDoubleMap> >(); |
75 | 75 |
|
76 | 76 |
checkConcept<ReadMap<B,double>, ComposeMap<DoubleMap,ReadMap<B,A> > >(); |
77 | 77 |
|
78 | 78 |
checkConcept<ReadMap<A,B>, FunctorMap<F, A, B> >(); |
79 | 79 |
|
80 | 80 |
checkConcept<ReadMap<A, bool>, NotMap<BoolMap> >(); |
81 | 81 |
checkConcept<ReadWriteMap<A, bool>, NotWriteMap<BoolWriteMap> >(); |
82 | 82 |
|
83 | 83 |
checkConcept<WriteMap<A, bool>, StoreBoolMap<A*> >(); |
84 | 84 |
checkConcept<WriteMap<A, bool>, BackInserterBoolMap<std::deque<A> > >(); |
85 | 85 |
checkConcept<WriteMap<A, bool>, FrontInserterBoolMap<std::deque<A> > >(); |
86 | 86 |
checkConcept<WriteMap<A, bool>, InserterBoolMap<std::set<A> > >(); |
87 | 87 |
checkConcept<WriteMap<A, bool>, FillBoolMap<WriteMap<A, B> > >(); |
88 | 88 |
checkConcept<WriteMap<A, bool>, SettingOrderBoolMap<WriteMap<A, int> > >(); |
89 | 89 |
|
90 | 90 |
int a; |
91 | 91 |
|
92 | 92 |
a=mapFunctor(constMap<A,int>(2))(A()); |
93 | 93 |
check(a==2,"Something is wrong with mapFunctor"); |
94 | 94 |
|
95 | 95 |
B b; |
96 | 96 |
b=functorMap(F())[A()]; |
97 | 97 |
|
98 | 98 |
a=functorMap(&func)[A()]; |
99 | 99 |
check(a==3,"Something is wrong with functorMap"); |
100 | 100 |
|
101 | 101 |
a=combineMap(constMap<B, int, 1>(), identityMap<B>(), &binc)[B()]; |
102 | 102 |
check(a==4,"Something is wrong with combineMap"); |
103 | 103 |
|
104 | 104 |
|
105 | 105 |
std::cout << __FILE__ ": All tests passed.\n"; |
106 | 106 |
|
107 | 107 |
return 0; |
108 | 108 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
#include <lemon/random.h> |
20 | 20 |
#include "test_tools.h" |
21 | 21 |
|
22 | 22 |
///\file \brief Test cases for random.h |
23 | 23 |
/// |
24 | 24 |
///\todo To be extended |
25 | 25 |
/// |
26 | 26 |
|
27 | 27 |
int main() |
28 | 28 |
{ |
29 | 29 |
double a=lemon::rnd(); |
30 | 30 |
check(a<1.0&&a>0.0,"This should be in [0,1)"); |
31 | 31 |
a=lemon::rnd.gauss(); |
32 | 32 |
a=lemon::rnd.gamma(3.45,0); |
33 | 33 |
a=lemon::rnd.gamma(4); |
34 | 34 |
//Does gamma work with integer k? |
35 | 35 |
a=lemon::rnd.gamma(4.0,0); |
36 | 36 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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_TEST_TEST_TOOLS_H |
20 | 20 |
#define LEMON_TEST_TEST_TOOLS_H |
21 | 21 |
|
22 | 22 |
#include <iostream> |
23 | 23 |
|
24 | 24 |
//! \ingroup misc |
25 | 25 |
//! \file |
26 | 26 |
//! \brief Some utilities to write test programs. |
27 | 27 |
|
28 | 28 |
|
29 | 29 |
///If \c rc is fail, writes an error message end exit. |
30 | 30 |
|
31 | 31 |
///If \c rc is fail, writes an error message end exit. |
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 this (and then exits). |
39 | 39 |
///\verbatim graph_test.cc:123: error: This is obviously false. \endverbatim |
40 | 40 |
/// |
41 | 41 |
///\todo It should be in \c error.h |
42 | 42 |
#define check(rc, msg) \ |
43 | 43 |
if(!(rc)) { \ |
44 | 44 |
std::cerr << __FILE__ ":" << __LINE__ << ": error: " << msg << std::endl; \ |
45 | 45 |
abort(); \ |
46 | 46 |
} else { } \ |
47 | 47 |
|
48 | 48 |
#endif |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
#include "test_tools.h" |
20 | 20 |
|
21 | 21 |
int main() |
22 | 22 |
{ |
23 | 23 |
check(false, "Don't panic. Failing is the right behaviour here."); |
24 | 24 |
return 0; |
25 | 25 |
} |
1 | 1 |
/* -*- C++ -*- |
2 | 2 |
* |
3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
4 | 4 |
* |
5 |
* Copyright (C) 2003- |
|
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 |
#include "test_tools.h" |
20 | 20 |
|
21 | 21 |
int main() |
22 | 22 |
{ |
23 | 23 |
check(true, "It should pass."); |
24 | 24 |
return 0; |
25 | 25 |
} |
0 comments (0 inline)