1 //=======================================================================
2 // Copyright 2007 Aaron Windsor
4 // Distributed under the Boost Software License, Version 1.0. (See
5 // accompanying file LICENSE_1_0.txt or copy at
6 // http://www.boost.org/LICENSE_1_0.txt)
7 //=======================================================================
8 #ifndef __MAKE_CONNECTED_HPP__
9 #define __MAKE_CONNECTED_HPP__
11 #include <boost/config.hpp>
12 #include <boost/next_prior.hpp>
13 #include <boost/tuple/tuple.hpp> //for tie
14 #include <boost/graph/connected_components.hpp>
15 #include <boost/property_map/property_map.hpp>
18 #include <boost/graph/planar_detail/add_edge_visitors.hpp>
19 #include <boost/graph/planar_detail/bucket_sort.hpp>
26 template <typename Graph,
27 typename VertexIndexMap,
28 typename AddEdgeVisitor
30 void make_connected(Graph& g, VertexIndexMap vm, AddEdgeVisitor& vis)
32 typedef typename graph_traits<Graph>::vertex_iterator vertex_iterator_t;
33 typedef typename graph_traits<Graph>::vertex_descriptor vertex_t;
34 typedef typename graph_traits<Graph>::vertices_size_type v_size_t;
35 typedef iterator_property_map< typename std::vector<v_size_t>::iterator,
37 > vertex_to_v_size_map_t;
39 std::vector<v_size_t> component_vector(num_vertices(g));
40 vertex_to_v_size_map_t component(component_vector.begin(), vm);
41 std::vector<vertex_t> vertices_by_component(num_vertices(g));
43 v_size_t num_components = connected_components(g, component);
45 if (num_components < 2)
48 vertex_iterator_t vi, vi_end;
49 boost::tie(vi,vi_end) = vertices(g);
50 std::copy(vi, vi_end, vertices_by_component.begin());
52 bucket_sort(vertices_by_component.begin(),
53 vertices_by_component.end(),
58 typedef typename std::vector<vertex_t>::iterator vec_of_vertices_itr_t;
60 vec_of_vertices_itr_t ci_end = vertices_by_component.end();
61 vec_of_vertices_itr_t ci_prev = vertices_by_component.begin();
62 if (ci_prev == ci_end)
65 for(vec_of_vertices_itr_t ci = boost::next(ci_prev);
66 ci != ci_end; ci_prev = ci, ++ci
69 if (component[*ci_prev] != component[*ci])
70 vis.visit_vertex_pair(*ci_prev, *ci, g);
78 template <typename Graph, typename VertexIndexMap>
79 inline void make_connected(Graph& g, VertexIndexMap vm)
81 default_add_edge_visitor vis;
82 make_connected(g, vm, vis);
88 template <typename Graph>
89 inline void make_connected(Graph& g)
91 make_connected(g, get(vertex_index,g));
99 #endif //__MAKE_CONNECTED_HPP__