X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?p=moebius2.git;a=blobdiff_plain;f=mgraph.h;h=63fd5dca82ae60e4a78d865bcbb3cdad71a04527;hp=85f1a3c8d0ef621e8417a2715c1af014df187e91;hb=HEAD;hpb=c01a4fe080a448234b39afb7ea7f7a3c9c229fb7

diff --git a/mgraph.h b/mgraph.h
index 85f1a3c..63fd5dc 100644
--- a/mgraph.h
+++ b/mgraph.h
@@ -89,14 +89,22 @@
 #define V6 6
 #define V3 3
 
-#define FOR_VERTEX(v) \
-  for ((v)=0; (v)<N; (v)++)
 
-#define FOR_VPEDGE(v,e) \
+/* Loop constructors are macros of the form
+ *    LOOP(v,zero,n, precomp)
+ * which work much like this one:
+ */
+#define INNER(v,zero,n, precomp)  \
+  for ((v)=(zero); precomp, (v)<(n); (v)++)
+
+#define FOR_VERTEX(v,loop) \
+  loop ((v), 0, N, NOTHING)
+
+#define FOR_VPEDGE(e) \
   for ((e)=0; (e)<V6; (e)++)
 
-int edge_end2(unsigned v1, int e);
-#define EDGE_END2 edge_end2
+extern short edge_end2_memo[][V6];
+#define EDGE_END2(v1,e) edge_end2_memo[v1][e]
 
 /* given    v1,e     s.t.  v2==EDGE_END2(v1,e) >= 0,
  * returns  eprime   s.t.  v1==EDGE_END2(v2,eprime) */
@@ -107,7 +115,7 @@ int edge_reverse(int v1, int e);
 #define RIM_VERTEX_P(v) (((v) & ~XMASK) == 0 || ((v) & ~XMASK) == (Y-1)*Y1)
 
 #define FOR_VEDGE_X(v1,e,v2,init,otherwise)	\
-  FOR_VPEDGE((v1),(e))				\
+  FOR_VPEDGE((e))				\
     if (((v2)= EDGE_END2((v1),(e)),		\
 	 (init),				\
 	 (v2)) < 0) { otherwise; } else
@@ -117,15 +125,22 @@ int edge_reverse(int v1, int e);
 #define FOR_VEDGE(v1,e,v2)			\
   FOR_VEDGE_X(v1,e,v2,NOTHING,NOTHING)
 
-#define FOR_EDGE(v1,e,v2)			\
-  FOR_VERTEX((v1))				\
+#define FOR_EDGE(v1,e,v2, loop)			\
+  FOR_VERTEX((v1), loop)			\
     FOR_VEDGE((v1),(e),(v2))
 
-#define FOR_RIM_VERTEX(vy,vx,v)					\
-  for ((vy)=0; (vy)<Y; (vy)+=Y-1)				\
-    for ((vx)=0; (v)= (vy)<<YSHIFT | (vx), (vx)<X; (vx)++)
+#define FOR_NEAR_RIM_VERTEX(vy,vx,v, disttorim, loop)		\
+  for ((vy)=(disttorim); (vy)<Y; (vy)+=Y-1-2*(disttorim))	\
+    loop ((vx), 0, X, (v)= (vy)<<YSHIFT | (vx))
+
+#define FOR_RIM_VERTEX(vy,vx,v, loop)		\
+  FOR_NEAR_RIM_VERTEX((vy),(vx),(v), 0, loop)
+
+int vertices_span_join_p(int v0, int v1);
 
 typedef double Vertices[N][D3];
 struct Vertices { Vertices a; };
 
+void mgraph_prepare(void);
+
 #endif /*MGRAPH_H*/