/*
* edge descriptor f = 00 | e | y | x
- * 2 YBITS XBITS
+ * 3 YBITS XBITS
*
- * e is 0,1 or 2. The edge is edge e out of vertex (x,y).
+ * e is 0..5. The edge is edge e out of vertex (x,y).
+ *
+ * BGL expects an undirected graph's edges to have two descriptors
+ * each, one in each direction.
*/
/*
#define VMASK (YMASK|XMASK)
#define ESHIFT (YBITS|XBITS)
+/* When we enumerate the edges we mainly just increment f.
+ * To avoid nonexistent edges
+ * we start with e=0 | y=1 | x=0
+ * and explicitly skip e=1 | y=0 | ...
+ * We go only up to e=2 | YMAX | XMAX
+ * ie we stop just before e=3 | y=0 | x=0
+ * so that we get each edge once rather than twice.
+ */
+
+#define F_ALL_MIN (1 << YSHIFT)
+#define F_ALL_MAX (3 << ESHIFT)
+#define F_ALL_SKIP_AT (1 << ESHIFT)
+#define F_ALL_SKIP_SET (1 << YSHIFT)
+
namespace boost {
+ // We make Layout a model of various BGL Graph concepts.
+ // This mainly means that graph_traits<Layout> has lots of stuff.
+
+ // First, some definitions used later:
+
struct layout_graph_traversal_category :
public virtual incidence_graph_tag,
- public virtual edge_list_graph_tag
- public virtual vertex_list_graph_tag { };
+ public virtual vertex_list_graph_tag,
+ public virtual edge_list_graph_tag { };
+ struct OutEdgeIncrable {
+ int f;
+ OutEdgeIncrable& operator++() { f += 1<<ESHIFT; return self; }
+ OutEdgeIncrable(int v, int e) : f(v | (e << ESHIFT)) { }
+ };
+ struct AnyEdgeIncrable {
+ int f;
+ AnyEdgeIncrable& operator++() {
+ f++;
+ if (f == F_ALL_SKIP_AT) f |= F_ALL_SKIP_SET
+ return self;
+ }
+ AnyEdgeIncrable(int _f) : f(_f) { }
+ AnyEdgeIncrable() : f(F_ALL_MIN) { }
+ };
+
struct graph_traits<Layout> {
- /* Concept Graph: */
+
+ // Concept Graph:
+
typedef int vertex_descriptor; /* vertex number, -1 => none */
typedef int edge_descriptor; /* see above */
typedef undirected_tag directed_category;
typedef disallow_parallel_ege edge_parallel_category;
typedef layout_graph_traversal_category traversal_category;
inline int null_vertex() { return -1; }
- }
- struct out_edge_iterator_policies {
- static void increment(int& f) { f += 1<<ESHIFT; }
- static void decrement(int& f) { f -= 1<<ESHIFT; }
- template <class Reference>
- static Reference dereference(type<Reference>, const int& f)
- { return const_cast<Reference>(f); }
-
- static bool equal(int x, int y) { return x == y; }
- }
- struct graph_traits<Layout> {
- /* Concept IncidenceGraph: */
- typedef iterator_adaptor<int, out_edge_iterator_policies,
- iterator<std::bidirectional_iterator_tag,int>
- > out_edge_iterator;
+ // Concept IncidenceGraph:
+
+ typedef counting_iterator<OutEdgeIncrable,
+ forward_iterator_tag> out_edge_iterator;
+ typedef int degree_size_type;
- inline int source(int f, const Layout& g) { return f&VMASK; }
- inline int target(int f, const Layout& g) { return EDGE_END2(f&VMASK, f>>ESHIFT); }
- inline std::pair<typename graph_traits<Layout>::out_edge_iterator,
- typename graph_traits<Layout>::out_edge_iterator>
- out_edges(int v, const Layout& g) {
- return std::make_pair(VE_MIN(v), VE_MAX(v));
+ inline int source(int f, const Layout&) { return f&VMASK; }
+ inline int target(int f, const Layout&) { return EDGE_END2(f&VMASK, f>>ESHIFT); }
+ inline std::pair<out_edge_iterator,out_edge_iterator>
+ out_edges(int v, const Layout&) {
+ return std::make_pair(out_edge_iterator(OutEdgeIncrable(v, VE_MIN(v))),
+ out_edge_iterator(OutEdgeIncrable(v, VE_MAX(v))));
}
-
- out_edge_iterator> {
-
+ out_degree(int v, const Layout&) { return VE_MAX(v) - VE_MIN(v); }
- /* Concept VertexListGraph: */
+ // Concept VertexListGraph:
typedef counting_iterator<int> vertex_iterator;
+ typedef unsigned vertices_size_type;
+ inline std::pair<vertex_iterator,vertex_iterator>
+ vertices(const Layout&) {
+ return std::make_pair(vertex_iterator(0), vertex_iterator(N));
+ }
+ inline unsigned num_vertices(const Layout&) { return N; }
+
+ // Concept EdgeListGraph:
+ typedef counting_iterator<AnyEdgeIncrable,
+ forward_iterator_tag> edge_iterator;
+ typedef unsigned edges_size_type;
+ inline std::pair<edge_iterator,edge_iterator>
+ edges(const Layout&) {
+ return std::make_pair(edge_iterator(AnyEdgeIncrable()),
+ edge_iterator(AnyEdgeIncrable(F_ALL_MAX)));
+ }
+ unsigned num_edges(const Layout&) {
+ return F_ALL_MAX - F_ALL_MIN - ((F_SKIP_SET|F_SKIP_AT) - F_SKIP_AT);
+ }
}