Skip to content Mastodon

Notes and Tones

Oranges in Many Dimensions

Peeling oranges in high dimensionality k ≫ 3 is hard. This post discusses a fundamental phenomenon in asymptotic geometry and explains at which rate you should peel an orange.

Suppose we have a special branch of hyper-orange existing in d ≥ 1 and we would like to peel such that 95% of pulp is preserved. We will denote the volume of an Lp norm in d dimensions, that is C = {x ∈ ℝd∣ ‖x‖p ≤ r}, as V ⁡ C = V ⁡ p,d. Because we first deal with balls we have p = 2 and set by default V ⁡ 2,d = V ⁡ d.

We first look at ratios of volume of balls radii scaled by φ such that V ⁡ d(φr)∕V ⁡ d(r) = 0.95. For d ∈{1,2,3} we know from middle school that V ⁡ 1(r) = 2r, V ⁡ 2(r) = πr2 and V ⁡ 3(r) = 4 3πr3, and hence φ1,2,3 = 0.951,2,3. For d > 3 we need some material from real analysis to work out the volume.

Volume of n-ball Exponential Function

One proof works with two observations: first scaling area identity An−1(r) = rn−1An−1(1) and evaluates the integral for rotational invariant function f(x1,… ⁡,xd) = exp ⁡ (−x12⋯ − xd2).

On the one side, we have

∫ ℝdf(x1,…,xd)dx1…dxd = ∫ ℝ exp ⁡ (−x12)dx 1… ⁡ ∫ ℝ exp ⁡ (−xd2)dx d = πd∕2

On the other side, isolating one variable and applying the area identity

∫ ℝdf(x1,…,xd)dx1…dxd = ∫ 0∞∫ Sn−1(r) exp ⁡ (−r2)dAdr = ∫ 0∞A n−1(1)rd−1 exp ⁡ (−r2)dr = An−1(1)∫ 0∞exp ⁡ (−r2)rd−1dr = 1 2An−1(r)Γ(d 2).

Solving for the d − 1 dimensional area we get Ad−1(1) = 2πd∕2 Γ(d∕2) and substituting into the volume formula, we get the desired V ⁡ d(R) = πd∕2Rd d 2 Γ(d∕2).

Consequence for Peeling Oranges

Now consider the volume of two balls V ⁡ d(1), V ⁡ d(1 − 𝜀) and compare their volume

V ⁡ d(1 − 𝜀) V ⁡ d(1) = (1 − 𝜀)d.

Then, for every constant 𝜀 the ratio tends to zero, and we are left with no pulp at all! Realizing that for 𝜀 ∼ d−1, the limit yields exponential function

V ⁡ d(1 − 𝜀) V ⁡ d(1) = (1 −M d )d → exp ⁡ (−M) for d →∞,

we should peel (−log ⁡ 0.95)∕d of the skin to get a fair share of pulp for high dimensions.

Consequence for Approximating Norms

Let’s pack the orange into a box. For simplicity we assume that the box has side length 1 and recognize that an d-cube is described with the maximum norm. The volume of such a box is V ⁡ ∞,d = 1d = 1, and comparing that to the volume of our orange

V ⁡ 2,d V ⁡ ∞,d = πd∕2 d 2Γ(d∕2)1

Express this in log-space and use log-gamma approximation log ⁡ Γ(z) ≈ (z −1 2)log ⁡ z − z + 1 2 log ⁡ (2π)

log ⁡ V ⁡ 2,d V ⁡ ∞,d = d 2log ⁡ π − log ⁡ d 2 − log ⁡ Γ(d 2) ∼ d 2log ⁡ 2πe d −1 2log ⁡ 2πd ∼ d 2log ⁡ 2πe d

which turns negative once d > 2πe and grows linear in the dimensionality. Hence, the fraction diminishes exponentially fast. That means for large d we can pack exponential many oranges into a single box!