A convergent finite difference method for optimal transport on the sphere

Document Type

Article

Publication Date

11-15-2021

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ère 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.

Identifier

85112487402 (Scopus)

Publication Title

Journal of Computational Physics

External Full Text Location

https://doi.org/10.1016/j.jcp.2021.110621

e-ISSN

10902716

ISSN

00219991

Volume

445

Grant

DMS-1751996

Fund Ref

National Science Foundation

This document is currently not available here.

Share

COinS