Consistent High Dimensional Rounding with Side Information

Abstract: In standard rounding, we want to map each value $X$ in a large continuous space (e.g., $R$) to a nearby point $P$ from a discrete subset (e.g., $Z$). This process seems to be inherently discontinuous in the sense that two consecutive noisy measurements $X_1$ and $X_2$ of the same value may be extremely close to each other and yet they can be rounded to different points $P_1\ne P_2$, which is undesirable in many applications. In this paper we show how to make the rounding process perfectly continuous in the sense that it maps any pair of sufficiently close measurements to the same point. We call such a process consistent rounding, and make it possible by allowing a small amount of information about the first measurement $X_1$ to be unidirectionally communicated to and used by the rounding process of $X_2$.The fault tolerance of a consistent rounding scheme is defined by the maximum distance between pairs of measurements which guarantees that they are always rounded to the same point, and our goal is to study the possible tradeoffs between the amount of information provided and the achievable fault tolerance for various types of spaces. When the measurements $X_i$ are arbitrary vectors in $R^d$, we show that communicating $\log_2(d+1)$ bits of information is both sufficient and necessary (in the worst case) in order to achieve consistent rounding for some positive fault tolerance, and when d=3 we obtain a tight upper and lower asymptotic bound of $(0.561+o(1))k^{1/3}$ on the achievable fault tolerance when we reveal $\log_2(k)$ bits of information about how $X_1$ was rounded. We analyze the problem by considering the possible colored tilings of the space with $k$ available colors, and obtain our upper and lower bounds with a variety of mathematical techniques including isoperimetric inequalities, the Brunn-Minkowski theorem, sphere packing bounds, and Čech cohomology.

\section{Introduction}

\paragraph{Rounding.} Whenever we want to digitally process or store an analog value, we have to perform some kind of rounding in order to represent this value (which is typically a number in $\mathbb{R}$) by a nearby point chosen from a discrete subset (such as the set of integers $\mathbb{Z}$). In the multidimensional generalization of this problem we are given a vector $X \in \mathbb{R}^d$, and want to represent it by a nearby point $P$ chosen from some discrete subset of $\mathbb{R}^d$ (which may be $\mathbb{Z}^d$, some lattice of points, or any other choice). Many rounding procedures had been proposed and studied in the literature (see, e.g.,~\cite{Balinsky96,CE82}), but all of them are inherently discontinuous in the sense that one can always find two inputs which are extremely close to each other (such as $0.49999$ and $0.50001$) which are mapped to different outputs (such as $0$ and $1$). In fact, any such rounding scheme partitions the space into a disjoint union of bounded subsets called {\em tiles}, in which each tile is rounded to a different representative point (e.g., its center), and the discontinuities happen along the tiles' boundaries. The algorithmic aspect of the problem of actually finding the tile that contains a given $X$ in a given partition is a classical problem in computational geometry, and had been studied extensively under the general name of {\em point location} (see~\cite{dBvKOS08,Snoeyink17}, and the references therein).

In many applications, such discontinuities are undesirable. One natural solution is to try to minimize the fraction of pairs $X_1, X_2$ with $distance (X_1,X_2)<\epsilon $ which are mapped to different tiles by considering {\em foam tilings} that minimize the total surface area of unit volume tiles (such a tiling is called `foam' since it emerges automatically in physical collections of soap bubbles). In a beautiful FOCS'2008 paper~\cite{KROW08} (which later appeared as a research highlight at CACM~\cite{KROW12}), Kindler, Rao, O'Donnell and Wigderson introduced a clever new construction of such tiles which they called {\em spherical cubes}. What makes these tiles special is that they have the $O(\sqrt{d})$ surface area of a ball and yet they can tile the whole $\mathbb{R}^d$ space without gaps.  Surprisingly, the main motivation for this construction came from the seemingly unrelated field of  computational hardness amplification, and it solved an interesting open problem posed by Lord Kelvin in $1887$.

A related class of schemes that aim at minimizing the proportion of close points mapped to different places is \emph{locality sensitive hashing} schemes (see, e.g.,~\cite{IMRV97}). However, they typically deal with situations in which the input is a vector and the output is a scalar, and thus they do not actually round inputs to nearby outputs in the same space.

\paragraph{Our approach -- consistent rounding.} In this paper we consider the more ambitious goal of completely eliminating all the discontinuities in the rounding process rather than reducing their fraction. We call any such scheme {\em consistent rounding}.\footnote{We note that in statistics, the term `consistent rounding' is used to denote a rounding that is consistent with some external constraints; see~\cite[p.~237]{WdW01}.} To make it possible, we have to slightly modify the model by thinking about $X_1$ and $X_2$ as two consecutive noisy measurements of the same $X$. When the first measurement $X_1$ is rounded, we allow it to produce a few bits of side information about how it was rounded and to provide them as an auxiliary input to the process that decides how to round $X_2$.

\medskip \noindent \emph{The one-dimensional case.} To demonstrate the basic idea, consider the one dimensional case in which $X_1$ and $X_2$ are real values which have to be consistently rounded to a nearby integer. $X_1$ is always rounded to the nearest integer, and it produces a single bit of side information which is whether it was rounded to an even or odd integer $P$. When $X_2$ is measured, it is rounded to the nearest integer which is of the same parity as $P$. To demonstrate this process, consider again the problematic inputs $X_1=0.49999$ and $X_2=0.50001$: $X_1$ is rounded to $0$, and $X_2$ is also rounded to $0$ since it is the closest even integer. In fact, $X_2$ could be anywhere between $-1$ and $1$ and still be consistently rounded to $0$, and thus the rounding scheme can tolerate any additive fault up to $0.5$. It is not difficult to show (see Section~\ref{ssec:BM}) that this is the highest possible fault tolerance of any one dimensional consistent rounding scheme, and that other natural schemes (such as providing one bit of side information about whether $X_1$ was rounded up or down) provide lower fault tolerance.

The way we think about the problem is to consider a colored tiling of the real line with two colors:
All the values in $[-0.5, 0.5)$ are colored by 1, all the values in $[0.5, 1.5)$ are colored by 2, all the values in $[1.5,2.5)$ are colored by 1, all the values in $[2.5,3.5)$ are colored by 2, etc. The side information provided about $X_1$ is the color of the tile in which it is located, and the way we process $X_2$ is to round it to the center of the closest tile which has the same color as that of $X_1$. The essential property of our partition is that the minimum distance between any two tiles which have the same color is $1$, and thus we can ``inflate'' all the tiles of a particular color to include any erroneously measured value $X_2$ up to a distance of $0.5$ away from the original tile, and still get nonoverlapping tiles which make it possible to uniquely associate such points with original tiles. This is reminiscent of how error correcting codes operate: To guarantee unique decoding, we ``inflate'' each codeword into a ball around it, and demand that all these balls will be disjoint. The error correction bound is thus half of the distance between the two closest codewords. The main difference is that in our case we inflate tiles rather than points, and try to maximize the smallest distance between tiles that have the same color (for each color separately).

\begin{figure}
	\centering
	\includegraphics[width=1.05\columnwidth]{Hexagons-joint.pdf}				\caption{A 3-colored hexagonal tiling of the plane, and a maximal non-intersecting inflation of the tiles colored 1}
	\label{fig:hexagons}
\end{figure}

\medskip \noindent \emph{The two-dimensional case.} To make this perspective clearer, consider the two dimensional plane. If we use the obvious checkerboard tiling by unit squares, we need at least 4 colors (and thus 2 bits of side information) to color the tiles  if we want any two tiles with the same color to have some positive distance between them. We can reduce the number of colors to 3 (and thus, provide only $\log_2(3)=1.58$ bits of side information) by considering the hexagonal partition of the plane depicted in the left part of Figure~\ref{fig:hexagons}. Given a two dimensional point $X_1$, we always round it to the center of the hexagon in which it is located, and given $X_2$ we round it to the center of the nearest hexagon which has the same color as the hexagon that contains $X_1$. To determine the fault tolerance of this consistent rounding scheme, we inflate all the hexagonal tiles of a particular color by the same amount until two of them touch each other for the first time, as depicted in the right part of Figure~\ref{fig:hexagons}. As it turns out, this scheme is not optimal: even though the hexagonal partition is very attractive since its circle-like tiles map each point to a nearby center, when we inflate the hexagons of a particular color we are stopped prematurely when their corners touch, leaving large gaps between the inflated hexagons. A 3-colored tiling with a higher fault tolerance will be described in Section~\ref{ssec:bricks}, and an asymptotically optimal tiling for a large number of colors will be described in Section~\ref{ssec:hcr}.

A fine point about consistent rounding schemes is how to break ties, and here we deal differently with $X_1$ and $X_2$. We want to be able to deal with any $X_1$, and thus we think about the tiles as being closed sets which include their boundaries. Therefore, points $X_1$ which are on the boundary between tiles can have more than one possible color. We allow such ties to be broken arbitrarily in the sense that $X_1$ can be rounded to the center of any one of the tiles that it belong to. However, when we think about $X_2$ we allow it to be at a distance of strictly less than some bound, and thus the inflated tiles (that contain all the possible $X_2$'s we are interested in) are open sets which have no intersections. Consequently, each $X_2$ can belong to at most one inflated tile, and is rounded to the center of that tile with no possible ties.

\paragraph{The minimal amount of side information that allows high-dimensional consistent rounding.} When we consider higher dimensional inputs $X \in \mathbb{R}^d$, we can provide one bit of side information about each one of its $d$ entries separately, but for large $d$ this requires a large number of side information bits. In our colored tiling formulation of the problem, it suffices to reveal the color of $X_1$ in order to consistently round $X_2$, and thus if we tile the space with $k$ colors, we need only $\log_2(k)$ bits of side information to specify this color. This naturally leads to the question of what is the minimum number of colors needed to tile $\mathbb{R}^d$ by bounded sized tiles so that all the tiles of the same color will have some positive distance from each other (note that two colors always suffice if we relax the problem by allowing an infinite series of nested boxes, or allow tiles of the same color to touch at corners). Surprisingly, we could not find in the literature any reference to this natural question. As we show in Section~\ref{sec:qualititave}, there can be no such colored tiling with $d$ colors, and as we show in Section~\ref{ssec:dim-reduce}, $d+1$ colors are sufficient. Consequently, $\log_2(d+1)$ bits of side information about $X_1$ are necessary and sufficient (in the worst case) in order to obtain a consistent rounding scheme with positive fault tolerance. To prove the negative result, we use techniques borrowed from algebraic topology (namely, the \v{C}ech cohomology and other cohomology theories), and to prove the positive result we provide an explicit construction of such a colored tiling.

\paragraph{The maximal fault tolerance that can be achieved with a given amount of side information.} In addition to aiming at minimizing the amount of side information, we study the question of maximizing the fault tolerance for a given amount of side information. In Section~\ref{sec:2dim} we study the special case of the two dimensional plane. In the negative direction, we obtain several upper bounds on the fault tolerance that can be achieved, using different techniques from geometry and analysis (including isoperimetry, the Brunn-Minkowski inequality and results on the circle packing problem). In the positive direction, we construct a variety of tiling schemes for various numbers of colors. In particular, while for three colors the hexagonal tiling scheme described above can tolerate an additive fault of up to $0.31$, we present a tiling with fault tolerance of $0.354$, and show that no 3-color tiling can have fault tolerance higher than $0.413$. When the number $k$ of colors tends to infinity, we derive the asymptotic bound of $(0.537 + o(1)) \sqrt{k}$ on the achievable fault tolerance, and prove its tightness by an explicit construction.

In Section~\ref{sec:3dim} we study the general case of $\mathbb{R}^d$, $d>2$. Like in the two-dimensional case, we obtain a number of lower and upper bounds on the fault tolerance. In particular, we show that the maximal fault tolerance achieved by a $(d+1)$-coloring of $\mathbb{R}^d$ is between $\Omega(1/d)$ and $O(\log d/\sqrt{d})$. In addition, we use the recent breakthrough results on the sphere packing problem~\cite{CKMRV17,Hales05,V17} to obtain tight asymptotic lower and upper bounds on the fault tolerance in dimensions $3,8,$ and $24$. Our bounds on fault tolerance are summarized in Table~\ref{tab:summary}.

\begin{table}[tb]
	\begin{tabular}{|c|c|c|c|c|}
		\hline
		Scenario & Lower Bound & Upper Bound & Techniques & Source \\
		& (LB) on FT  & (UB) on FT  & & \\
		\hline \hline
		3 colors & 0.354 & 0.413 & Brunn-Minkowski ineq.~(UB),
		& Sec.~\ref{ssec:BM} (UB), \\
		in $\mathbb{R}^2$ & & & Brick wall const.~(LB) & Sec.~\ref{ssec:bricks} (LB) \\
		\hline
		4 colors  & 0.5 & 0.564 & Brunn-Minkowski ineq.~(UB),  & Sec.~\ref{ssec:BM} (UB), \\
		in $\mathbb{R}^2$ & & & Brick wall const.~(LB) & Sec.~\ref{ssec:bricks} (LB) \\
		\hline
		$k$ colors & $(0.537+o(1))\sqrt{k}$ & $(0.537+o(1))\sqrt{k}$ & Circle packing~(UB), & Sec.~\ref{ssec:CP} (UB),
		\\
		in $\mathbb{R}^2$ & & & HCR const.~(LB) &	Sec.~\ref{ssec:hcr} (LB) \\
		\hline
		4 colors & 0.175 & 0.365 & Brunn-Minkowski ineq.~(UB), & Sec.~\ref{ssec:BM-high} (UB), \\
		in $\mathbb{R}^3$ & & & Bricks and balloons const.~(LB) & Sec.~\ref{ssec:bb} (LB) \\
		\hline
		$d+1$ colors & $\Omega(1/d)$ & $O(\log d/\sqrt{d})$ & Brunn-Minkowski ineq.~(UB), & Sec.~\ref{ssec:BM-high} (UB), \\
		in $\mathbb{R}^d$ & & & Dimension reducing const.~(LB) &  Sec.~\ref{ssec:dim-reduce} (LB) \\
		\hline
		$k$ colors & $(0.561+o(1))k^{1/3}$ & $(0.561+o(1))k^{1/3}$ & Sphere packing~(UB), & Sec.~\ref{ssec:SP-high} (UB), \\
		in $\mathbb{R}^3$ & & & CPB const.~(LB) & Sec.~\ref{ssec:CPB} (LB) \\
		\hline
		$k$ colors & $(0.707+o(1))k^{1/8}$ & $(0.707+o(1))k^{1/8}$ & Sphere packing~(UB), & Sec.~\ref{ssec:SP-high} (UB), \\
		in $\mathbb{R}^8$ & & & CPB const.~(LB) & Sec.~\ref{ssec:CPB} (LB) \\
		\hline
		$k$ colors & $(1+o(1))k^{1/24}$ & $(1+o(1))k^{1/24}$ & Sphere packing~(UB), & Sec.~\ref{ssec:SP-high} (UB), \\
		in $\mathbb{R}^{24}$ & & & CPB const.~(LB) & Sec.~\ref{ssec:CPB} (LB) \\
		\hline
		\multicolumn{5}{l}{FT -- fault tolerance,  LB -- lower bound, UB -- upper bound, Const. -- construction}\\
		\multicolumn{5}{l}{HCR -- honeycomb of rectangles, CPB -- close packing of boxes}
	\end{tabular}
	\caption{Summary of our lower and upper bounds on the fault tolerance, for different values of the dimension $d$ and the number of colors $k$}
	\label{tab:summary}
\end{table}

\paragraph{Applications.} The problem of consistent rounding has numerous applications. For example, consider the problem of biometric identification, in which multiple measurements of the same fingerprint are similar but not identical. It would be very helpful if all these slight variants could be represented by the same rounded point $P$, and we can achieve this by storing a small amount of side information in the biometric database (as part of the initial acquisition of the biometric data of a new employee).

Another example is the problem of developing a contact tracing app for the COVID-19 pandemic, where we want to record all the cases in which two smart phones were at roughly the same place at roughly the same time. We can do this by measuring in each phone the GPS location, the local time, and perhaps other parameters such as the ambient noise level or the existence of sudden accelerations (to determine that the two phones were in the same moving car). When someone tests positive for COVID-19, the health authority can release a list of his measurements, but in order to keep the patient's privacy, it wants to apply a cryptographically strong hash function to each measurement before publishing it. Since we expect these measurements to be slightly different for the infected and exposed persons, the health authority can publish the small amount of side information along with the consistently rounded and then hashed measurements.

Finally, in industrial control systems we can collect analog measurements from thousands of sensors (temperature, pressure, flow, etc.) and the use of a consistent rounding scheme could be ideal in order to check the repeatability of the process in spite of slight variations.

\paragraph{Related work.} A line of study that seems related to our work is \emph{fuzzy} constructions that were widely studied in the cryptographic literature. For example,~\cite{DORS08} introduces the notions of {\em fuzzy extractors} and {\em secure sketches}, which enable two parties to secretly reach a consensus value from multiple noisy measurements of some high entropy source (a recent survey of such techniques can be found in~\cite{FRS20}). However, such schemes concentrate on the aspects of cryptographic security (which we do not consider), and produce sketches whose size depends on the number of possible inputs (which is meaningless for real valued inputs). In this sense our consistent rounding scheme can be viewed as an exceptionally efficient reconciliation process, since it can produce for each million entry vector of arbitrarily large real numbers a $20$ bit ``sketch'' (in the form of its color side information), and process this information with trivial point location algorithms.

An actual application in which side information is used during an ``error reconciliation'' process can be found in the post quantum key exchange scheme {\em New Hope} by Alkim et.al.~\cite{ADPS16}. Their approach is limited to lattice points, and  requires side information whose size is linear in the number of dimensions, while our approach requires only a logarithmic number of bits.

Another natural idea (which is used in many fuzzy cryptographic constructions such as {\em fuzzy commitment}~\cite{JW99}) is to use an error correction scheme to map inputs to nearby codewords. However, most error correction schemes cannot be applied to real valued inputs. In addition, even if we could perfectly tile the continuous input space with balls surrounding each codeword, this would not solve the problem of consistent rounding near the boundary between adjacent balls. Finally, this boundary problem  becomes worse as the dimension grows, since in high dimensions almost all the points are likely to be near a boundary. Consequently, these code based cryptographic solutions manage to solve a variety of interesting related tasks, but not the consistent rounding problem we consider in this paper.

\paragraph{Open problems.} While we fully solved the question of minimizing the amount of side information required for fault tolerance, several questions remain open regarding the maximal tolerance rate that can be achieved for a given amount of side information. In particular, for dimensions $2,3,8,24$ we determined the exact asymptotic fault tolerance rate when $\log_2 k$ bits of information are allowed, using a connection to the \emph{densest sphere packing} problem. When only very few bits of information are allowed, the situation is much less clear. For example, we do not even know whether the brick wall constructions we present in Section~\ref{ssec:bricks} have the highest fault tolerance among rounding schemes in $\mathbb{R}^2$ with $\log_2 3$ and $\log_2 4$ bits of side information. It will be interesting to obtain new upper bounds via different techniques or new lower bound constructions.

\section{Our Setting}
\label{sec:notation}

In this section we present the basic setting that will be assumed throughout the paper.

\paragraph{Tiling.} We study tilings of $\mathbb{R}^d$, where each tile is \emph{connected}, \emph{closed} and \emph{bounded}, and the tiles intersect only in their boundaries. We further assume that the tiling is \emph{locally finite}, meaning that the number of tiles that intersect any bounded ball $B(0,r)$ is finite. In some of the results we make additional assumptions on the tiles; such assumptions are stated explicitly.

\paragraph{Coloring.} In our tiling, each tile is colored in one of $k$ colors. We note that the natural assumption that the tiles are unicolored is made for the sake of simplicity; all our results hold also if the tiles contain several colors, provided that the color classes inside each tile satisfy the assumptions we made on the tiles.

\paragraph{Fault tolerance and inflation.} In order to compute the fault tolerance of a given tiling (with respect to the $L_2$ distance), we consider all tiles of the same color and inflate them (i.e., replace the tile $T$ by the set $T'=\{y: \exists x \in T, |x-y|<r\}$ for some $r>0$) until they touch each other. Clearly, the fault tolerance is the maximal $r$ for which such a non-intersecting inflation is possible.

\medskip \noindent We denote the minimal distance between two same-colored points in different tiles by $t$, and so, the fault tolerance is $t/2$.

\paragraph{Normalization.} The $d$-dimensional volume of a figure $T \subset \mathbb{R}^d$ is denoted by $\lambda(T)$. We normalize the tiling by assuming that \emph{the volume of each tile is bounded by $1$} (like in the 1-dimensional case presented in the introduction, where all tiles are segments of length $1$). We make the natural assumption that any inflated tile $T'$ satisfies $\lambda(\bar{T}')=\lambda(T')$, where $\bar{T}'$ is the closure of $T'$.

We choose this normalization since it complies with using the Brunn-Minkowski theorem (see Section~\ref{ssec:BM}) and makes computations more convenient.
We note that our results can be easily translated to results with respect to
other natural normalizations. For example, the assumption that \emph{the radius of each tile is at most $1$} (which means that the $L_2$ distance between each $X \in \mathbb{R}^d$ and its rounded value is at most $1$) implies that the volume of each tile is at most the volume of the unit ball in $\mathbb{R}^d$, which allows translating all our fault tolerance upper bounds to this `bounded radius' setting. As for the lower bounds, they come from explicit constructions whose fault tolerance can be recomputed with respect to any other natural normalization.

\paragraph{Alternative metrics.} Instead of the $L_2$ distance (which is probably the most natural distance metric), one may consider the fault tolerance problem with respect to other distance metrics.

For example, in order to measure the fault tolerance with respect to the $L_\infty$ distance, one has to inflate each tile $T$ into $T''=\{y: \exists x \in T, \max_i |x_i-y_i|<r\}$. It turns out that the problem is much easier with respect to this metric. Indeed, assume that the number of colors is $k=m^d$ for some $2 \leq m \in \mathbb{N}$. Consider a periodic tiling in which each basic unit is a cube with side length $m$ that consists of $m^d$ unit-cube tiles, each colored in a different color. This tiling achieves fault tolerance of $(m-1)/2$. On the other hand, it is easy to see that the argument via the Brunn-Minkowski theorem presented in Section~\ref{ssec:BM} implies that any tiling of $\mathbb{R}^d$ with $k=m^d$ colors achieves fault tolerance of at most $(m-1)/2$ with respect to the $L_{\infty}$ metric, and thus, the cubic tiling we described is optimal. Lower and upper bounds for other values of $k$ can be obtained by variants of this construction and the corresponding upper bound proof.

The difference between the metrics comes from the fact that in `cubic' inflation (which is done with respect to the $L_{\infty}$ metric), non-intersecting inflations of cubic tiles fill the entire space, while in inflation by balls (which we have with respect to the $L_2$ metric), large gaps are left between the inflated tiles.

\section{The Minimal Number of Colors Required for Fault Tolerance}
\label{sec:qualititave}

In this section we prove that the minimal number of colors required for tolerating any rate $t>0$ of faults in a tiling of $\mathbb{R}^d$ is $d+1$. We provide two proofs, under different additional natural assumptions on the tiles. The first proof assumes that the tiles are polytopes and uses an inductive argument. The second proof assumes that the tiles and their non-empty intersections are contractible and uses an algebraic-topologic argument. The lower bound $d+1$ on the number of required colors is tight; a matching construction for any $d$ is presented in Section~\ref{ssec:dim-reduce}.

\subsection{Lower bound for polytopic tiles}
\label{sec:sub:polytopic}

In this section we prove the following.
\begin{proposition}\label{Prop:Crumpled}
	For any $m>0$, the following holds. Let $T_1,T_2,\ldots$ be a colored tiling of $\mathbb{R}^d$ in which the tiles are polytopes and each tile is contained in a box with side length $m$. If the fault tolerance of the tiling is $>0$, then the number of colors is at least $d+1$.
\end{proposition}

\paragraph{Informal proof.} We show by induction that there is a point $x$ that is included in at least $d+1$ tiles. This implies that if the tiling uses only $d$ colors, there must be two same-colored tiles that intersect in $x$, and thus, no fault tolerance is possible.

For $d=1$, as each single tile is included in a segment of length $m$, for any $M>m$, the entire segment $[-M,M]$ cannot be covered by a single tile. Hence, the segment must contain a point in which two tiles intersect.

For the induction step, we consider the restriction of the tiling to a large box $[-M,M]^d$. Then we look at the $d$'th coordinate, consider the union of all tiles that contain a point in the `lower' half-box, and take the `upper boundary' of this union. This gives a `crumpled' variant of the equator hyperplane. By the induction hypothesis, the restriction of the tiling to this `crumpled hyperplane' contains a point $v$ that is included in at least $d$ tiles. However, there is also at least one tile in the `upper half-box' that includes $v$, and thus, $v$ is included in at least $d+1$ tiles. This completes the proof by induction.

\paragraph{Formal proof.} Making the argument described above formal is somewhat cumbersome, and requires an auxiliary notion, somewhat close to the notion of a \emph{polytopal complex}~\cite[p.~127]{Ziegler95}.

\paragraph{Notation.} Let $k,d,m,M \in \mathbb{N}$, be such that $d \geq k$ and $M>m$. A \emph{crumpled $(k,d)$-flat} $F$ is a union of a finite collection of $k$-dimensional polytopes $P_1,P_2,\ldots,P_{\ell} \in \mathbb{R}^d$ that satisfies the following conditions:
\begin{enumerate}
	\item $F$ is contained in the slab $[-M,M]^k \times [0,m]^{d-k}$.

	\item The projection of $F$ onto the hyperplane $H_{k}=\mathbb{R}^k \times \{0\}^{d-k}$ is $[-M,M]^k \times \{0\}^{d-k}$.

	\item $F$ is connected.

	\item Each $P_i$ is contained in a box with side length $m$.

	\item The intersection of any $P_i,P_j \in F$ is either empty or a polytope of dimension~$\leq k-1$.
\end{enumerate}

\begin{figure}
	\centering
	\includegraphics[width=0.5\columnwidth]{Crumpled.pdf}
	\caption{An illustration of the inductive process of Lemma~\ref{Lem:Aux3}. A polygonal tiling of $[-3,3]^2$ is depicted in ordinary lines, and the `crumpled equator' to which we reduce in the inductive process (that has to be included in the slab $[-3,3]\times[0,2]$ and whose projection on the $x$-axis must be $[-3,3]$) is depicted by a bold line. The bold line is a crumpled $(1,2)$-flat with the parameters $M=3,m=2$.}
	\label{fig:crumpled}
\end{figure}

For example, a crumpled $(1,d)$-flat is a polygonal line included in the slab $[-M,M] \times [0,m]^{d-1}$, whose projection on the $x$-axis is $[-M,M]$. A crumpled $(1,2)$-flat with the parameters $M=3,m=2$ is demonstrated in Figure~\ref{fig:crumpled} (in boldface).

\begin{lemma}\label{Lem:Aux3}
	Let $k,d,m,M$ be as above. Any crumpled $(k,d)$-flat $F$ contains a point that belongs to at least $k+1$ polytopes.
\end{lemma}

\begin{proof}
	The proof goes by induction on $k$. For $k=1$, observe that $F$ cannot consist of a single polytope. Indeed, as the projection of $F$ on the $x$-axis is $[-M,M]$, there are two points $v,v' \in F$ with $|v-v'|\geq 2M$. These points cannot belong to the same polytope $P_i$, since each $P_i$ is included in a segment of length $m$. As $F$ is connected, it must contain an intersection point of two $P_i$'s.

	For the induction step, consider a crumpled $(k,d)$-flat $F$. Let $F'$ be the union of all polytopes in $F$ whose intersection with the set $\{x \in \mathbb{R}^d: x_k \leq 0\}$ is non-empty, and let
	\[
	F'' = \partial(F') \cap \{x \in \mathbb{R}^d:x_k \geq 0\}.
	\]
	(That is, we look at the $k$'th coordinate and take to $F'$ all polytopes in $F$ that contain a point in the `lower' half-space. Then, we define $F''$ to be the intersection of the boundary of $F'$ with the `upper' half-space. We view $F''$ as a `crumpled equator flat' of $F$). This process, in the case $d=2$, $k=1$, is demonstrated in Figure~\ref{fig:crumpled}, where the tiling $F$ is depicted in ordinary lines and $F''$ is depicted by a bold line.

	We view $F''$ as the union of the polytopes $P'_i = P_i \cap F''$. We claim that $F''$ is a crumpled $(k-1,d)$-flat. Indeed:
	\begin{enumerate}
		\item $F''$ is contained in the slab $[-M,M]^{k-1} \times [0,m]^{d-k+1}$, since each $v \in \partial(F')$ belongs to some $P_i$ that contains also a point $v'$ with $v'_k \leq 0$. As each $P_i$ is included in a box of side length $m$, we have $v'_k \leq m$. (The condition in all other coordinates clearly holds.)

		\item For any fixed value $(x_1,\ldots,x_{k-1}) \in [-M,M]^{k-1}$, the projection of the set  $\{v \in F: (v_1,\ldots,v_{k-1})=(x_1,\ldots,x_{k-1})\}$ on the $k$'th axis is $[-M,M]$. In particular, $F'$ contains a point whose first $k$ coordinates are $(x_1,\ldots,x_{k-1},0)$. We can take such a point, and move from it in the positive direction of the $k$'th axis until we reach a point on the boundary of $F'$. The resulting point is an element in $F''$ whose first $k-1$ coordinates are $x_1,\ldots,x_{k-1}$. Hence, the projection of $F''$ on $H_{k-1}$ is $[-M,M]^{k-1} \times \{0\}^{d-k+1}$.

		\item $F'$ is connected, and hence, $F'' = \partial(F') \cap \{x \in \mathbb{R}^d:x_k \geq 0\}$ is connected as well.

		\item Each polytope in $F''$ is contained in a box of side length $m$, since it is part of a polytope that was contained in $F'$.

		\item The intersection of any $P'_i,P'_j \in F''$ is either empty or a polytope, since $P'_i,P'_j$ are parts of faces of polytopes (from $F$) created by intersection with other polytopes or with the halfspace $\{x \in \mathbb{R}^d: x_k \geq 0\}$. The dimension of any such intersection is at most $k-2$ since the transformation from $H'$ to its boundary reduces all dimensions by 1.
	\end{enumerate}

	By the induction hypothesis, $H''$ contains an intersection point $v$ of $k$ polytopes $P'_i$. This means that the corresponding polytopes $P_i$ contain $v$ as well.

	\medskip Furthermore, as $v \in \partial(H')$ and $0 \leq v_k \leq m$, there exists $P_j \in H \setminus H'$ that contains $v$. (Here, we use the assumption that for any sufficiently small $\epsilon>0$, $H$ contains some $w$ with $(w_1,\ldots,w_k)=(v_1,\ldots,v_{k-1},v_k+\epsilon)$ and the assumption that $F$ consists of a finite collection of polytopes.) Therefore, $v$ belongs to at least $k+1$ polytopes of $H$, completing the proof by induction.
\end{proof}

Now we are ready to prove Proposition~\ref{Prop:Crumpled}.

\medskip

\begin{proof}
	Let $T_1,T_2,\ldots$ be a tiling of $\mathbb{R}^d$ that satisfies the assumptions of the proposition. Let $T'_i = T_i \cap [-M,M]^d$, and let $F=\cup T'_i$. It is easy to see that $F$ is a crumpled $(d,d)$-flat. Hence, by the lemma, $F$ contains a point $v$ that belongs to $d+1$ polytopes. The point $v$ is included in at least $d+1$ tiles $T_i$, proving the assertion.

\end{proof}

\paragraph{Remark.}	We note that if one assumes that the tiles are \emph{convex} polytopes, then a much simpler argument shows that at least $d+1$ colors are needed. Indeed, start with a point in some tile, and move in some direction until you reach the boundary of the tile. The points on this boundary belong to two tiles. Move along this boundary face until you reach a point on its boundary (i.e., now we reduce from dimension $d-1$ to $d-2$.) The points on this boundary already belong to at least three tiles. We can continue reducing dimension, until we reach dimension $0$ with a point that belongs to $d+1$ tiles. This argument fails in nonconvex polytopes since the corners of a box nested inside a bigger box touch only two tiles.

\subsection{Lower bound for contractible tiles}
\label{sec:sub:contractible}

Recall that a set in $\mathbb{R}^d$ is called \emph{contractible} if is can be continuously shrunk to a point in $\mathbb{R}^d$. (The formal definition is that the identity map is homotopic to some constant map.)

Informally, in this section we prove that if the tiles and their non-empty intersections are finite unions of contractible sets, then at least $d+1$ colors are required for fault tolerance. Due to the possibility of pathologies, the formal statement is a bit more cumbersome:
\begin{proposition}\label{Prop:Contractible}
Let $T_1,T_2,\ldots$ be a colored tiling of $\mathbb{R}^d$ in which the tiles and all their non-empty intersections are disjoint unions of finitely many contractible sets. Assume that the tiling is locally finite and that all $T_i$'s are bounded. In addition, assume that each $T_i$ has an open neighborhood $U_i$ such that for any set of indices $I,$
\[\bigcap_{i\in I}U_i\neq\emptyset\Leftrightarrow\bigcap_{i\in I}T_i\neq\emptyset,\]
and the $U_i$'s and their non-empty intersections are disjoint unions of finitely many contractibles.

Then the number of colors is at least $d+1$.
\end{proposition}

A similar method proves an analogous statement for colored tilings of the sphere $\mathbb{S}^d$:
\begin{proposition}\label{Prop:Contractible_Sn}
	Let $T_1,T_2,\ldots,T_N$ be a colored tiling of $\mathbb{S}^d$ in which the tiles and all their non-empty intersections are disjoint unions of finitely many contractible sets. Assume that each $T_i$ has an open neighborhood $U_i$ such that for any set of indices $I,$
	\[\bigcap_{i\in I}U_i\neq\emptyset\Leftrightarrow\bigcap_{i\in I}T_i\neq\emptyset,\]
	and the $U_i$'s and their non-empty intersections are disjoint unions of finitely many contractibles.

	Then the number of colors is at least $d+1$.
\end{proposition}

\begin{remark}
	We stress that for most natural tilings the additional assumption on the existence of the neighborhoods $U_i$ follows from the existence of $T_i$'s with the corresponding properties. However, there are topological pathologies in which this is not the case.
\end{remark}

\iffalse
\paragraph{Reduction to a topologic formulation.}\footnote{Nathan: I suggest to erase the reduction, in all "natural cases" given $T_i$ your suggestion of $U_i$ will work, but there are pathologies. } We observe that Proposition~\ref{Prop:Contractible} is implied by the following topologic proposition.
\begin{proposition}\label{Prop:Contractible-Top}
	Let $U_1,U_2,\ldots$ be an open covering of $\mathbb{R}^d$, where the $U_i$'s and all their non-empty intersections are unions of contractible sets. Then there exists a point that belongs to at least $d+1$ $U_i$'s.
\end{proposition}
To see that Proposition~\ref{Prop:Contractible} follows from Proposition~\ref{Prop:Contractible-Top}, consider a tiling $T_1,T_2,\ldots$ that satisfies the assumptions of Proposition~\ref{Prop:Contractible}. Let $U_i$ be an $(\epsilon/2)$-neighborhood of $T_i$, that is,
\[
U_i = \{x \in \mathbb{R}^d: \exists y \in T_i, |x-y|<\epsilon/2\}.
\]
It is clear that $U_1,U_2,\ldots$ is an open cover of $\mathbb{R}^d$, and that (assuming $\epsilon$ is sufficiently small) all $U_i$'s and their non-empty intersections are unions of contractible sets. Hence, by Proposition~\ref{Prop:Contractible-Top}, there exists a point $x$ that belongs to $d+1$ of the $U_i$'s. Consequently, there exist $d+1$ points $x_1,\ldots,x_{d+1}$ in different tiles such that $\forall i,j: |x_i-x_j|<\epsilon$. Therefore, in order to have distance of at least $\epsilon$ between two same-colored points in different tiles, at least $d+1$ colors are needed.
\fi

\medskip \noindent The proof of Propositions~\ref{Prop:Contractible} and~\ref{Prop:Contractible_Sn} uses the notion of \v{C}ech cohomology and classical results regarding its properties. For the ease of reading, we begin with an intuitive explanation of the proof ideas, and then present the formal proof.

\paragraph{Intuitive proof.}
The $d$'th (singular) cohomology group is a topological invariant of a manifold which roughly counts ``non trivial holes'' of dimension $d$. A classical result asserts that the $d$'th cohomology group of a $d$-dimensional compact oriented manifold like $\mathbb{S}^d$ is $\mathbb{R}$. (This corresponds to the intuitive understanding that $\mathbb{S}^d$ has one $d$-dimensional hole.) The de-Rham cohomology and the \v{C}ech cohomology are analytic and algebro-geometric/combinatorial invariants, that in many cases agree with their topological cousin. In particular, the $d$'th de-Rham and \v{C}ech cohomologies of $\mathbb{S}^d$ are equal to $\mathbb{R}$ as well.

The $d$'th \v{C}ech cohomology with respect to an open cover of the manifold depends on properties of intersections of $d+1$ sets in that cover. In general, it depends on the sets which form the cover, however, it is known that if these sets and their non-empty intersections are finite disjoint unions of contractibles, then the cohomology groups remain the same, independently of the cover. In particular, if the $d$'th \v{C}ech cohomology with respect to such a cover is non trivial, then there must be $d+1$ sets with a non-empty intersection.

Hence, for our cover $U_1,U_2,\ldots$, we know that its $d$'th \v{C}ech cohomology is $\mathbb{R}$. This readily completes the proof of the proposition for $\mathbb{S}^d$, as this implies that there must be a point that belongs to at least $d+1$ of the $U_i$'s.

The proof in $\mathbb{R}^d$ works in essentially the same way, with cohomology groups replaced by cohomology groups with compact support.

\paragraph{Formal proof.}
For the proof we recall the notion of \v{C}ech cohomology with values in the constant sheaf $\underline{\Re}$, and describe the slightly less standard concept of \v{C}ech cohomology with compact support.

\paragraph{Definitions.} Let $S$ be either $\Re^{d}$ or a compact manifold such as $\mathbb{S}^{d}.$
Let $\mathcal{U}=\{U_1, U_2,\ldots\}$ be an open cover of $S.$ If $S$ is compact, we assume the collection to be finite. If $S$ is $\Re^{d}$, we assume it to be locally finite and assume in addition that each $U_i$ is bounded.
\begin{itemize}
	\item A \emph{$q-$simplex} $\sigma=\{U_{i_0},\ldots, U_{i_k}\}$ of ${\mathcal {U}}$ is an ordered collection of $q+1$ different sets chosen from $\mathcal {U},$ such that \[\bigcap_{k=0}^{q}U_{i_k}\neq\emptyset.\]

	\item For a $q-$simplex $\sigma =(U_{i_k})_{k\in \{0,\ldots ,q\}}$, the $j$'th \emph{partial boundary} is the $(q-1)$-simplex
	\[\partial_{j}\sigma :=(U_{i_k})_{k\in \{0,\ldots ,q\}\setminus \{j\}},\]
	obtained by removing the $j$'th set from $\sigma$.

	\item A \emph{$q-$cochain} of $\mathcal{U}$ is a function which associates to any $q-$simplex a real number. The $q-$cochains form a vector space denoted by $C^{q}(\mathcal{U},\underline{\Re}),$ with operations
	\[(\lambda f +\mu g)(\sigma)=\lambda f(\sigma)+\mu g(\sigma), \quad \mbox{where} \quad \lambda,\mu\in\Re,~f,g\in C^{q}(\mathcal{U},\underline{\Re}),~\sigma~\mbox{ is a $q-$simplex}.\]
	Similarly, we define $C^{q}_c(\mathcal{U},\underline{\Re}),$ as the vector space of $q-$cochains \emph{with compact support}, meaning those cochains which assign $0$ to all $q-$simplices, except for finitely many.

	\item There is a \emph{differential map} $\delta_{q}:C^{q}(\mathcal{U},\underline{\Re})\to C^{q+1}(\mathcal{U},\underline{\Re})$ whose application to $f\in C^{q}(\mathcal{U},\underline{\Re})$ is the $(q+1)-$cochain $\delta_q(f)$ whose value at a $(q+1)-$simplex $\sigma$ is
	\[(\delta_{q}f)(\sigma)=\sum_{j=0}^{q+1}(-1)^{j}f(\partial_j\sigma).\] The restriction of $\delta_q$ to $C^{q}_c(\mathcal{U},\underline{\Re})$ maps it to $C^{q+1}_c(\mathcal{U},\underline{\Re})$.

	\item  It can be easily seen that
	\[\delta_{q+1}\circ\delta_q=0.\]

	\item The \emph{$q$'th \v{C}ech cohomology group (with compact support)} of $S$ with respect to the cover $\mathcal{U}$ and values in $\underline{\Re}$ is
	\[\check{{H}}^{q}(\mathcal{U},\underline{\Re}):=\mbox{Ker}(\delta_q)/\mbox{Image}(\delta_{q-1}),\]
	\[(\check{{H}}^{q}_c(\mathcal{U},\underline{\Re}):=\mbox{Ker}(\delta_q|_{C^{q}_c(\mathcal{U},\underline{\Re})})/\mbox{Image}(\delta_{q-1}|_{C^{q-1}_c(\mathcal{U},\underline{\Re})})).\]

	\item A cover (by open sets) is \emph{good} if all its sets as well as their multiple intersections are either empty or contractible. It is \emph{almost good} if all non empty intersections are unions of finitely many disjoint contractible components.
\end{itemize}

\paragraph{Classical results we use.} The first result we use is the following:
\begin{theorem}\label{theorem:cech_equals_dR}
If $S$ is a compact smooth orientable manifold (such as $\mathbb{S}^{d}$), and $\mathcal{U}$ is a good or an almost good finite cover, then
\[\check{H}^{i}(\mathcal{U},\underline{\Re})\simeq {H}^{i}_{dR}(S),\] where the right hand side is the standard de-Rham cohomology group.

\medskip \noindent
Similarly, if $S=\Re^{d}$ and $\mathcal{U}$ is a locally finite good or almost good cover whose sets are bounded, then
\[\check{H}^{i}_c(\mathcal{U},\underline{\Re})\simeq {H}^{i}_{dR,c}(S),\] where the right hand side is the $i$'th de-Rham cohomology group with compact support.
\end{theorem}

For further reading about de-Rham cohomology, with or without compact support, we refer the reader to \cite[Sec.~1]{BT}. For further reading about the \v{C}ech cohomology, we refer to \cite[Sec.~8]{BT}. In particular, Theorem \ref{theorem:cech_equals_dR}, for the compact case and good covers is Theorem~8.9 there. The passage to almost good covers is straightforward: In the paragraph which precedes the proof, it is explained that the obstructions to the isomorphism between \v{C}ech and de-Rham cohomologies are given by products of the $i$'th de-Rham cohomology groups, for $i\geq 1,$ of the different intersections $\bigcap_{k=0}^{q} U_{i_k}.$ Since those intersections are disjoint unions of contractibles, their higher cohomology groups vanish, hence there is no obstruction to the isomorphism.

Regarding the case $S=\Re^{d},$ the proof in \cite[Sec.~8]{BT} requires a few small changes: In the statement of Proposition 8.5 there, one needs to replace the de-Rham complex of the manifold with the de-Rham complex with compact support, and the direct product with direct sum.
The maps $r,\delta$ which appear there will still be well defined by our local finiteness assumption on the cover, and the assumption that $U_i$'s are bounded. The proof requires no change.
Then, the double complex in the definition of Proposition~8.8 should also be defined using direct sum rather than direct product, but again there is no change in the proof. Given these changes in definitions, the proof of Theorem~8.9 (also for the almost good case) is unchanged.

The second standard result, which is a consequence of Poincar\'e duality, is the following:
\begin{theorem}\label{theorem:PD}
For a compact smooth oriented manifold $S$ of dimension $d$ (such as $\mathbb{S}^{d}$),
\[{{H}}^{d}_{dR}(S)\simeq\Re.\]
Similarly, for $S=\Re^{d},$
\[{H}^{d}_{dR,c}(\Re^{d})\simeq \Re.\]
\end{theorem}
See, for example, \cite[Sec.~7]{BT} for the compact case, and~\cite[Sec.~4]{BT} for $\Re^{d}.$

\medskip Theorems~\ref{theorem:cech_equals_dR} and~\ref{theorem:PD} yield:
\begin{corollary}\label{cor:non_vanish}
If $S$ is a compact smooth orientable manifold (such as $\mathbb{S}^{d}$), and $\mathcal{U}$ is a good or an almost good finite cover, then
\[\check{{H}}^{d}(\mathcal{U},\underline{\Re})=\Re.\]
Similarly, if $S=\Re^{d}$ and $\mathcal{U}$ is a locally finite good or almost good cover whose sets are bounded, then
\[\check{{H}}^{d}_c(\mathcal{U},\underline{\Re})=\Re.\]
\end{corollary}

\paragraph{Proof of Propositions~\ref{Prop:Contractible} and~\ref{Prop:Contractible_Sn}.} We show that there must exist $d+1$ $T_i$'s whose intersection is non-empty. This clearly implies that for any fault tolerance, at least $d+1$ colors are needed.

\medskip Assume on the contrary that any $(d+1)$-intersection of the $T_i$'s is empty. Let $U_i$ be as in the statement of the theorem. Then by definition, they form an almost good cover. All intersections of at least $d+1$ $U_i$'s are empty  by our assumptions. Therefore, there are no $d-$simplices, and so $C^{d}(\mathcal{U},\underline{\Re})=0.$ Thus, in the compact case, $\check{{H}}^{d}(\mathcal{U},\underline{\Re})=0.$ But on the other hand, by Corollary \ref{cor:non_vanish}, \[\check{{H}}^{d}(\mathcal{U},\underline{\Re})\simeq\Re,\] a contradiction. In the case of $\Re^{d}$ the same argument works, with $\check{{H}}^{d}_c(\mathcal{U},\underline{\Re})$ in place of $\check{{H}}^{d}(\mathcal{U},\underline{\Re})$.

\section{Fault-Tolerant Rounding in the Plane}
\label{sec:2dim}

In this section we consider tilings of the plane by tiles $T_1,T_2,\ldots$ of area at most $1$. Each point in the plane is colored in one of $k \geq 3$ colors, and our goal is to maximize the minimal distance $t$ between two points of the same color that belong to different tiles. (The maximum is taken over all possible tilings that satisfy the mild regularity conditions stated in Section~\ref{sec:notation} and over all possible colorings.) Clearly, the fault tolerance of a rounding scheme based on such a colored tiling is $t/2$.

We present three upper bounds on $t$, using the Brunn-Minkowski inequality, the circle packing problem, and the Minkowski-Steiner formula, respectively. In the other direction, we prove lower bounds on $t$ by a series of explicit tilings. In particular, we obtain a tight asymptotic estimate for $t$, namely, $t=(\beta+o(1))\sqrt{k}$, where $\beta=\sqrt{2/\sqrt{3}} \approx 1.074$. Our main upper and lower bounds are summarized in Table~\ref{tab:summary}.

\subsection{Upper bounds on the fault tolerance}
\label{sec:sub:upper}

The basic idea behind our upper bound proofs is as follows. Assume we have a colored tiling of the plane, with minimal distance $t$. Pick a single color -- say, black -- and consider all black tiles inside a large square $S$.
We obtain a new collection of tiles $T'_1,T'_2,\ldots,T'_m$ that covers part of $S$. Now, the assumption that the minimal distance between two same-colored points in different tiles is $t$ implies that if we inflate each black tile $T'_i$ into the set
\begin{equation}\label{Eq:Blowup1}
T''_i =\{x: \exists y \in T'_i, |x-y|<t/2\},
\end{equation}
then the inflations $T''_i$ are pairwise disjoint. Hence, the sum of their areas essentially cannot exceed the area of the large square, and this allows bounding $t$ from above.

\subsubsection{An upper bound using the Brunn-Minkowski inequality}
\label{ssec:BM}

The inflations $T''_i$ have a convenient representation in terms of the \emph{Minkowski sum} of sets in $\mathbb{R}^2$.
\begin{definition}
	For two sets $A,B \subset \mathbb{R}^d$, the Minkowski sum of $A,B$ is $A+B=\{a+b:a \in A, b \in B\}$.
\end{definition}
In terms of this definition, we have
\begin{equation}\label{Eq:BM1}
T''_i = T'_i + B(0,t/2),
\end{equation}
where $B(0,t/2)$ is a disk of radius $t/2$ around the origin. This allows us to lower bound the area of each $T''_i$, using the classical Brunn-Minkowski (BM) inequality (see, e.g.,~\cite{Busemann58}). Recall the inequality asserts the following.
\begin{theorem}[Brunn-Minkowski]
	Let $A,B$ be compact sets in $\mathbb{R}^d$. Then
	\[
	\lambda(A+B)^{1/d} \geq \lambda(A)^{1/d}+\lambda(B)^{1/d},
	\]
	where $\lambda(X)$ is the volume of $X$ (formally, the $d$-dimensional Lebesgue measure of $X$).

	\medskip \noindent The inequality is tight if and only if $A,B$ are positive homothetic, i.e., are equal up to translation and dilation by a positive factor.
\end{theorem}

\begin{proposition}\label{Prop:BM}
	Let $T_1,T_2,\ldots$ be a $k$-colored tiling of the plane, with tiles of area $\leq 1$ and minimal distance $t$. Then
	\[
	t \leq \left(\frac{2}{\sqrt{\pi}}+o(1) \right) \cdot (\sqrt{k}-1) \approx 1.128(\sqrt{k}-1).
	\]
\end{proposition}

\begin{proof}
	Consider a square $S$ such that $\lambda(S)=n^2$ (for some `large' $n$). By the pigeonhole principle, there exists a color (say, black) that covers at least $n^2/k$ of the area of $S$. Look at the black tiles whose intersection with $S$ is non-empty, and denote their intersections with $S$ by $T'_1,T'_2,\ldots,T'_m$. Hence, we have $m$ `black' subsets of $S$, each of area at most 1, whose total area is at least $n^2/k$.

	For each $T'_i$, define $T''_i = T'_i + B(0,t/2)$. By assumption, the regions $T''_i$ are disjoint. Furthermore, they are included in $S+B(0,t/2)$ whose area is less than $(n+t)^2$. By taking a sufficiently large $n$, we may assume
	\begin{equation}\label{Eq:BM2}
	\sum_i \lambda(T''_i) \leq (1+o(1))n^2.
	\end{equation}
	By the Brunn-Minkowski inequality, we have
	\[
	\forall i: \sqrt{\lambda(T''_i)} \geq \sqrt{\lambda(T'_i)} + \sqrt{\pi} \cdot \frac{t}{2},
	\]
	and thus,
	\[
	\forall i: \lambda(T''_i) \geq \lambda(T'_i) + \sqrt{\pi} \cdot t \cdot \sqrt{\lambda(T'_i)} + \pi t^2/4.
	\]
	Summing over $i$ and using~(\ref{Eq:BM2}), we get
	\begin{equation}\label{Eq:BM3}
	(1+o(1))n^2 \geq \sum_{i=1}^m \lambda(T''_i) \geq \sum_{i=1}^m \lambda(T'_i) + \sqrt{\pi} \cdot t \cdot \sum_{i=1}^m \sqrt{\lambda(T'_i)} + m \pi t^2/4.
	\end{equation}
	Note that we have $\sum_i \lambda(T'_i) \geq n^2/k$ (by assumption), $\sum_i \sqrt{\lambda(T'_i)} \geq n^2/k$ (since $\forall 0 \leq x \leq 1: \sqrt{x} \geq x$), and $m \geq n^2/k$ (since $\forall i: \lambda(T'_i) \leq 1$). Hence, we get
	\[
	(1+o(1))n^2 \geq \frac{n^2}{k} \cdot (1+\sqrt{\pi} \cdot t + \pi t^2/4) = \frac{n^2}{k} \cdot (1+\sqrt{\pi}t/2)^2.
	\]
	Dividing by $n^2$ and rearranging, we get
	$(1+o(1))(\sqrt{k}-1) \geq \sqrt{\pi}t/2$,
	and hence,
	\[
	t \geq \left(\frac{2}{\sqrt{\pi}}+o(1) \right) \cdot (\sqrt{k}-1) \approx 1.128(\sqrt{k}-1),
	\]
	as asserted.
\end{proof}

\paragraph{Discussion.} Substituting specific values of $k$, the upper bounds that follow from Proposition~\ref{Prop:BM} are $t \leq 0.826$ for $k=3$, $t \leq 1.128$ for $k=4$ and $t \leq 1.128\sqrt{k}$ for a large $k$.

These upper bounds are lossy in two ways. One source of loss is the application of the Brunn-Minkowski inequality. Here, the inequality is tight if the tiles are \emph{disks}, and the farther they are from disks, the larger is the loss. Another source of loss is the space left between the inflations, that is not taken into account in the proof.

Interestingly, there is a dichotomy between these two sources of loss. As follows from the circle packing problem, when the tiles are disks (and so, there is no loss in the BM inequality), the space between the inflations (and so, the loss of the second type) is relatively large. The space between the inflations can be made smaller if the tiles are taken to be polygons with a few vertices. However, this comes at the expense of increased loss in the BM inequality, as is demonstrated in Section~\ref{ssec:Steiner}.

\paragraph{Optimality of our 1-dimensional rounding scheme.} The argument described above gives an easy proof of the optimality of the 1-dimensional rounding scheme presented in the introduction. Indeed, consider a 2-colored tiling of the line and look at the segment $I = [-n,n]$ for some large $n$. By the pigeonhole principle, we may assume that black tiles cover at least half of $I$. By the 1-dimensional Brunn-Minkowski inequality, for each black tile $T'_i \subset I$ and the corresponding inflation $T''_i=T'_i + (-t/2,t/2)$, we have $\lambda(T''_i) \geq \lambda(T'_i)+t$. As $\forall i: \lambda(T'_i) \leq 1$, there are at least $n$ tiles. Since the $T''_i$'s are pairwise disjoint and included in $[-n-1,n+1]$, we obtain
\[
2n+2 \geq \sum_i \lambda(T''_i) \geq \sum_i \lambda(T'_i)+ \sum_i t \geq n + nt,
\]
and thus, $t \leq (n+2)/n$. By letting $n$ tend to infinity, we obtain $t \leq 1$, implying that the fault tolerance of any two-colored tiling of the line is at most $1/2$.

\subsubsection{An upper bound using the Circle Packing problem}
\label{ssec:CP}

This upper bound is interesting mainly for a large number $k$ of colors, where the size and form of the original tiles is negligible with respect to the size of the inflations. In this setting, we can represent each tile by a single point from it, and consider our problem as an instance of the well-known \emph{circle packing} problem (see, e.g.,~\cite{CS13}), as we know that circles of radius $t/2$ around the points must be disjoint. We use the following classical result of Fejes-T{\'{o}}th~\cite{FT50}:
\begin{theorem}[Fejes-T{\'{o}}th]\label{Thm:FT}
	The maximal number of points that can be placed in a square of area $a$, such that the minimal distance between two points is $1$, is at most $2a/\sqrt{3}$.
\end{theorem}
By rescaling, the theorem implies:
\begin{corollary}\label{Cor:CP}
	The maximal number of points that can be placed in a square of area $a$, such that the minimal distance between two points is $t$, is at most $\frac{2a}{\sqrt{3}t^2}$.
\end{corollary}

\begin{proposition}\label{Prop:CP}
	Let $T_1,T_2,\ldots$ be a tiling of the plane in $k$ colors, with tiles of area $\leq 1$ and minimal distance $t$. Then
	\[
	t \leq \left(\sqrt{\frac{2}{\sqrt{3}}}+o(1) \right) \cdot \sqrt{k} \approx 1.074\sqrt{k}.
	\]
\end{proposition}

\begin{proof}
	Like in the proof of Proposition~\ref{Prop:BM}, we begin with a tiling $T_1,T_2,\ldots$ with minimal distance $t$, consider the intersection with a square $S$ of area $n^2$, pick a single color (say, black) whose intersection with $S$ has area of at least $n^2/k$, and obtain disjoint black tiles $T'_1,\ldots,T'_m$ inside a square $S$ of area $(1+o(1))n^2$, whose total area is at least $n^2/k$. We pick a single point $x_i$ from each tile $T'_i$. By assumption, the minimal distance between the $x_i$'s is at least $t$. By Corollary~\ref{Cor:CP}, this implies
	\[
	m \leq (1+o(1)) \cdot \frac{2n^2}{\sqrt{3}t^2}.
	\]
	On the other hand, $m \geq n^2/k$, since the total area of the black tiles is at least $n^2/k$ and the area of each tile is at most $1$. Therefore,
	\[
	\frac{n^2}{k} \leq (1+o(1)) \cdot \frac{2n^2}{\sqrt{3}t^2},
	\]
	and consequently, $t \leq \left(\sqrt{\frac{2}{\sqrt{3}}}+o(1) \right) \cdot \sqrt{k}$, as asserted.
\end{proof}

\paragraph{Discussion.} For small values of $k$, this upper bound is rather weak -- e.g., for $k=3$ we obtain $t \leq 1.861$ (compared to $t \leq 0.826$ via the BM inequality). This happens, since we completely neglect the tiles and consider only the area added in the inflation process. For a large $k$, when the area is dominated by the inflation step, the upper bound $t \leq (1.074+o(1))\sqrt{k}$ is tight, as will be shown in Section~\ref{ssec:hcr}.

\subsubsection{An upper bound using the Minkowski-Steiner formula}
\label{ssec:Steiner}

This upper bound is applicable under the additional condition that the tiles $T_i$ are convex. It uses another well-known result, called the \emph{Minkowski-Steiner formula}.
\begin{theorem}[Minkowski-Steiner]\label{Thm:MS}
	Let $A$ be a convex compact set in $\mathbb{R}^2$. Then for any circle $B(0,r)$,
	\begin{equation}\label{Eq:MS1}
	\lambda(A+B(0,r)) = \lambda(A)+\ell(\partial A)r + \pi r^2,
	\end{equation}
	where $\ell(\partial(A))$ is the (1-dimensional) length of the boundary of $A$.
\end{theorem}

\medskip \noindent By substituting~(\ref{Eq:MS1}) into the proof of Proposition~\ref{Prop:BM}, we obtain
\begin{equation}\label{Eq:MS2}
(1+o(1))n^2 \geq \sum_{i=1}^m \lambda(T''_i) = \sum_{i=1}^m \lambda(T'_i) + \frac{t}{2} \cdot \sum_{i=1}^m \ell(\partial(T'_i)) + m \pi t^2/4.
\end{equation}
To proceed, we may bound the length of the boundary of each $T'_i$ in terms of the area $\lambda(T'_i)$. If we do not make further assumptions on the $T'_i$'s, then the best possible bound of this type is the \emph{isoperimetric inequality}, which asserts
\[
\ell(\partial(A)) \geq \sqrt{4\pi}\sqrt{\lambda(A)},
\]
for any region $A$ in the plane bounded by a closed curve. Plugging this into~(\ref{Eq:MS2}) yields exactly the same upper bound as Proposition~\ref{Prop:BM} (which comes by no surprize, as the tightness case of the isoperimetric inequality is circles, just like the tightness case in the proof of Proposition~\ref{Prop:BM}).

However, if we further assume that each tile is a convex polygon, then we can obtain an improved bound, as function of the number of vertices in each such polygon. We use another classical result, going back to an Ancient Greec mathematician:
\begin{theorem}[Zenodorus]
	Among all polygons on $n$ vertices with the same area, the perimeter is minimized for the regular $n$-gon.
\end{theorem}
Recall that the perimeter of a regular $l$-gon $P_l$ of area $\lambda(P_l)$ is
\[
\ell(\partial(P_l)) = \frac{2\sqrt{l}}{\sqrt{\cot(\pi/l)}} \cdot \sqrt{\lambda(P_l)}.
\]
Plugging this into~(\ref{Eq:MS2}), we obtain the following.
\begin{proposition}\label{Prop:MS}
	Let $T_1,T_2,\ldots$ be a tiling of the plane in $k$ colors, with tiles of area $\leq 1$ and minimal distance $t$. Assume than all tiles are convex polygons with at most $l$ vertices. Then
	\[
	t \leq (1+o(1)) \cdot \frac{-\alpha_l+\sqrt{\alpha_l^2-4(1-k)\pi}}{\pi},
	\qquad \mbox{where} \qquad \alpha_l = \frac{2\sqrt{l}}{\sqrt{\cot(\pi/l)}}.
	\]
\end{proposition}

\begin{proof}
	First, we follow the proof sketch presented above until Equation~(\ref{Eq:MS2}), which asserts
	\[
	(1+o(1))n^2 \geq \sum_{i=1}^m \lambda(T''_i) = \sum_{i=1}^m \lambda(T'_i) + \frac{t}{2} \cdot \sum_{i=1}^m \ell(\partial(T'_i)) + m \pi t^2/4.
	\]
	Using the theorem of Zenodorus, we obtain
	\[
	(1+o(1))n^2 \geq \sum_{i=1}^m \lambda(T'_i) + \frac{t}{2} \cdot \alpha_l \cdot \sum_{i=1}^m \sqrt{\lambda(T'_i)} + m \pi t^2/4,
	\]
	where $\alpha_l=\frac{2\sqrt{l}}{\sqrt{\cot(\pi/l)}}$. As $m\geq n^2/k$ and $\sum_i \sqrt{\lambda(T'_i)} \geq \sum_i \lambda(T'_i) \geq n^2/k$, this implies
	\[
	(1+o(1))n^2 \geq \frac{n^2}{k} \cdot \left(1 + \frac{t}{2} \cdot \alpha_l+\pi t^2/4\right),
	\]
	and consequently, $(1+o(1))k \geq 1+\frac{\alpha_l}{2} t + \frac{\pi}{4}t^2$. Solving the quadratic inequality, we obtain
	\[
	t \leq (1+o(1))\cdot \frac{-\alpha_l+\sqrt{\alpha_l^2-4(1-k)\pi}}{\pi},
	\]
	as asserted.
\end{proof}

\paragraph{Discussion.} For small values of $l$, the upper bound obtained in Proposition~\ref{Prop:MS} is rather strong. For example, for $l=3$ we obtain $t \leq 0.707$ for 3 colors and $t \leq 0.985$ for four colors. As the constructions we obtain below, based on rectangular tiles, satisfy $t=1/\sqrt{2}$ for 3 colors and $t=1$ for four colors, Proposition~\ref{Prop:MS} shows that these constructions cannot be improved using triangular tiles.

As $l$ increases, the bound of Proposition~\ref{Prop:MS} becomes weaker and approaches the bound of Proposition~\ref{Prop:BM}, since circles (for which Proposition~\ref{Prop:BM} is tight) can be approximated to any precision by convex polygons with a sufficiently large number of vertices.

\subsection{Lower bounds on the fault tolerance}
\label{sec:sub:lower}

In this subsection we present two series of explicit tilings. The first series provides the best lower bounds on $t$ for a small number of colors we aware of, and the second series shows the tightness of the asymptotic upper bound $t \leq (1.074+o(1))\sqrt{k}$.

\subsubsection{The brick wall construction}
\label{ssec:bricks}

In the brick wall construction with $k$ colors, demonstrated in Figure~\ref{fig:brick_wall}, each tile is a rectangle with side lengths $\sqrt{(k-2)/2}$ and $\sqrt{2/(k-2)}$ (and so, the area of each tile is $1$). The construction is periodic, where the basic unit is two rows of adjacent rectangles, colored in a round robin fashion. For an even $k$, the second row is placed exactly below the first row, and the sequence of colors is shifted by $k/2$. For an odd $k$, the second row is indented by half a brick (making the construction look like a brick wall), and the sequence of colors is shifted by $(k+1)/2$.

\begin{figure}
	\centering
	\includegraphics[width=1.0\columnwidth]{Rectangles.pdf}				\caption{The brick wall construction for $k=3,4,5$ colors. The ratio between the width and the height of each tile is $2:k-2$.}
	\label{fig:brick_wall}
\end{figure}

It is easy to see that the minimal distance between two same-colored points in different tiles is $\sqrt{(k-2)/2}$. (This distance is attained both in the vertical and in the horizontal directions. Having the same minimal distance in both directions is the optimization that dictates the side lengths of the bricks.)

In particular, we obtain the lower bounds $t \geq 1/\sqrt{2}$ for 3 colors, $t \geq 1$ for 4 colors, and $t \geq \sqrt{3/2}$ for 5 colors.

\subsubsection{The honeycomb of rectangles construction}
\label{ssec:hcr}

In the honeycomb of rectangles construction with $k=m^2$ colors, each tile is a rectangle with side lengths $a$ and $1/a$, where
\begin{equation}\label{Eq:Honey_rect1}
a=\left(\frac{m^2-2m+1}{\frac{3}{4}m^2-m} \right)^{1/4}.
\end{equation}
The construction is periodic, where the basic unit is composed as follows. First, we construct a basic `large rectangle', which is an $m$-by-$m$ square block of tiles, using all the $k=m^2$ colors (in arbitrary order). Then, the basic unit of the construction is two `fat rows' of adjacent large rectangles, where the second row is indented by half a large rectangle. The coloring of each large rectangle is the same. The construction, for $k=16$, is demonstrated in Figure~\ref{fig:honey_rect}.

\begin{figure}
	\centering
	\includegraphics[width=.45\columnwidth]{Grid.pdf}				\caption{The honeycomb of rectangles construction for $k=16$ colors. The boundaries of the basic `large rectangles' are depicted in bold. The placement of the tiles in each single color corresponds to the honeycomb lattice.}
	\label{fig:honey_rect}
\end{figure}

Note that the set of tiles colored in some single color forms the shape of a honeycomb lattice (including the centers of the hexagons). This is why we call the construction `honeycomb of rectangles'.

It is easy to see that the minimal horizontal and diagonal distances between two same-colored tiles are
\[
(m-1)a \qquad \mbox{and} \qquad \sqrt{\left(\frac{m-1}{a}\right)^2+\left(\left(\frac{m}{2}-1\right)a\right)^2},
\]
respectively. By choosing $a$ such that the two distances are equal, we obtain~(\ref{Eq:Honey_rect1}), and the asymptotic lower bound
\[
t\geq (1+o(1)) \left(\frac{4}{3}\right)^{1/4} \cdot m \approx 1.074 \sqrt{k},
\]
that matches the upper bound obtained above.

\section{Fault-Tolerant Rounding in Higher Dimensions}
\label{sec:3dim}

In this section we consider tilings of $\mathbb{R}^d$, for $d \geq 3$, by tiles $T_1,T_2,\ldots$ of volume at most $1$. As was shown in Section~\ref{sec:qualititave}, the minimal number of colors that allows for any fault tolerance is $d+1$.

First, we obtain two upper bounds on the fault tolerance, by generalizing the bounds via the Brunn-Minkowski inequality and via the circle packing problem presented in Section~\ref{sec:sub:upper}. In the other direction, we present an explicit tiling with $d+1$ colors that attains a minimal distance of $t=\frac{1}{1+2(d-1)\sqrt{2}}$ (and hence, fault tolerance of $t/2$), a specific 4-color tiling of $\mathbb{R}^3$ that attains a larger minimal distance, and a $k$-color tiling of $\mathbb{R}^3$ that attains the optimal asymptotic fault tolerance.

\subsection{Upper bounds on the fault tolerance in $\mathbb{R}^d$}
\label{sec:sub:upper-high}

We present two upper bounds. The first -- via the Brunn-Minkowski inequality -- is most effective for a small number of colors. The second -- via sphere packing -- is most effective in dimensions $3,8,$ and $24$, for which the sphere packing problem was solved completely.

\subsubsection{An upper bound using the Brunn-Minkowski inequality}
\label{ssec:BM-high}

Our bound is a straightforward generalization of Proposition~\ref{Prop:BM}. We start with a large cube $S$ of volume $n^d$, find a color whose intersection with the cube has volume at least $n^d/k$, consider the tiles $T'_1,\ldots,T'_m$ in that color, and use the fact that by assumption, inflations of these tiles by Minkowski sum with the ball $B(0,t/2)$ are pairwise disjoint.
\begin{proposition}\label{Prop:BM-higher}
	Let $T_1,T_2,\ldots$ be a tiling of $\mathbb{R}^d$ in $k$ colors, with tiles of volume $\leq 1$ and minimal distance $t$. Then
	\[
	t \leq \left(\frac{2\Gamma(d/2+1)^{1/d}}{\sqrt{\pi}} +o(1) \right) \cdot (k^{1/d}-1),
	\]
	where $\Gamma(\cdot)$ is the Gamma function.
\end{proposition}

\begin{proof}
	The proof is almost identical to the proof of Proposition~\ref{Prop:BM}. Instead of~(\ref{Eq:BM3}), we obtain
	\[
	(1+o(1))n^d \geq \sum_{i=1}^m \sum_{j=0}^d {{d}\choose{j}} \lambda(T'_i)^{j/d} (b_{t/2})^{1-j/d},
	\]
	where $b_{t/2}$ is the volume of the $d$-dimensional ball $B(0,t/2)$. As $0 \leq \lambda(T'_i) \leq 1$, we have $\sum_i \lambda(T'_i)^{j/d} \geq \sum_i \lambda(T'_i) \geq n^d/k$, and hence we obtain
	\[
	(1+o(1))n^d \geq \frac{n^d}{k} \cdot \left(1+b_{t/2}^{1/d} \right)^d.
	\]
	This implies
	\[
	(1+o(1))(k^{1/d}-1) \geq b_{t/2}^{1/d} = \frac{\pi^{1/2}}{\Gamma(d/2+1)^{1/d}} \cdot \frac{t}{2},
	\]
	and consequently,
	\[
	t \leq \left(\frac{2\Gamma(d/2+1)^{1/d}}{\sqrt{\pi}} +o(1) \right) \cdot (k^{1/d}-1),
	\]
	as asserted.
\end{proof}

\paragraph{Discussion.} For a large number $k \gg d$ of colors, Proposition~\ref{Prop:BM-higher} gives the upper bound $t \leq (\sqrt{2/\pi e}+o(1))k^{1/d}$ (see~(\ref{Eq:Aux5.1}) below). This bound is not far from being tight. Indeed, its dependence on $k$ is correct, as it can be easily matched by a periodic cubic tiling, in which each tile is a cube with side length 1 and the basic unit is a large cube with side length $k^{1/d}$ that contains each color in exactly one tile (in the same order). Moreover, even regarding the coefficient of $k^{1/d}$, the optimal asymptotic upper bounds for $d=3,8,24$ which we obtain below via the sphere packing problem, improve over this bound by only a small factor.

\medskip \noindent To estimate the upper bound on $t$ for a large $d$ and $k=d+1$ colors, note that
\begin{equation}\label{Eq:Aux5.1}
(d+1)^{1/d}-1 = (1+o(1))\frac{\ln(d)}{d}, \qquad \mbox{and} \qquad \Gamma(d/2+1)^{1/d} = \left(\frac{1}{\sqrt{2e}}+o(1)\right)\sqrt{d}.
\end{equation}
Therefore, the bound we obtain in this case is
\[
t \leq \left(\sqrt{\frac{2}{\pi e}}+o(1) \right)\frac{\ln d}{\sqrt{d}},
\]
which implies that the fault tolerance decreases to zero as $d$ tends to infinity. For comparison, the lower bound we obtain in Section~\ref{ssec:dim-reduce} is $t \geq \Omega(1/d)$.

\medskip \noindent For the `smallest' case $d=3,k=4$, in which the bound of Proposition~\ref{Prop:BM-higher} is much stronger than the bound we obtain via sphere packing, the bound is
\[
t \leq \frac{2\Gamma(2.5)^{1/3}}{\sqrt{\pi}} \cdot (4^{1/3}-1) \approx 0.729.
\]
For comparison, the best construction we have in this setting satisfies $t=0.350$.

\paragraph{The possibility of obtaining an upper bound via the Minkowski-Steiner formula.} In Section~\ref{ssec:Steiner} we obtained an upper bound on $t$, under the additional assumption that the tiles are convex polygons with a bounded number of vertices, using the Minkowski-Steiner formula. This formula has a higher-dimensional analogue:
\begin{equation}\label{Eq:MS-higher}
\lambda(A+B(0,r)) = \lambda(A)+\lambda_{d-1}(\partial A)r + \sum_{j=2}^{d-1} \lambda_j(A) r^j + \frac{2\pi^{d/2}}{d\Gamma(d/2)} r^d,
\end{equation}
where $\lambda_{d-1}(\partial(A))$ is the $(d-1)$-dimensional volume of the boundary of $A$, and $\lambda_j(A)$ are continuous functions of $A$ called \emph{mixed volumes}.

In order to obtain an upper bound on the fault tolerance using~(\ref{Eq:MS-higher}), one has to bound the mixed volumes $\lambda_j(A)$ from below in terms of $\lambda(A)$, which seems to be a challenging task.

\subsubsection{An upper bound using the Sphere Packing problem}
\label{ssec:SP-high}

In order to generalize the argument of Proposition~\ref{Prop:CP} to $\mathbb{R}^d$, we need a generalization to $d>2$ of Theorem~\ref{Thm:FT}, namely, an answer to the following well-known question:
\begin{question}\label{Q:Packing}
	What is the maximal number of points that can be placed in a $d$-dimensional cube of volume $v$, such that the minimal distance between two points is 1?
\end{question}
However, as far as we know, this question is open for all $d>2$.

\medskip \noindent Instead, we obtain an upper bound via a reduction to the classical \emph{sphere packing} problem, which asks for the maximal possible \emph{density} of a set of non-intersecting congruent spheres in $\mathbb{R}^d$.
\begin{definition}
	The density of a sphere packing (i.e., collection of pairwise disjoint congruent spheres) $P= \cup P_i$ in $\mathbb{R}^d$ is
	\[
	\limsup_{r \rightarrow \infty}  \frac{\lambda \left(B(0,r) \cap \bigcup P \right)}{\lambda(B(0,r))}.
	\]
\end{definition}
Intuitively, this measures the fraction of the volume of a large ball covered by the spheres of the packing.
\begin{notation}
	Denote the maximal density of a sphere packing in $\mathbb{R}^d$ by $\delta_d$, and the volume of the unit ball $B(0,1) \subset \mathbb{R}^d$ by $v_d = \pi^{d/2}/\Gamma(d/2+1)$.
\end{notation}

\begin{proposition}\label{Prop:SP-higher}
Let $T_1,T_2,\ldots$ be a tiling of $\mathbb{R}^d$ in $k$ colors, with tiles of volume $\leq 1$ and minimal distance $t$. Then
\[
t \leq \left(2(\delta_d/v_d)^{1/d}+o(1) \right) \cdot k^{1/d} = \left(\frac{2\Gamma(d/2+1)^{1/d} \cdot \delta_d^{1/d}}{\sqrt{\pi}}+o(1) \right) \cdot k^{1/d},
\]
where $\Gamma(\cdot)$ is the Gamma function.
\end{proposition}
Note that the asymptotic upper bound of Proposition~\ref{Prop:SP-higher} is stronger than the asymptotic upper bound that follows from Proposition~\ref{Prop:BM-higher} by the constant factor $(\delta_d)^{1/d}$.

\medskip

\begin{proof}
	Let $T$ be a $k$-colored tiling of $\mathbb{R}^d$ that satisfies the assumptions of the proposition, and consider the sequence of balls  $\{B(0,n)\}_{n=1,2,3,\ldots}$. By the pigeonhole principle, there exists a color (say, black) such that for each $n_{\ell}$ in an infinite subsequence $\{n_\ell\}_{\ell=1,2,\ldots}$, the intersection of the black tiles with the ball $B(0,n_\ell)$ has volume of at least
	\[
	\lambda(B(0,n_\ell))/k = \frac{n_\ell^d \cdot v_d}{ k}.
	\]
	As the volume of each tile is at most $1$, we know that for each $n_\ell$, the number of black tiles that intersect $B(0,n_\ell)$ is at least $n_\ell^d v_d / k$.

	\medskip \noindent Pick some value $n_{\ell}$, denote the intersections of black tiles with $B(0,n_{\ell})$ by $T'_1,T'_2,\ldots$, and take one point $x_i$ from each tile $T'_i$. As the minimal distance between two black points in different tiles is $t$, balls of radius $t/2$ around the points $x_i$ are pairwise disjoint. Hence, their total volume is at least
	\[
	\frac{n_\ell^d \cdot v_d}{k} \cdot (t/2)^d \cdot v_d.
	\]
	On the other hand, each such ball is contained in the ball $B(0,(1+o(1))n_{\ell})$ (since its volume is at most $1$, and it contains a point in $B(0,n_\ell)$). This implies that for a sufficiently large $\ell$, the total volume of these balls must be smaller than $\delta_d \cdot \lambda(B(0,(1+o(1))n_{\ell}))$, as otherwise, the infinite collection of the balls $B(x_i,t/2)$ (where for each ball $B(0,n_\ell)$ we select $x_i$'s in the way described above, respecting the $x_i$'s selected for smaller values of $n_\ell$) would be a sphere packing of $\mathbb{R}^d$ whose density is larger than $\delta_d$.	Therefore, for a sufficiently large $n_\ell$, we have
	\[
	\frac{n_\ell^d \cdot v_d}{k} \cdot (t/2)^d \cdot v_d \leq (1+o(1))\delta_d \cdot n_\ell^d v_d,
	\]
	or equivalently,
	\[
	t \leq (2+o(1)) (\delta_d/v_d)^{1/d} \cdot k^{1/d},
	\]
	as asserted.
\end{proof}

\paragraph{Discussion.} In the two last decades, there has been a tremendous progress in the research of the sphere packing problem. In 2005, Hales~\cite{Hales05} solved the problem for $d=3$, proving a 17'th century conjecture of Kepler. Three years ago, in a beautiful short paper, Viazovska~\cite{V17} solved the problem for $d=8$, and shortly after, Cohn, Kumar, Miller, Radchenko, and Viazovska~\cite{CKMRV17} used Viazovska's method along with other tools to solve the problem for $d=24$. For other dimensions, the problem is still open. We can use the results of~\cite{CKMRV17,Hales05,V17} to obtain tight upper bounds on $t$ in dimensions $3,8,$ and $24$.
\begin{itemize}
	\item For $d=3$, Hales~\cite{Hales05} showed that $\delta_3=\frac{\pi}{3\sqrt{2}} \approx 0.740$. Hence, we obtain the bound $t\leq (2^{1/6}+o(1))k^{1/3} \approx (1.122+o(1))k^{1/3}$.

	\item For $d=8$, Viazovska~\cite{V17} showed that $\delta_8=\frac{\pi^4}{2^4 4!} \approx 0.254$. Hence, we obtain the bound $t\leq (\sqrt{2}+o(1))k^{1/8} \approx (1.414+o(1))k^{1/8}$.

	\item For $d=24$, Cohn et al.~\cite{CKMRV17} showed that $\delta_{24}=\frac{\pi^{12}}{12!} \approx 0.0019$. Hence, we obtain the bound $t\leq (2+o(1))k^{1/24}$.
\end{itemize}
As we show in Section~\ref{sec:sub:lower-high}, these bounds are asymptotically tight.

Using the same method, we can transform any upper bound for the sphere packing problem to an upper bound on the fault tolerance of a rounding scheme in the corresponding dimension. A list of such conjectured bounds for $d \leq 10$ can be found in~\cite{CS95}.

\subsection{Lower bounds on the fault tolerance in $\mathbb{R}^d$}
\label{sec:sub:lower-high}

In this section we present three explicit constructions that yield lower bounds on $t$. The first construction considers $d+1$ colors in $\mathbb{R}^d$. The second construction is specific to 4-colored tiling of $\mathbb{R}^3$ but yields a larger lower bound. The third construction shows that the asymptotic upper bound $t \leq (2^{1/6}+o(1))k^{1/3}$ in $\mathbb{R}^3$ is tight.

\subsubsection{The dimension reducing construction}
\label{ssec:dim-reduce}

We exemplify the construction in $\mathbb{R}^3$, but it will be apparent how to generalize it to higher dimensions.

\paragraph{Informal description of the construction.} Informally, the construction works as follows.

\medskip \noindent \textbf{Step~1:} We begin with dividing $\mathbb{R}^3$ into cubes with side length $a$ (to be determined below), depicted in Figure~\ref{fig:dim_reduce}(a). Then, we make the `interior' of each cube into a tile and give all these tiles the color~$1$. In order to keep a minimal distance of $t$ between two points colored~1 in different tiles, we must leave a neighborhood of width $t/2$ in each side of each facet of the basic cube. Hence, we are left with `fattened' walls of total width $t$. The part of such a wall included in a basic cube is shown in Figure~\ref{fig:dim_reduce}(b).

\medskip \noindent \textbf{Step~2:} We make the `interior' of each wall (i.e., fattened facet) into a tile and give all these tiles the color~$2$, as is demonstrated in Figure~\ref{fig:dim_reduce}(b). (Note that each tile contains points from two adjacent basic cubes.) In order to keep a minimal distance of $t$ between two points colored~2 in different tiles, we must leave a neighborhood of $t/\sqrt{2}$ near each edge (i.e., intersection of facets), as is shown in Figure~\ref{fig:dim_reduce}(c). Hence, we are left with a `fattened skeleton'.

\medskip \noindent \textbf{Step~3:} We make the `interior' of each edge of the skeleton into a tile and give all these tiles the color~$3$, as is shown in Figure~\ref{fig:dim_reduce}(d). (Note that each tile contains points from four adjacent basic cubes.) In order to keep a minimal distance of $t$ between two points colored~3 in different tiles, we must leave an additional neighborhood of $t/\sqrt{2}$ near each vertex (i.e., intersection of edges), as is demonstrated in Figure~\ref{fig:dim_reduce}(d). Hence, we are left with neighborhoods of the corners (a.k.a. vertices).

\medskip \noindent \textbf{Step~4:} We make the neighborhoods of the corners into tiles and give them the color~$4$. (Note that each tile contains points from 8 adjacent basic cubes.) We have to make sure that the distance between each two such `fattened corners' is at least $t$, and this requirement dictates the choice of $t$, as is explained below.

\begin{figure}
	\centering
	\includegraphics[width=0.9\columnwidth]{3D-Mess}
	\caption{The dimension reducing construction. Part~(a) shows the basic cubes we start with. Part~(b) shows a fattened wall of width $t/2$, and its interior that gets the color~2. Part~(c) shows where the minimal distance between two tiles colored~2 is attained. Part~(d) presents (in full lines) the interiors of the fattened edges that get the color~3 and shows where the minimal distance between two such tiles is attained.}
	\label{fig:dim_reduce}
\end{figure}

It is clear that the resulting 4-colored tiling can be generalized into a $(d+1)$-colored tiling of $\mathbb{R}^d$.

\paragraph{Formal definition of the construction.} For the sake of formality, we give the exact definitions of the tiles below.

\medskip \noindent For $(x,y,z) \in \mathbb{R}^3$, let $f(x,y,z) \in [0,a/2]^3$ be a motonone non-decreasing ordering of $\{\bar{x},\bar{y},\bar{z}\}$, where $\bar{b}=\min\{|b-an|:n \in \mathbb{Z}\}$. That is, we measure the minimal distance to a wall in each coordinate and arrange these minimal distances in a non-decreasing order. For example, for any $a \geq 1$, the point $(2.4a,7.8a,3.3a)$ is mapped by $f$ to $(0.2a,0.3a,0.4a)$.
\begin{itemize}
	\item The tiles colored~1 consist of all points $(x,y,z)$ in a basic cube such that $f(x,y,z)_1 > t/2$. (These are exactly the points that are at least $t/2$-far from each wall -- the `interior' of the cube.)

	\item The tiles colored~2 consist of points $(x,y,z)$ in a basic cube such that $f(x,y,z)_1 \leq t/2$ and $f(x,y,z)_2 > t/2+t/\sqrt{2}$. (These are the points that are close to a wall only in a single coordinate -- the `interiors' of the fattened facets.) Note that each basic cube contains six such tiles, and each such tile contains points of two adjacent basic cubes.

	\item The tiles colored~3 consist of points $(x,y,z)$ in a basic cube such that $f(x,y,z)_1 \leq t/2$, $f(x,y,z)_2 \leq t/2+t/\sqrt{2}$, and $f(x,y,z)_3 > t/2+2t/\sqrt{2}$. (These are the points that are close to a wall in two coordinates -- the `interiors' of the fattened edges.) Note that each basic cube contains 12 such tiles (one for each edge), and each such tile contains points of four adjacent basic cubes.

	\item The tiles colored~4 consist of the rest of the points (i.e., the fattened neighborhoods of the vertices). Note that each basic cube contains 8 such tiles (one for each vertex), and each such tile contains points of eight adjacent basic cubes.
\end{itemize}

\paragraph{The choice of $t$.} The construction makes sure that for $i=1,2,3$, the distance between two points colored~$i$ in different tiles is at least $t$. In order to guarantee the same condition for the color~4 as well, we have to choose $t$ such that the distance between two `neighborhoods of corners' will be at least $t$. By the construction, this amounts to the inequality
\[
a-2 \left(\frac{t}{2}+\frac{2t}{\sqrt{2}} \right) \geq t,
\]
or equivalently, $t \leq \frac{a}{2+2\sqrt{2}}$. Similarly, in $\mathbb{R}^d$ we obtain $t \leq \frac{a}{2+2(d-1)\sqrt{2}}$.

\paragraph{The choice of $a$.} An easy computation shows that the largest tiles are those colored~1. Hence, in order to make the volume of all tiles $\leq 1$, we have to choose $a$ such that the volume of each such tile is~1. By construction, this amounts to the equality
\[
1 = (a-t)^3= a^3 \cdot \left(1- \frac{1}{2+2\sqrt{2}}\right)^3,
\]
or equivalently, $a=1+\frac{1}{1+2\sqrt{2}} \approx 1.261$. Similarly, in $\mathbb{R}^d$ we obtain $a=1+\frac{1}{1+2(d-1)\sqrt{2}}$.

\paragraph{The lower bound on $t$.} Summarizing the above, the value of $t$ we obtain in $\mathbb{R}^3$ is
\[
t = \frac{a}{2+2\sqrt{2}} = \frac{1}{1+2\sqrt{2}} \approx 0.261.
\]
In $\mathbb{R}^d$, we obtain $t=\frac{1}{1+2(d-1)\sqrt{2}}$. This gives lower bounds of $0.131$ and $\frac{1}{2+4(d-1)\sqrt{2}}$ on the fault tolerance in $\mathbb{R}^3$ and $\mathbb{R}^d$, respectively. The lower bound in $\mathbb{R}^3$ is superseded by the `bricks and balloons' construction described below. The lower bound in $\mathbb{R}^d$, which is $\Omega(1/d)$, is not very far from the upper bound we obtained in Section~\ref{ssec:BM-high} -- namely, $O(\log d / \sqrt{d})$.

\subsubsection{The bricks and balloons construction}
\label{ssec:bb}

\begin{figure}
	\centering
	\includegraphics[width=1.0\columnwidth]{2D-4Col.pdf}				\caption{The bricks and balloons construction}
	\label{fig:b&b}
\end{figure}

The bricks and balloons (B\&B) construction, demonstrated in Figure~\ref{fig:b&b}, is a periodic tiling of $\mathbb{R}^3$, colored in 4 colors. In order to present the tiling, we need an auxiliary notation.

\paragraph{Notation.} Consider the slab $D=\mathbb{R}\times \mathbb{R} \times [z_1,z_2] \subset \mathbb{R}^3$. We say that a tiling $T_1,T_2,\ldots$ of $D$ is a \emph{fattened plane tiling} if there exists a tiling $T'_1,T'_2,\ldots$ of the plane such that $\forall i: T_i = T'_i \times [z_1,z_2]$.

\paragraph{The layers of B\&B.} The B\&B construction is based on two fattened plane tilings:
\begin{itemize}
	\item The \emph{brick wall layer} -- a fattening of a brick wall tiling of the plane. The underlying plane tiling is a periodic tiling, in which the basic unit consists of four rows of adjacent rectangles, where the even rows are indented by half a brick, making the construction look like a brick wall. In the odd rows, the colors 1,2 are used alternately, and in the even rows, the colors 3,4 are used alternately. Furthermore, the colors in the third and fourth rows are shifted by one, so that a rectangle colored $i$ is not placed right under another rectangle colored $i$, see Figure~\ref{fig:b&b}. The side lengths of the rectangles and the width of the slab will be specified below.

	\item The \emph{balloon layer} -- a fattening of another periodic tiling of the plane. This tiling is similar in its high-level structure to the brick wall layer, but instead of rectangles, the tiles are `balloons' consisting of a regular octagon and a square built on one of its edges. The basic unit of the tiling is four columns of adjacent balloons, where the even columns are shifted in such a way that the balloons tile the plane.  In the odd columns, the colors 1,2 are used alternately, and in the even columns, the colors 3,4 are used alternately. Furthermore, the colors in the third and fourth columns are shifted by one, so that a balloon colored $i$ is not placed near another balloon colored $i$ (see Figure~\ref{fig:b&b}). The side lengths of the balloons and the width of the slab will be specified below.
\end{itemize}

\paragraph{The structure of B\&B.} The B\&B tiling is periodic, where the basic unit consists of a brick wall layer and a balloon layer, placed one on top of the other in the way presented in Figure~\ref{fig:b&b}. (Note that once the two layers are placed, the coloring of one layer fully determines the coloring of the other.)

\paragraph{The side lengths and the slab heights.} Denote the side length of the balloon by $2a$. It is easy to see that the side lengths of the rectangles are $(2+\sqrt{2})a$ and $(4+2\sqrt{2})a$. Furthermore, the area of both a balloon and a rectangle is equal to $(12+8\sqrt{2})a^2$. Hence, in order to make all tiles equal-volume, we choose the height of all slabs (both in the brick wall layers and in the balloon layers) to be the same value $z$.

\medskip To find the minimal distance between two same-colored tiles, we consider several cases:
\begin{itemize}
	\item Two same-colored balloons in the same layer: The minimal distance between two such balloons is $\sqrt{8+4\sqrt{2}}a \approx 3.696a$.

	\item Two same-colored bricks in the same layer: Here, the minimal distance is $(2+\sqrt{2})a$.

	\item A balloon and a same-colored brick at an adjacent layer: Here, the minimal distance is $a$.

	\item A balloon/brick and a same colored balloon/brick two layers apart: Here, the minimal distance is $z$.
\end{itemize}
Hence, the minimal distance between two same-colored tiles is $\min\{a,z\}$, and in order to optimize it we choose $z=a$.

\paragraph{The lower bound on $t$.} Since $z=a$, the volume of each tile is $(12+8\sqrt{2})a^2 \cdot a = (12+8\sqrt{2})a^3$. Thus, in order to make the volume of all tiles equal to $1$, we fix
\[
a = \left(\frac{1}{12+8\sqrt{2}} \right)^{1/3} \approx 0.350.
\]
Therefore, this construction satisfies $t = a \approx 0.35$, and hence, provides a lower bound of $t/2 \approx 0.175$ on the fault tolerance. This improves significantly over the fault tolerance of the dimension reducing construction (i.e.,~$0.131$), but is still far from the upper bound obtained in Section~\ref{ssec:BM-high} (namely, $0.365$).

\subsubsection{The close packing of boxes construction}
\label{ssec:CPB}

\paragraph{Motivation.} This construction, a $k$-colored tiling of $\mathbb{R}^3$ where $k \gg 3$, is a natural generalization to $\mathbb{R}^3$ of the `honeycomb of rectangles' construction presented in Section~\ref{ssec:hcr}. The idea behind the construction is to choose an optimal sphere packing in $\mathbb{R}^3$, and construct a fattened plane tiling in which the tiles in each color are placed at the centers of the spheres of the packing. (Recall that as the number of colors is large, the size of each tile is negligible with respect to the size of its inflation, and hence, we can treat the tiles as single points.)

We use the classical HCP lattice (one of the most common \emph{closed packings}, see~\cite{CS13}), which corresponds to a periodic sphere packing in which the basic unit is two hexagonal layers of spheres, where in the top layer, each sphere is placed on top in the hollow between three spheres in the bottom layer. The coordinates of the centers of these spheres are:
\[
(r,r,r), (3r,r,r), (5r,r,r),\ldots,(2r,r+\sqrt{3}r,r),(4r,r+\sqrt{3}r,r), (6r,r+\sqrt{3}r,r), \ldots
\]
for the bottom layer, and
\[
(2r,r+\frac{\sqrt{3}}{3}r,r+\frac{2\sqrt{6}}{3}r), (4r,r+\frac{\sqrt{3}}{3}r,r+\frac{2\sqrt{6}}{3}r), \ldots,
(r,r+\frac{4\sqrt{3}}{3}r,r+\frac{2\sqrt{6}}{3}r),
(3r,r+\frac{4\sqrt{3}}{3}r,r+\frac{2\sqrt{6}}{3}r), \ldots
\]
for the top layer.

\paragraph{The structure of the tiling.} Assume that the number of colors is $k=m^3$. Each tile is a box with side lengths $(a,b,c)$ to be determined below, and the basic unit is a `large box', which is an $m\times m \times m$ cubic block of tiles, using all the $k=m^3$ colors (in arbitrary order). Then, the basic unit of the construction is a two-layer fattened plane tiling, in which each layer is a fattened copy of the `honeycomb of rectangles' tiling. The upper layer is shifted by $\frac{m}{2}a$ in the $x$-coordinate and by $\frac{\sqrt{3}m}{6}a$ in the $y$-coordinate, so that the corners of the large boxes lie in the coordinates of the sphere centers described above (for $r=\frac{m}{2}a$). A quick calculation shows that in order to make this possible, the proportion $(a:b:c)$ should be \emph{approximately} $(2:\sqrt{3}:2\sqrt{6}/3)$ (where we neglect the size of each tile with respect to the size of the inflation, which can be absorbed in an $1+o(1)$ multiplicative factor in the final value of $t$). The volume of each tile is clearly $abc$. In order to make the volumes of all tiles equal to $1$, we need
\[
a \cdot \frac{\sqrt{3}}{2}a \cdot \frac{\sqrt{6}}{3}a = 1,
\]
and thus, $a=(6/\sqrt{18})^{1/3} = 2^{1/6} \approx 1.122$. Hence, the side lengths of each tile are
\[
(2^{1/6}, 3^{1/2}/2^{5/6}, 2^{2/3}/3^{1/2}) \approx (1.122, 0.972, 0.916),
\]
and the minimal distance between two same-colored points in different tiles is
\[
(m-1)a = (2^{1/6}+o(1))k^{1/d},
\]
which matches the upper bound proved in Section~\ref{ssec:BM-high}.

\paragraph{Generalization to higher dimensions.} A similar tiling can be constructed to match any \emph{lattice sphere packing}, assuming the number of colors $k$ is sufficiently large with respect to $d$. Hence, any dense lattice sphere packing can be translated into a lower bound on the asymptotic fault tolerance of rounding schemes in the corresponding dimension. In particular, as the \emph{E8 lattice} and the \emph{Leech lattice} which attain the maximal possible density of sphere packings in dimension 8 and 24 (respectively) are lattice packings, they can be translated to box tilings showing that the asymptotic upper bounds on the fault tolerance in $\mathbb{R}^8$ and $\mathbb{R}^{24}$ proved in Section~\ref{ssec:BM-high} are tight.

\section{Summary}
\label{sec:summary}

In this paper we demonstrated that the natural problem of how to round points in $\mathbb{R}^d$ with optimal fault tolerance using a tiny amount of side information is tied to deep problems on space partitions and sphere packings.  We proposed a new way to look at the problem as colored tilings of the space, and related it to a generalized type of error correcting code in which we inflate tiles into Minkowski sums rather than points into balls. We obtained a number of tight asymptotic bounds for a large number of colors, as well as lower and upper bounds for small numbers of colors, but finding the optimal solutions for particular numbers of colors and dimensions is still wide open.

\section*{Acknowledgements}

The authors are grateful to Stephen D.~Miller for valuable suggestions
regarding the sphere packing problem.