Delone set

The red points form part of an ε-net for the Euclidean plane, where ε is the radius of the large yellow disks. The blue disks of half the radius are disjoint, and the yellow disks together cover the whole plane, satisfying the two definitional requirements on an ε-net.

In the mathematical theory of metric spaces, ε-nets, ε-packings, ε-coverings, uniformly discrete sets, relatively dense sets, and Delone sets (named after Boris Delone) are several closely related definitions of well-spaced sets of points, and the packing radius and covering radius of these sets measure how well-spaced they are. These sets have applications in coding theory, approximation algorithms, and the theory of quasicrystals.

Definitions

If (M,d) is a metric space, and X is a subset of M, then the packing radius of X is half of the infimum of distances between distinct members of X. If the packing radius is r, then open balls of radius r centered at the points of X will all be disjoint from each other, and each open ball centered at one of the points of X with radius 2r will be disjoint from the rest of X. The covering radius of X is the infimum of the numbers r such that every point of M is within distance r of at least one point in X; that is, it is the smallest radius such that closed balls of that radius centered at the points of X have all of M as their union. An ε-packing is a set X of packing radius ≥ ε/2 (equivalently, minimum distance ≥ ε), an ε-covering is a set X of covering radius ≤ ε, and an ε-net is a set that is both an ε-packing and an ε-covering. A set is uniformly discrete if it has a nonzero packing radius, and relatively dense if it has a finite covering radius. A Delone set is a set that is both uniformly discrete and relatively dense; thus, every ε-net is Delone, but not vice versa.[1][2]

Construction of ε-nets

As the most restrictive of the definitions above, ε-nets are at least as difficult to construct as ε-packings, ε-coverings, and Delone sets. However, whenever the points of M have a well-ordering, transfinite induction shows that it is possible to construct an ε-net N, by including in N every point for which the infimum of distances to the set of earlier points in the ordering is at least ε. For finite sets of points in a Euclidean space of bounded dimension, each point may be tested in constant time by mapping it to a grid of cells of diameter ε, and using a hash table to test which nearby cells already contain points of N; thus, in this case, an ε-net can be constructed in linear time.[3][4]

For more general finite or compact metric spaces, an alternative algorithm of Teo Gonzalez can be used to construct a finite ε-net. This algorithm initializes the net N to be empty, and then repeatedly adds to N the farthest point in M from N, breaking ties arbitrarily and stopping when all points of M are within distance ε of N.[5] In spaces of bounded doubling dimension, Gonzalez' algorithm can be implemented in O(n log n) time for point sets with a polynomial ratio between their farthest and closest distances, and approximated in the same time bound for arbitrary point sets.[6]

Applications

Coding theory

In the theory of error-correcting codes, the metric space containing a block code C consists of strings of a fixed length, say n, taken over an alphabet of size q (can be thought of as vectors), with the Hamming metric. This space is denoted by \scriptstyle \mathcal{A}_q^n. The covering radius and packing radius of this metric space are related to the code's ability to correct errors.

Approximation algorithms

Har-Peled & Raichel (2013) describe an algorithmic paradigm that they call "net and prune" for designing approximation algorithms for certain types of geometric optimization problems defined on sets of points in Euclidean spaces. An algorithm of this type works by performing the following steps:

  1. Choose a random point p from the point set, find its nearest neighbor q, and set ε to the distance between p and q.
  2. Test whether ε is (approximately) bigger than or smaller than the optimal solution value (using a technique specific to the particular optimization problem being solved)
    • If it is bigger, remove from the input the points whose closest neighbor is farther than ε
    • If it is smaller, construct an ε-net N, and remove from the input the points that are not in N.

In both cases, the expected number of remaining points decreases by a constant factor, so the time is dominated by the testing step. As they show, this paradigm can be used to construct fast approximation algorithms for k-center clustering, finding a pair of points with median distance, and several related problems.

A hierarchical system of nets, called a net-tree, may be used in spaces of bounded doubling dimension to construct well-separated pair decompositions, geometric spanners, and approximate nearest neighbors.[6][7]

Crystallography

For points in Euclidean space, a set X is a Meyer set if it is relatively dense and its difference set X − X is uniformly discrete. Equivalently, X is a Meyer set if both X and X − X are Delone. Meyer sets are named after Yves Meyer, who introduced them (with a different but equivalent definition based on harmonic analysis) as a mathematical model for quasicrystals. They include the point sets of lattices, Penrose tilings, and the Minkowski sums of these sets with finite sets.[8]

References

  1. ↑ Clarkson, Kenneth L. (2006), "Building triangulations using ε-nets", STOC'06: Proceedings of the 38th Annual ACM Symposium on Theory of Computing, New York: ACM, pp. 326–335, doi:10.1145/1132516.1132564, MR 2277158.
  2. ↑ Some sources use "ε-net" for what is here called an "ε-covering"; see, e.g. Sutherland, W.A. (1975), Introduction to metric and topological spaces, Oxford University Press, p. 110, ISBN 0-19-853161-3, Zbl 0304.54002.
  3. ↑ Har-Peled, S. (2004), "Clustering motion", Discrete and Computational Geometry 31 (4): 545–565, doi:10.1007/s00454-004-2822-7, MR 2053498.
  4. ↑ Har-Peled, S.; Raichel, B. (2013), "Net and prune: A linear time algorithm for Euclidean distance problems", STOC'13: Proceedings of the 45th Annual ACM Symposium on Theory of Computing, pp. 605–614.
  5. ↑ Gonzalez, T. F. (1985), "Clustering to minimize the maximum intercluster distance", Theoretical Computer Science 38 (2-3): 293–306, doi:10.1016/0304-3975(85)90224-5, MR 807927.
  6. ↑ 6.0 6.1 Har-Peled, S.; Mendel, M. (2006), "Fast construction of nets in low-dimensional metrics, and their applications", SIAM Journal on Computing 35 (5): 1148–1184, doi:10.1137/S0097539704446281, MR 2217141.
  7. ↑ Krauthgamer, Robert; Lee, James R. (2004), "Navigating nets: simple algorithms for proximity search", Proceedings of the 15th Annual ACM-SIAM Symposium on Discrete Algorithms (SODA '04), Philadelphia, PA, USA: Society for Industrial and Applied Mathematics, pp. 798–807, ISBN 0-89871-558-X.
  8. ↑ Moody, Robert V. (1997), "Meyer sets and their duals", The Mathematics of Long-Range Aperiodic Order (Waterloo, ON, 1995), NATO Advanced Science Institutes Series C: Mathematical and Physical Sciences 489, Dordrecht: Kluwer Academic Publishers, pp. 403–441, MR 1460032.

External links