Pith. sign in

REVIEW 4 major objections 4 minor 66 references

GEORCE: A Fast New Control Algorithm for Computing Geodesics

T0 review · 4 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read The paper claims that GEORCE, a control-theoretic reformulation of the discretized geodesic energy, converges globally and locally quadratically while using no Hessians, and reports much faster and more accurate geodesics than ADAM, BFGS…

desk verdict Novel fixed-metric control solver for geodesics with strong empirical results, but the quadratic local convergence proof has a real gap and should not be taken as established. read the letter →

arxiv 2505.05961 v3 pith:SWYBQYKQ submitted 2025-05-09 math.DG stat.CO

classification math.DGstat.CO MSC 53C2249Q9965D15
keywords geodesicsRiemannianmanifoldsFinsleroptimalcontroldiscreteenergyminimizationglobalconvergencequadraticinformationgeometry
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

GEORCE is a new algorithm for computing geodesics—locally shortest paths—on Riemannian and Finsler manifolds. The paper's central claim is that discretizing the geodesic energy as a discrete optimal control problem makes each iteration an explicit linear solve, because the cost becomes convex in the velocity variables once the metric matrices are frozen. The authors prove global convergence to a local minimum and local quadratic convergence near it, matching Newton's rate without forming a Hessian. They report that GEORCE is much faster and finds shorter curves than ADAM, BFGS, and sparse Newton methods on spheres, ellipsoids, hyperbolic spaces, statistical manifolds, and learned generative-model manifolds.

What carries the argument

The central mechanism is the discrete control reformulation: the curve is represented by states $x_t$ and velocity controls $u_t$ linked by $x_{t+1}=x_t+u_t$, with endpoints fixed. Each GEORCE iteration freezes the positive-definite metric matrices $G_t=G(x_t^{(i)})$ and the derivative vectors at the current iterate, so the Pontryagin necessary conditions become a linear system whose closed-form solution (Proposition 2) is a sum of inverses $G_t^{-1}$. The soft line search (Armijo backtracking with decay $\rho=0.5$) converts each descent direction into a globally convergent step. The cost per iteration is $O(Td^3)$, linear in the number of grid points $T$ and cubic in the manifold dimension $d$, and only first-order derivatives of the metric are needed.

What would settle it

Run GEORCE with an exact line search on a strongly convex metric near a known geodesic—for example, a small perturbation of a great-circle geodesic on $S^2$—and record the accepted steps $\alpha^*$ as the iterates approach the optimum. If $\alpha^*$ is not 1 in an entire neighborhood of the optimum, or if the ratio $\|z^{(i+1)}-z^*\|/\|z^{(i)}-z^*\|^2$ fails to stay bounded, then the assumptions of Proposition 8 do not hold and the quadratic-rate claim would not hold on that example.

Watch

Extended reading notes

Core claim

The discovery is that the geodesic boundary-value problem, discretized as minimising $\sum_{t=0}^{T-1} (x_{t+1}-x_t)^\top G(x_t)(x_{t+1}-x_t)$, can be rewritten as a discrete control problem in which the objective is strictly convex in the controls once the metric is frozen. The time-discrete Pontryagin maximum principle yields first-order necessary conditions, and fixing $G_t=G(x_t^{(i)})$ and the metric-derivative terms at the current iterate turns those conditions into a closed-form update built from a sum of inverse metric matrices. The resulting scheme, GEORCE, is proven to converge from any feasible straight-line initialization to a local minimum and, under local strict convexity and acceptance of the full step, to converge quadratically near the minimum. The same construction, with one additional velocity-derivative term $\zeta_t$, extends to Finsler manifolds.

Load-bearing premise

The load-bearing premise is that near the optimum the line search always accepts the full GEORCE step ($\alpha^* = 1$); the paper does not prove this, and converting the proven quadratic decrease of the energy into quadratic convergence of the iterates needs an unstated quadratic-growth condition on the energy.

Editorial extensions

If this is right

  • A geodesic boundary-value problem can be solved by a sequence of closed-form linear solves, each costing $O(Td^3)$, with no Christoffel symbols, Hessians, or second-order metric information.
  • Because the same update with one extra term handles Finsler metrics, asymmetric shortest-path problems such as navigation under a force field inherit the same convergence guarantees and practical speed.
  • On the benchmark manifolds, GEORCE reaches the gradient tolerance in a handful of iterations and returns shorter computed lengths than ADAM, BFGS, and sparse Newton at both $T=50$ and $T=100$.
  • On Fisher-Rao statistical manifolds and on variational-autoencoder latent spaces, the few-iteration behavior makes geodesic computation practical at higher dimension and grid size.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural reading is that GEORCE is a metric-aware quasi-Newton method: each iteration solves the exact quadratic model obtained by freezing $G$, which explains Newton-like local speed when the frozen model approximates the true energy well.
  • Replacing the exact inverse metric sum with a cheaper surrogate (low-rank, diagonal-plus-sampling, or iterative solve) could reduce the $O(d^3)$ per-iteration cost; the paper names inverse evaluation as the main limitation but does not pursue such surrogates.
  • Because the first control $u_0$ estimates the logarithmic map, GEORCE could plug directly into downstream tasks such as Fréchet means, geodesic regression, or parallel transport without extra machinery.
  • The Randers-type benchmarks suggest an immediate test in time-optimal path planning with asymmetric costs, where reversing the force field should reverse the optimal geodesic and expose the Finsler distance asymmetry.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The paper introduces GEORCE, an algorithm for computing geodesics on Riemannian and Finsler manifolds by rewriting the discretized geodesic energy (Eq. 5) as a discrete optimal control problem (Eq. 6), freezing the metric and its derivative at the current iterate, solving the resulting linear system explicitly (Eq. 11), and combining this with a line search. The authors claim global convergence and local quadratic convergence for the discretized energy, an extension to Finsler manifolds, and substantial empirical speedups and lower final energies relative to ADAM, BFGS, and sparse Newton methods across spheres, ellipsoids, tori, hyperbolic space, SPD matrices, Fisher-Rao statistical manifolds, and VAE-learned manifolds. The paper includes pseudocode for the Riemannian and Finsler versions, derivations for a sparse Newton baseline, extensive tables for T=50 and T=100, and a public code repository.

Significance. If the convergence claims are made rigorous, GEORCE would be a genuinely useful contribution: it is a first-order method with an O(Td^3) per-iteration cost that scales linearly in the number of grid points, and the benchmarks cover a broad and relevant set of manifolds, including information-geometric and learned-manifold settings. The paper is also commendable for shipping open-source code, using no fitted parameters in the algorithm itself, and comparing against a wide range of optimization baselines on CPU and GPU. However, the theoretical results currently outrun the proofs: the quadratic-convergence theorem has a substantive logical gap, the Finsler proof is deferred to a sketch, and the global-convergence proof conflates convergence of the energy with convergence of the iterates. The empirical evidence is suggestive but cannot substitute for those proofs, so the central theoretical claims are not yet established as stated.

major comments (4)
  1. [Section 3.2, Proposition 8 and Eqs. (15)-(16)] The proof of local quadratic convergence is not complete. Eq. (15) and the surrounding argument establish, at most, a bound on the energy difference E(z^{(i+1)}) - E(z*) of order O(||z^{(i)} - z*||^2). To convert such an energy bound into the iterate bound ||z^{(i+1)} - z*|| <= c||z^{(i)} - z*||^2 in Eq. (16), one needs a quadratic lower bound E(z) - E(z*) >= m||z - z*||^2 near z*. Assumption 6 only postulates local strict convexity, which does not imply quadratic growth; for example, f(z) = z^4 is locally strictly convex at z* = 0 but has no quadratic lower bound there. The final step of the proof uses Lipschitz continuity to write E(z^{(i+1)}) - E(z*) = O(||z^{(i+1)} - z*||), which is an upper bound on the energy gap, not the required lower comparison. Without an additional strong-convexity or quadratic-growth hypothesis, the claimed Q-quadratic convergence does not follow.
  2. [Section 3.2, Proposition 8 and Algorithm 1] The proposition assumes "locally alpha* = 1, i.e. no line search," but this is not proved for either the exact line search or the Armijo/backtracking procedure described in Algorithm 1. If, in a neighborhood of z*, the accepted step size is bounded away from 1, the update is a damped iteration and one would generally expect only linear convergence. The paragraph after the proof gives a heuristic justification, not a mathematical argument. This hypothesis is load-bearing for the quadratic-convergence claim and must either be proved or be made an explicit assumption with the theorem restated accordingly.
  3. [Section 3.2, Proposition 5] The global-convergence proof shows that the scalar energy sequence E^{(i)} is decreasing and bounded below, hence convergent, but it does not establish convergence of the iterates z^{(i)}. The sentence "denote the convergence point (x-hat, u-hat)" presupposes the desired conclusion. Without compactness of the feasible set or a separate argument that every cluster point is stationary and that the energy limit is attained, an energy-convergent sequence need not have a well-defined limit point. The subsequent paragraph about jumping between two equal-energy local minima is not a substitute for a standard global-convergence argument. Thus the claim that GEORCE "will converge to a (local) minimum" is underproved as written.
  4. [Appendix A.2] The Finsler extension is deferred with the statement that the proof is "completely similar" to the Riemannian case. This is not sufficient because the Finsler version inherits both gaps identified above: the unproved assumption that alpha* = 1 locally and the missing quadratic-growth condition. In addition, because the fundamental tensor G(x,u) now depends on the control variable, the Taylor expansion and the identity in Eq. (A5) need to be verified carefully with the velocity-dependent metric; the sketch does not carry the proof. The abstract's claim of quadratic local convergence for Finsler manifolds is therefore not currently supported.
minor comments (4)
  1. [Appendix B.2, Algorithms 3 and 4] Both algorithms contain malformed conditional statements ("if then ...") and the switching condition between the gradient step and the Newton step is not fully specified; the pseudocode should be corrected and clarified.
  2. [Section 2, paragraph on Levi-Civita connection] The text states that existence and uniqueness of the Levi-Civita connection are "guaranteed by completeness of the metric." The Levi-Civita connection exists for every smooth Riemannian metric without a completeness assumption; this sentence should be corrected.
  3. [Section 5 and Tables 2-7] The empirical accuracy metric is the same discrete energy functional in Eq. (5) that all compared algorithms minimize. Phrases such as "more accurate" should be framed as "lower discrete energy and shorter computed length" rather than accuracy against an independent ground truth, unless a reference solution is available for the specific manifold.
  4. [Figure 2] The caption and text mention dimensions n = 2,3,5,10,20, but the horizontal axis of the figure is labeled only up to 10; please align the description with the plotted dimensions.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: GEORCE's derivation is self-contained, and its convergence claims are conditional theorems rather than reductions to fitted inputs or self-citations.

full rationale

The derivation of GEORCE is a genuine construction: the discretized energy is reformulated as a control problem, the optimality conditions are stated, and an iterative scheme is obtained by fixing the metric and its derivative. No parameter is fitted to data and then renamed as a prediction. Proposition 5 is a standard descent-with-line-search global convergence argument, and Proposition 8 assumes local strict convexity and alpha* = 1 and derives a quadratic rate; its conclusion is not identical to its hypotheses. The Finsler extension is deferred as 'completely similar', which is a proof gap rather than a circular reduction. The empirical benchmarks compare methods on the same discretized energy objective, which is standard optimization benchmarking and does not make the comparison circular; no fitted value is relabeled as an independent ground-truth result. Self-citations are used for background and for learned-manifold context, but they are not load-bearing for the main convergence or derivation claims. Therefore no specific circular step can be exhibited from the paper's equations or citations.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The algorithm introduces no fitted parameters or invented entities. It rests on standard differential geometry assumptions plus a specific regularity condition for the quadratic convergence proof.

assumptions (5)
  • domain assumption The manifold is given in a local chart and the entire candidate curve remains within that chart.
    Stated in Section 3: 'We assume access to a local chart of the manifold and assume that the entire candidate curve of the geodesic can be represented within the same chart.' Without this, the linear-system update may not be defined.
  • domain assumption End points are not conjugate points and the discretization is fine enough that critical points of the discrete energy are local or global minima.
    Stated in Section 3 bullet list; needed for the limit of GEORCE to be a true local minimum of the energy.
  • standard math The metric matrix G(x_t) is symmetric positive definite at every point used.
    Definition of a Riemannian metric; used to invert G_t in Proposition 2.
  • ad hoc to paper Close to the solution, the line search step alpha* equals 1 and E is locally strictly convex.
    Assumption 6 and Proposition 8; this is load-bearing for the quadratic convergence claim and is not proven for the line search.
  • standard math The discretized energy functional is smooth, so first-order Taylor expansion yields O(||Delta z||^2) terms.
    Used throughout the convergence proofs; smoothness of the metric implies smoothness of E.

how reviews work

0 comments
Cite this review

Pith. "Pith review of GEORCE: A Fast New Control Algorithm for Computing Geodesics." pith.science (2026). https://pith.science/paper/SWYBQYKQ

@misc{pith2026250505961,
  author       = {Pith},
  title        = {Pith review of: GEORCE: A Fast New Control Algorithm for Computing Geodesics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SWYBQYKQ}},
  note         = {Machine review of arXiv:2505.05961}
}
read the original abstract

Computing geodesics for Riemannian manifolds is a difficult task that often relies on numerical approximations. However, these approximations tend to be either numerically unstable, have slow convergence, or scale poorly with manifold dimension and number of grid points. We introduce a new algorithm called GEORCE that computes geodesics in a local chart via a transformation into a discrete control problem. We show that GEORCE has global convergence and quadratic local convergence. In addition, we show that it extends to Finsler manifolds. For both Finslerian and Riemannian manifolds, we thoroughly benchmark GEORCE against several alternative optimization algorithms and show empirically that it has a much faster and more accurate performance for a variety of manifolds, including key manifolds from information theory and manifolds that are learned using generative models.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

66 extracted references · 52 canonical work pages

  1. [1]

    & Fletcher, T.Riemannian Geometric Statistics in Medical Image Analysis(Academic Press, 2019)

    Pennec, X., Sommer, S. & Fletcher, T.Riemannian Geometric Statistics in Medical Image Analysis(Academic Press, 2019)

  2. [2]

    & Thirion, J.-P.Landmark-based registration using fea- tures identified through differential geometry, 499–513 (Academic Press, Inc., USA, 2000)

    Pennec, X., Ayache, N. & Thirion, J.-P.Landmark-based registration using fea- tures identified through differential geometry, 499–513 (Academic Press, Inc., USA, 2000)

  3. [3]

    & Faugeras, O

    Lenglet, C., Rousson, M., Deriche, R. & Faugeras, O. Statistics on Multi- variate Normal Distributions: A Geometric Approach and its Application to Diffusion Tensor MRI.Rapport de recherche, Institut National de Recherche en Informatique et en Automatique (INRIA)(2004)

  4. [4]

    Watson, Z., Townshend, R. J. L. & Dror, R. O. Protein Design Using Diffusion Models.bioRxiv(2022). URL https://www.biorxiv.org/content/10.1101/2022. 12.09.519842v1

  5. [5]

    & Others

    Deftelsen, N. & Others. Protein Structure Prediction Using Rosetta.Methods in Molecular Biology2021, 1–30 (2019)

  6. [6]

    & Dunbrack, R

    Shapovalov, M. & Dunbrack, R. A Smoothed Backbone-Dependent Rotamer Library for Proteins Derived from Adaptive Kernel Density Estimates and Regressions.Structure (London, England : 1993)19, 844–58 (2011)

  7. [7]

    A.et al.Reactive Motion Generation on Learned Riemannian Man- ifolds.Int

    Hsieh, M. A.et al.Reactive Motion Generation on Learned Riemannian Man- ifolds.Int. J. Rob. Res.42, 729–754 (2023). URL https://doi.org/10.1177/ 02783649231193046

  8. [8]

    Simeonov, A.et al.Neural Descriptor Fields: SE(3)-Equivariant Object Repre- sentations for Manipulation (2022)

Show all 66 references
  1. [9]

    & Ramos, F

    Senanayake, R. & Ramos, F. T. Directional Grid Maps: Modeling Multimodal Angular Uncertainty in Dynamic Environments.2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS)3241–3248 (2018). URL https://api.semanticscholar.org/CorpusID:52158241

  2. [10]

    & Hirche, S

    Feiten, W., Lang, M. & Hirche, S. Rigid Motion Estimation Using Mixtures of Projected Gaussians 1465–1472 (2013)

  3. [11]

    Arvanitidis, G., Hansen, L. K. & Hauberg, S. Latent Space Oddity: on the Cur- vature of Deep Generative Models (2018). URL https://openreview.net/forum? id=SJzRZ-WCZ

  4. [12]

    & Thomas Fletcher, P

    Shao, H., Kumar, A. & Thomas Fletcher, P. The Riemannian Geometry of Deep Generative Models (2018). 34

  5. [13]

    D.et al.Oh, A

    Bortoli, V. D.et al.Oh, A. H., Agarwal, A., Belgrave, D. & Cho, K. (eds) Riemannian Score-Based Generative Modelling. (eds Oh, A. H., Agarwal, A., Belgrave, D. & Cho, K.)Advances in Neural Information Processing Systems (2022). URL https://openreview.net/forum?id=oDRQGo8I7P

  6. [14]

    & Hwang, S

    Jo, J. & Hwang, S. J. Generative Modeling on Manifolds through Mixture of Riemannian Diffusion Processes (2024)

  7. [15]

    & Nickel, M

    Mathieu, E. & Nickel, M. Riemannian Continuous Normalizing Flows (2020)

  8. [16]

    & Lipman, Y

    Rozen, N., Grover, A., Nickel, M. & Lipman, Y. Moser Flow: Divergence-based Generative Modeling on Manifolds (2021)

  9. [17]

    & Hauberg, S

    Arvanitidis, G., Gonz´ alez-Duque, M., Pouplin, A., Kalatzis, D. & Hauberg, S. Camps-Valls, G., Ruiz, F. J. R. & Valera, I. (eds)Pulling back information geometry. (eds Camps-Valls, G., Ruiz, F. J. R. & Valera, I.)Proceedings of The 25th International Conference on Artificial ...

  10. [18]

    An Elementary Introduction to Information Geometry.Entropy22, 1100 (2020)

    Nielsen, F. An Elementary Introduction to Information Geometry.Entropy22, 1100 (2020). URL http://dx.doi.org/10.3390/e22101100

  11. [19]

    & Liu, G

    Cheng, C., Li, J., Peng, J. & Liu, G. Globerson, A.et al.(eds)Categorical Flow Matching on Statistical Manifolds. (eds Globerson, A.et al.)Advances in Neural Information Processing Systems, Vol. 37, 54787–54819 (Curran Asso- ciates, Inc., 2024). URL https://proceedings.neurips...

  12. [20]

    & Dragan, A

    Myers, V., Ellis, E., Levine, S., Eysenbach, B. & Dragan, A. Learning to Assist Humans without Inferring Rewards (2024). URL https://openreview.net/forum? id=WCnJmb7cv1

  13. [21]

    (eds Globerson, A.et al.)Advances in Neural Information Processing Systems, Vol

    Davis, O.et al.Globerson, A.et al.(eds)Fisher Flow Matching for Gen- erative Modeling over Discrete Data. (eds Globerson, A.et al.)Advances in Neural Information Processing Systems, Vol. 37, 139054–139084 (Curran Asso- ciates, Inc., 2024). URL https://proceedings.neurips.cc/pa...

  14. [22]

    Hauberg, S.Differential Geometry for Generative Modeling(2025)

  15. [23]

    URL https://books.google.dk/books?id=uXJQQgAACAAJ

    do Carmo, M.Riemannian GeometryMathematics (Boston, Mass.) (Birkh¨ auser, 1992). URL https://books.google.dk/books?id=uXJQQgAACAAJ

  16. [24]

    S.et al.StochMan.GitHub

    Detlefsen, N. S.et al.StochMan.GitHub. Note: https://github.com/MachineLearningLifeScience/stochman/(2021). 35

  17. [25]

    A Global algorithm for geodesics.Journal of the Australian Mathematical Society

    Noakes, L. A Global algorithm for geodesics.Journal of the Australian Mathematical Society. Series A. Pure Mathematics and Statistics65, 37–50 (1998)

  18. [26]

    & Lawrence, N

    Tosi, A., Hauberg, S., Vellido, A. & Lawrence, N. D. Metrics for probabilistic geometries (2014)

  19. [27]

    C1 Isometric Imbeddings.Annals of Mathematics60, 383–396 (1954)

    Nash, J. C1 Isometric Imbeddings.Annals of Mathematics60, 383–396 (1954). URL http://www.jstor.org/stable/1969840

  20. [28]

    Differentiable Manifolds.Annals of Mathematics37, 645–680 (1936)

    Whitney, H. Differentiable Manifolds.Annals of Mathematics37, 645–680 (1936). URL http://www.jstor.org/stable/1968482

  21. [29]

    URL https://books.google.dk/books? id=4Mh6zgEACAAJ

    Ohta, S.Comparison Finsler GeometrySpringer Monographs in Mathematics (Springer International Publishing, 2021). URL https://books.google.dk/books? id=4Mh6zgEACAAJ

  22. [30]

    K., Meneghetti, F

    Miyamoto, H. K., Meneghetti, F. C. C., Pinele, J. & Costa, S. I. R. On Closed- Form Expressions for the Fisher-Rao Distance.Information Geometry7(2024)

  23. [31]

    & Rangarajan, A

    Peter, A. & Rangarajan, A. Shape analysis using the Fisher-Rao Riemannian metric: Unifying shape representation and deformation.In: IEEE ISBI. IEEE 2006, 1164–1167 (2006)

  24. [32]

    Broyden, C. G. The Convergence of a Class of Double-rank Minimization Algo- rithms 1. General Considerations.IMA Journal of Applied Mathematics6, 76–90 (1970). URL https://doi.org/10.1093/imamat/6.1.76

  25. [33]

    A new approach to variable metric algorithms.The Computer Journal13, 317–322 (1970)

    Fletcher, R. A new approach to variable metric algorithms.The Computer Journal13, 317–322 (1970). URL https://doi.org/10.1093/comjnl/13.3.317

  26. [34]

    URL https://api.semanticscholar.org/CorpusID:790344

    A family of variable-metric methods derived by variational means, author=Donald Goldfarb.Mathematics of Computation24, 23–26 (1970). URL https://api.semanticscholar.org/CorpusID:790344

  27. [35]

    & Kettler, P

    Shanno, D. & Kettler, P. Optimal Conditioning of Quasi-Newton Methods. Mathematics of Computation - Math. Comput.24, 657–657 (1970)

  28. [36]

    & Ye, Y.Linear and Nonlinear ProgrammingInternational Series in Operations Research & Management Science (Springer US, 2008)

    Luenberger, D. & Ye, Y.Linear and Nonlinear ProgrammingInternational Series in Operations Research & Management Science (Springer US, 2008). URL https: //books.google.dk/books?id=-pD62uvi9lgC

  29. [37]

    Kaya, C. Y. & Noakes, J. L. Leapfrog for Optimal Control.SIAM Journal on Numerical Analysis46, 2795–2817 (2008). URL https://doi.org/10.1137/ 060675034

  30. [38]

    & Hauberg, S

    Hennig, P. & Hauberg, S. Kaski, S. & Corander, J. (eds)Probabilistic Solutions to Differential Equations and their Application to Riemannian Statistics. (eds Kaski, 36 S. & Corander, J.)Proceedings of the Seventeenth International Conference on Artificial Intelligence and Stat...

  31. [39]

    & Schober, M

    Arvanitidis, G., Hauberg, S., Hennig, P. & Schober, M. Chaudhuri, K. & Sugiyama, M. (eds)Fast and Robust Shortest Paths on Manifolds Learned from Data. (eds Chaudhuri, K. & Sugiyama, M.)Proceedings of the Twenty-Second International Conference on Artificial Intelligence and St...

  32. [40]

    Kingma, D. P. & Ba, J. Adam: A Method for Stochastic Optimization (2014). URL http://arxiv.org/abs/1412.6980. Cite arxiv:1412.6980Comment: Published as a conference paper at the 3rd International Conference for Learning Representations, San Diego, 2015

  33. [41]

    & Qin, Y

    Crane, K., Livesu, M., Puppo, E. & Qin, Y. A Survey of Algorithms for Geodesic Paths and Distances.CoRRabs/2007.10430(2020). URL https://arxiv.org/ abs/2007.10430

  34. [42]

    Surazhsky, V., Surazhsky, T., Kirsanov, D., Gortler, S. J. & Hoppe, H. Fast exact and approximate geodesics on meshes.ACM Trans. Graph.24, 553–560 (2005). URL https://doi.org/10.1145/1073204.1073228

  35. [43]

    & Han, Y

    Chen, J. & Han, Y. Shortest paths on a polyhedron.Proceedings of the Sixth Annual Symposium on Computational Geometry360–369 (1990). URL https: //doi.org/10.1145/98524.98601

  36. [44]

    & Wang, G.-J

    Xin, S.-Q. & Wang, G.-J. Improving Chen and Han’s Algorithm on the Discrete Geodesic Problem.ACM Trans. Graph.28(2009)

  37. [45]

    & Zhang, J

    Qin, Y., Han, X., Yu, H., Yu, Y. & Zhang, J. Fast and exact discrete geodesic com- putation based on triangle-oriented wavefront propagation.ACM Trans. Graph. 35(2016). URL https://doi.org/10.1145/2897824.2925930

  38. [46]

    & Rozo, L

    Beik-Mohammadi, H., Hauberg, S., Arvanitidis, G., Neumann, G. & Rozo, L. Learning Riemannian Manifolds for Geodesic Motion Skills (2021)

  39. [47]

    URL https://books.google.dk/ books?id=fCh2SAtWIdwC

    Kirk, D.Optimal Control Theory: An IntroductionDover Books on Electrical Engineering Series (Dover Publications, 2004). URL https://books.google.dk/ books?id=fCh2SAtWIdwC

  40. [48]

    Boyd, S. P. & Vandenberghe, L.Convex Optimization(Cambridge University Press, Cambridge, UK, 2004). 37

  41. [49]

    Convergence Conditions for Ascent Methods.SIAM Review11, 226–235 (1969)

    Wolfe, P. Convergence Conditions for Ascent Methods.SIAM Review11, 226–235 (1969). URL https://doi.org/10.1137/1011036

  42. [50]

    Minimization of functions having Lipschitz continuous first partial derivatives.Pacific Journal of Mathematics16, 1–3 (1966)

    Armijo, L. Minimization of functions having Lipschitz continuous first partial derivatives.Pacific Journal of Mathematics16, 1–3 (1966)

  43. [51]

    & Golestanian, R

    Piro, L., Tang, E. & Golestanian, R. Optimal navigation strategies for microswim- mers on curved manifolds.Physical Review Research3(2021). URL http: //dx.doi.org/10.1103/PhysRevResearch.3.023125

  44. [52]

    Kingma, D. P. & Welling, M. Auto-Encoding Variational Bayes (2014). http://arxiv.org/abs/1312.6114v10

  45. [53]

    & Tang, X

    Liu, Z., Luo, P., Wang, X. & Tang, X. Deep Learning Face Attributes in the Wild. Proceedings of International Conference on Computer Vision (ICCV)(2015)

  46. [54]

    Elliptic Problems in Linear Difference Equations over a Network (1949)

    Thomas, L. Elliptic Problems in Linear Difference Equations over a Network (1949)

  47. [55]

    & Zhang, E

    Noakes, L. & Zhang, E. Finding geodesics joining given points.Advances in Computational Mathematics48(2022)

  48. [56]

    The mnist database of handwritten digit images for machine learning research.IEEE Signal Processing Magazine29, 141–142 (2012)

    Deng, L. The mnist database of handwritten digit images for machine learning research.IEEE Signal Processing Magazine29, 141–142 (2012)

  49. [57]

    URL http://ufldl.stanford.edu/housenumbers/nips2011 housenumbers.pdf

    Netzer, Y.et al.Reading Digits in Natural Images with Unsupervised Fea- ture Learning.NIPS Workshop on Deep Learning and Unsupervised Feature Learning 2011(2011). URL http://ufldl.stanford.edu/housenumbers/nips2011 housenumbers.pdf

  50. [58]

    An Overview of Gradient Descent OptimizationAlgorithms (2017)

    Ruder, S. An Overview of Gradient Descent OptimizationAlgorithms (2017). URL https://arxiv.org/abs/1609.04747. arXiv:1609.04747

  51. [59]

    & Wright, S.Numerical OptimizationSpringer Series in Operations Research and Financial Engineering (Springer New York, 2006)

    Nocedal, J. & Wright, S.Numerical OptimizationSpringer Series in Operations Research and Financial Engineering (Springer New York, 2006). URL https: //books.google.dk/books?id=VbHYoSyelFcC

  52. [60]

    Virtanen, P.et al.SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python.Nature Methods17, 261–272 (2020)

  53. [61]

    & Prince, P

    Dormand, J. & Prince, P. A Family of Embedded Runge-Kutta Formulae.Journal of Computational and Applied Mathematics6, 19–26 (1980). URL https://www. sciencedirect.com/science/article/pii/0771050X80900133

  54. [62]

    & Shampine, L

    Bogacki, P. & Shampine, L. A 3(2) Pair of Runge - Kutta Formulas.Applied Mathematics Letters2, 321–325 (1989). URL https://www.sciencedirect.com/ science/article/pii/0893965989900797. 38

  55. [63]

    & Wanner, G.Solving Ordinary Differential Equations I: Nonstiff ProblemsVol

    Hairer, E., Norsett, S. & Wanner, G.Solving Ordinary Differential Equations I: Nonstiff ProblemsVol. 8 (1993)

  56. [64]

    & Wanner, G.Solving Ordinary Differential Equations II

    Hairer, E. & Wanner, G.Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic ProblemsVol. 14 (1996)

  57. [65]

    Shampine, L. F. & Reichelt, M. W. The MATLAB ODE Suite.SIAM Jour- nal on Scientific Computing18, 1–22 (1997). URL https://doi.org/10.1137/ S1064827594276424

  58. [66]

    Hindmarsh, A. C. Stepleman, R. S. (ed.)ODEPACK, a Systematized Collection of ODE Solvers. (ed.Stepleman, R. S.)Scientific Computing, 55–64 (North-Holland, Amsterdam, 1983). 39 Appendix A Proofs and derivations A.1 Proof of update scheme for Finslerian geometry Consider the fol...

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.