REVIEW 3 minor 42 references
Semi-discrete convex order and Laguerre tessellation fitting
T0 review · 0 major / 3 minor · reviewed 2026-06-30 · grok-4.3
Pith's one-line read Reconstructing a Laguerre tessellation from cell barycenters reduces to a Wasserstein projection onto discrete measures dominated in convex order by an absolutely continuous measure.
desk verdict The paper recasts Laguerre fitting as a Wasserstein projection onto convex-order dominated measures, which is a clean geometric move but rests on an AC assumption whose practical impact needs checking. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The set of discrete measures dominated in convex order by an absolutely continuous measure, with the Wasserstein projection onto this set serving as the approximation device for Laguerre reconstruction.
What would settle it
A concrete counter-example in which the Wasserstein projection onto the convex-order set produces cell volumes that deviate substantially from the prescribed volumes would show the approximation does not work.
Extended reading notes
Core claim
The reconstruction problem of finding a Laguerre tessellation with prescribed cell volumes from the barycenters of its cells admits a geometric interpretation as finding a discrete measure dominated in convex order by an absolutely continuous measure. The problem can therefore be solved approximately by computing the Wasserstein projection onto the set of all such discrete measures.
Load-bearing premise
The target measure must be absolutely continuous so that the convex-order domination relation is well-defined and the projection supplies a useful approximation.
Editorial extensions
If this is right
- The exact Laguerre reconstruction problem is replaced by a tractable convex-order projection that can be computed numerically.
- The same projection procedure yields a Laguerre tessellation fit even when the supplied barycenters do not come from any Laguerre tessellation.
- The method directly supplies a practical tool for fitting convex partitions to experimental data such as EBSD images in materials science.
Reading between the lines
- The convex-order viewpoint may extend to other semi-discrete fitting problems where cell volumes and centers must be matched simultaneously.
- Iterative refinement around the projection could convert the approximate solution into an exact one when the data are consistent.
- Numerical schemes for Wasserstein projection on this set could be reused as subroutines in related optimal-transport discretizations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies reconstruction of a Laguerre tessellation from prescribed cell volumes and barycenters. It establishes a geometric equivalence between this problem and Wasserstein projection onto the set of discrete measures dominated in convex order by a given absolutely continuous measure, shows that the reconstruction can be solved approximately via this projection, and extends the method to fitting Laguerre tessellations to arbitrary barycenter data. A concrete application to fitting a tessellation to an EBSD image of steel microstructure is presented.
Significance. If the claimed geometric link and approximation result hold with controlled error, the work supplies a new optimal-transport route to a class of inverse problems that arise in computational geometry, imaging, and materials science. The explicit reduction to a Wasserstein projection onto a convex-order constrained set is a clean conceptual contribution; the EBSD example demonstrates immediate applicability.
minor comments (3)
- The abstract and introduction state that the reconstruction is solved 'approximately' by the Wasserstein projection, but the precise sense of approximation (e.g., in which metric, under what quantitative error bound) is not made explicit in the opening paragraphs; a short clarifying sentence would help readers.
- Notation for the convex-order domination relation and the admissible set of discrete measures should be introduced once, early, and used consistently; occasional re-definition of symbols across sections slows reading.
- In the materials-science application, the precise preprocessing steps that turn the EBSD image into a point cloud of barycenters and target volumes are only sketched; a short algorithmic box or pseudocode would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. No major comments were provided in the report, so we have no specific points requiring response or revision at this stage. We will proceed with minor polishing as appropriate for the final version.
Circularity Check
No circularity: derivation uses standard convex order and Wasserstein properties
full rationale
The paper's central result interprets Laguerre reconstruction as an approximate Wasserstein projection onto the convex-order dominated set. This rests on established properties of convex order (for absolutely continuous measures) and the Wasserstein metric, which are external to the paper and not defined or fitted inside it. No equations reduce a prediction to a fitted input by construction, no uniqueness theorem is imported via self-citation, and the absolute-continuity hypothesis is stated explicitly as an enabling assumption rather than smuggled in. The EBSD application is a downstream use case, not part of the derivation chain. The argument is therefore self-contained against external benchmarks.
Assumptions & free parameters
assumptions (1)
- standard math Standard properties of convex order and Wasserstein distance between measures
Cite this review
Pith. "Pith review of Semi-discrete convex order and Laguerre tessellation fitting." pith.science (2026). https://pith.science/paper/UHNSSFMB
@misc{pith2026260629913,
author = {Pith},
title = {Pith review of: Semi-discrete convex order and Laguerre tessellation fitting},
year = {2026},
howpublished = {\url{https://pith.science/paper/UHNSSFMB}},
note = {Machine review of arXiv:2606.29913}
}
read the original abstract
Laguerre tessellations offer an efficient way to parameterize a large class of convex partitions of Euclidean space using only a set of points and scalar weights. For this reason, they have become popular in computational geometry, imaging and numerical analysis, both as a modeling and a discretization tool. In this paper we study the problem of reconstructing a Laguerre tessellation with prescribed cell volumes from the barycenters of its cells. We establish a geometric interpretation of this problem in terms of the set of discrete measures dominated in convex order by an absolutely continuous measure. In particular, we show that the reconstruction problem can be solved approximately by computing a Wasserstein projection onto this set. More generally, our method can also be applied to fit a Laguerre tessellation to an arbitrary set of barycenters. We give a concrete application of this in materials science, of fitting a Laguerre tessellation to an electron backscatter diffraction (EBSD) image of a steel.
Figures
Figures from the paper (13 more)
Reference graph
Works this paper leans on
-
[1]
Kadets-type theorems for partitions of a convex body
Arseniy Akopyan and Roman Karasev. Kadets-type theorems for partitions of a convex body. Discrete & Computational Geometry, 48:766–776, 2012. SEMI-DISCRETE CONVEX ORDER AND LAGUERRE TESSELLATIONS 31 Figure 12.Generators of true Laguerre tessellation in Figure 11 (left) and generators for reconstruction fort= 100 (subgradient descent, center, and Frank-Wol...
2012
-
[2]
Aur´ elien Alfonsi, Jacopo Corbetta, and Benjamin Jourdain. Sampling of one-dimensional proba- bility measures in the convex order and computation of robust option price bounds.International Journal of Theoretical and Applied Finance, 22(03):1950002, 2019
2019
-
[3]
Sampling of probability measures in the convex order by Wasserstein projection.Annales de l’Institut Henri Poincar´ e - Probabilit´ es et Statistiques, 56(3):1706–1729, 2020
Aur´ elien Alfonsi, Jacopo Corbetta, and Benjamin Jourdain. Sampling of probability measures in the convex order by Wasserstein projection.Annales de l’Institut Henri Poincar´ e - Probabilit´ es et Statistiques, 56(3):1706–1729, 2020
2020
-
[4]
A criterion for the affine equivalence of cell complexes inR d and convex polyhedra inR d+1.Discrete & Computational Geometry, 2(1):49–64, 1987
Franz Aurenhammer. A criterion for the affine equivalence of cell complexes inR d and convex polyhedra inR d+1.Discrete & Computational Geometry, 2(1):49–64, 1987
1987
-
[5]
Power diagrams: properties, algorithms and applications.SIAM Journal on Computing, 16(1):78–96, 1987
Franz Aurenhammer. Power diagrams: properties, algorithms and applications.SIAM Journal on Computing, 16(1):78–96, 1987
1987
-
[6]
Minkowski-type theorems and least- squares clustering.Algorithmica, 20(1):61–76, 1998
Franz Aurenhammer, Friedrich Hoffmann, and Boris Aronov. Minkowski-type theorems and least- squares clustering.Algorithmica, 20(1):61–76, 1998. 32 D. P. BOURNE, T. O. GALLOU ¨ET, Q. M ´ERIGOT, AND A. NATALE Figure 14.Convergence of the scheme for the data represented in Fig- ure 11 andt= 100, shown in terms ofG N(Yk)−G N(YK) for the subgra- dient scheme (...
1998
-
[7]
Bauschke and Walaa M
Heinz H. Bauschke and Walaa M. Moursi.An Introduction to Convexity, Optimization, and Al- gorithms. SIAM, 2023
2023
-
[8]
SIAM, 2017
Amir Beck.First-order methods in optimization. SIAM, 2017
2017
Show all 42 references
-
[9]
Gerald A. Beer. The Hausdorff metric and convergence in measure.Michigan Mathematical Jour- nal, 21(1):63–64, 1974
1974
-
[10]
Generalized Voronoi tessellation as a model of two-dimensional cell tissue dynamics.Bulletin of Mathematical Biology, 72:1696–1731, 2010
Martin Bock, Amit Kumar Tyagi, Jan-Ulrich Kreft, and Wolfgang Alt. Generalized Voronoi tessellation as a model of two-dimensional cell tissue dynamics.Bulletin of Mathematical Biology, 72:1696–1731, 2010
2010
-
[11]
Bourne, Piet J
David P. Bourne, Piet J. J. Kok, Steven M. Roper, and Wil D. T. Spanjer. Laguerre tessellations and polycrystalline microstructures: a fast algorithm for generating grains of given volumes. Philosophical Magazine, 100(21):2677–2707, 2020
2020
-
[12]
Bourne, Mason Pearce, and Steven M
David P. Bourne, Mason Pearce, and Steven M. Roper. Inverting Laguerre tessellations: Recov- ering tessellations from the volumes and centroids of their cells using optimal transport.ESAIM: Mathematical Modelling and Numerical Analysis, 59(5):841–871, 2025
2025
-
[13]
Bourne and Steven M
David P. Bourne and Steven M. Roper. Centroidal power diagrams, Lloyd’s algorithm, and ap- plications to optimal location problems.SIAM Journal on Numerical Analysis, 53(6):2545–2569, 2015
2015
-
[14]
On optimal weighted balanced clusterings: Gravity bodies and power diagrams.SIAM Journal on Discrete Mathematics, 26(2):415–434, 2012
Andreas Brieden and Peter Gritzmann. On optimal weighted balanced clusterings: Gravity bodies and power diagrams.SIAM Journal on Discrete Mathematics, 26(2):415–434, 2012
2012
-
[15]
Constrained clustering via diagrams: A unified theory and its application to electoral district design.European Journal of Operational Research, 263(1):18–34, 2017
Andreas Brieden, Peter Gritzmann, and Fabian Klemm. Constrained clustering via diagrams: A unified theory and its application to electoral district design.European Journal of Operational Research, 263(1):18–34, 2017
2017
-
[16]
Convex optimization: Algorithms and complexity.Foundations and Trends® in Machine Learning, 8(3-4):231–357, 2015
S´ ebastien Bubeck. Convex optimization: Algorithms and complexity.Foundations and Trends® in Machine Learning, 8(3-4):231–357, 2015
2015
-
[17]
Krzysztof J. Ciosmak. Applications of Strassen’s theorem and Choquet theory to optimal transport problems, to uniformly convex functions and to uniformly smooth functions.Nonlinear Analysis, 232:113267, 2023. 34 D. P. BOURNE, T. O. GALLOU ¨ET, Q. M ´ERIGOT, AND A. NATALE
2023
-
[18]
Optimal transport and the geometry ofL 1(Rd).Pro- ceedings of the American Mathematical Society, 142(10):3585–3596, 2014
Ivar Ekeland and Walter Schachermayer. Optimal transport and the geometry ofL 1(Rd).Pro- ceedings of the American Mathematical Society, 142(10):3585–3596, 2014
2014
-
[19]
An algorithm for quadratic programming.Naval research logistics quarterly, 3(1-2):95–110, 1956
Marguerite Frank, Philip Wolfe, et al. An algorithm for quadratic programming.Naval research logistics quarterly, 3(1-2):95–110, 1956
1956
-
[20]
Gallou¨ et, Andrea Natale, and Gabriele Todeschi
Thomas O. Gallou¨ et, Andrea Natale, and Gabriele Todeschi. Metric extrapolation in the Wasser- stein space.Calculus of Variations and Partial Differential Equations, 64(5):147, 2025
2025
-
[21]
On a mixture of Brenier and Strassen theorems.Proceedings of the London Mathematical Society, 120(3):434–463, 2020
Nathael Gozlan and Nicolas Juillet. On a mixture of Brenier and Strassen theorems.Proceedings of the London Mathematical Society, 120(3):434–463, 2020
2020
-
[22]
Kantorovich duality for general transport costs and applications.Journal of Functional Analysis, 273(11):3327–3405, 2017
Nathael Gozlan, Cyril Roberto, Paul-Marie Samson, and Prasad Tetali. Kantorovich duality for general transport costs and applications.Journal of Functional Analysis, 273(11):3327–3405, 2017
2017
-
[23]
Constrained clustering, diagrams, coresets, and their applications.Geometry, Analysis and Convexity, page 87, 2026
Peter Gritzmann. Constrained clustering, diagrams, coresets, and their applications.Geometry, Analysis and Convexity, page 87, 2026
2026
-
[24]
Recursively-regular subdivisions and applications.Journal of Com- putational Geometry, 7(1), 2016
Rafel Jaume and G¨ unter Rote. Recursively-regular subdivisions and applications.Journal of Com- putational Geometry, 7(1), 2016
2016
-
[25]
Statistical inference of convex order by Wasserstein projection.arXiv:2406.02840, 2024
Jakwang Kim, Young-Heon Kim, Yuanlong Ruan, and Andrew Warren. Statistical inference of convex order by Wasserstein projection.arXiv:2406.02840, 2024
2024
-
[26]
Backward and forward Wasserstein projections in stochas- tic order.Journal of Functional Analysis, 286(2):110201, 2024
Young-Heon Kim and Yuanlong Ruan. Backward and forward Wasserstein projections in stochas- tic order.Journal of Functional Analysis, 286(2):110201, 2024
2024
-
[27]
Lagrangian discretiza- tion of crowd motion and linear diffusion.SIAM Journal on Numerical Analysis, 58(4):2093–2118, 2020
Hugo Leclerc, Quentin M´ erigot, Filippo Santambrogio, and Federico Stra. Lagrangian discretiza- tion of crowd motion and linear diffusion.SIAM Journal on Numerical Analysis, 58(4):2093–2118, 2020
-
[28]
Emerson Le´ on and G¨ unter M. Ziegler. Spaces of convex n-partitions. InNew Trends in Intuitive Geometry, pages 279–306. Springer, 2018
2018
-
[29]
On the use of Laguerre tessellations for representations of 3d grain structures
Allan Lyckegaard, Erik Mejdal Lauridsen, Wolfgang Ludwig, Richard Warren Fonda, and Hen- ning Friis Poulsen. On the use of Laguerre tessellations for representations of 3d grain structures. Advanced Engineering Materials, 13(3):165–170, 2011
2011
-
[30]
Non-asymptotic convergence bounds for Wasserstein approximation using point clouds
Quentin M´ erigot, Filippo Santambrogio, and Cl´ ement Sarrazin. Non-asymptotic convergence bounds for Wasserstein approximation using point clouds. In M. Ranzato, A. Beygelzimer, Y. Dauphin, P.S. Liang, and J. Wortman Vaughan, editors,Proceedings of Advances in Neu- ral Infor...
2021
-
[31]
Mirrors, lenses and Monge–Amp` ere equations.European Mathematical Society Magazine, (120):16–28, 2021
Quentin M´ erigot and Boris Thibert. Mirrors, lenses and Monge–Amp` ere equations.European Mathematical Society Magazine, (120):16–28, 2021
2021
-
[32]
Optimal transport: discretization and algorithms
Quentin M´ erigot and Boris Thibert. Optimal transport: discretization and algorithms. In A. Bonito and R. H. Nochetto, editors,Geometric Partial Differential Equations - Part II, vol- ume 22 ofHandbook of Numerical Analysis, pages 133–212. North-Holland, 2021
2021
-
[33]
Initialization procedures for discrete and semi-discrete optimal transport
Jocelyn Meyron. Initialization procedures for discrete and semi-discrete optimal transport. Computer-Aided Design, 115:13–22, 2019
2019
-
[34]
Reconstruction of grains in polycrystalline materials from incomplete data using Laguerre tessellations.Microscopy and Microanalysis, 25(3):743–752, 2019
Lukas Petrich, Jakub Stanˇ ek, Mingyan Wang, Daniel Westhoff, Ludˇ ek Heller, PetrˇSittner, Carl E Krill III, Viktor Beneˇ s, and Volker Schmidt. Reconstruction of grains in polycrystalline materials from incomplete data using Laguerre tessellations.Microscopy and Microanalysi...
2019
-
[35]
On the sequential convergence of Lloyd’s algorithms.Mathematics of Operations Research, 51(2):1120–1138, 2025
Leo Portales, Elsa Cazelles, and Edouard Pauwels. On the sequential convergence of Lloyd’s algorithms.Mathematics of Operations Research, 51(2):1120–1138, 2025
2025
-
[36]
Romain Quey and Lo¨ ıc Renversade. Optimal polyhedral description of 3D polycrystals: Method and application to statistical and synchrotron X-ray diffraction data.Computer Methods in Applied Mechanics and Engineering, 330:308–333, 2018
2018
-
[37]
Birkh¨ auser, 2015
Filippo Santambrogio.Optimal Transport for Applied Mathematicians. Birkh¨ auser, 2015
2015
-
[38]
The existence of probability measures with given marginals.The Annals of Math- ematical Statistics, 36(2):423–439, 1965
Volker Strassen. The existence of probability measures with given marginals.The Annals of Math- ematical Statistics, 36(2):423–439, 1965
1965
-
[39]
¨Uber exponierte punkte abgeschlossener punktmengen.Fundamenta Mathe- maticae, 24(1):139–143, 1935
Stefan Straszewicz. ¨Uber exponierte punkte abgeschlossener punktmengen.Fundamenta Mathe- maticae, 24(1):139–143, 1935. SEMI-DISCRETE CONVEX ORDER AND LAGUERRE TESSELLATIONS 35
1935
-
[40]
John F. Toland. A duality principle for non-convex optimisation and the calculus of variations. Archive for Rational Mechanics and Analysis, 71:41–61, 1979
1979
-
[41]
Springer, 2009
C´ edric Villani.Optimal transport: old and new, volume 338. Springer, 2009
2009
-
[42]
Centroidal power diagrams with capacity constraints: Computation, applications, and extension
Shi-Qing Xin, Bruno L´ evy, Zhonggui Chen, Lei Chu, Yaohui Yu, Changhe Tu, and Wenping Wang. Centroidal power diagrams with capacity constraints: Computation, applications, and extension. ACM Transactions on Graphics, 35(6):244:1–244:12, 2016. David P. Bourne, Maxwell Institut...
2016
Reviewed June 30, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.