Honeycomb Networks Topological Properties and Communication Algorithms
Honeycomb Networks Topological Properties and Communication Algorithms
Abstract—The honeycomb mesh, based on hexagonal plane tessellation, is considered as a multiprocessor interconnection
network. A honeycomb mesh network with n nodes has degree 3 and diameter ª 1.63 n - 1, which is 25 percent smaller degree
and 18.5 percent smaller diameter than the mesh-connected computer with approximately the same number of nodes. Vertex and
edge symmetric honeycomb torus network is obtained by adding wraparound edges to the honeycomb mesh. The network cost,
defined as the product of degree and diameter, is better for honeycomb networks than for the two other families based on square
(mesh-connected computers and tori) and triangular (hexagonal meshes and tori) tessellations. A convenient addressing scheme for
nodes is introduced which provides simple computation of shortest paths and the diameter. Simple and optimal (in the number of
required communication steps) routing, broadcasting, and semigroup computation algorithms are developed. The average distance
in honeycomb torus with n nodes is proved to be approximately 0.54 n . In addition to honeycomb meshes bounded by a regular
hexagon, we consider also honeycomb networks with rhombus and rectangle as the bounding polygons.
Index Terms—Computer architecture, interconnection topology, routing, broadcasting, mesh-connected computer, hexagonal
tessellation.
—————————— ✦ ——————————
1 INTRODUCTION
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio
STOJMENOVIC: HONEYCOMB NETWORKS: TOPOLOGICAL PROPERTIES AND COMMUNICATION ALGORITHMS 1037
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio
1038 IEEE TRANSACTIONS ON PARALLEL AND DISTRIBUTED SYSTEMS, VOL. 8, NO. 10, OCTOBER 1997
TABLE 1
COMPARISON OF NETWORKS
network degree diameter cost bisection width
mesh-connected computer 4 2 n 8 n n
hexagonal mesh 6 1.16 n 6.93 n 2.31 n
honeycomb mesh 3 1.63 n 4.90 n 0.82 n
torus 4 n 4 n 2 n
hexagonal torus 6 0.58 n 3.46 n 4.61 n
honeycomb torus 3 0.81 n 2.45 n 2.04 n
honeycomb rhombic mesh 3 2.83 n 8.49 n 0.71 n
honeycomb square mesh 3 2 n 6 n 0.5 n
honeycomb rhombic torus 3 1.06 n 3.18 n 1.41 n
honeycomb square torus 3 n 3 n n
2
hypercube log n log n log n n/2
cube-connected cycles 3 O(log n) O(log n) O(n/log n)
butterfly 4 O(log n) O(log n) O(n/log n)
deBruijn 4 O(log n) O(log n) O(n/log n)
Fig. 4 indicates black and white nodes in HM3. Edges going + w = 2 for each white one. Two nodes are connected
from a black node to a white node have vector coordinates by an edge iff they differ in exactly one coordinate.
+ + +
x = (1, 0, 0), y = (0, 1, 0), or z = (0, 0, 1), while edges going
out of white vertices toward black ones have vector coordi-
- - -
nates x = (-1, 0, 0), y = (0, -1, 0), or z = (0, 0, -1). Note
+ - + - + -
that x and x , y and y , z and z are parallel to x-, y-, and
z-axes, respectively. Once edge vectors are defined, con-
sider what are coordinates of nodes. The z-coordinate of all
nodes reachable from a given node A by traveling along
edges (±1, 0, 0) or (0, ±1, 0) (i.e., parallel to x- or y-axis) is
the same as the z-coordinate of A. All nodes with a fixed
(say) z-coordinate belong to a zigzag chain. The z-
coordinate will change only by following an edge parallel
+ -
to z-axis, i.e., edge z or z . If two central chains are given z-
coordinates 1 and 0, the six zigzag chains with respect to z-
axis are z = -2, z = -1, z = 0, z = 1, z = 2, and z = 3; they are
drawn in bold in Fig. 4. Similar zigzag chains can be de-
fined for x- and y-coordinates. For instance, the coordinates
of node A in Fig. 4 are (1, 2, -1). If each zigzag chain with a
fixed z-coordinate is drawn as a vertical line, the honey-
comb mesh receives the shape of bricks in a wall network Fig. 4. Z-coordinates of nodes.
(Fig. 5). The network drawn in this way clearly shows that
the number of edges parallel to z-axis in a shortest path
between two nodes is greater or equal to the difference in
their z coordinate. It easily follows that the coordinates of
all nodes are integers (x, y, z) such that -t + 1 ≤ x, y, z ≤ t.
Moreover, we can prove the following lemma.
LEMMA 1. Nodes of HMt can be coded by integer triples (u, v, w)
such that - t + 1 ≤ u, v, w ≤ t, and 1 ≤ u + v + w ≤ 2. Two
nodes (u¢, v¢, w¢) and (u¢¢, v¢¢, w¢¢) are connected by an
edge iff |u¢ - u¢¢| + |v¢ - v¢¢| + |w¢ - w¢¢| = 1.
PROOF. Edges going from a black node to a white node join
(u, v, w) to either (u + 1, v, w), (u, v + 1, w), or (u, v,
w + 1) and, thus, the sum u + v + w increases by one.
Analogously, edges joining from a white node to a black
one join (u, v, w) to one of (u - 1, v, w), (u, v - 1, w), or
(u, v, w - 1), thus decreasing the sum u + v + w by
one. Since u + v + w alternatively increases (decreases)
by one for each white (black) node, it obviously fol-
lows that u + v + w = 1 for each black node and u + v Fig. 5. Brick drawing of HM3.
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio
STOJMENOVIC: HONEYCOMB NETWORKS: TOPOLOGICAL PROPERTIES AND COMMUNICATION ALGORITHMS 1039
Using mathematical induction, we can prove that the edges parallel to z-axis. In short, the routing algorithm
distance between nodes (u¢, v¢, w¢) and (u¢¢, v¢¢, w¢¢) of a checks at each current node which of the edge directions x,
y, or z (in this order) would reduce the distance to the des-
hexagonal mesh is |u¢ - u¢¢| + |v¢ - v¢¢| + |w¢ - w¢¢|, and
tination, and will send the message on that edge. At least
that the diameter of HMt is 4t - 1 (nodes (t, -t + 1, 0) and one of the edge directions would lead to a node closer to
(-t + 1, t, 1) are at distance |2t - 1| + |2t - 1| + 1 = 4t - 1). the destination. The algorithm is simple, and details and
Thus, honeycomb mesh network with n nodes has degree code are omitted. The algorithm is optimal in the number of
three and diameter 4 n6 - 1 ª 1.63 n - 1 and its network communication steps, and asymptotically optimal in the
n
number of computation steps (i.e., has constant number of
cost is 12 6
- 3 = 2 6 n - 3 ª 4.90 n . The number of instructions between any two communication steps). The
2
nodes of a hexagonal mesh of size t is n = 3t - 3t + 1 [3], shortest path between two nodes can be explicitly described
(with constant number of parameters) on the basis of their
while its diameter is 2t - 2. The degree of nodes is six.
2 3
addresses and consists of two zigzag chains.
Therefore, n ª 3t and t ª 3
n . The network cost is there- The broadcasting process should satisfy some desirable
fore ª 6 2e 3
j
n = 4 3 n ª 6.93 n . The diameter of a properties. A node should not receive the message simulta-
3
neously from more than one of its neighbors. Moreover, it
square mesh-connected computer with n nodes is 2 n - 2 is desirable that a node receives the message exactly once
and its degree is four. Thus, the cost of a mesh-connected for the whole duration of the broadcasting process. This is
computer is ª 8 n . achieved by a broadcasting algorithm that follows the same
Therefore, honeycomb meshes have 25 percent smaller paths generated by the corresponding routing scheme. The
degree and 18.5 percent smaller diameter than the mesh- source will send the message to all three neighbors. Each
connected computer with the same number of nodes. This other node will examine the incoming message and for-
is an important advantage of honeycomb meshes over ward it to zero, one, or two of its neighbors (it will never be
square ones. The advantage can be explained as follows. It sent back to the neighbor the message was received from).
is clear from Fig. 5, which also shows that honeycomb mesh The algorithm counts the number of times the message
can be embedded onto the square grid, that the distance used edge parallel to the x-, y-, and z-axes, respectively. The
between two nodes remains virtually the same if a quarter counters count edges in positive directions by positive
of edges is eliminated from the grid. Further, by making a numbers (i.e., by increasing their value by one) and edges
hexagonal rather than square boundary, the boundary is in negative directions (opposite to axis orientations) by
closer to a circular shape, thus reducing its diameter. negative numbers (i.e., by decreasing their value), and do
Let us compare the network costs of the networks based not allow forwarding messages in opposite directions. The
on regular plane tessellations. The hexagonal mesh has ap- number of communication steps between a source and any
proximately 13 percent smaller cost than a mesh-connected node is equal to the distance between them. The details of
computer with the same number of nodes. Honeycomb the algorithm are omitted. In the one-port communication
mesh, however, has smaller cost than either of these two model (in which nodes can forward the message to only
networks that are used in practice for designing multiproc- one neighbor at a time), the number of communication
essors. The cost of a honeycomb network is approximately steps necessary for the message to reach each particular
29 percent smaller than the cost of a hexagonal mesh and destination is equal to the length of the shortest path be-
about 39 percent smaller than the cost of a mesh-connected tween the source and a given node plus the delay that cor-
computer. The latter is a significant reduction in the cost responds to messages being sent on other edges (leading to
compared to (arguably) the most popular network on the other destinations) of intermediate nodes before it being
market. sent on the edge of the shortest path. We can show that the
delay is, at most, three communication steps. If all port
communication model (in which each node is allowed to
3 ROUTING, BROADCASTING, AND SEMIGROUP
send the same message to all or some of its neighbors at
COMPUTATION IN HONEYCOMB MESHES once) is used, then the number of communication steps
The problem of finding a path from a source to destination between the source and any node is equal to the distance
and forwarding a message along the path is known as the between them.
routing problem. The broadcasting task is to send a mes- The semigroup computation problem is to compute x1 ©
sage from a source to all other nodes. Routing and broad- x2 © ... © xn, where © is any binary associative and com-
casting are the basic data communication problems for mutative operation (e.g., the sum, product, minimum,
every network. In this section, we will describe routing and maximum, or, and, etc.) and the data are distributed one
broadcasting algorithms for the honeycomb mesh net- per processor on a given network. The problem can be
works. solved by running the broadcasting algorithm from any
Suppose that source = (u¢, v¢, w¢) is the source node node in the reverse order, such that data are being collected
while the destination node is node dest = (u¢¢, v¢¢, w¢¢). Let (each node on the path applies “©” to reduce the informa-
Dx = u¢¢ - u¢, Dy = v¢¢ - v¢, and Dz = w¢¢ - w¢. The vector tion forwarded to the next node to a single message).
(Dx, Dy, Dz) is the “translation” vector for the message. The Broadcasting can be then applied to send the result to other
shortest path between the two nodes consists of |Dx| edges nodes. The performance is the same as that of the broad-
parallel to x-axis, |Dy| edges parallel to y-axis, and |Dz| casting algorithm.
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio
1040 IEEE TRANSACTIONS ON PARALLEL AND DISTRIBUTED SYSTEMS, VOL. 8, NO. 10, OCTOBER 1997
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio
STOJMENOVIC: HONEYCOMB NETWORKS: TOPOLOGICAL PROPERTIES AND COMMUNICATION ALGORITHMS 1041
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio
1042 IEEE TRANSACTIONS ON PARALLEL AND DISTRIBUTED SYSTEMS, VOL. 8, NO. 10, OCTOBER 1997
REFERENCES
[1] S.G. Akl, Parallel Computation: Models and Methods. Prentice Hall,
1996.
[2] R. Ben-Natan and A. Barak, “Parallel Contractions of Grids for Task
Assignment to Processor Networks,” Networks, vol. 22, pp. 539-562,
1992.
[3] M.S. Chen, K.G. Shin, and D.D. Kandlur, “Addressing, Routing,
and Broadcasting in Hexagonal Mesh Multiprocessors,” IEEE
Trans. Computers, vol. 39, no. 1, pp. 10-18, Jan. 1990.
[4] F.T. Leighton, Introduction to Parallel Algorithms and Architectures:
Arrays, Trees, and Hypercubes. Morgan Kaufmann, 1992.
[5] D. Milutinovic, V. Milutinovic, and B. Soucek, “The Honeycomb
Architecture,” Computer, vol. 20, no. 4, pp. 81-83, Apr. 1987.
[6] V. Milutinovic, “Mapping of Neural Networks on the Honey-
comb Architecture,” Proc. IEEE, vol. 77, no. 12, pp. 1,875-1,878,
1989.
[7] G. Rote, “On the Connection Between Hexagonal and Unidirec-
tional Rectangular Systolic Arrays,” Lecture Notes in Computer Sci-
ence, vol. 227, pp. 70-83. Springer-Verlag, 1986.
[8] B. Robic and J. Silc, “High-Performance Computing on a Honey-
comb Architecture,” Proc. Second Int’l ACPC Parallel Computation
Conf., pp. 22-28. Austria: Springer-Verlag, 1993.
[9] K.G. Shin, “HARTS: A Distributed Real-Time Architecture,” Com-
Fig. 9. Honeycomb rhombic torus. puter, vol. 24, no. 5, pp. 25-35, May 1991.
[10] I. Stojmenovic, “Direct Interconnection Networks,” Parallel and
Distributed Computing Handbook, A.Y. Zomaya, ed., pp. 537-567.
McGraw-Hill, 1996.
[11] I. Stojmenovic, “Honeycomb Networks,” Proc. Math. Foundations
of Computer Science MFCS ’95, Lecture Notes in Computer Science,
vol. 969, pp. 267-276. Prague: Springer-Verlag, 1995.
[12] H.Y. Youn and J.Y. Lee, “An Efficient Dictionary Machine Using
Hexagonal Processor Arrays,” IEEE Trans. Parallel and Distributed
Systems, vol. 7, no. 3, pp. 266-273, Mar. 1996.
ACKNOWLEDGMENTS
The author is grateful to the referees whose careful reading
has greatly improved the clarity and conciseness of the
presentation, and to Professor Michael A. Palis for his
timely and professional handling of my submission. This
research is partially supported by NSERC.
d licensed use limited to: AMRITA VISHWA VIDYAPEETHAM AMRITA SCHOOL OF ENGINEERING. Downloaded on September 12,2022 at 07:50:44 UTC from IEEE Xplore. Restrictio