Wasserstein Mirror Gradient Flow as the limit of the Sinkhorn Algorithm

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Deb, Nabarun, Kim, Young-Heon, Pal, Soumik, Schiebinger, Geoffrey
Formato: Preprint
Publicado: 2023
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866914486843080704
author Deb, Nabarun
Kim, Young-Heon
Pal, Soumik
Schiebinger, Geoffrey
author_facet Deb, Nabarun
Kim, Young-Heon
Pal, Soumik
Schiebinger, Geoffrey
contents We prove that the sequence of marginals obtained from the iterations of the Sinkhorn algorithm or the iterative proportional fitting procedure (IPFP) on joint densities, converges to an absolutely continuous curve on the $2$-Wasserstein space, as the regularization parameter $\varepsilon$ goes to zero and the number of iterations is scaled as $1/\varepsilon$ (and other technical assumptions). This limit, which we call the Sinkhorn flow, is an example of a Wasserstein mirror gradient flow, a concept we introduce here inspired by the well-known Euclidean mirror gradient flows. In the case of Sinkhorn, the gradient is that of the relative entropy functional with respect to one of the marginals and the mirror is half of the squared Wasserstein distance functional from the other marginal. Interestingly, the norm of the velocity field of this flow can be interpreted as the metric derivative with respect to the linearized optimal transport (LOT) distance. An equivalent description of this flow is provided by the parabolic Monge-Ampère PDE whose connection to the Sinkhorn algorithm was noticed by Berman (2020). We derive conditions for exponential convergence for this limiting flow. We also construct a Mckean-Vlasov diffusion whose marginal distributions follow the Sinkhorn flow.
format Preprint
id arxiv_https___arxiv_org_abs_2307_16421
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Wasserstein Mirror Gradient Flow as the limit of the Sinkhorn Algorithm
Deb, Nabarun
Kim, Young-Heon
Pal, Soumik
Schiebinger, Geoffrey
Probability
Analysis of PDEs
Machine Learning
49N99, 49Q22, 60J60
We prove that the sequence of marginals obtained from the iterations of the Sinkhorn algorithm or the iterative proportional fitting procedure (IPFP) on joint densities, converges to an absolutely continuous curve on the $2$-Wasserstein space, as the regularization parameter $\varepsilon$ goes to zero and the number of iterations is scaled as $1/\varepsilon$ (and other technical assumptions). This limit, which we call the Sinkhorn flow, is an example of a Wasserstein mirror gradient flow, a concept we introduce here inspired by the well-known Euclidean mirror gradient flows. In the case of Sinkhorn, the gradient is that of the relative entropy functional with respect to one of the marginals and the mirror is half of the squared Wasserstein distance functional from the other marginal. Interestingly, the norm of the velocity field of this flow can be interpreted as the metric derivative with respect to the linearized optimal transport (LOT) distance. An equivalent description of this flow is provided by the parabolic Monge-Ampère PDE whose connection to the Sinkhorn algorithm was noticed by Berman (2020). We derive conditions for exponential convergence for this limiting flow. We also construct a Mckean-Vlasov diffusion whose marginal distributions follow the Sinkhorn flow.
title Wasserstein Mirror Gradient Flow as the limit of the Sinkhorn Algorithm
topic Probability
Analysis of PDEs
Machine Learning
49N99, 49Q22, 60J60
url https://arxiv.org/abs/2307.16421