REVIEW 1 major objections 5 minor 17 references
Gromov-Wasserstein Bound between Reeb and Mapper Graphs
T0 review · 1 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read The Mapper graph converges to the Reeb graph in expected Gromov-Wasserstein distance at rate $n^{-\nu/(d+\alpha)}$, with $\nu = \min\{1/2, d/(p(d+1))\}$.
desk verdict The measure-aware framework is a genuine step forward, but the main rate theorem is unproven: the proof wrongly assumes a small-ball lower bound that Assumption 2 does not give. 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 load-bearing objects are the two metric measure spaces $(R,d_H,m\circ\pi^{-1})$ and $(M_n,d_H,m_n\circ\Delta^{-1})$, together with the p-Gromov-Wasserstein distance between them. The Hausdorff distance on compact subsets of the manifold makes the Reeb graph a metric space, and pushing the volume (or empirical) measure forward under the quotient map makes it an mm-space; the Mapper is handled by representing each simplex as a subset of the sample and pushing the empirical measure forward. The proof rides on an approximation-estimation decomposition, $GW_p(R,M_n)\le W_p(m\circ\pi^{-1},m_n\circ\pi^{-1}) + \left(n^{-1}\sum_i d_H([x_i]_{\sim_f},\Delta(x_i))^p\right)^{1/p}$. Estimation is handled by bounding the covering number of the Reeb graph (using the Bishop-Cheeger-Gromov volume comparison and Morse-level-set structure) and invoking the Weed-Bach empirical-measure rates; approximation is handled by the $(a,b)$-standard assumption plus the Chazal et al. Hausdorff rate, by a modulus-of-continuity argument on the cover, and by a coarea-formula estimate of the volume of preimages of refined cover elements. The mechanism is that each of the two terms imposes the same optimal resolution scale, $r(n)\sim n^{1/(d+\alpha)}$, and their exponents combine into the final rate.
What would settle it
Take the circle with a smooth probability density that vanishes at one point, e.g., proportional to $\sin^2(\theta)$; at that point $m(B(x,r))$ decays like $r^3$ rather than $r^1$, so no $(a,1)$-standard constant exists and the Hausdorff bound used in Proposition 5.4 has no justification. Simulating the Mapper on such a distribution and comparing the empirical rate of $E[GW_p]$ against $n^{-\nu/(d+\alpha)}$ would directly test whether the main theorem survives this failure.
Extended reading notes
Core claim
The central claim is that the Mapper graph converges to its Reeb target in expected p-Gromov-Wasserstein distance at an explicit algebraic rate, and that the rate splits cleanly into an estimation term and an approximation term. The estimation term is the Wasserstein distance between the pushed-forward sampling measure and its empirical version on the Reeb graph, which is controlled by a covering-number bound $N_\varepsilon(R)\lesssim (1/\varepsilon)^{2d}$ and the Weed-Bach empirical measure theorem. The approximation term measures how far, in Hausdorff distance, the Mapper simplex containing a sampled point lies from the Reeb element containing it; it is controlled by the Hausdorff sample-to-manifold rate and by a volume bound on preimages of refined cover elements, $\max_J \mathrm{Vol}(f^{-1}(J))\lesssim W(r)^{d/(d+1)}$, obtained from the coarea formula. The decisive contribution is that the proof gets the same dimensional trade-off from both sides: a resolution that balances them yields exponent $\nu/(d+\alpha)$.
Load-bearing premise
The rate rests on the sampling measure being $(a,d)$-standard—every ball of radius $r$ carrying at least a constant multiple of $r^d$ of the total mass—since that assumption drives the Hausdorff sample-to-manifold bound; absolute continuity with a bounded density and full support, as imposed by Assumption 2, does not by itself guarantee it when the density vanishes at a point.
Editorial extensions
If this is right
- With resolution $r(n)\sim n^{1/(d+\alpha)}$, the expected p-Gromov-Wasserstein distance between Mapper and Reeb graphs is $O(n^{-\nu/(d+\alpha)})$, so the two graphs are statistically distinguishable at an explicit rate as $n$ grows.
- Because $\hat{GW}_p\le GW_p$ for the alternative formulation of Mémoli, the same rate automatically holds for the computationally popular $\hat{GW}_p$ distance.
- The estimation error term converges at the empirical-measure rate for the Reeb graph, whose covering number grows only like $(1/\varepsilon)^{2d}$; this shows the Reeb graph is a low-complexity target for measure estimation.
- The approximation error reflects the volume of preimages of cover elements: with maximal width $W(r)$, it decays like $W(r)^{d/(d+1)}$, which is what forces the resolution to grow as $n^{1/(d+\alpha)}$.
- By Remark 1, when $d/(p(d+1))\le 1/2$ the volume bound is saturated (parabola example), so the stated dependence on $d$ and $p$ cannot be improved in that regime.
Reading between the lines
- Editorial inference: if the rate holds in expected GW distance, it should also yield high-probability statements by Markov's inequality, and the exponential tail in the approximation bound suggests nearly sure convergence along carefully chosen resolutions; the paper only states the expectation.
- Editorial inference: the proof's dependence on the $(a,d)$-standard property means that densities that vanish somewhere on the manifold may genuinely slow the rate; checking this with a density proportional to a power of distance from a point would separate the paper's conclusion from its weakest assumption.
- Editorial inference: because the GW distance is sensitive to node masses, the framework can separate Mapper graphs that are identical as simplicial complexes, as the torus experiments illustrate; this suggests practical use for comparing data sets with different sampling distributions even when topological summaries coincide.
- Editorial inference: the same approximation-estimation strategy should extend to Mapper complexes of higher dimension or to Reeb spaces with multivariate filters, provided a covering-number control and a coarea-type volume estimate are available.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper treats Reeb and Mapper graphs of a Morse function on a compact Riemannian manifold as metric measure spaces: the Reeb graph is equipped with the Hausdorff metric on its level-set components and the pushforward of the sampling measure, while the Mapper graph carries the empirical pushforward measure. The authors prove geometric facts about the Hausdorff metric on the Reeb graph, derive covering-number bounds, and propose an approximation-estimation decomposition for the p-Gromov-Wasserstein distance between the Reeb graph and a random Mapper graph. Their main result (Theorem 5.1) claims that under Assumptions 1 and 2, with Mapper resolution r(n) ~ n^{1/(d+alpha)}, the expected GW distance decays as n^{-nu/(d+alpha)} with nu = min{1/2, d/(p(d+1))}. The proof combines a Wasserstein empirical-measure bound, a bound on the Hausdorff distance between the sample and the manifold, and a volume estimate for preimages of cover elements.
Significance. The paper addresses a real gap in TDA-inference: existing Mapper-to-Reeb convergence results use combinatorial or interleaving distances and ignore the sampling measure, whereas here the measure is encoded in the objects themselves and the GW distance gives a natural measure-aware comparison. The approximation-estimation decomposition in Proposition 5.1 is a clean and potentially reusable idea, and the paper uses external tools such as Weed-Bach, Chazal et al., and Morse theory rather than circular reasoning. The public code for computing GW distances between Mapper graphs is a useful addition. The significance is conditional, however, because the main theorem's stated assumptions are insufficient; the flaw is localized and repairable, but the hypothesis set must be changed.
major comments (1)
- [Section 5.4.4, Proposition 5.4] The assertion that 'the measure m satisfies the (a, b)-standard assumption for b = d ... by Assumption 2' is false. Assumption 2 gives only an upper bound on the density and full support; it does not provide a uniform lower bound on m(B(x,r))/r^d. On a flat torus with density proportional to d(p,x)^2 near some point p, one has m(B(p,r)) ~ c r^{d+2}, so no a>0 satisfies m(B(p,r)) >= a r^d for all small r, despite the measure being absolutely continuous with bounded density and full support. This invalidates the invocation of Theorem 5.4 and leaves the indicator term 1_{dH(M,X_n)>lambda} in the proof of Proposition 5.4 uncontrolled; consequently Theorem 5.1 is not established under the stated assumptions. The proof is repairable by explicitly assuming the (a,d)-standard condition or a positive lower bound on the density, but the theorem statement and the discussion around Definition 8 must be changed accordingly.
minor comments (5)
- [Section 5.4.1, Lemma 5.3] In the proof of Lemma 5.3, the sentence 'The same argument made in the proof of Lemma 5.3' should refer to Lemma 5.2.
- [Section 5.4.1, Lemma 5.3] The step 'We can also immediately see that dH(Xn cap S, S_lambda) <= lambda' is valid, but it deserves one sentence of justification: for p in S_lambda the ball B(p,lambda) remains inside the connected component S, so any sample point within lambda of p lies in S and hence in Delta Delta(xi).
- [Section 1 and Theorem 5.1] The introduction states a convergence rate 'for any p > 0' while Theorem 5.1 assumes p >= 1; the two statements should be harmonized.
- [Section 4] The map denoted Delta Delta (the simplex-assignment map) is not given a proper name; define it explicitly, e.g., as the assignment map sending each sample point to the unique Mapper simplex that contains it.
- [Section 5.4.4, Proposition 5.4] The constant a appears in the exponent of the final bound without being introduced in the theorem hypotheses; after the suggested repair, state explicitly which assumption supplies a and how it enters the rate.
Circularity Check
No significant circularity: the convergence-rate derivation is built from external statistical and geometric results, with no fitted parameter or self-referential construction.
full rationale
The paper's main theorem 5.1 is derived through an explicit approximation-estimation decomposition (Proposition 5.1). The estimation term is controlled by covering-number estimates for the Reeb graph (Proposition 5.2) combined with the external Weed--Bach empirical-measure bound (Proposition 5.3, Corollary 5.1.1). The approximation term is controlled by direct geometric bounds on Mapper graph elements (Lemmas 5.2 and 5.3), a volume estimate for cover preimages (Lemma 5.5), and an external Hausdorff concentration bound from Chazal et al. (Theorem 5.4). None of these ingredients is defined in terms of the target Gromov--Wasserstein distance, no parameter is fitted to data to produce the rate, and the resolution choice r(n) ~ n^{1/(d+alpha)} is a design choice, not a fitted input. The paper does contain a substantive mathematical gap in Proposition 5.4: Assumption 2 only gives an upper density bound plus full support, which does not imply the (a,d)-standard lower bound needed to apply Theorem 5.4, so the stated proof is incomplete. This is a correctness or assumptions gap, not circularity: the missing condition is an additional hypothesis, not a hidden restatement of the conclusion. Self-citations such as [CM22] and [CMO18] appear only as context or related work and are not load-bearing in the proof chain. Accordingly, no circular step can be exhibited from the paper's own equations.
Assumptions & free parameters
free parameters (1)
- a (small-ball constant) =
unstated
assumptions (7)
- standard math Morse theory background: Morse Lemma, gradient flow, cylindrical level-set structure, coarea formula, Bishop-Gromov volume comparison.
- standard math Weed-Bach Proposition 5.3: empirical measure Wasserstein convergence under covering-number conditions.
- standard math Chazal et al. Theorem 2: Hausdorff convergence rate of a sample to its support under (a,b)-standard measures.
- ad hoc to paper The sampling measure is (a,d)-standard for some a>0.
- domain assumption Perfect clustering oracle for Mapper connected components.
- domain assumption Assumption 1: Ricci curvature lower bound Ric >= (d-1)k.
- domain assumption Assumption 2: m absolutely continuous with respect to Vol, bounded density, fully supported.
Cite this review
Pith. "Pith review of Gromov-Wasserstein Bound between Reeb and Mapper Graphs." pith.science (2026). https://pith.science/paper/VCYQNDLX
@misc{pith2026250602810,
author = {Pith},
title = {Pith review of: Gromov-Wasserstein Bound between Reeb and Mapper Graphs},
year = {2026},
howpublished = {\url{https://pith.science/paper/VCYQNDLX}},
note = {Machine review of arXiv:2506.02810}
}
read the original abstract
Since its introduction as a computable approximation of the Reeb graph, the Mapper graph has become one of the most popular tools from topological data analysis for performing data visualization and inference. However, finding an appropriate metric (that is, a tractable metric with theoretical guarantees) for comparing Reeb and Mapper graphs, in order to, e.g., quantify the rate of convergence of the Mapper graph to the Reeb graph, is a difficult problem. While several metrics have been proposed in the literature, none is able to incorporate measure information, when data points are sampled according to an underlying probability measure. The resulting Reeb and Mapper graphs are therefore purely deterministic and combinatorial, and substantial effort is thus required to ensure their statistical validity. In this article, we handle this issue by treating Reeb and Mapper graphs as metric measure spaces. This allows us to use Gromov-Wasserstein metrics to compare these graphs directly in order to better incorporate the probability measures that data points are sampled from. Then, we describe the geometry that arises from this perspective, and we derive rates of convergence of the Mapper graph to the Reeb graph in this context. Finally, we showcase the usefulness of such metrics for Reeb and Mapper graphs in a few numerical experiments.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
‘Estimating the Reach of a Manifold’
[Aam+19] Eddie Aamari et al. ‘Estimating the Reach of a Manifold’. In: Electronic Journal of Statistics (2019). [ACL11] Ery Arias-Castro, Guangliang Chen and Gilad Lerman. ‘Spectral clustering based on local linear approximations’. In: Electronic Journal of Statistics 5 (2011), pp. 1537–
work page 2019
-
[5]
isbn: 978-3-03868-004-8. doi: 10 . 2312 / 3dor . 20161084. url: https://doi.org/10.2312/3dor.20161084. [BG14] Emmanuel Boissard and Thibaut Le Gouic. ‘On the mean speed of convergence of empirical and occupation measures in Wasserstein distance’. In: Annales de l’Institut Henri Poincar´ e, Probabilit´ es et Statistiques50.2 (2014), pp. 539–563. doi: 10.12...
arXiv 2014
-
[10]
Covariance alignment: from maximum likelihood estimation to Gromov-Wasserstein
42 [HRS23] Yanjun Han, Philippe Rigollet and George Stepaniants. ‘Covariance alignment: from maximum likelihood estimation to Gromov-Wasserstein’. In:arXiv preprint arXiv:2311.13595 (2023). [JPH21] Brian B Joseph, Trami Pham and Christopher Hastings. ‘Topological Data Analysis in Conjunction with Traditional Machine Learning Techniques to Predict Future M...
work page Pith review arXiv 2023
-
[12]
‘Finding the homology of sub- manifolds with high confidence from random samples’
[NSW08] Partha Niyogi, Stephen Smale and Shmuel Weinberger. ‘Finding the homology of sub- manifolds with high confidence from random samples’. In: Discrete & Computational Geometry 39 (2008), pp. 419–441. [Oul] Ziyad Oulhaj. Mapper Gromov-Wasserstein distance examples . url: https://github. com/ZiyadOulhaj/Mapper-MMS. [Pet06] P Petersen. ‘Riemannian geome...
arXiv 2008
-
[16]
‘Duality and Sample Complexity for the Gromov-Wasserstein Distance’
[Zha+23] Zhengxin Zhang et al. ‘Duality and Sample Complexity for the Gromov-Wasserstein Distance’. In: NeurIPS 2023 Workshop Optimal Transport and Machine Learning
work page 2023
-
[51]
Schloss Dagstuhl – Leibniz-Zentrum f¨ ur Informatik, 2016, 53:1–53:16.isbn: 978-3-95977-009-5
Leibniz In- ternational Proceedings in Informatics (LIPIcs). Schloss Dagstuhl – Leibniz-Zentrum f¨ ur Informatik, 2016, 53:1–53:16.isbn: 978-3-95977-009-5. doi: 10.4230/LIPIcs.SoCG. 2016.53. [Nai+18] Gregory Naitzat et al. ‘M-Boost: profiling and refining deep neural networks with topological data analysis’. In: KDD Workshop on Interactive Data Exploratio...
-
[131]
‘Spectral Gromov-Wasserstein distances for shape matching’
[M´ em09] Facundo M´ emoli. ‘Spectral Gromov-Wasserstein distances for shape matching’. In: 2009 IEEE 12th International Conference on Computer Vision Workshops, ICCV Workshops. IEEE. 2009, pp. 256–263. [M´ em11] Facundo M´ emoli. ‘Gromov–Wasserstein distances and the metric approach to object matching’. In: Foundations of computational mathematics 11 (20...
work page 2011
-
[293]
Schloss Dagstuhl – Leibniz-Zentrum f¨ ur Inform- atik, 2024, 80:1–80:18. isbn: 978-3-95977-316-4. doi: 10.4230/LIPIcs.SoCG.2024.80. [Wan20] Ziqi Wang. ‘Exploration of Topological Data Analysis In 3D Printing’. In: 2020 In- ternational Conference on Information Science, Parallel and Distributed Systems (IS- PDS). IEEE. 2020, pp. 150–153. [WB19] Jonathan We...
Show all 17 references
-
[1587]
‘Exposition and interpretation of the topology of neural networks’
[BC18] Rickard Br¨ uel-Gabrielsson and Gunnar Carlsson. ‘Exposition and interpretation of the topology of neural networks’. In: CoRR. arXiv:1810.03234,
-
[1963]
‘Experiments on Fraud Detection use case with QML and TDA Mapper’
[MR21] Satanik Mitra and Kameshwar Rao JV. ‘Experiments on Fraud Detection use case with QML and TDA Mapper’. In: 2021 IEEE International Conference on Quantum Computing and Engineering (QCE) . 2021, pp. 471–472. doi: 10.1109/QCE52317.2021. 00083. [MW16] Elizabeth Munch and Be...
2021
-
[1986]
‘Probabilistic convergence and stability of random mapper graphs’
[Bro+21] Adam Brown et al. ‘Probabilistic convergence and stability of random mapper graphs’. In: Journal of Applied and Computational Topology 5.1 (2021), pp. 99–140. [Cha+14] Fr´ ed´ eric Chazal et al. ‘Convergence rates for persistence diagram estimation in to- pological da...
2021
-
[1996]
‘Topological methods for the analysis of high dimensional data sets and 3d object recognition.’ In: PBG@ Euro- graphics 2 (2007), pp
[SMC+07] Gurjeet Singh, Facundo M´ emoli, Gunnar E Carlsson et al. ‘Topological methods for the analysis of high dimensional data sets and 3d object recognition.’ In: PBG@ Euro- graphics 2 (2007), pp. 091–100. [Stu06] Karl-Theodor Sturm. ‘On the geometry of metric measure spac...
2007
-
[2016]
by Alfredo Ferreira, Andrea Giachetti and Daniela Giorgi
Ed. by Alfredo Ferreira, Andrea Giachetti and Daniela Giorgi. ISSN: 1997-0471. The Euro- graphics Association,
1997
-
[2018]
‘Computing the homology of basic semialgebraic sets in weak exponential time’
[BCL18] Peter B¨ urgisser, Felipe Cucker and Pierre Lairez. ‘Computing the homology of basic semialgebraic sets in weak exponential time’. In: Journal of the ACM (JACM) 66.1 (2018), pp. 1–30. [BFL16] Ulrich Bauer, Barbara Di Fabio and Claudia Landi. ‘An Edit Distance for Reeb ...
2018
-
[2021]
‘Computing the Gromov-Wasserstein distance between two surface meshes using optimal transport’
[KDO23] Patrice Koehl, Marc Delarue and Henri Orland. ‘Computing the Gromov-Wasserstein distance between two surface meshes using optimal transport’. In: Algorithms 16.3 (2023), p
2023
-
[2023]
‘Gromov–Wasserstein distances: Entropic regularization, dual- ity and sample complexity’
[Zha+24] Zhengxin Zhang et al. ‘Gromov–Wasserstein distances: Entropic regularization, dual- ity and sample complexity’. In: The Annals of Statistics 52.4 (2024), pp. 1616–1645. 44
2024
-
[6941]
‘Topographical transcriptome mapping of the mouse medial gan- glionic eminence by spatially resolved RNA-seq’
[Zec+14] Sabrina Zechel et al. ‘Topographical transcriptome mapping of the mouse medial gan- glionic eminence by spatially resolved RNA-seq’. In: Genome biology 15 (2014), pp. 1–
2014
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.