Pith. sign in

REVIEW

Minimum enclosing Bregman balls are exactly power-distance MEBs on dual weighted points, so ordinary Frank–Wolfe geometry applies.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Left Bregman MEBs equal power MEBs on dual Laguerre points; Frank-Wolfe power approximation recovers the 2005 Bregman algorithm, and Bregman liftings equal paraboloid liftings.

T0 review reviewed 2026-07-31 challenge →

load-bearing objection Central Prop 1 equating Bregman MEBs to power MEBs is algebraically false as stated; Voronoi/lifting parts mostly survive.

arxiv 2607.24197 v2 pith:PQH2TADW submitted 2026-07-27 cs.IT cs.CGmath.IT

Minimum enclosing Bregman balls made easy

classification cs.IT cs.CGmath.IT
keywords Bregman divergenceminimum enclosing ballpower distanceFrank–Wolfe algorithmFenchel–Young divergenceBregman Voronoi diagramLaguerre diagramconvex duality
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper shows that finding the smallest ball that covers a finite set of parameters under a Bregman divergence is the same problem as finding a minimum enclosing ball under the classical power distance on a carefully weighted dual point set. Because of that identity, the familiar Frank–Wolfe (1+ε)-approximation used for ordinary MEBs immediately yields a matching algorithm for Bregman MEBs; when rewritten in dual gradient coordinates it recovers the earlier Bregman Badøiu–Clarkson procedure. The same rewriting also explains why Bregman Voronoi diagrams are simply power diagrams clipped to the parameter domain, and why the usual paraboloid lifting of computational geometry is interchangeable with the Bregman potential lifting. A sympathetic reader cares because every existing Euclidean MEB tool—coresets, diagrams, exact LP-type solvers—transfers at once to the information-geometric setting without new analysis.

Core claim

The left Bregman MEB circumcenter of a finite parameter set T equals the power-MEB circumcenter of the corresponding Laguerre weighted points whose sites are the dual gradients η_i = ∇F(θ_i) and whose weights are ∥η_i∥² + 2F*(η_i). Consequently the Bregman Badøiu–Clarkson algorithm is identical to Frank–Wolfe on that weighted set, and Bregman potential liftings coincide with ordinary paraboloid liftings of the same weighted points.

What carries the argument

The mixed-parameterized rewrite of a Bregman divergence as a weighted squared Euclidean (Fenchel–Young) distance: Y_F(θ:η') = ½∥θ−η'∥² − ω_F(θ) − ω_{F*}(η'). This single algebraic identity converts every left Bregman ball query into a power-distance query and every Bregman bisector into a radical hyperplane.

Load-bearing premise

The generator must be Legendre type so the gradient map is a global bijection and the dual weights are well-defined on the whole ambient space.

What would settle it

Pick any concrete Legendre generator (e.g., negative entropy) and a small point set whose Bregman MEB center is already known by exact LP-type or by symmetry; construct the dual weighted points and solve the ordinary power MEB; the two centers must coincide to machine precision, otherwise the claimed equivalence fails.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Any core-set or Frank–Wolfe guarantee proved for power MEBs immediately supplies a (1+ε)-approximation and core-set size bound for Bregman MEBs.
  • Bregman MEB circumcenters lie on the farthest Bregman Voronoi diagram, which is just the farthest power diagram of the dual weighted sites.
  • Bregman Voronoi diagrams can be computed with any off-the-shelf power-diagram code after the dual weighting map is applied.
  • The same reduction extends verbatim to duo-Bregman MEBs (different generators per site), again yielding an ordinary power MEB.
  • Chernoff points between exponential-family densities appear as special two-point Bregman MEB centers and are therefore power-MEB centers.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Once Bregman MEBs are power MEBs, hardware or GPU libraries already tuned for ordinary MEBs become drop-in solvers for information-geometric centering tasks such as minimax priors or KL balls.
  • The dual-weight construction suggests a practical numerical test: monitor condition of the dual map near the domain boundary to decide when the Legendre assumption is about to break.
  • Because the reduction is purely algebraic, the same identity should convert other Bregman facility-location problems (1-median, k-center) into weighted Euclidean problems that inherit existing approximation schemes.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Circularity Check

1 steps flagged

Algebraic equivalences from definitions; only minor non-load-bearing self-citation of prior BVD=PD result. No circular derivation of the central MEB claim.

specific steps
  1. self citation load bearing [§4.2, Proposition 3]
    "Proposition 3 ([8]) The left Bregman Voronoi diagram of a set {θ_i} of n parameters amounts to the power diagram of a set {q̂_i} of n corresponding weighted points clipped to the domain Θ: Vor_F({θ_i}) = Vor_σ({q̂_i}) ∩ Θ, where q̂_i = (η_i, ∥η_i∥² + 2(F(θ_i) − ⟨θ_i, η_i⟩)) with η_i = ∇F(θ_i)."

    The BVD≡PD equivalence is imported wholesale from the authors' own 2010 paper [8] rather than re-derived. This is only mildly circular and not load-bearing for the paper's main MEB claim (Prop. 1), which is attempted from the mixed-parameterization rewrite independently of [8]. Standard self-citation of prior geometry; does not force the MEB center identity.

full rationale

The paper's central chain (Prop. 1: left Bregman MEB circumcenter equals power-MEB circumcenter of Laguerre-weighted dual points; Prop. 2: BregmanBC coincides with FWPowerMEB in dual space; Prop. 4: Bregman potential lifting matches shifted paraboloid lifting) is presented as direct algebraic rewriting of the definitions of Bregman divergence, Fenchel-Young divergence, mixed weighted squared Euclidean form (Eq. 1), and power distance. There are no fitted parameters, no empirical predictions, and no uniqueness theorem used to forbid alternatives. Self-citations (Nock & Nielsen 2005 for BregmanBC; Boissonnat–Nielsen–Nock 2010 for BVD as PD) point to independently published, externally checkable geometric/algorithmic results and are not used as unexamined premises that force the new MEB identity. Prop. 3 explicitly attributes the BVD=PD fact to [8] rather than re-deriving it as novel. Any algebraic sign/structural slips in the §2 chain are correctness defects, not circularity: wrong expansion is not reduction-to-input-by-construction. Score 1 only for the routine restatement of the authors' prior BVD result; the MEB equivalence claim itself is not circular.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The paper rests on standard convex-analysis facts (Legendre-type generators, Fenchel conjugacy, power distance) plus the algorithmic template of Frank-Wolfe / Badǐiu-Clarkson. No numerical parameters are fitted; the only modeling choice is the Legendre-type hypothesis needed for duality. No new physical or geometric entities are postulated.

axioms (3)
  • domain assumption F is a smooth strictly convex Legendre-type generator on an open convex domain (gradient norm diverges at the boundary).
    Required from §2 onward so that η= abla F is bijective and the dual Bregman / Fenchel-Young identities hold globally.
  • standard math Power distance σ(x,p̂)=∥x-p∥²-w and its bisectors / diagrams are the classical Laguerre objects.
    Used throughout Propositions 1–3; taken from Aurenhammer 1987 and standard computational-geometry references.
  • ad hoc to paper Frank-Wolfe curvature of the power-MEB objective is finite and bounded by 8R*.
    Stated without full derivation in Appendix A; needed for the O(1/ε) iteration claim.

reviewed 2026-07-31 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Minimum enclosing Bregman balls made easy." pith.science (2026). https://pith.science/paper/PQH2TADW

@misc{pith2026260724197,
  author       = {Pith},
  title        = {Pith review of: Minimum enclosing Bregman balls made easy},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PQH2TADW}},
  note         = {Machine review of arXiv:2607.24197}
}
Share X Bluesky LinkedIn Reddit HN
abstract

In this work, we revisit the problem of computing minimum enclosing Bregman balls (Bregman MEBs) of finite sets of parameters. First, we show that Bregman MEBs are equivalent to MEBs of corresponding weighted point sets with respect to the power distance. We then report an efficient Frank--Wolfe $(1+\epsilon)$-approximation algorithm for computing power MEBs, for any $\epsilon>0$. This power MEB approximation algorithm coincides with the Bregman MEB approximation algorithm of Nock and Nielsen (2005) when expressed in the dual gradient space. Finally, we show that the Bregman potential lifting transforms used to construct Bregman Voronoi diagrams can be reinterpreted as the classical paraboloid lifting transform applied to corresponding weighted point sets. In particular, Bregman MEB circumcenters lie on the farthest Bregman Voronoi diagrams or equivalently on the corresponding farthest power diagrams.

Figures

Figures reproduced from arXiv: 2607.24197 by Frank Nielsen.

Figure 1
Figure 1. Figure 1: Bad˘oiu–Clarkson’s approximation of the Euclidean MEB after [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The right KL minimax center of a set of n = 8 univariate Gaussians (right column) is computed equivalently as a left Bregman circumcenter of their corresponding natural parameters θ = ( µ σ2 , − 1 2σ2 ). The top row shows an example where the Bregman MEB is defined by two support distributions while the bottom row displays an example with the Bregman MEB defined by three support distributions. two Gaussian… view at source ↗
Figure 3
Figure 3. Figure 3: Laguerre geometry: The power bisector between two weighted points with non-negative [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Top: Power bisector between (a) two n = 2 weighted points with non-negative weights corresponding to squared radius, (b) weights shifted by some amount W > 0 showing that power bisectors coincide with (a), and (c) vanishing weights yielding the ordinary Euclidean bisector. Bottom: (a) Power diagram (PD) for n = 8 points with some empty Voronoi cells, (b) PD with weights shifted by the same amount W > 0 ill… view at source ↗
Figure 5
Figure 5. Figure 5: Bregman Voronoi diagrams (eKLD) of n = 8 points with equivalent power diagrams obtained for different values of λ. The power diagram (PD) of a weighted point set {pˆi} is the Voronoi diagram with respect to the power distance where the Voronoi cell of ˆpi is defined by Vorσ(ˆpi) = n x ∈ R d : σ(x, pˆi) ≤ σ(x, pˆj ) ∀j ∈ [n] o . Notice that some cells Vor(ˆpi) = {x ∈ R d : σ(x, pˆi) ≤ σ(x, pˆj ) ∀i ∈ [n]} o… view at source ↗
Figure 6
Figure 6. Figure 6: Correspondence between (left) Bregman Voronoi diagrams (BVDs) and power diagrams [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The supporting hyperplane HF θ of the potential graph F at ˆθ = (θ, F(θ)) is parallel to the lifted hyperplane HP p at p = 1 2 η of the paraboloid graph P. By shifting vertically HP p by the amounts w = 1 4 ∥η∥ 2 − ⟨θ, η⟩ + F(θ), the hyperplanes match: HP pˆ = HF θ . Proposition 4 The tangent hyperplane at ˆθ = (θ, F(θ)) to F coincides with the tangent hyperplane at (η,⟨η, η⟩) of the paraboloid P shifted v… view at source ↗
Figure 8
Figure 8. Figure 8: Left: Equivalence of hyperplanes HP pˆ = HF θ illustrated for the case of the Shannon negentropy Bregman generator: F(θ) = θ log θ − θ. Right: Equivalent of hyperplanes H F0 pˆ = HF θ for the case of F0(θ) = − log θ (Burg negentropy) and F(θ) = θ log θ − θ. = ⟨θ − θ ′ , ∇F(θ) − ∇F(θ ′ )⟩ (11) for a d-dimensional Bregman generator F(θ) is generally not a Bregman divergence (except for quadratic Bregman gene… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

This paper was first reviewed by grok-4.5 on July 31, 2026.