English

A Convergent Finite Difference Method for Optimal Transport on the Sphere

Numerical Analysis 2021-05-11 v1 Numerical Analysis

Abstract

We introduce a convergent finite difference method for solving the optimal transportation problem on the sphere. The method applies to both the traditional squared geodesic cost (arising in mesh generation) and a logarithmic cost (arising in the reflector antenna design problem). At each point on the sphere, we replace the surface PDE with a Generated Jacobian equation posed on the local tangent plane using geodesic normal coordinates. The discretization is inspired by recent monotone methods for the Monge-Amp\`ere equation, but requires significant adaptations in order to correctly handle the mix of gradient and Hessian terms appearing inside the nonlinear determinant operator, as well as the singular logarithmic cost function. Numerical results demonstrate the success of this method on a wide range of challenging problems involving both the squared geodesic and the logarithmic cost functions.

Keywords

Cite

@article{arxiv.2105.03500,
  title  = {A Convergent Finite Difference Method for Optimal Transport on the Sphere},
  author = {Brittany Froese Hamfeldt and Axel G. R. Turnquist},
  journal= {arXiv preprint arXiv:2105.03500},
  year   = {2021}
}

Comments

34 pages, 21 figures

R2 v1 2026-06-24T01:53:28.936Z