A convergence rate for the entropic JKO scheme

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Baradat, Aymeric, Cherf, Sofiane
Formato: Preprint
Publicado: 2026
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866910116685545472
author Baradat, Aymeric
Cherf, Sofiane
author_facet Baradat, Aymeric
Cherf, Sofiane
contents The so-called JKO scheme, named after Jordan, Kinderlehrer and Otto, provides a variational way to construct discrete time approximations of certain partial differential equations (PDEs) appearing as gradient flows in the space of probability measures equipped with the Wasserstein metric. The method consists of an implicit Euler scheme, which can be implemented numerically. Yet, in practice, evaluating the Wasserstein distance can be numerically expensive. To address this problem, a common strategy introduced by Peyré in 2015 and which has been shown to produce faster computations, is to replace the Wasserstein distance with its entropic regularization, also known as the Schrödinger cost. In 2026, the first author, Hraivoronska and Santambrogio, proved that if the regularization parameter $\varepsilon$ is proportional to the time step $τ$, that is, $\varepsilon = ατ$ for some $α> 0$, then as $τ\to 0$, this change results in adding to the limiting PDE the additional linear diffusion term $\fracα{2} Δρ$. Our goal in this article is to provide a convergence rate under convexity assumptions between the entropic JKO scheme and the solution of the initial PDE as both $α$ and $τ$ tend to zero. This will appear as a consequence of a new bound between the classical and entropic JKO schemes.
format Preprint
id arxiv_https___arxiv_org_abs_2604_08283
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle A convergence rate for the entropic JKO scheme
Baradat, Aymeric
Cherf, Sofiane
Analysis of PDEs
Numerical Analysis
49Q22
The so-called JKO scheme, named after Jordan, Kinderlehrer and Otto, provides a variational way to construct discrete time approximations of certain partial differential equations (PDEs) appearing as gradient flows in the space of probability measures equipped with the Wasserstein metric. The method consists of an implicit Euler scheme, which can be implemented numerically. Yet, in practice, evaluating the Wasserstein distance can be numerically expensive. To address this problem, a common strategy introduced by Peyré in 2015 and which has been shown to produce faster computations, is to replace the Wasserstein distance with its entropic regularization, also known as the Schrödinger cost. In 2026, the first author, Hraivoronska and Santambrogio, proved that if the regularization parameter $\varepsilon$ is proportional to the time step $τ$, that is, $\varepsilon = ατ$ for some $α> 0$, then as $τ\to 0$, this change results in adding to the limiting PDE the additional linear diffusion term $\fracα{2} Δρ$. Our goal in this article is to provide a convergence rate under convexity assumptions between the entropic JKO scheme and the solution of the initial PDE as both $α$ and $τ$ tend to zero. This will appear as a consequence of a new bound between the classical and entropic JKO schemes.
title A convergence rate for the entropic JKO scheme
topic Analysis of PDEs
Numerical Analysis
49Q22
url https://arxiv.org/abs/2604.08283