The family of ellipsoids we will consider is defined by Bonnard, Caillau and Rifford have shown that for smaller than , there exists a point on the equator whose injectivity domain is non-convex. This fact implies the existence of two absolutely continuous probability measure with densities in such that the optimal transport plan between and is discontinuous.
To obtain a qualitative description of this phenomena, we introduce two non-vanishing densities with equal total mass on , whose formula in Cartesian coordinates is given by:
with . Previous parameters have been selected to define non-vanishing densities which are not too concentrated. Even in this regular situation, we will see that optimal transport maps may have very different qualitative behaviors depending on the eccentricity of the ellipsoid.
We construct a quantization by a sum of 10000 unit Dirac masses of the two measures of densities and with respect to the two dimensional Hausdorff measure on $S_1$ using Lloyd's algorithm. Our goal is to compare here the optimal transport between the two discrete measures on the sphere and between two analogous measures defined on the oblate ellipsoid , with . By analogous we mean that we deduce our two measures on by the naive bi-Lipschitz transformation
To illustrate the difference between the two optimal transports, we computed for every sampling point of the source densities the geodesic distance between its image and the center of the Gaussian part of the source density, i.e. the point . More formally, this is equivalent to compute the pullback of the geodesic distance function to by the optimal transport. These functions are rescaled to have the same infinity norm, and their values are reported in the following figures, with the same color map. One can observe that the function is smooth along the equator in the spherical case ( ) whereas in the oblate case ( ) some points close to the equator can be transported to points whose geodesic distance to is very different. This qualitative behavior, observed here at a discrete level, is the first numerical illustration of the discontinuity of optimal transport plans between two non-vanishing densities supported on a convex manifold.