A convergence framework for optimal transport on the sphere
share
We consider a PDE approach to numerically solving the optimal transportation problem on the sphere. We focus on both the traditional squared geodesic cost and a logarithmic cost, which arises in the reflector antenna design problem. At each point on the sphere, we replace the surface PDE with a generalized Monge-Ampère type equation posed on the tangent plane using normal coordinates. The resulting nonlinear PDE can then be approximated by any consistent, monotone scheme for generalized Monge-Ampère type equations on the plane. By augmenting this discretization with an additional term that constrains the solution gradient, we obtain a strong form of stability. A modification of the Barles-Souganidis convergence framework then establishes convergence to the mean-zero solution of the original PDE.
READ FULL TEXT