REVIEW 3 major objections 5 minor 52 references
Convergence of Finite Element Methods for Ricci Flow
T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read Two-dimensional Ricci flow can be approximated by finite elements with proven error bounds of order h^{q+1}+h^{r+1} for both the metric and the Gaussian curvature.
desk verdict First convergence proof for a finite element discretization of Ricci flow, with the main theorem solid but the optimal-rate theorem only sketched and requiring a fuller proof before acceptance. 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 central object is the coupled system: the metric evolves as a solution-driven evolution with the Gaussian curvature as driving term, while the curvature satisfies a parabolic equation depending on the metric. This system is discretized with Regge finite elements for the symmetric (0,2)-tensor metric and Lagrange finite elements for the scalar curvature. The analysis uses metric-dependent projection operators for the metric equation, a Ritz projection for the curvature, and a matrix-vector mass/stiffness formulation that converts the discrete system into a differential-algebraic equation. Stability is obtained by estimating consistency defects, applying norm-equivalence lemmas for the dis
What would settle it
Run the scheme with r at least 1 and q at least 1 on a smooth initial metric for which the exact Ricci flow solution is known analytically, using a negligible time step, and check whether the L^2 errors of the metric and curvature decay like h^(q+1)+h^(r+1) over successive mesh refinements; a systematically slower rate would contradict Theorem 3.2.
Extended reading notes
Core claim
The central claim is that the coupled formulation (2.3) of the two-dimensional Ricci flow, metric evolving under Gaussian curvature and curvature solving a metric-dependent parabolic equation, admits convergent finite element spatial discretizations. With Regge elements of degree r and Lagrange elements of degree q, the semidiscrete error satisfies the bound stated in Theorem 3.1, and the optimal L^2 bound of Theorem 3.2 holds for r at least 1 and q at least 1. The proof controls consistency defects through metric-dependent projections and a matrix-vector formulation, closing with Grönwall stability estimates.
Load-bearing premise
The paper assumes that an exact Ricci flow solution (g,kappa) exists on the whole time interval [0,T] and is sufficiently smooth, with the metric uniformly positive definite; if that smooth existence fails, for example at a finite-time singularity, the error estimates no longer apply.
Editorial extensions
If this is right
- For every r at least 0 and q at least 1 the spatial semi-discretization converges, so choosing higher polynomial degrees yields higher-order accuracy in both the metric and the Gaussian curvature.
- For r at least 1 and q at least 1 the optimal L^2 rate h^{q+1}+h^{r+1} holds, meaning low-order elements already give O(h^2) accuracy.
- The discrete solutions satisfy the Gauss-Bonnet identity exactly, and the L^2-projection variant conserves area, preserving two geometric invariants of the continuum flow.
- The error estimates hold uniformly in time up to any T on which the exact solution stays smooth and uniformly positive definite, so the method is stable over the whole smooth existence window.
- The proof covers both proposed schemes, with the metric-dependent-projection variant valid for all r at least 0 and the L^2-projection variant requiring r at least 1 and q at least 1.
Reading between the lines
- Beyond the paper: the same consistency-defect and matrix-vector machinery may extend to other intrinsic flows whose curvature evolution is parabolic, such as Calabi flow, 3D Yamabe flow, or 3D Ricci flow, since the proof identifies that parabolic structure as the enabling mechanism.
- Beyond the paper: the numerical observation that the (r,q)=(0,1) case gives O(h^2) for curvature, one order better than Theorem 3.1, suggests the logarithmic losses in the theorem are an artifact of the proof and may be removable for lowest-order elements.
- Beyond the paper: the proposed embedding postprocessor could be tested independently by measuring how closely the discrete surface velocity satisfies the constraint (6.1) on examples with known isometric embeddings.
- Beyond the paper: the error estimates are conditional on smooth existence; for flows that form a finite-time singularity, the method may still produce useful approximations but the stated high-order rates should be expected to degrade.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes spatial semi-discretizations of a coupled reformulation of two-dimensional Ricci flow in which the metric evolves by the Gauss curvature and the curvature satisfies a metric-dependent parabolic equation (system (2.3)). The proposed schemes (3.3) and (3.4) use Regge elements of degree r for the metric and Lagrange elements of degree q for the curvature. The main convergence result, Theorem 3.1, states an L^p metric error plus L^2 curvature error of order (ln(1/h))^{\bar q} h^{q+1} + ln(1/h) h^{r+1} for scheme (3.3), under a smoothness assumption on the exact solution. Theorem 3.2 states an optimal L^2 rate h^{q+1}+h^{r+1} for r,q≥1 for both schemes. The paper also proves discrete Gauss-Bonnet and area conservation properties, proposes a post-processing embedding algorithm, and reports numerical experiments. Theorem 3.1 is proved in detail via consistency-defect estimates, a bootstrap stability argument, and Gronwall estimates; Theorem 3.2 is not proved, and the text explicitly says only that it follows by a similar approach with remarks on the differences.
Significance. If fully substantiated, this would be the first rigorous convergence proof for a finite element discretization of two-dimensional Ricci flow, and the techniques—solution-driven metric evolution, Regge interpolation, and a matrix-vector stability framework—are promising for other intrinsic curvature flows. The detailed proof of Theorem 3.1, together with the technical appendices on projection stability and Ritz projection time derivatives, is a substantial positive contribution. However, Theorem 3.2 is advertised in the abstract and introduction as an optimal L^2 error estimate, yet the proof is missing; the stability estimates on which it relies are only asserted in Remarks 4.4, 4.7, and 4.10. The central claim of the paper is therefore only partially established.
major comments (3)
- [Section 3.2 and Section 5] Theorem 3.2 is not proved. The text immediately after Theorem 3.2 states that it 'can be proved using a very similar approach' and only remarks are given; the proof of Theorem 3.2 in Section 5 simply invokes the unproved stability estimates (4.20a), (4.20b), (4.26a), and (4.26b) from Remarks 4.7 and 4.10. In particular, the L^2 analogue of Proposition 4.9 is not derived, and the constants in the bootstrap are not tracked. Since Theorem 3.2 is one of the two headline results and is explicitly advertised in the abstract and introduction as an optimal L^2 rate for both schemes, this omission is load-bearing. The paper should supply the full proof, or, if the proof cannot be completed, the theorem should be rephrased as a conjecture and removed from the list of established results.
- [Remark 4.7, Eq. (4.18)] The estimate (4.18) for the difference of L^2 projections is only sketched, and the sketch relies on bounding the term \|\partial_t g_h^* - d_{L2,g}\|_{L^\infty(M_h)} without stating how this is controlled. From (4.12a) one only has an L^2 bound on d_{L2,g}; the L^\infty bound presumably follows from an inverse estimate for r,q≥1, but this step is not given. More generally, the L^2 stability estimates (4.20) and (4.26) are asserted without the full derivation of the analogue of Proposition 4.9. Because these are the exact estimates needed to obtain the advertised rate h^{q+1}+h^{r+1}, this is not a purely technical presentation issue.
- [Section 3.2, hypotheses of Theorems 3.1 and 3.2] The theorems assume the exact solution is 'sufficiently smooth' without specifying the required Sobolev classes or norms. The proof uses W^{r+1,\infty}-type bounds on the metric, time-derivative bounds, and W^{1,\infty} bounds on the curvature, but the precise hypotheses are not stated. For unnormalized Ricci flow, finite-time singularities can occur, so the assumption is not vacuous even for smooth initial data. The paper should state a precise regularity hypothesis (e.g., g∈L^\infty(0,T;W^{r+2,\infty}(M)), \partial_t g∈L^\infty(0,T;W^{r+1,\infty}(M)), κ∈L^\infty(0,T;W^{1,\infty}(M)), etc.) and, ideally, indicate when this is satisfied by the normalized flow.
minor comments (5)
- [Abstract and Section 3.2] The abstract and introduction state that the proposed method preserves area conservation. Theorem 3.3 proves area conservation only for scheme (3.4); scheme (3.3) is not shown to have this property. Please qualify the claim accordingly.
- [Appendix B, Lemma B.1] The norm-equivalences in Lemma B.1 are labelled (C.1) and (C.4), but the appendix is labelled B; the equation numbers should be (B.1) and (B.4).
- [Lemma 3.5] The proof of Theorem 3.1 relies on the Ritz projection estimates of Lemma 3.5, whose proof invokes assumptions A1–A4 of [17] with the statement 'one can verify'. The verification is not shown. This is likely routine, but a short discussion of why the assumptions hold for the metric g_h^* would help the reader.
- [Section 2.2] The normalized Ricci flow constant \bar κ is defined, but the paper later also compares \bar κ_h to \bar κ in consistency estimates without displaying a proof of the bound \|\bar κ_h - \bar κ\|_{L^\infty} ≲ h^{q+1}+h^{r+1}. This estimate is plausible from the approximation properties of the initial data but should be stated explicitly with a brief justification.
- [Section 6] The embedding PDE (6.5) is presented without a convergence analysis, and Remark 6.2 notes that well-posedness is only established for planar domains. This is acceptable as a post-processing tool, but the text should more clearly separate this computational heuristic from the rigorously analyzed results in Sections 3–5.
Circularity Check
No significant circularity: consistency-plus-stability proof; self-citations are not load-bearing.
full rationale
The error analysis follows a standard consistency/stability paradigm. The discrete error functions e_g, e_kappa are defined against Regge interpolation and Ritz projection, and the error equations (4.3a)/(4.3b) and (4.2) are derived by subtracting the discrete scheme from the consistency-defect equations (4.1); nothing in these equations defines the output in terms of the input. The final rates in Theorem 3.1 and 3.2 are obtained by combining the consistency defect estimates (Lemma 4.5), which are genuine interpolation/projection approximation bounds, with stability estimates (Prop. 4.6, Prop. 4.9, and the L2 analogues sketched in Remarks 4.4, 4.7, 4.10) via Gronwall and a bootstrap time interval. No parameter is fitted to data and no prediction is an identity in disguise. Citations to the first author's earlier work [25,26,27] and to [9] supply the discretization and auxiliary geometric approximation lemmas, but the convergence conclusion is not an input to those cited results; in particular, the structure preservation (Theorem 3.3) is re-proved here. The one honest caveat is that Theorem 3.2 is not proved in full: Section 3.2 states 'Since Theorem 3.2 can be proved using a very similar approach, we only provide several remarks to clarify the key differences in the proof (see Remark 4.4, 4.7, 4.10).' That is a proof gap (and the regularity of the exact solution is assumed), but a missing proof is not circularity. Therefore no circular step is exhibited.
Assumptions & free parameters
assumptions (5)
- domain assumption Existence of a sufficiently smooth exact solution (g,κ) to the coupled system on [0,T] with g uniformly positive definite.
- domain assumption The triangulated surface M_h is a quasi-uniform approximation of M with a bijective lift a: M_h -> M for sufficiently small h.
- standard math The assumptions A1-A4 of Demlow [17, Section 3.1] are satisfied for the Ritz projection with respect to the time-dependent metric g*_h(t).
- standard math Regge interpolation approximation properties from [37, Section 2.3.2] and the metric perturbation estimates from Gawlik [26, Lemmas 4.5-4.6] hold.
- standard math Lemma 2.1: the Gaussian curvature of a metric evolving under Ricci flow satisfies the parabolic evolution equation (2.2).
Cite this review
Pith. "Pith review of Convergence of Finite Element Methods for Ricci Flow." pith.science (2026). https://pith.science/paper/XTUIZQ5N
@misc{pith2026260717051,
author = {Pith},
title = {Pith review of: Convergence of Finite Element Methods for Ricci Flow},
year = {2026},
howpublished = {\url{https://pith.science/paper/XTUIZQ5N}},
note = {Machine review of arXiv:2607.17051}
}
read the original abstract
The convergence of a finite element discretization for the two-dimensional Ricci flow is proved. In this method, the Ricci flow on a two-dimensional surface is formulated into solution-driven metric evolution, with the metric evolution driven by the Gauss curvature. The Gauss curvature satisfies a parabolic equation that in turn depends on the metric, thereby enhancing the parabolic structure of the problem. The solution-driven metric evolution formulation is discretized by the finite element method, and the convergence of finite element approximations is proved by adapting the matrix-vector formulation developed in the literature initially for studying solution-driven surface evolution in extrinsic curvature flow. In addition to its convergence, the proposed method also preserves important geometric structures of the Ricci flow at the discrete level, such as area conservation and the Gauss-Bonnet theorem. Extensive numerical experiments are presented to demonstrate the convergence of the proposed method as well as the simulation of Ricci flow.
Figures
Reference graph
Works this paper leans on
-
[1]
Bai and B
G. Bai and B. Li. Convergence of a stabilized parametric finite element method of the Barrett–Garcke– N¨ urnberg type for curve shortening flow.Mathematics of Computation, 2024
2024
-
[2]
B¨ ansch, P
E. B¨ ansch, P. Morin, and R. H. Nochetto. A finite element method for surface diffusion: the parametric case.Journal of Computational Physics, 203(1):321–343, 2005
2005
-
[3]
Bao and Y
W. Bao and Y. Li. A structure-preserving parametric finite element method for geometric flows with anisotropic surface energy.Numerische Mathematik, 156(2):609–639, 2024
2024
-
[4]
Bao and Q
W. Bao and Q. Zhao. A structure-preserving parametric finite element method for surface diffusion. SIAM Journal on Numerical Analysis, 59(5):2775–2799, 2021
2021
-
[5]
J. W. Barrett, K. Deckelnick, and V. Styles. Numerical analysis for a system coupling curve evolution to reaction diffusion on the curve.SIAM Journal on Numerical Analysis, 55(2):1080–1100, 2017
2017
-
[6]
J. W. Barrett, H. Garcke, and R. N¨ urnberg. A parametric finite element method for fourth order geometric evolution equations.Journal of Computational Physics, 222(1):441–467, 2007
2007
-
[7]
J. W. Barrett, H. Garcke, and R. N¨ urnberg. On the parametric finite element approximation of evolving hypersurfaces inR 3.Journal of Computational Physics, 227(9):4281–4307, 2008
2008
-
[8]
S. Bartels. A simple scheme for the approximation of the elastic flow of inextensible curves.IMA Journal of Numerical Analysis, 33(4):1115–1125, 2013
2013
Show all 52 references
-
[9]
Berchenko-Kogan and E
Y. Berchenko-Kogan and E. S. Gawlik. Finite element approximation of the Levi-Civita connection and its curvature in two dimensions.Foundations of Computational Mathematics, 24(2):587–637, 2024
2024
-
[10]
Bonito, R
A. Bonito, R. H. Nochetto, and M. S. Pauletti. Parametric FEM for geometric biomembranes.Journal of Computational Physics, 229(9):3171–3188, 2010
2010
-
[11]
Chow and F
B. Chow and F. Luo. Combinatorial Ricci flows on surfaces.Journal of Differential Geometry, 63(1):97– 129, 2003
2003
-
[12]
S. H. Christiansen. On the linearization of Regge calculus.Numerische Mathematik, 119:613–640, 2011
2011
-
[13]
S. H. Christiansen and K. Hu. Finite element systems for vector bundles: elasticity and curvature. Foundations of Computational Mathematics, 23(2):545–596, 2023
2023
-
[14]
S. H. Christiansen, K. Hu, and T. Lin. Extended Regge complex for linearized Riemann-Cartan geometry and cohomology.arXiv preprint arXiv:2312.11709, 2023
2023 arXiv
-
[15]
Deckelnick and G
K. Deckelnick and G. Dziuk. On the approximation of the curve shortening flow.Pitman Research Notes in Mathematics Series, pages 100–100, 1995
1995
-
[16]
Deckelnick and V
K. Deckelnick and V. Styles. Stability and error analysis for a diffuse interface approach to an advection– diffusion equation on a moving surface.Numerische mathematik, 139:709–741, 2018
2018
-
[17]
A. Demlow. Higher-order finite element methods and pointwise error estimates for elliptic problems on surfaces.SIAM Journal on Numerical Analysis, 47(2):805–827, 2009
2009
-
[18]
Duan and B
B. Duan and B. Li. New artificial tangential motions for parametric finite element approximation of surface evolution.SIAM Journal on Scientific Computing, 46(1):A587–A608, 2024. 24
2024
-
[19]
G. Dziuk. An algorithm for evolutionary surfaces.Numerische Mathematik, 58(1):603–611, 1990
1990
-
[20]
Dziuk and C
G. Dziuk and C. M. Elliott. Finite element methods for surface PDEs.Acta Numerica, 22:289–396, 2013
2013
-
[21]
C. M. Elliott, H. Garcke, and B. Kov´ acs. Numerical analysis for the interaction of mean curvature flow and diffusion on closed surfaces.Numerische Mathematik, 151(4):873–925, 2022
2022
-
[22]
H. Fritz. Numerical Ricci–DeTurck flow.Numerische Mathematik, 131(2):241–271, 2015
2015
-
[23]
Gao and W
H. Gao and W. Sun. Optimal analysis of non-uniform Galerkin-mixed finite element approximations to the Ginzburg–Landau equations in superconductivity.SIAM Journal on Numerical Analysis, 61(2):929– 951, 2023
2023
-
[24]
Garfinkle and J
D. Garfinkle and J. Isenberg. Numerical studies of the behavior of Ricci flow.Contemporary Mathematics, 367:103–114, 2005
2005
-
[25]
E. S. Gawlik. Finite element methods for geometric evolution equations. InInternational Conference on Geometric Science of Information, pages 532–540. Springer, 2019
2019
-
[26]
E. S. Gawlik. High-order approximation of Gaussian curvature with Regge finite elements.SIAM Journal on Numerical Analysis, 58(3):1801–1821, 2020
2020
-
[27]
Gopalakrishnan, M
J. Gopalakrishnan, M. Neunteufel, J. Sch¨ oberl, and M. Wardetzky. Analysis of curvature approximations via covariant curl and incompatibility for Regge metrics.The SMAI Journal of computational mathe- matics, 9, 2023
2023
-
[28]
X. D. Gu and S.-T. Yau.Computational conformal geometry. International Press, Somerville, MA; Higher Education Press, Beijing, 2008
2008
-
[29]
R. S. Hamilton. Three-manifolds with positive Ricci curvature.Journal of Differential geometry, 17(2):255–306, 1982
1982
-
[30]
Hu and B
J. Hu and B. Li. Evolving finite element methods with an artificial tangential velocity for mean curvature flow and Willmore flow.Numerische Mathematik, 152(1):127–181, 2022
2022
-
[31]
M. Jin, J. Kim, F. Luo, and X. Gu. Discrete surface Ricci flow.IEEE Transactions on Visualization and Computer Graphics, 14(5):1030–1043, 2008
2008
-
[32]
M. Jin, J. Kim, F. Luo, S. Lee, and X. Gu. Conformal surface parameterization using Euclidean Ricci flow.Technique report, May, 2006
2006
-
[33]
Kov´ acs, B
B. Kov´ acs, B. Li, and C. Lubich. A convergent evolving finite element algorithm for mean curvature flow of closed surfaces.Numerische Mathematik, 143:797–853, 2019
2019
-
[34]
Kov´ acs, B
B. Kov´ acs, B. Li, and C. Lubich. A convergent evolving finite element algorithm for Willmore flow of closed surfaces.Numerische Mathematik, 149(3):595–643, 2021
2021
-
[35]
Kov´ acs, B
B. Kov´ acs, B. Li, C. Lubich, and C. A. Power Guerra. Convergence of finite elements on an evolving surface driven by diffusion on the surface.Numerische Mathematik, 137:643–689, 2017
2017
-
[36]
B. Li. Convergence of Dziuk’s linearly implicit parametric finite element method for curve shortening flow.SIAM Journal on Numerical Analysis, 58(4):2315–2333, 2020
2020
-
[37]
Li.Regge finite elements with applications in solid mechanics and relativity
L. Li.Regge finite elements with applications in solid mechanics and relativity. PhD thesis, University of Minnesota, 2018
2018
-
[38]
C. M. Elliott and H. Fritz. On approximations of the curve shortening flow and of the mean curvature flow based on the DeTurck trick.IMA Journal of Numerical Analysis, 37(2):543–603, 2017
2017
-
[39]
Perelman
G. Perelman. The entropy formula for the Ricci flow and its geometric applications.arXiv preprint math/0211159, 2002
2002 arXiv
-
[40]
Perelman
G. Perelman. Finite extinction time for the solutions to the Ricci flow on certain three-manifolds.arXiv preprint math/0307245, 2003
2003 arXiv
-
[41]
Perelman
G. Perelman. Ricci flow with surgery on three-manifolds.arXiv preprint math/0303109, 2003
2003 arXiv
-
[42]
Pozzi and B
P. Pozzi and B. Stinner. Curve shortening flow coupled to lateral diffusion.Numerische Mathematik, 135:1171–1205, 2017
2017
-
[43]
Pozzi and B
P. Pozzi and B. Stinner. On motion by curvature of a network with a triple junction.The SMAI journal of computational mathematics, 7:27–55, 2021
2021
-
[44]
T. Regge. General relativity without coordinates.Il Nuovo Cimento (1955-1965), 19:558–571, 1961
1955
-
[45]
J. H. Rubinstein and R. Sinclair. Visualizing Ricci flow of manifolds of revolution.Experimental Mathe- matics, 14(3):285–298, 2005
2005
-
[46]
Sch¨ oberl
J. Sch¨ oberl. C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014
2014
-
[47]
E. Sonn, E. Saucan, E. Appelboim, and Y. Y. Zeevi. Ricci flow for image processing. In2014 IEEE 28th Convention of Electrical & Electronics Engineers in Israel (IEEEI), pages 1–5. IEEE, 2014
2014
-
[48]
W. Sun. Analysis of lowest-order characteristics-mixed FEMs for incompressible miscible flow in porous media.SIAM Journal on Numerical Analysis, 59(4):1875–1895, 2021
2021
-
[49]
J. C. Taft.Intrinsic geometric flows on manifolds of revolution. The University of Arizona, 2010
2010
-
[50]
Y. Wang, J. Shi, X. Yin, X. Gu, T. F. Chan, S.-T. Yau, A. W. Toga, and P. M. Thompson. Brain surface conformal parameterization with the Ricci flow.IEEE transactions on medical imaging, 31(2):251–264, 2011. 25
2011
-
[51]
Ye and J
C. Ye and J. Cui. Convergence of Dziuk’s fully discrete linearly implicit scheme for curve shortening flow. SIAM Journal on Numerical Analysis, 59(6):2823–2842, 2021
2021
-
[52]
W. Zeng, X. Yin, Y. Zeng, Y. Lai, X. Gu, and D. Samaras. 3D face matching and registration based on hyperbolic Ricci flow. In2008 IEEE Computer Society Conference on Computer Vision and Pattern Recognition Workshops, pages 1–8. IEEE, 2008. 26 Supplementary material This supple...
2008
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.