1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586 |
- //=======================================================================
- // Copyright 2002 Indiana University.
- // Authors: Andrew Lumsdaine, Lie-Quan Lee, Jeremy G. Siek
- //
- // Distributed under the Boost Software License, Version 1.0. (See
- // accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- //=======================================================================
- #ifndef BOOST_CREATE_CONDENSATION_GRAPH_HPP
- #define BOOST_CREATE_CONDENSATION_GRAPH_HPP
- #include <boost/graph/graph_traits.hpp>
- #include <boost/property_map/property_map.hpp>
- namespace boost
- {
- template < typename Graph, typename ComponentLists, typename ComponentNumberMap,
- typename CondensationGraph, typename EdgeMultiplicityMap >
- void create_condensation_graph(const Graph& g, const ComponentLists& components,
- ComponentNumberMap component_number, CondensationGraph& cg,
- EdgeMultiplicityMap edge_mult_map)
- {
- typedef typename graph_traits< Graph >::vertex_descriptor vertex;
- typedef typename graph_traits< Graph >::vertices_size_type size_type;
- typedef
- typename graph_traits< CondensationGraph >::vertex_descriptor cg_vertex;
- std::vector< cg_vertex > to_cg_vertex(components.size());
- for (size_type s = 0; s < components.size(); ++s)
- to_cg_vertex[s] = add_vertex(cg);
- for (size_type si = 0; si < components.size(); ++si)
- {
- cg_vertex s = to_cg_vertex[si];
- std::vector< cg_vertex > adj;
- for (size_type i = 0; i < components[si].size(); ++i)
- {
- vertex u = components[s][i];
- typename graph_traits< Graph >::adjacency_iterator v, v_end;
- for (boost::tie(v, v_end) = adjacent_vertices(u, g); v != v_end;
- ++v)
- {
- cg_vertex t = to_cg_vertex[component_number[*v]];
- if (s != t) // Avoid loops in the condensation graph
- adj.push_back(t);
- }
- }
- std::sort(adj.begin(), adj.end());
- if (!adj.empty())
- {
- size_type i = 0;
- cg_vertex t = adj[i];
- typename graph_traits< CondensationGraph >::edge_descriptor e;
- bool inserted;
- boost::tie(e, inserted) = add_edge(s, t, cg);
- put(edge_mult_map, e, 1);
- ++i;
- while (i < adj.size())
- {
- if (adj[i] == t)
- put(edge_mult_map, e, get(edge_mult_map, e) + 1);
- else
- {
- t = adj[i];
- boost::tie(e, inserted) = add_edge(s, t, cg);
- put(edge_mult_map, e, 1);
- }
- ++i;
- }
- }
- }
- }
- template < typename Graph, typename ComponentLists, typename ComponentNumberMap,
- typename CondensationGraph >
- void create_condensation_graph(const Graph& g, const ComponentLists& components,
- ComponentNumberMap component_number, CondensationGraph& cg)
- {
- create_condensation_graph(
- g, components, component_number, cg, dummy_property_map());
- }
- } // namespace boost
- #endif // BOOST_CREATE_CONDENSATION_GRAPH_HPP
|