Benchmarking Numerical Algorithms for Harmonic Maps into the Sphere

Research output: Preprint/documentation/reportPreprint

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 a stationary point of the Dirichlet energy. We observe that while the nonconforming and the conforming discretizations both show similar behavior, the second-order trust-region method needs less iterations than the solver based on gradient flow.

Details

Original languageEnglish
PublisherarXiv
Number of pages22
Volume2209.13665
Publication statusPublished - 27 Sept 2022
No renderer: customAssociatesEventsRenderPortal,dk.atira.pure.api.shared.model.researchoutput.WorkingPaper

External IDs

ArXiv 2209.13665
ORCID /0000-0003-1093-6374/work/142660180