Conditional Stability of the Euler Method on Riemannian Manifolds
Fuente:
arXiv
Salvato in:
| Autori principali: | , , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2025
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866912889166626816 |
|---|---|
| author | Ghirardelli, Marta Owren, Brynjulf Celledoni, Elena |
| author_facet | Ghirardelli, Marta Owren, Brynjulf Celledoni, Elena |
| contents | We derive nonlinear stability results for numerical integrators on Riemannian manifolds, by imposing conditions on the ODE vector field and the step size that makes the numerical solution non-expansive whenever the exact solution is non-expansive over the same time step. Our model case is a geodesic version of the explicit Euler method. Precise bounds are obtained in the case of Riemannian manifolds of constant sectional curvature. The approach is based on a cocoercivity property of the vector field adapted to manifolds from Euclidean space. It allows us to compare the new results to the corresponding well-known results in flat spaces, and in general we find that a non-zero curvature will deteriorate the stability region of the geodesic Euler method. The step size bounds depend on the distance traveled over a step from the initial point. Numerical examples for spheres and hyperbolic 2-space confirm that the bounds are tight. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2503_09434 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Conditional Stability of the Euler Method on Riemannian Manifolds Ghirardelli, Marta Owren, Brynjulf Celledoni, Elena Numerical Analysis Mathematical Physics We derive nonlinear stability results for numerical integrators on Riemannian manifolds, by imposing conditions on the ODE vector field and the step size that makes the numerical solution non-expansive whenever the exact solution is non-expansive over the same time step. Our model case is a geodesic version of the explicit Euler method. Precise bounds are obtained in the case of Riemannian manifolds of constant sectional curvature. The approach is based on a cocoercivity property of the vector field adapted to manifolds from Euclidean space. It allows us to compare the new results to the corresponding well-known results in flat spaces, and in general we find that a non-zero curvature will deteriorate the stability region of the geodesic Euler method. The step size bounds depend on the distance traveled over a step from the initial point. Numerical examples for spheres and hyperbolic 2-space confirm that the bounds are tight. |
| title | Conditional Stability of the Euler Method on Riemannian Manifolds |
| topic | Numerical Analysis Mathematical Physics |
| url | https://arxiv.org/abs/2503.09434 |