Benchmarking Numerical Algorithms for Harmonic Maps into the Sphere
Research output: Contribution to journal › Research article › Contributed › peer-review
Contributors
Abstract
We numerically benchmark methods for computing harmonic maps into the unit sphere, with particular focus on harmonic maps with singularities. For the discretization we compare two different approaches, both based on Lagrange finite elements. While the first method enforces the unit-length constraint only at the Lagrange nodes, the other one adds a pointwise projection to fulfill the constraint everywhere. For the solution of the resulting algebraic problems we compare a nonconforming gradient flow with a Riemannian trust-region method. Both are energy-decreasing and can be shown to converge globally to stationary points of the discretized Dirichlet energy. We observe that while the nonconforming and the conforming discretizations both show similar behavior for smooth problems, the nonconforming discretization handles singularities better. On the solver side, the second-order trust-region method converges after few steps, whereas the number of gradient-flow steps increases proportionally to the inverse grid element diameter.
Details
Original language | English |
---|---|
Number of pages | 29 |
Journal | Journal of Numerical Mathematics |
Publication status | Accepted/In press - 7 Feb 2025 |
Peer-reviewed | Yes |
External IDs
ArXiv | 2209.13665 |
---|---|
ORCID | /0000-0003-1093-6374/work/142660180 |