... | ... |
@@ -167,104 +167,114 @@ |
167 | 167 |
typename enable_if<MapInputIndicator<In>, void>::type> |
168 | 168 |
{ |
169 | 169 |
typedef typename In::Value Value; |
170 | 170 |
}; |
171 | 171 |
|
172 | 172 |
template <typename Graph, typename In, typename Out, |
173 | 173 |
typename InEnable = void> |
174 | 174 |
struct KruskalInputSelector {}; |
175 | 175 |
|
176 | 176 |
template <typename Graph, typename In, typename Out, |
177 | 177 |
typename InEnable = void> |
178 | 178 |
struct KruskalOutputSelector {}; |
179 | 179 |
|
180 | 180 |
template <typename Graph, typename In, typename Out> |
181 | 181 |
struct KruskalInputSelector<Graph, In, Out, |
182 | 182 |
typename enable_if<SequenceInputIndicator<In>, void>::type > |
183 | 183 |
{ |
184 | 184 |
typedef typename In::value_type::second_type Value; |
185 | 185 |
|
186 | 186 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
187 | 187 |
return KruskalOutputSelector<Graph, In, Out>:: |
188 | 188 |
kruskal(graph, in, out); |
189 | 189 |
} |
190 | 190 |
|
191 | 191 |
}; |
192 | 192 |
|
193 | 193 |
template <typename Graph, typename In, typename Out> |
194 | 194 |
struct KruskalInputSelector<Graph, In, Out, |
195 | 195 |
typename enable_if<MapInputIndicator<In>, void>::type > |
196 | 196 |
{ |
197 | 197 |
typedef typename In::Value Value; |
198 | 198 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
199 | 199 |
typedef typename In::Key MapArc; |
200 | 200 |
typedef typename In::Value Value; |
201 | 201 |
typedef typename ItemSetTraits<Graph, MapArc>::ItemIt MapArcIt; |
202 | 202 |
typedef std::vector<std::pair<MapArc, Value> > Sequence; |
203 | 203 |
Sequence seq; |
204 | 204 |
|
205 | 205 |
for (MapArcIt it(graph); it != INVALID; ++it) { |
206 | 206 |
seq.push_back(std::make_pair(it, in[it])); |
207 | 207 |
} |
208 | 208 |
|
209 | 209 |
std::sort(seq.begin(), seq.end(), PairComp<Sequence>()); |
210 | 210 |
return KruskalOutputSelector<Graph, Sequence, Out>:: |
211 | 211 |
kruskal(graph, seq, out); |
212 | 212 |
} |
213 | 213 |
}; |
214 | 214 |
|
215 |
template <typename T> |
|
216 |
struct RemoveConst { |
|
217 |
typedef T type; |
|
218 |
}; |
|
219 |
|
|
220 |
template <typename T> |
|
221 |
struct RemoveConst<const T> { |
|
222 |
typedef T type; |
|
223 |
}; |
|
224 |
|
|
215 | 225 |
template <typename Graph, typename In, typename Out> |
216 | 226 |
struct KruskalOutputSelector<Graph, In, Out, |
217 | 227 |
typename enable_if<SequenceOutputIndicator<Out>, void>::type > |
218 | 228 |
{ |
219 | 229 |
typedef typename In::value_type::second_type Value; |
220 | 230 |
|
221 | 231 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
222 |
typedef StoreBoolMap<Out> Map; |
|
232 |
typedef StoreBoolMap<typename RemoveConst<Out>::type> Map; |
|
223 | 233 |
Map map(out); |
224 | 234 |
return _kruskal_bits::kruskal(graph, in, map); |
225 | 235 |
} |
226 | 236 |
|
227 | 237 |
}; |
228 | 238 |
|
229 | 239 |
template <typename Graph, typename In, typename Out> |
230 | 240 |
struct KruskalOutputSelector<Graph, In, Out, |
231 | 241 |
typename enable_if<MapOutputIndicator<Out>, void>::type > |
232 | 242 |
{ |
233 | 243 |
typedef typename In::value_type::second_type Value; |
234 | 244 |
|
235 | 245 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
236 | 246 |
return _kruskal_bits::kruskal(graph, in, out); |
237 | 247 |
} |
238 | 248 |
}; |
239 | 249 |
|
240 | 250 |
} |
241 | 251 |
|
242 | 252 |
/// \ingroup spantree |
243 | 253 |
/// |
244 | 254 |
/// \brief Kruskal's algorithm to find a minimum cost tree of a graph. |
245 | 255 |
/// |
246 | 256 |
/// This function runs Kruskal's algorithm to find a minimum cost tree. |
247 | 257 |
/// Due to some C++ hacking, it accepts various input and output types. |
248 | 258 |
/// |
249 | 259 |
/// \param g The graph the algorithm runs on. |
250 | 260 |
/// It can be either \ref concepts::Digraph "directed" or |
251 | 261 |
/// \ref concepts::Graph "undirected". |
252 | 262 |
/// If the graph is directed, the algorithm consider it to be |
253 | 263 |
/// undirected by disregarding the direction of the arcs. |
254 | 264 |
/// |
255 | 265 |
/// \param in This object is used to describe the arc costs. It can be one |
256 | 266 |
/// of the following choices. |
257 | 267 |
/// - An STL compatible 'Forward Container' with |
258 | 268 |
/// <tt>std::pair<GR::Edge,X></tt> or |
259 | 269 |
/// <tt>std::pair<GR::Arc,X></tt> as its <tt>value_type</tt>, where |
260 | 270 |
/// \c X is the type of the costs. The pairs indicates the arcs |
261 | 271 |
/// along with the assigned cost. <em>They must be in a |
262 | 272 |
/// cost-ascending order.</em> |
263 | 273 |
/// - Any readable Arc map. The values of the map indicate the arc costs. |
264 | 274 |
/// |
265 | 275 |
/// \retval out Here we also have a choise. |
266 | 276 |
/// - It can be a writable \c bool arc map. After running the |
267 | 277 |
/// algorithm this will contain the found minimum cost spanning |
268 | 278 |
/// tree: the value of an arc will be set to \c true if it belongs |
269 | 279 |
/// to the tree, otherwise it will be set to \c false. The value of |
270 | 280 |
/// each arc will be set exactly once. |
... | ... |
@@ -40,155 +40,155 @@ |
40 | 40 |
|
41 | 41 |
class F { |
42 | 42 |
public: |
43 | 43 |
typedef A argument_type; |
44 | 44 |
typedef B result_type; |
45 | 45 |
|
46 | 46 |
B operator()(const A&) const { return B(); } |
47 | 47 |
private: |
48 | 48 |
F& operator=(const F&); |
49 | 49 |
}; |
50 | 50 |
|
51 | 51 |
int func(A) { return 3; } |
52 | 52 |
|
53 | 53 |
int binc(int a, B) { return a+1; } |
54 | 54 |
|
55 | 55 |
typedef ReadMap<A, double> DoubleMap; |
56 | 56 |
typedef ReadWriteMap<A, double> DoubleWriteMap; |
57 | 57 |
typedef ReferenceMap<A, double, double&, const double&> DoubleRefMap; |
58 | 58 |
|
59 | 59 |
typedef ReadMap<A, bool> BoolMap; |
60 | 60 |
typedef ReadWriteMap<A, bool> BoolWriteMap; |
61 | 61 |
typedef ReferenceMap<A, bool, bool&, const bool&> BoolRefMap; |
62 | 62 |
|
63 | 63 |
int main() |
64 | 64 |
{ |
65 | 65 |
// Map concepts |
66 | 66 |
checkConcept<ReadMap<A,B>, ReadMap<A,B> >(); |
67 | 67 |
checkConcept<ReadMap<A,C>, ReadMap<A,C> >(); |
68 | 68 |
checkConcept<WriteMap<A,B>, WriteMap<A,B> >(); |
69 | 69 |
checkConcept<WriteMap<A,C>, WriteMap<A,C> >(); |
70 | 70 |
checkConcept<ReadWriteMap<A,B>, ReadWriteMap<A,B> >(); |
71 | 71 |
checkConcept<ReadWriteMap<A,C>, ReadWriteMap<A,C> >(); |
72 | 72 |
checkConcept<ReferenceMap<A,B,B&,const B&>, ReferenceMap<A,B,B&,const B&> >(); |
73 | 73 |
checkConcept<ReferenceMap<A,C,C&,const C&>, ReferenceMap<A,C,C&,const C&> >(); |
74 | 74 |
|
75 | 75 |
// NullMap |
76 | 76 |
{ |
77 | 77 |
checkConcept<ReadWriteMap<A,B>, NullMap<A,B> >(); |
78 | 78 |
NullMap<A,B> map1; |
79 | 79 |
NullMap<A,B> map2 = map1; |
80 | 80 |
map1 = nullMap<A,B>(); |
81 | 81 |
} |
82 | 82 |
|
83 | 83 |
// ConstMap |
84 | 84 |
{ |
85 | 85 |
checkConcept<ReadWriteMap<A,B>, ConstMap<A,B> >(); |
86 | 86 |
checkConcept<ReadWriteMap<A,C>, ConstMap<A,C> >(); |
87 | 87 |
ConstMap<A,B> map1; |
88 |
ConstMap<A,B> map2 |
|
88 |
ConstMap<A,B> map2 = B(); |
|
89 | 89 |
ConstMap<A,B> map3 = map1; |
90 | 90 |
map1 = constMap<A>(B()); |
91 | 91 |
map1 = constMap<A,B>(); |
92 | 92 |
map1.setAll(B()); |
93 | 93 |
ConstMap<A,C> map4(C(1)); |
94 | 94 |
ConstMap<A,C> map5 = map4; |
95 | 95 |
map4 = constMap<A>(C(2)); |
96 | 96 |
map4.setAll(C(3)); |
97 | 97 |
|
98 | 98 |
checkConcept<ReadWriteMap<A,int>, ConstMap<A,int> >(); |
99 | 99 |
check(constMap<A>(10)[A()] == 10, "Something is wrong with ConstMap"); |
100 | 100 |
|
101 | 101 |
checkConcept<ReadWriteMap<A,int>, ConstMap<A,Const<int,10> > >(); |
102 | 102 |
ConstMap<A,Const<int,10> > map6; |
103 | 103 |
ConstMap<A,Const<int,10> > map7 = map6; |
104 | 104 |
map6 = constMap<A,int,10>(); |
105 | 105 |
map7 = constMap<A,Const<int,10> >(); |
106 | 106 |
check(map6[A()] == 10 && map7[A()] == 10, "Something is wrong with ConstMap"); |
107 | 107 |
} |
108 | 108 |
|
109 | 109 |
// IdentityMap |
110 | 110 |
{ |
111 | 111 |
checkConcept<ReadMap<A,A>, IdentityMap<A> >(); |
112 | 112 |
IdentityMap<A> map1; |
113 | 113 |
IdentityMap<A> map2 = map1; |
114 | 114 |
map1 = identityMap<A>(); |
115 | 115 |
|
116 | 116 |
checkConcept<ReadMap<double,double>, IdentityMap<double> >(); |
117 | 117 |
check(identityMap<double>()[1.0] == 1.0 && identityMap<double>()[3.14] == 3.14, |
118 | 118 |
"Something is wrong with IdentityMap"); |
119 | 119 |
} |
120 | 120 |
|
121 | 121 |
// RangeMap |
122 | 122 |
{ |
123 | 123 |
checkConcept<ReferenceMap<int,B,B&,const B&>, RangeMap<B> >(); |
124 | 124 |
RangeMap<B> map1; |
125 | 125 |
RangeMap<B> map2(10); |
126 | 126 |
RangeMap<B> map3(10,B()); |
127 | 127 |
RangeMap<B> map4 = map1; |
128 | 128 |
RangeMap<B> map5 = rangeMap<B>(); |
129 | 129 |
RangeMap<B> map6 = rangeMap<B>(10); |
130 | 130 |
RangeMap<B> map7 = rangeMap(10,B()); |
131 | 131 |
|
132 | 132 |
checkConcept< ReferenceMap<int, double, double&, const double&>, |
133 | 133 |
RangeMap<double> >(); |
134 | 134 |
std::vector<double> v(10, 0); |
135 | 135 |
v[5] = 100; |
136 | 136 |
RangeMap<double> map8(v); |
137 | 137 |
RangeMap<double> map9 = rangeMap(v); |
138 | 138 |
check(map9.size() == 10 && map9[2] == 0 && map9[5] == 100, |
139 | 139 |
"Something is wrong with RangeMap"); |
140 | 140 |
} |
141 | 141 |
|
142 | 142 |
// SparseMap |
143 | 143 |
{ |
144 | 144 |
checkConcept<ReferenceMap<A,B,B&,const B&>, SparseMap<A,B> >(); |
145 | 145 |
SparseMap<A,B> map1; |
146 |
SparseMap<A,B> map2 |
|
146 |
SparseMap<A,B> map2 = B(); |
|
147 | 147 |
SparseMap<A,B> map3 = sparseMap<A,B>(); |
148 | 148 |
SparseMap<A,B> map4 = sparseMap<A>(B()); |
149 | 149 |
|
150 | 150 |
checkConcept< ReferenceMap<double, int, int&, const int&>, |
151 | 151 |
SparseMap<double, int> >(); |
152 | 152 |
std::map<double, int> m; |
153 | 153 |
SparseMap<double, int> map5(m); |
154 | 154 |
SparseMap<double, int> map6(m,10); |
155 | 155 |
SparseMap<double, int> map7 = sparseMap(m); |
156 | 156 |
SparseMap<double, int> map8 = sparseMap(m,10); |
157 | 157 |
|
158 | 158 |
check(map5[1.0] == 0 && map5[3.14] == 0 && map6[1.0] == 10 && map6[3.14] == 10, |
159 | 159 |
"Something is wrong with SparseMap"); |
160 | 160 |
map5[1.0] = map6[3.14] = 100; |
161 | 161 |
check(map5[1.0] == 100 && map5[3.14] == 0 && map6[1.0] == 10 && map6[3.14] == 100, |
162 | 162 |
"Something is wrong with SparseMap"); |
163 | 163 |
} |
164 | 164 |
|
165 | 165 |
// ComposeMap |
166 | 166 |
{ |
167 | 167 |
typedef ComposeMap<DoubleMap, ReadMap<B,A> > CompMap; |
168 | 168 |
checkConcept<ReadMap<B,double>, CompMap>(); |
169 | 169 |
CompMap map1(DoubleMap(),ReadMap<B,A>()); |
170 | 170 |
CompMap map2 = composeMap(DoubleMap(), ReadMap<B,A>()); |
171 | 171 |
|
172 | 172 |
SparseMap<double, bool> m1(false); m1[3.14] = true; |
173 | 173 |
RangeMap<double> m2(2); m2[0] = 3.0; m2[1] = 3.14; |
174 | 174 |
check(!composeMap(m1,m2)[0] && composeMap(m1,m2)[1], "Something is wrong with ComposeMap") |
175 | 175 |
} |
176 | 176 |
|
177 | 177 |
// CombineMap |
178 | 178 |
{ |
179 | 179 |
typedef CombineMap<DoubleMap, DoubleMap, std::plus<double> > CombMap; |
180 | 180 |
checkConcept<ReadMap<A,double>, CombMap>(); |
181 | 181 |
CombMap map1(DoubleMap(), DoubleMap()); |
182 | 182 |
CombMap map2 = combineMap(DoubleMap(), DoubleMap(), std::plus<double>()); |
183 | 183 |
|
184 | 184 |
check(combineMap(constMap<B,int,2>(), identityMap<B>(), &binc)[B()] == 3, |
185 | 185 |
"Something is wrong with CombineMap"); |
186 | 186 |
} |
187 | 187 |
|
188 | 188 |
// FunctorToMap, MapToFunctor |
189 | 189 |
{ |
190 | 190 |
checkConcept<ReadMap<A,B>, FunctorToMap<F,A,B> >(); |
191 | 191 |
checkConcept<ReadMap<A,B>, FunctorToMap<F> >(); |
192 | 192 |
FunctorToMap<F> map1; |
193 | 193 |
FunctorToMap<F> map2(F()); |
194 | 194 |
B b = functorToMap(F())[A()]; |
0 comments (0 inline)