FOR_VERTEX((v1), loop) \
FOR_VEDGE((v1),(e),(v2))
-#define FOR_RIM_VERTEX(vy,vx,v, loop) \
- for ((vy)=0; (vy)<Y; (vy)+=Y-1) \
+#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];