diff -r 000000000000 -r e4d67989cc36 stdcpp/tsrc/Boost_test/graph/src/king_ordering.cpp --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/stdcpp/tsrc/Boost_test/graph/src/king_ordering.cpp Tue Feb 02 02:01:42 2010 +0200 @@ -0,0 +1,143 @@ +//======================================================================= +// Copyright 1997, 1998, 1999, 2000 University of Notre Dame. +// Authors: Andrew Lumsdaine, Lie-Quan Lee, Jeremy G. Siek +// Doug Gregor, D. Kevin McGrath +// +// 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) +//======================================================================= +/* + * © Portions copyright (c) 2006-2007 Nokia Corporation. All rights reserved. +*/ + +#include +#include +#include +#include +#include +#include +#include +#ifdef __SYMBIAN32__ +#include "std_log_result.h" +#define LOG_FILENAME_LINE __FILE__, __LINE__ +#endif +/* + Sample Output + original bandwidth: 8 + Reverse Cuthill-McKee ordering starting at: 6 + 8 3 0 9 2 5 1 4 7 6 + bandwidth: 4 + Reverse Cuthill-McKee ordering starting at: 0 + 9 1 4 6 7 2 8 5 3 0 + bandwidth: 4 + Reverse Cuthill-McKee ordering: + 0 8 5 7 3 6 4 2 1 9 + bandwidth: 4 + */ +int main(int , char* []) +{ + using namespace boost; + using namespace std; + typedef adjacency_list > > Graph; + typedef graph_traits::vertex_descriptor Vertex; + typedef graph_traits::vertices_size_type size_type; + + typedef std::pair Pair; + Pair edges[14] = { Pair(0,3), //a-d + Pair(0,5), //a-f + Pair(1,2), //b-c + Pair(1,4), //b-e + Pair(1,6), //b-g + Pair(1,9), //b-j + Pair(2,3), //c-d + Pair(2,4), //c-e + Pair(3,5), //d-f + Pair(3,8), //d-i + Pair(4,6), //e-g + Pair(5,6), //f-g + Pair(5,7), //f-h + Pair(6,7) }; //g-h + + Graph G(10); + for (int i = 0; i < 14; ++i) + add_edge(edges[i].first, edges[i].second, G); + + graph_traits::vertex_iterator ui, ui_end; + + property_map::type deg = get(vertex_degree, G); + for (boost::tie(ui, ui_end) = vertices(G); ui != ui_end; ++ui) + deg[*ui] = degree(*ui, G); + + property_map::type + index_map = get(vertex_index, G); + + std::cout << "original bandwidth: " << bandwidth(G) << std::endl; + + std::vector inv_perm(num_vertices(G)); + std::vector perm(num_vertices(G)); + { + Vertex s = vertex(6, G); + //king_ordering + king_ordering(G, s, inv_perm.rbegin(), get(vertex_color, G), + get(vertex_degree, G), get(vertex_index, G)); + cout << "King ordering starting at: " << s << endl; + cout << " "; + for (std::vector::const_iterator i = inv_perm.begin(); + i != inv_perm.end(); ++i) + cout << index_map[*i] << " "; + cout << endl; + + for (size_type c = 0; c != inv_perm.size(); ++c) + perm[index_map[inv_perm[c]]] = c; + std::cout << " bandwidth: " + << bandwidth(G, make_iterator_property_map(&perm[0], index_map, perm[0])) + << std::endl; + } + { + Vertex s = vertex(0, G); + //king_ordering + king_ordering(G, s, inv_perm.rbegin(), get(vertex_color, G), + get(vertex_degree, G), get(vertex_index, G)); + cout << "King ordering starting at: " << s << endl; + cout << " "; + for (std::vector::const_iterator i=inv_perm.begin(); + i != inv_perm.end(); ++i) + cout << index_map[*i] << " "; + cout << endl; + + for (size_type c = 0; c != inv_perm.size(); ++c) + perm[index_map[inv_perm[c]]] = c; + std::cout << " bandwidth: " + << bandwidth(G, make_iterator_property_map(&perm[0], index_map, perm[0])) + << std::endl; + } + + { + //king_ordering + king_ordering(G, inv_perm.rbegin()); + + cout << "King ordering:" << endl; + cout << " "; + for (std::vector::const_iterator i=inv_perm.begin(); + i != inv_perm.end(); ++i) + cout << index_map[*i] << " "; + cout << endl; + + for (size_type c = 0; c != inv_perm.size(); ++c) + perm[index_map[inv_perm[c]]] = c; + std::cout << " bandwidth: " + << bandwidth(G, make_iterator_property_map(&perm[0], index_map, perm[0])) + << std::endl; + } + +#ifdef __SYMBIAN32__ +std_log(LOG_FILENAME_LINE,"[End Test Case ]"); + + testResultXml("king_ordering"); + close_log_file(); +#endif + return 0; +}