2024Activity reportProject-TeamMOKAPLAN
RNSR: 201321083P- Research center Inria Paris Centre
- In partnership with:Université Paris-Dauphine, CNRS
- Team name: Advances in Numerical Calculus of Variations
- In collaboration with:CEREMADE
- Domain:Applied Mathematics, Computation and Simulation
- Theme:Numerical schemes and simulations
Keywords
Computer Science and Digital Science
- A5.3. Image processing and analysis
- A5.9. Signal processing
- A6.1.1. Continuous Modeling (PDE, ODE)
- A6.2.1. Numerical analysis of PDE and ODE
- A6.2.6. Optimization
- A6.3.1. Inverse problems
- A8.2.3. Calculus of variations
- A8.12. Optimal transport
- A9. Artificial intelligence
Other Research Topics and Application Domains
- B9.5.2. Mathematics
- B9.5.3. Physics
- B9.6.3. Economy, Finance
1 Team members, visitors, external collaborators
Research Scientists
- Vincent Duval [Team leader, INRIA, Senior Researcher, HDR]
- Jean-David Benamou [INRIA, Senior Researcher, HDR]
- Antonin Chambolle [CNRS, Senior Researcher, HDR]
- Flavien Leger [INRIA, Researcher]
- Irene Waldspurger [CNRS, Researcher, Université Paris Dauphine]
Faculty Members
- Guillaume Carlier [DAUPHINE PSL, Professor, HDR]
- Christian Leonard [UNIV PARIS NANTERRE, Professor Delegation, from Sep 2024, HDR]
- Paul Pegon [DAUPHINE PSL, Associate Professor]
Post-Doctoral Fellows
- Annette Dumas [DAUPHINE PSL, Post-Doctoral Fellow, from Oct 2024]
- Adrien Vacher [INRIA, Post-Doctoral Fellow, until Aug 2024]
PhD Students
- Guillaume Chazareix [BNP PARIBAS, CIFRE]
- Hugo Malamut [DAUPHINE PSL]
- Joao Pinto Anastacio Machado [DAUPHINE PSL, until Sep 2024]
- Faniriana Rakoto Endor [CNRS, Université Paris Dauphine]
- Maxime Sylvestre [DAUPHINE PSL]
- Louis Tocquec [UNIV PARIS SACLAY, from Oct 2024]
Interns and Apprentices
- Majid Arthaud [INRIA, Intern, from May 2024 until Oct 2024]
Administrative Assistant
- Derya Gok [INRIA]
Visiting Scientists
- Luigi De Pascale [UNIV FLORENCE, from May 2024 until Sep 2024]
- Jacob Francis [IMPERIAL COLLEGE LDN, from Mar 2024 until May 2024]
2 Overall objectives
2.1 Introduction
The last decade has witnessed a remarkable convergence between several sub-domains of the calculus of variations, namely optimal transport (and its many generalizations), infinite dimensional geometry of diffeomorphisms groups and inverse problems in imaging (in particular sparsity-based regularization). This convergence is due to (i) the mathematical objects manipulated in these problems, namely sparse measures (e.g. coupling in transport, edge location in imaging, displacement fields for diffeomorphisms) and (ii) the use of similar numerical tools from non-smooth optimization and geometric discretization schemes. Optimal Transportation, diffeomorphisms and sparsity-based methods are powerful modeling tools, that impact a rapidly expanding list of scientific applications and call for efficient numerical strategies. Our research program shows the important part played by the team members in the development of these numerical methods and their application to challenging problems.
2.2 Static Optimal Transport and Generalizations
Optimal Transport, Old and New.
Optimal Mass Transportation is a mathematical research topic which started two centuries ago with Monge's work on the “Théorie des déblais et des remblais" (see 95). This engineering problem consists in minimizing the transport cost between two given mass densities. In the 40's, Kantorovich 104 introduced a powerful linear relaxation and introduced its dual formulation. The Monge-Kantorovich problem became a specialized research topic in optimization and Kantorovich obtained the 1975 Nobel prize in economics for his contributions to resource allocations problems. Since the seminal discoveries of Brenier in the 90's 57, Optimal Transportation has received renewed attention from mathematical analysts and the Fields Medal awarded in 2010 to C. Villani, who gave important contributions to Optimal Transportation and wrote the modern reference monographs 129, 130, arrived at a culminating moment for this theory. Optimal Mass Transportation is today a mature area of mathematical analysis with a constantly growing range of applications. Optimal Transportation has also received a lot of attention from probabilists (see for instance the recent survey 108 for an overview of the Schrödinger problem which is a stochastic variant of the Benamou-Brenier dynamical formulation of optimal transport). The development of numerical methods for Optimal Transportation and Optimal Transportation related problems is a difficult topic and comparatively underdeveloped. This research field has experienced a surge of activity in the last five years, with important contributions of the Mokaplan group (see the list of important publications of the team). We describe below a few of recent and less recent Optimal Transportation concepts and methods which are connected to the future activities of Mokaplan :
Brenier's theorem 60 characterizes the unique optimal map as the gradient of a convex potential. As such Optimal Transportation may be interpreted as an infinite dimensional optimisation problem under “convexity constraint": i.e. the solution of this infinite dimensional optimisation problem is a convex potential. This connects Optimal Transportation to “convexity constrained" non-linear variational problems such as, for instance, Newton's problem of the body of minimal resistance. The value function of the optimal transport problem is also known to define a distance between source and target densities called the Wasserstein distance which plays a key role in many applications such as image processing.
Monge-Ampère Methods.
A formal substitution of the optimal transport map as the gradient of a convex potential in the mass conservation constraint (a Jacobian equation) gives a non-linear Monge-Ampère equation. Caffarelli 64 used this result to extend the regularity theory for the Monge-Ampère equation. In the last ten years, it also motivated new research on numerical solvers for non-linear degenerate Elliptic equations 8711250 49 and the references therein. Geometric approaches based on Laguerre diagrams and discrete data 115 have also been developed. Monge-Ampère based Optimal Transportation solvers have recently given the first linear cost computations of Optimal Transportation (smooth) maps.
Generalizations of OT.
In recent years, the classical Optimal Transportation problem has been extended in several directions. First, different ground costs measuring the “physical" displacement have been considered. In particular, well posedness for a large class of convex and concave costs has been established by McCann and Gangbo 94. Optimal Transportation techniques have been applied for example to a Coulomb ground cost in Quantum chemistry in relation with Density Functional theory 81. Given the densities of electrons Optimal Transportation models the potential energy and their relative positions. For more than more than 2 electrons (and therefore more than 2 densities) the natural extension of Optimal Transportation is the so called Multi-marginal Optimal Transport (see 120 and the references therein). Another instance of multi-marginal Optimal Transportation arises in the so-called Wasserstein barycenter problem between an arbitrary number of densities 34. An interesting overview of this emerging new field of optimal transport and its applications can be found in the recent survey of Ghoussoub and Pass 119.

Example of color transfer between two images, computed using the method developed in 47, see also 124. The image framed in red and blue are the input images. Top and middle row: adjusted image where the color of the transported histogram has been imposed. Bottom row: geodesic (displacement) interpolation between the histogram of the chrominance of the image.
Numerical Applications of Optimal Transportation.
Optimal transport has found many applications, starting from its relation with several physical models such as the semi-geostrophic equations in meteorology 99, 84, 83, 45, 111, mesh adaptation 110, the reconstruction of the early mass distribution of the Universe 92, 58 in Astrophysics, and the numerical optimisation of reflectors following the Optimal Transportation interpretation of Oliker 63 and Wang 131. Extensions of OT such as multi-marginal transport has potential applications in Density Functional Theory , Generalized solution of Euler equations 59 (DFT) and in statistics and finance 43, 93 .... Recently, there has been a spread of interest in applications of OT methods in imaging sciences 53, statistics 51 and machine learning 85. This is largely due to the emergence of fast numerical schemes to approximate the transportation distance and its generalizations, see for instance 47. Figure 1 shows an example of application of OT to color transfer. Figure 2 shows an example of application in computer graphics to interpolate between input shapes.

Example of barycenter between shapes computed using optimal transport barycenters of the uniform densities inside the 3 extremal shapes, computed as detailed in 124. Note that the barycenters are not in general uniform distributions, and we display them as the surface defined by a suitable level-set of the density.
2.3 Diffeomorphisms and Dynamical Transport
Dynamical transport.
While the optimal transport problem, in its original formulation, is a static problem (no time evolution is considered), it makes sense in many applications to rather consider time evolution. This is relevant for instance in applications to fluid dynamics or in medical images to perform registration of organs and model tumor growth.
In this perspective, the optimal transport in Euclidean space corresponds to an evolution where each particule of mass evolves in straight line. This interpretation corresponds to the Computational Fluid Dynamic (CFD) formulation proposed by Brenier and Benamou in 44. These solutions are time curves in the space of densities and geodesics for the Wasserstein distance. The CFD formulation relaxes the non-linear mass conservation constraint into a time dependent continuity equation, the cost function remains convex but is highly non smooth. A remarkable feature of this dynamical formulation is that it can be re-cast as a convex but non smooth optimization problem. This convex dynamical formulation finds many non-trivial extensions and applications, see for instance 46. The CFD formulation also appears to be a limit case of Mean Fields games (MFGs), a large class of economic models introduced by Lasry and Lions 106 leading to a system coupling an Hamilton-Jacobi with a Fokker-Planck equation. In contrast, the Monge case where the ground cost is the euclidan distance leads to a static system of PDEs 54.

Examples of displacement interpolation (geodesic for optimal transport) according to a non-Euclidean Riemannian metric (the mass is constrained to move inside a maze) between to input Gaussian distributions. Note that the maze is dynamic: its topology change over time, the mass being “trapped” at time
Gradient Flows for the Wasserstein Distance.
Another extension is, instead of considering geodesic for transportation metric (i.e. minimizing the Wasserstein distance to a target measure), to make the density evolve in order to minimize some functional. Computing the steepest descent direction with respect to the Wasserstein distance defines a so-called Wasserstein gradient flow, also known as JKO gradient flows after its authors 103. This is a popular tool to study a large class of non-linear diffusion equations. Two interesting examples are the Keller-Segel system for chemotaxis 102, 73 and a model of congested crowd motion proposed by Maury, Santambrogio and Roudneff-Chupin 114. From the numerical point of view, these schemes are understood to be the natural analogue of implicit scheme for linear parabolic equations. The resolution is however costly as it involves taking the derivative in the Wasserstein sense of the relevant energy, which in turn requires the resolution of a large scale convex but non-smooth minimization.
Geodesic on infinite dimensional Riemannian spaces.
To tackle more complicated warping problems, such as those encountered in medical image analysis, one unfortunately has to drop the convexity of the functional involved in defining the gradient flow. This gradient flow can either be understood as defining a geodesic on the (infinite dimensional) group of diffeomorphisms 42, or on a (infinite dimensional) space of curves or surfaces 132. The de-facto standard to define, analyze and compute these geodesics is the “Large Deformation Diffeomorphic Metric Mapping” (LDDMM) framework of Trouvé, Younes, Holm and co-authors 42, 98. While in the CFD formulation of optimal transport, the metric on infinitesimal deformations is just the
2.4 Sparsity in Imaging
Sparse regularization.
Beside image warping and registration in medical image analysis, a key problem in nearly all imaging applications is the reconstruction of high quality data from low resolution observations. This field, commonly referred to as “inverse problems”, is very often concerned with the precise location of features such as point sources (modeled as Dirac masses) or sharp contours of objects (modeled as gradients being Dirac masses along curves). The underlying intuition behind these ideas is the so-called sparsity model (either of the data itself, its gradient, or other more complicated representations such as wavelets, curvelets, bandlets 113 and learned representation 133).
The huge interest in these ideas started mostly from the introduction of convex methods to serve as proxy for these sparse regularizations. The most well known is the
Regularization over measure spaces.
However, the theoretical analysis of sparse reconstructions involving real-life acquisition operators (such as those found in seismic imaging, neuro-imaging, astro-physical imaging, etc.) is still mostly an open problem. A recent research direction, triggered by a paper of Candès and Fernandez-Granda 67, is to study directly the infinite dimensional problem of reconstruction of sparse measures (i.e. sum of Dirac masses) using the total variation of measures (not to be mistaken for the total variation of 2-D functions). Several works 66, 89, 88 have used this framework to provide theoretical performance guarantees by basically studying how the distance between neighboring spikes impacts noise stability.
![]() |
![]() |
![]() |
![]() |
Segmentation input | output | Zooming input | output |
Two example of application of the total variation regularization of functions. Left: image segmentation into homogeneous color regions. Right: image zooming (increasing the number of pixels while keeping the edges sharp).
Low complexity regularization and partial smoothness.
In image processing, one of the most popular methods is the total variation regularization 123, 61. It favors low-complexity images that are piecewise constant, see Figure 4 for some examples on how to solve some image processing problems. Beside applications in image processing, sparsity-related ideas also had a deep impact in statistics 125 and machine learning 38. As a typical example, for applications to recommendation systems, it makes sense to consider sparsity of the singular values of matrices, which can be relaxed using the so-called nuclear norm (a.k.a. trace norm) 37. The underlying methodology is to make use of low-complexity regularization models, which turns out to be equivalent to the use of partly-smooth regularization functionals 109, 128 enforcing the solution to belong to a low-dimensional manifold.
2.5 Mokaplan unified point of view
The dynamical formulation of optimal transport creates a link between optimal transport and geodesics on diffeomorphisms groups. This formal link has at least two strong implications that Mokaplan will elaborate on: (i) the development of novel models that bridge the gap between these two fields ; (ii) the introduction of novel fast numerical solvers based on ideas from both non-smooth optimization techniques and Bregman metrics.
In a similar line of ideas, we believe a unified approach is needed to tackle both sparse regularization in imaging and various generalized OT problems. Both require to solve related non-smooth and large scale optimization problems. Ideas from proximal optimization has proved crucial to address problems in both fields (see for instance 44, 122). Transportation metrics are also the correct way to compare and regularize variational problems that arise in image processing (see for instance the Radon inversion method proposed in 47) and machine learning (see 85).
3 Research program
Since its creation, the Mokaplan team has made important contributions in Optimal Transport both on the theoretical and the numerical side, together with applications such as fluid mechanics, the simulation biological systems, machine learning. We have also contributed to to the field of inverse problems in signal and image processing (super-resolution, nonconvex low rank matrix recovery). In 2022, the team was renewed with the following research program which broadens our spectrum and addresses exciting new problems.
3.1 OT and related variational problems solvers encore et toujours
Participants: Flavien Léger, Jean-David Benamou, Guillaume Carlier, Thomas Gallouët, Guillaume Chazareix , Adrien Vacher , Paul Pegon.
-
Asymptotic analysis of entropic OT
for a small entropic parameter is well understood for regular data on compact manifolds and standard quadratic ground cost 79, the team will extend this study to more general settings and also establish rigorous asymptotic estimates for the transports maps. This is important to provide a sound theoretical background to efficient and useful debiasing approaches like Sinkhorn Divergences 90. Guillaume Carlier , Paul Pegon and Luca Tamanini are investigating speed of convergence and quantitative stability results under general conditions on the cost (so that optimal maps may not be continuous or even fail to exist). Some sharp bounds have already been obtained, the next challenging goal is to extend the Laplace method to a nonsmooth setting and understand what entropic OT really selects when there are several optimal OT plans.
- High dimensional - Curse of dimensionality
-
Back-and-forth
The back-and-forth method 101, 100 is a state-of-the-art solver to compute optimal transport with convex costs and 2-Wasserstein gradient flows on grids. Based on simple but new ideas it has great potential to be useful for related problems. We plan to investigate: OT on point clouds in low dimension, the principal-agent problem in economics and more generally optimization under convex constraints 105, 117.
-
Transport and diffusion
The diffusion induced by the entropic regularization is fixed and now well understood. For recent variations of the OT problem (Martingale OT, Weak OT see 39) the diffusion becomes an explicit constraint or the control itself 96. The entropic regularisation of these problems can then be understood as metric/ground cost learning 69 (see also 121) and offers a tractable numerical method.
- Wasserstein Hamiltonian systems
-
Nonlinear fourth-order diffusion equations
such as thin-films or (the more involved) DLSS quantum drift equations are WGF. Such WGF are challenging both in terms of mathematical analysis (lack of maximum principle...) and of numerics. They are currently investigated by Jean-David Benamou , Guillaume Carlier in collaboration with Daniel Matthes. Note also that Mokaplan already contributed to a related topic through the TV-JKO scheme 71.
-
Lagrangian approaches for fluid mechanics
More generally we want to extend the design and implementation of Lagrangian numerical scheme for a large class of problem coming from fluids mechanics (WHS or WGF) using semi-discrete OT or entropic regularization. We will also take a special attention to link this approaches with problems in machine/statistical learning. To achieve this part of the project we will join forces with colleagues in Orsay University: Y. Brenier, H. Leclerc, Q. Mérigot, L. Nenna.
-
optimal transportis a variant of OT where we want to minimize the maximal displacement of the transport plan, instead of the average distance. Following the seminal work of 77, and more recent developments 86, Guillaume Carlier , Paul Pegon and Luigi De Pascale are working on the description of restrictable solutions (which are cyclically -monotone) through some potential maps, in the spirit of Mange-Kantorovich potentials provided by a duality theory. Some progress has been made to partially describe cyclically quasi-motonone maps (related in some sense to cyclically -monotone maps), through quasi-convex potentials.
3.2 Application of OT numerics to non-variational and non convex problems
Participants: Flavien Léger, Guillaume Carlier, Jean-David Benamou.
-
Market design
Z-mappings form a theory of non-variational problems initiated in the '70s but that has been for the most part overlooked by mathematicians. We are developing a new theory of the algorithms associated with convergent regular splitting of Z-mappings. Various well-established algorithms for matching models can be grouped under this point of view (Sinkhorn, Gale–Shapley, Bertsekas' auction) and this new perspective has the potential to unlock new convergence results, rates and accelerated methods.
-
Non Convex inverse problems
The PhD 134 provided a first exploration of Unbalanced Sinkhorn Divergence in this context. Given enough resources, a branch of PySit, a public domain software to test misfit functions in the context of Seismic imaging will be created and will allow to test other signal processing strategies in Full Waveform Inversion. Likewise the numerical method tested for 1D reflectors in 48 coule be develloped further (in particular in 2D).
-
Equilibrium and transport
Equilibrium in labor markets can often be expressed in terms of the Kantorovich duality. In the context of urban modelling or spatial pricing, this observation can be fruitfully used to compute equilibrium prices or densities as fixed points of operators involving OT, this was used in 41 and 40. Quentin Petit, Guillaume Carlier and Yves Achdou are currently developing a (non-variational) new semi-discrete model for the structure of cities with applications to tele-working.
-
Non-convex Principal-Agent problems
Guillaume Carlier , Xavier Dupuis, Jean-Charles Rochet and John Thanassoulis are developing a new saddle-point approach to non-convex multidimensional screening problems arising in regulation (Barron-Myerson) and taxation (Mirrlees).
3.3 Inverse problems with structured priors
Participants: Irène Waldspurger, Antonin Chambolle, Vincent Duval, Faniriana Rakoto Endor, Annette Dumas.
-
Off-the-grid reconstruction of complex objects
Whereas, very recently, some methods were proposed for the reconstruction of curves and piecewise constant images on a continuous domain (56 and 72), those are mostly proofs of concept, and there is still some work to make them competitive in real applications. As they are much more complex than point source reconstruction methods, there is room for improvements (parametrization, introduction of several atoms...). In particular, we are currently working on an improvement of the algorithm 56 for inverse problems in imaging which involve Optimal Transport as a regularizer (see 126 for preliminary results). Moreover, we need to better understand their convergence and the robustness of such methods, using sensitivity analysis.
-
Correctness guarantees for Burer-Monteiro methods
Burer-Monteiro methods work well in practice and are therefore widely used, but existing correctness guarantees 55 hold under unrealistic assumptions only. In the long term, we aim at proposing new guarantees, which would be slightly weaker but would hold in settings more relevant to practice. A first step is to understand the “average” behavior of Burer-Monteiro methods, when applied to random problems, and could be the subject of a PhD thesis.
3.4 Geometric variational problems, and their interactions with transport
Participants: Vincent Duval, Paul Pegon, Antonin Chambolle, Joao-Miguel Machado .
-
Approximation of measures with geometric constraints
Optimal Transport is a powerful tool to compare and approximate densities, but its interaction with geometric constraints is still not well understood. In applications such as optimal design of structures, one aims at approximating an optimal pattern while taking into account fabrication constraints 52. In Magnetic Resonance Imaging (MRI), one tries to sample the Fourier transform of the unknown image according to an optimal density but the acquisition device can only proceed along curves with bounded speed and bounded curvature 107. Our goal is to understand how OT interacts with energy terms which involve, e.g. the length, the perimeter or the curvature of the support... We want to understand the regularity of the solutions and to quantify the approximation error. Moreover, we want to design numerical methods for the resolution of such problems, with guaranteed performance.
-
Discretization of singular measures
Beyond the (B)Lasso and the total variation (possibly off-the-grid), numerically solving branched transportation problems requires the ability to faithfully discretize and represent 1-dimensional structures in the space. The research program of A. Chambolle consists in part in developing the numerical analysis of variational problems involving singular measures, such as lower-dimensional currents or free surfaces. We will explore both phase-field methods (with P. Pegon, V. Duval) 74, 118 which easily represent non-convex problems, but lack precision, and (with V. Duval) precise discretizations of convex problems, based either on finite elements (and relying to the FEM discrete exterior calculus 36, cf 75 for the case of the total variation), or on finite differences and possibly a clever design of dual constraints as studied in 80, 76 again for the total variation.
-
Transport problems with metric optimization
In urban planning models, one looks at building a network (of roads, metro or train lines, etc.) so as to minimize a transport cost between two distributions, penalized by the cost for building the network, usually its length. A typical transport cost is Monge cost
with a metric which is modified as a fraction of the euclidean metric on the network . We would like to consider general problems involving a construction cost to generate a conductance field (having in mind 1-dimensional integral of some function of ), and a transport cost depending on this conductance field. The afore-mentioned case studied in 62 falls into this category, as well as classical branched transport. The biologically-inspired network evolution model of 97 seems to provide such an energy in the vanishing diffusivity limit, with a cost for building a 1-dimensional permeability tensor and an congested transport cost with associated resistivity metric ; such a cost seems particularly relevant to model urban planning. Finally, we would like to design numerical methods to solve such problems, taking advantage of the separable structure of the whole cost.
4 Application domains
4.1 Natural Sciences
FreeForm Optics, Fluid Mechanics (Incompressible Euler, Semi-Geostrophic equations), Quantum Chemistry (Density Functional Theory), Statistical Physics (Schroedinger problem), Porous Media.
4.2 Signal Processing and inverse problems
Full Waveform Inversion (Geophysics), Super-resolution microscopy (Biology), Satellite imaging (Meteorology)
4.3 Social Sciences
Mean-field games, spatial economics, principal-agent models, taxation, nonlinear pricing.
5 Highlights of the year
We organized a workshop for the 10th anniversary of Mokaplan (MOKA10) at Ecole des Mines (Fontainebleau) on June, 4th to 6th. It gathered approximately 45 people who are close collaborators or former members of the team. It was a huge success.
5.1 Awards
Frontiers of Science 2024 prize (Beijing)for: “Iterative Bregman projections for regularized transportation problems”, SIAM Journal on Scientific Computing (2015) Authors: Jean-David Benamou, Guillaume Carlier, Marco Cuturi, Luca Nenna, Gabriel Peyré
6 New results
6.1 From entropic transport to martingale transport, and applications to model calibration
Participants: Jean-David Benamou, Guillaume Chazareix, Gregoire Loeper.
We propose a discrete time formulation of the semi martingale optimal transport problem based on multi-marginal entropic transport. This approach offers a new way to formulate and solve numerically the calibration problem proposed by Guo et al. 2022, using a multi-marginal extension of Sinkhorn algorithm as in Benamou, Carlier, and Nenna 2019; Carlier et al. 2017; Benamou et al. 2019. In the limit when the time step goes to zero we recover, as detailed in the companion paper Benamou et al. 2024, a semi-martingale process, solution to a semi-martingale optimal transport problem, with a cost function involving the so-called specific entropy introduced in Gantert 1991, see also Föllmer 2022 and Backhoff-Veraguas and Unterberger 2023.
6.2 Entropic Semi-Martingale Optimal Transport
Participants: Jean-David Benamou, Guillaume Chazareix, Marc Hoffman, Gregoire Loeper, François-Xavier Vialard.
23 Entropic Optimal Transport (EOT), also known as the Schrodinger problem, aims to find random processes with given initial and final marginals, minimizing the relative entropy (RE) with respect to a reference measure. Both processes (the reference and the controlled one) necessarily share the same diffusion coeffcients to ensure finiteness of the RE. This initially suggests that controlled-diffusion Semi-Martingale Optimal Transport (SMOT) problems may be incompatible with entropic regularization. However, when the process is observed at discrete times, forming a Markov chain, the RE remains finite even with variable diffusion coefficients. In this case, discrete semi-martingales can emerge as solutions to multi-marginal EOT problems. For smooth semi-martingales, the scaled limit of the relative entropy of their time discretizations converges to the specific relative entropy", a convex functional of the variance. This observation leads to an entropy regularized time discretization of the continuous SMOT problems, enabling the computation of discrete approximations via a multi-marginal Sinkhorn algorithm. We prove convergence of the time-discrete entropic problem to the continuous case, provide an implementation, and present numerical experiments supporting the theoretical results.
6.3 A geometric Laplace method
Participants: Flavien Léger, François-Xavier Vialard.
A classical tool for approximating integrals is the Laplace method. The first-order, as well as the higher-order Laplace formula is most often written in coordinates without any geometrical interpretation. In 9, motivated by a situation arising, among others, in optimal transport, we give a geometric formulation of the first-order term of the Laplace method. The central tool is the Kim–McCann Riemannian metric which was introduced in the field of optimal transportation. Our main result expresses the first-order term with standard geometric objects such as volume forms, Laplacians, covariant derivatives and scalar curvatures of two different metrics arising naturally in the Kim–McCann framework. Passing by, we give an explicitly quantified version of the Laplace formula, as well as examples of applications.
6.4 Convergence rate of general entropic optimal transport costs
Participants: Luca Nenna, Paul Pegon.
In 20 we improved and extended the results of 70 to the multi-marginal setting, investigating the convergence rate of the multi-marginal optimal entropic cost
6.5 A geometric approach to apriori estimates for optimal transport maps
Participants: Simon Brendle, Flavien Leger, Robert J. McCann, Cale Rankin.
A key inequality which underpins the regularity theory of optimal transport for costs satisfying the Ma–Trudinger–Wang condition is the Pogorelov second derivative bound. This translates to an apriori interior C1 estimate for smooth optimal maps. Here we give a new derivation of this estimate which relies in part on Kim, McCann and Warren's observation that the graph of an optimal map becomes a volume maximizing spacelike submanifold when the product of the source and target domains is endowed with a suitable pseudo-Riemannian geometry that combines both the marginal densities and the cost.
6.6 Gradient descent with a general cost
Participants: Flavien Leger, Pierre-Cyril Aubin.
We present a new class of gradient-type optimization methods that extends vanilla gradient descent, mirror descent, Riemannian gradient descent, and natural gradient descent. Our approach involves constructing a surrogate for the objective function in a systematic manner, based on a chosen cost function. This surrogate is then minimized using an alternating minimization scheme. Using optimal transport theory we establish convergence rates based on generalized notions of smoothness and convexity. We provide local versions of these two notions when the cost satisfies a condition known as nonnegative cross-curvature. In particular our framework provides the first global rates for natural gradient descent and the standard Newton's method.
6.7 Variational approximation of H-masses
Participants: Paul Pegon, Antonin Monteil.
In 18 we consider first order local minimization problems of the form
6.8 From geodesic extrapolation to a variational BDF2 scheme for Wasserstein gradient flows
Participants: Thomas Gallouët, Andrea Natale, Gabriele Todeschi.
We introduce a time discretization for Wasserstein gradient flows based on the classical Backward Differentiation Formula of order two. The main building block of the scheme is the notion of geodesic extrapolation in the Wasserstein space, which in general is not uniquely defined. We propose several possible definitions for such an operation, and we prove convergence of the resulting scheme to the limit PDE, in the case of the Fokker-Planck equation. For a specific choice of extrapolation we also prove a more general result, that is convergence towards EVI flows. Finally, we propose a variational finite volume discretization of the scheme which numerically achieves second order accuracy in both space and time.
6.9 Entropic approximation of infinity optimal transport problems
Participants: Camilla Brizzi, Guillaume Carlier, Luigi De-Pascale.
We propose an entropic approximation approach for optimal transportation problems with a supremal cost.
We establish
6.10 Quantitative Stability of the Pushforward Operation by an Optimal Transport Map
Participants: Guillaume Carlier, Alex Delalande, Quentin Mérigot.
We study the quantitative stability of the mapping that to a measure associates its pushforward measure by a fixed (non-smooth) optimal transport map. We exhibit a tight Hölder-behavior for this operation under minimal assumptions. Our proof essentially relies on a new bound that quantifies the size of the singular sets of a convex and Lipschitz continuous function on a bounded domain.
6.11 Quantitative Stability of Barycenters in the Wasserstein Space
Participants: Guillaume Carlier, Alex Delalande, Quentin Mérigot.
Wasserstein barycenters define averages of probability measures in a geometrically meaningful way. Their use is increasingly popular in applied fields, such as image, geometry or language processing. In these fields however, the probability measures of interest are often not accessible in their entirety and the practitioner may have to deal with statistical or computational approximations instead. In this article, we quantify the effect of such approximations on the corresponding barycenters. We show that Wasserstein barycenters depend in a Hölder-continuous way on their marginals under relatively mild assumptions. Our proof relies on recent estimates that quantify the strong convexity of the dual quadratic optimal transport problem and a new result that allows to control the modulus of continuity of the push-forward operation under a (not necessarily smooth) optimal transport map.
6.12 Wasserstein medians: robustness, PDE characterization and numerics
Participants: Guillaume Carlier, Enis Chenchene, Katharina Eichinger.
We investigate the notion of Wasserstein median as an alternative to the Wasserstein barycenter, which has become popular but may be sensitive to outliers. In terms of robustness to corrupted data, we indeed show that Wasserstein medians have a breakdown point of approximately 1 2. We give explicit constructions of Wasserstein medians in dimension one which enable us to obtain
6.13 Projected gradient descent accumulates at Bouligand stationary points
Participants: Guillaume Olikier, Irène Waldspurger.
In 33, we consider the problem of minimizing a continuously differentiable function on a nonempty closed subset of a Euclidean vector space using projected gradient descent (PGD). Without further assumptions on the objective function, PGD (as all standard algorithms) is not expected to find a global, or even local minimizer. Instead, it finds a so-called stationary point. Up to now, it was known that the accumulation points of sequences generated by PGD were Mordukhovich stationary, which is a relatively weak notion of stationarity. In the article, the accumulation points are proven to satisfy the stronger definition of Bouligand stationarity, and even proximal stationarity if the gradient is locally Lipschitz continuous. These are the strongest stationarity properties that can be expected for the considered problem.
6.14 Benign landscape for Burer-Monteiro factorizations of MaxCut-type semidefinite programs
Participants: Faniriana Rakoto Endor, Irène Waldspurger.
We study the Burer-Monteiro factorization, which is a well-known heuristic to reduce the computational cost of solving semidefinite programs (SDP) in the case where the solution is a priori known to be low rank. This factorization reduces the dimension of the SDP at the cost of its convexity, therefore possibly introducing spurious second-order critical points which could trap the optimization algorithm and prevent it from finding the desired minimizer.
In 31, for MaxCut-type SDP, we give a sharp condition on the conditioning of the associated Laplacian matrix under which any second-order critical point of the non-convex problem is a global minimizer. By applying our theorem, we improve on recent results about the correctness of the Burer-Monteiro factorization on
6.15 Nonnegative cross-curvature in infinite dimensions: synthetic definition and spaces of measures
Participants: Flavien Léger, Gabriele Todeschi, François-Xavier Vialard.
Nonnegative cross-curvature (NNCC) is a geometric property of a cost function defined on a product space that originates in optimal transportation and the Ma-Trudinger-Wang theory. Motivated by applications in optimization, gradient flows and mechanism design, we propose a variational formulation of nonnegative cross-curvature on c-convex domains applicable to infinite dimensions and nonsmooth settings. The resulting class of NNCC spaces is closed under Gromov-Hausdorff convergence and for this class, we extend many properties of classical nonnegative cross-curvature: stability under generalized Riemannian submersions, characterization in terms of the convexity of certain sets of c-concave functions, and in the metric case, it is a subclass of positively curved spaces in the sense of Alexandrov. One of our main results is that Wasserstein spaces of probability measures inherit the NNCC property from their base space. Additional examples of NNCC costs include the Bures-Wasserstein and Fisher-Rao squared distances, the Hellinger-Kantorovich squared distance (in some cases), the relative entropy on probability measures, and the 2-Gromov-Wasserstein squared distance on metric measure spaces.
6.16 Dynamical programming for off-the-grid dynamic inverse problems
Participants: Vincent Duval, Robert Tovey.
In 14, we consider off-the-grid algorithms for the reconstruction of sparse measures from time-varying data. In particular, the reconstruction is a finite collection of Dirac measures whose locations and masses vary continuously in time. Recent work showed that this decomposition was possible by minimising a convex variational model which combined a quadratic data fidelity with dynamical Optimal Transport. We generalise this framework and propose new numerical methods which leverage efficient classical algorithms for computing shortest paths on directed acyclic graphs. Our theoretical analysis confirms that these methods converge to globally optimal reconstructions. Numerically, we show new examples for unbalanced Optimal Transport penalties, and for balanced examples we are 100 times faster in comparison to the previously known method.
6.17 Inclusion and estimates for the jumps of minimizers in variational denoising
Participants: Antonin Chambolle, Michał Łasica.
We study in 13 the stability and inclusion of the jump set of minimizers of convex denoising functionals, such as the celebrated "Rudin-Osher-Fatemi" functional, for scalar or vectorial signals. We show that under mild regularity assumptions on the data fidelity term and the regularizer, the jump set of the minimizer is essentially a subset of the original jump set. Moreover, we give an estimate on the magnitude of jumps in terms of the data. This extends old results, in particular of the first author (with V. Caselles and M. Novaga) and of T. Valkonen, to much more general cases. We also consider the case where the original datum has unbounded variation, and define a notion of its jump set which, again, must contain the jump set of the solution.
6.18 L1-Gradient Flow of Convex Functionals
Participants: Antonin Chambolle, Matteo Novaga.
In 12, we are interested in the gradient flow of a general first order convex functional with respect to the L¹-topology. By means of an implicit minimization scheme, we show existence of a global limit solution, which satisfies an energy-dissipation estimate, and solves a non-linear and non-local gradient flow equation, under the assumption of strong convexity of the energy. Under a monotonicity assumption we can also prove uniqueness of the limit solution, even though this remains an open question in full generality. We also consider a geometric evolution corresponding to the L¹-gradient flow of the anisotropic perimeter. When the initial set is convex, we show that the limit solution is monotone for the inclusion, convex and unique until it reaches the Cheeger set of the initial datum. Eventually, we show with some examples that uniqueness cannot be expected in general in the geometric case.
6.19 Convergent plug-and-play with proximal denoiser and unconstrained regularization parameter
Participants: Antonin Chambolle, Samuel Hurault, Arthur Leclaire, Nicolas Papadakis.
In 17, we present new proofs of convergence for Plug-and-Play (PnP) algorithms. PnP methods are efficient iterative algorithms for solving image inverse problems where regularization is performed by plugging a pre-trained denoiser in a proximal algorithm, such as Proximal Gradient Descent (PGD) or Douglas-Rachford Splitting (DRS). Recent research has explored convergence by incorporating a denoiser that writes exactly as a proximal operator. However, the corresponding PnP algorithm has then to be run with stepsize equal to 1. The stepsize condition for nonconvex convergence of the proximal algorithm in use then translates to restrictive conditions on the regularization parameter of the inverse problem. This can severely degrade the restoration capacity of the algorithm. In this paper, we present two remedies for this limitation. First, we provide a novel convergence proof for PnP-DRS that does not impose any restrictions on the regularization parameter. Second, we examine a relaxed version of the PGD algorithm that converges across a broader range of regularization parameters. Our experimental study, conducted on deblurring and super-resolution experiments, demonstrate that both of these solutions enhance the accuracy of image restoration.
6.20 A Cahn–Hilliard–Willmore phase field model for non-oriented interfaces
Participants: Antonin Chambolle, Elie Bretin, Simon Masnou.
In 26, we investigate a new phase field model for representing non-oriented interfaces, approximating their area and simulating their area-minimizing flow. Our contribution is related to the approach proposed in arXiv:2105.09627 that involves ad hoc neural networks. We show here that, instead of neural networks, similar results can be obtained using a more standard variational approach that combines a Cahn-Hilliard-type functional involving an appropriate non-smooth potential and a Willmore-type stabilization energy. We show some properties of this phase field model in dimension 1 and, for radially symmetric functions, in arbitrary dimension. We propose a simple numerical scheme to approximate its L2-gradient flow. We illustrate numerically that the new flow approximates fairly well the mean curvature flow of codimension 1 or 2 interfaces in dimensions 2 and 3.
6.21 Phase-field approximation for 1-dimensional shape optimization problems
Participants: Joao-Miguel Machado.
In 32, we propose an unified framework for the phase field approximation of 1-dimensional shape optimization problems with connectedness constraints in any dimension. In particular, we focus on the average distance minimizers problem and the Wasserstein-
7 Bilateral contracts and grants with industry
Participants: Jean-David Benamou, Gregoire Loeper.
CIFRE PhD thesis scholarship (Guillaume Chazareix ) with BNP. Main supervisor Jean-David Benamou , co-supervision with Guillaume Carlier (Inria Mokaplan) and Gregoire Loeper (BNP). This contract is handled by Dauphine University.
8 Partnerships and cooperations
8.1 International research visitors
8.1.1 Visits of international scientists
Other international visits to the team
Participants: Guillaume Carlier, Jean-David Benamou, Paul Pegon.
Luigi De Pascale
-
Status
Associate Professor
-
Institution of origin:
Universita degli Studi Firenze
-
Country:
Italie
-
Dates:
06/03–06/07, 06/26–07-10, 09/09–09/13 (the visit has been split because of the relocation of our offices)
-
Context of the visit:
Luigi De Pascale is a long term collaborator of our team. He visited us to work on ongoing projects with Guillaume Carlier and Paul Pegon on
-type optimal transport. -
Mobility program/type of mobility:
research stay
Jacob Francis
-
Status
PhD student
-
Institution of origin:
Imperial College London
-
Country:
Royaume-Uni
-
Dates:
From March, 4th to May, 10th
-
Context of the visit:
This is part of an ongoing project on the use of Entropic Optimal Transport in the context of the semi-geostrophic approximation of rotating fluids. The research also involves Hugo Malamut , Colin Cotter (ICL) and Jean-David Benamou
-
Mobility program/type of mobility:
research stay
8.1.2 Visits to international teams
Research stays abroad
Jean-David Benamou
-
Visited institution:
Imperial College
-
Country:
Royaume-Uni
-
Dates:
Nov. 12th - Dec 11th
-
Context of the visit:
This is part of an ongoing project on the use of Entropic Optimal Transport in the context of the semi-geostrophic approximation of rotating fluids. The research also involves Hugo MalamutColin Cotter and Jacob Francis (ICL).
-
Mobility program/type of mobility:
CNRS Imperial Abraham de Moivre Fellowship.
8.2 National initiatives
-
ANR CIPRESSI
(2019-2024) is a JCJC grant (149k€) carried by Vincent Duval . Its aim is to develop off-the-grid methods for inverse problems involving the reconstruction of complex objects.
-
PDE AI
(2023-2027) Antonin Chambolle is the main coordinator of the PDE-AI project, funded by the PEPR IA (France 2030, ANR) and gathering 10 groups throughout France working on PDEs and nonlinear analysis for artificial intelligence.
-
ANR GOTA
is a JCJC grant (253k€) carried by Luca Nenna (PI), Paul Pegon and Maxime Laborde , dealing with some generalizations and applications of Optimal Transport theory with a particular focus on three main topics: multi-marginal optimal transport, urban planning and ulti-population models, and multi-marginal entropic optimal transport.
9 Dissemination
9.1 Promoting scientific activities
9.1.1 Scientific events: organisation
General chair, scientific chair
Jean-David Benamou organizer of : Stochastic Control, Optimal Transport, Calibration applications in Finance, workshop Carry le Rouet.
Member of the organizing committees
Antonin Chambolle is a co-organizer (member of the organizing and scientific committee) of “PGMO days 2024”, EDF Saclay, 19–20 Nov. 2024
Vincent Duval
was a member of the organizing committee of the Workhop Julia & Optimization Days 2024 in Toulouse in October 2024 (
Guillaume Carlier coorganized with Eric Séré a two days conference in september in honor of Ivar Ekeland at Dauphine, a two days workshop with Alessio Figalli, Quentin Mérigot and Filippo Santambrogio on Particle systems and dynamics, optimization and learning at the Lagrange Center in march and is a coorganizer of the (monthly) séminaire parisien d'optimisation.
Paul Pegon and Vincent Duval coorganized the workshop for the 10th anniversary of Mokaplan.
9.1.2 Scientific events: selection
Reviewer
- Antonin Chambolle is a reviewer for AISTATS 2025 and SSVM 2025
9.1.3 Journal
Member of the editorial boards
Antonin Chambolle is in the editorial board of
- Interfaces and Free Boundaries (EMS Press, one of the 4 co-editors in chief)
- Journal of the European Math. Society (JEMS, EMS Press)
- M2AN Mathematical Modeling and Numerical Analysis
- AMO Applied Math. Optim (Springer)
- IPI Inverse problems and imaging (AIMS)
- JMIV Journal of Math. Imaging and Vision (Springer)
- A special issue of JOTA (Journal of Optimization Theory and Applications, Springer, special issue on optimization for energy).
Vincent Duval is associate editor for the Journal of Mathematical Imaging and Vision journal (JMIV).
Guillaume Carlier i s associate editor of Journal of Mathematical Analysis and Applications, Applied Mathematics and Optimization, Journal of Dynamics and Games.
Reviewer - reviewing activities
- Irène Waldspurger : Mathematical Programming.
- Vincent Duval : Computational Optimization and Applications, SIAM Journal on Imaging Sciences, Mathematical reviews (x2).
9.1.4 Invited talks
- Flavien Léger : Journée SIGMA – MODE, INRIA Paris, January 2024 ;
- Flavien Léger : Learning and Optimization in Luminy, CIRM, June 2024 ;
- Flavien Léger : Math+Econ+Code seminar, Paris, October 2024 ;
- Flavien Léger : Séminaire du MIA Paris-Saclay, November 2024 ;
- Flavien Léger : Optimal Transportation and Applications, Pisa, December 2024 ;
- Flavien Léger : Online seminar on Geometry and Statistics, December 2024 ;
- Irène Waldspurger : journée thématique Low-rank approximation and optimization, at the Institut de mathématiques de Marseille ;
- Irène Waldspurger : session on Negative eigenvalues and nonconvex optimization, at the SIAM conference on Applied Linear Algebra (Paris) ;
- Irène Waldspurger : Women in Machine Learning and Data Science event ;
- Guillaume Carlier : Colloquium Créteil, Workshop on Optimal Transport from Theory to Applications, Berlin, TSE Madstat Seminar, Toulouse, GT Calva à Orsay, Seminar on Optimization, University of Vienna, Nonlinear Analysis seminar, ISTA, Vienna, Optimal Transportation and Applications Centro di Ricerca Matematica Ennio De Giorgi, Pisa.
- Paul Pegon : Analysis Seminar, École Polytechnique (Palaiseau), March 2024 ;
- Paul Pegon : GT Optimal Transport and Optimisation, Université Aix-Marseille and Université de Toulon, April 2024 ;
- Paul Pegon : Mini-symposium at CANUM 2024, Université de La Rochelle, May 2024 ;
- Jean-David Benamou : Numerical methods for optimal transport problems, mean field games and multi-agent dynamics, Universidad Federico Santa María in Valparaíso, Chile.
- Jean-David Benamou : MFO Workshop 2406 on Optimal Transportation (Oberwolfach).
- Jean-David Benamou : MOKA10 workshop (ten years of Mokaplan), Fontainebleau .
- Jean-David Benamou : Optimal transport workshop, IES CNRS, Cargese.
- Vincent Duval : Journée SIGMA – MODE, INRIA Paris, January 2024 ;
- Antonin Chambolle was invited, and gave a seminar a the MFO Oberwolfach workshop #2407, “Interfaces, Free Boundaries and Geometric Partial Differential Equations”, 11–16 Feb. 2024 (talk on “Discrete-to-continuum crystalline curvature flow”, 12 Feb. 2024).
-
Antonin Chambolle
was invited to give a talk at the (French) mini-workshop
“a
tour of optimal design”, Autrans, 6–8 March 2024, in honor of Grégoire Allaire, Eric Bonnetier, François Jouve. - Antonin Chambolle gave a plenary talk at the conference “ICIPE 24” (“Internat. Conf. on Inverse Problems in Engineering”), Buzios, Brazil, 23–27 June 2024, on “Qualitative Properties of Minimizers of Total-Variation Regularized Problems” (26 June).
- Antonin Chambolle gave an invited talk at the “IFIP TC7 Meeting on System Modeling and Optimization” in Hamburg, 12–16 Aug. 2024, “Qualitative Properties of Minimizers of Total-Variation Regularized Problems” (13 Aug.).
- Antonin Chambolle gave an invited talk at “VARANA 2024”, Erice, Italie, on “Discrete-to-continuum crystalline curvature flow”, 6 Sept.
- Antonin Chambolle participated to the “Conférence pour les 50 ans du CMAP”, 11–13 Sept. 2024 and gave a talk on “Discrete-to-continuum crystalline curvature flow” (13 Sept.).
- Antonin Chambolle gave a talk at the workshop “Mathematical Materials Science: Defects and Polycrystals”, Edinburgh, 13–19 Oct. 2024, on “Discrete-to-continuum crystalline curvature flow” (15 Oct.).
9.1.5 Research administration
- Antonin Chambolle is a member of the scientific council of Université Paris-Dauphine. He is also a member of the scientfic council and of the board of the PGMO “Programme Gaspard Monge pour l'Optimisation et la Recherche Opérationnelle”.
- Irène Waldspurger : hiring committee for a maîtresse de conférence position at INSA Rennes.
- Vincent Duval is "délégué scientifique adjoint" of the Inria Paris center.
- Vincent Duval was vice-president of the CRCN-ISFP selection committee of the Inria Paris center, and a member of the selection committee for the DR2 Inria selection.
- Guillaume Carlier is in the scientific board of the PGMO
- Paul Pegon was a member of a hiring committee for a Maître de conférences position at Université de Toulon
9.2 Teaching - Supervision - Juries
9.2.1 Teaching
- Master: Antonin Chambolle Optimisation Continue, 24h, niveau M2, Université Paris Dauphine-PSL, FR
- Master : Vincent Duval , Problèmes Inverses, 22,5 h équivalent TD, niveau M1, Université PSL/Mines ParisTech, FR
- Licence : Irène Waldspurger , Algorithmique et programmation, 34,5 h équivalent TD, niveau L2, Université Paris-Dauphine, FR
- Master : Irène Waldspurger , Optimization for Machine Learning, 6h, niveau M2, Université PSL/ENS, FR
- Master : Irène Waldspurger , Introduction à la géométrie différentielle et aux équations différentielles, 29,25 h équivalent TD, niveau M1, Université Paris Dauphine, FR
- Master : Irène Waldspurger , Non-convex inverse problems, 27 h d'équivalent TD, niveau M2, Université Paris Dauphine, FR
- Licence : Guillaume Carlier , algebre 1, L1 78h, Dauphine, FR
- Master/PhD : Guillaume Carlier Optimal transport, University of Vienna, 21h
- Guillaume Carlier : Licence Algèbre 1, Dauphine 70h, M2 Masef: Variational and transport problems in economics, 18h
- Licence : Paul Pegon , Analysis 1, L1, 146h eq TD, FR
- Licence : Paul Pegon , Analysis 4, L1, 50h eq TD, FR
9.2.2 Supervision
- L3 memoir completed: Enzo Regna (Univ Paris-Dauphine) supervised by Paul Pegon .
- M2 internship completed: Majid Arthaud (ENPC), co-supervised by Antonin Chambolle and Vincent Duval .
- M2 internship completed: Paul Caucheteux (master MASH, Dauphine) supervised by Irène Waldspurger .
- M2 internship completed: Lorenzo Braglia (Università di Padova), co-supervised by Paul Pegon and Antonin Chambolle .
- M2 internship completed: Louis Tocquec ( ENS Paris-Saclay), co-supervised by Paul Pegon and Luca Nenna .
- Phd completed: Joao Pinto Anastacio Machado , Transport optimal et structures géométriques, defended on 17/09/2024, co-supervised by Antonin Chambolle and Vincent Duval .
- PhD in progress: Faniriana Rakoto Endor , Pourquoi l'heuristique de Burer-Monteiro fonctionne-t-elle si bien?, started on 01/10/2023, co-supervised by Irène Waldspurger and Antonin Chambolle .
- Phd in progress: Hugo Malamut , Régularisation Entropique et Transport Optimal Généralisé, started on 1/09/2022, co-supervised by Jean-David Benamou and Guillaume Carlier .
- Phd in progress: Maxime Sylvestre , On Hybrid methods fot Optimal Transport, started on 01/09/2022, co-supervised by Guillaume Carlier and Alfred Galichon.
- Phd in progress: Guillaume Chazareix , Non Linear Parabolic equations and Volatility Calibration, started on 1/08/2021, co-supervised by Jean-David Benamou and Grégoire Loeper.
- PhD in progress: Louis Tocquec , On entropic optimal transport and some generalizations, started on 01/10/2024, co-supervised by Paul Pegon and Luca Nenna .
- Postdoc completed: Adrien Vacher , supervised by Flavien Léger .
9.2.3 Juries
- Vincent Duval was a jury member for the PhD of Théo Bertrand (Université Paris-Dauphine PSL)
- Irène Waldspurger was a jury member for the PhD of Jean-Jacques Godème (ENSICAEN).
- Irène Waldspurger was a jury member for the PhD of Pierre-Jean Bénard (Institut Mathématique de Bordeaux).
- Guillaume Carlier was the coordinator of the HDR of Thomas Gallouët and a referee for the PhD's of Jules Candau-Tilh (Lille) and Giacomo Greco (Eindhoven, Netherlands).
9.3 Popularization
9.3.1 Productions (articles, videos, podcasts, serious games, ...)
Vincent Duval has worked together with a cartoonist (Julien Joliclerc) for the writing of a popularization cartoon on his ANR JCJC project, for the ANR Medianum project.
9.3.2 Participation in Live events
Flavien Léger gave the “demi-heure de science” at inria Paris, October 2024.
10 Scientific production
10.1 Major publications
- 1 inproceedingsMirror Descent with Relative Smoothness in Measure Spaces, with application to Sinkhorn and EM.NeurIPS 2022 - Thirty-sixth Conference on Neural Information Processing SystemsNew Orleans, United States2022HAL
- 2 articleIterative Bregman Projections for Regularized Transportation Problems.SIAM Journal on Scientific Computing2372015, A1111-A1138HALDOI
- 3 articleSecond order models for optimal transport and cubic splines on the Wasserstein space.Foundations of Computational MathematicsOctober 2019HALDOI
- 4 articleOn Representer Theorems and Convex Regularization.SIAM Journal on Optimization292May 2019, 1260–1281HALDOI
- 5 articleA variational finite volume scheme for Wasserstein gradient flows.Numerische Mathematik14632020, pp 437 - 480HALDOI
- 6 articleConvergence of Entropic Schemes for Optimal Transport and Gradient Flows.SIAM Journal on Mathematical Analysis492April 2017HALDOI
- 7 miscConvergence rate of general entropic optimal transport costs.June 2022HAL
- 8 miscGradient descent with a general cost.December 2023HAL
- 9 miscA geometric Laplace method.December 2022HALback to text
- 10 articlePhase retrieval with random Gaussian sensing vectors by alternating projections.IEEE Transactions on Information Theory6452018, 3301-3312HAL
- 11 articleRank optimality for the Burer-Monteiro factorization.SIAM Journal on Optimization3032020, 2577-2602HALDOI
10.2 Publications of the year
International journals
- 12 articleL1-Gradient Flow of Convex Functionals.SIAM Journal on Mathematical Analysis565September 2024, 5747-5781HALDOIback to text
- 13 articleInclusion and estimates for the jumps of minimizers in variational denoising.SIAM Journal on Imaging Sciences173September 2024, 1844-1878HALDOIback to text
- 14 articleDynamical programming for off-the-grid dynamic inverse problems.ESAIM: Control, Optimisation and Calculus of Variations302024, 7HALDOIback to text
- 15 articleEfficient preconditioners for solving dynamical optimal transport via interior point methods.SIAM Journal on Scientific Computing4632024HALDOI
- 16 articleFrom geodesic extrapolation to a variational BDF2 scheme for Wasserstein gradient flows.Mathematics of Computation932024, 2769-2810HALDOI
- 17 articleConvergent plug-and-play with proximal denoiser and unconstrained regularization parameter.Journal of Mathematical Imaging and Vision2024. In press. HALback to text
- 18 articleMass concentration in rescaled first order integral functionals.Journal de l'École polytechnique — Mathématiques11February 2024, 42HALDOIback to text
- 19 articleNear-optimal estimation of smooth transport maps with kernel sums-of-squares.SIAM Journal on Mathematics of Data Science2024. In press. HAL
- 20 articleConvergence rate of entropy-regularized multi-marginal optimal transport costs.Canadian Journal of Mathematics = Journal Canadien de MathématiquesMarch 2024HALDOIback to text
- 21 articleModel-based Clustering with Missing Not At Random Data.Statistics and ComputingJune 2024HALDOI
- 22 articleProduct of Gaussian Mixture Diffusion Models.Journal of Mathematical Imaging and Vision664March 2024, 504-528HALDOI
Reports & preprints
- 23 miscEntropic Semi-Martingale Optimal Transport.2024HALDOIback to text
- 24 miscFrom entropic transport to martingale transport, and applications to model calibration.June 2024HALback to text
- 25 miscOn the global convergence of Wasserstein gradient flow of the Coulomb discrepancy.January 2024HAL
- 26 miscA Cahn--Hilliard--Willmore phase field model for non-oriented interfaces.2024HALback to text
- 27 miscDisplacement smoothness of entropic optimal transport.March 2024HAL
- 28 miscWell-posedness and convergence of entropic approximation of semi-geostrophic equations.April 2024HAL
- 29 miscDiscrete-to-continuous crystalline curvature flows.March 2024HAL
- 30 miscImage = Cartoon+Texture: How Yves Meyer's "Oscillating patterns in image processing and in some nonlinear evolution equations" ended up in a computer vision model.2024, 75-94HAL
- 31 miscBenign landscape for Burer-Monteiro factorizations of MaxCut-type semidefinite programs.November 2024HALback to text
- 32 miscPhase-field approximation for 1-dimensional shape optimization problems.June 2024HALback to text
- 33 miscProjected gradient descent accumulates at Bouligand stationary points.March 2024HALback to text
10.3 Cited publications
- 34 articleBarycenters in the Wasserstein space.SIAM J. Math. Anal.4322011, 904--924back to text
- 35 articleHamiltonian ODEs in the Wasserstein space of probability measures.Communications on Pure and Applied Mathematics: A Journal Issued by the Courant Institute of Mathematical Sciences6112008, 18--53back to text
- 36 articleFinite element exterior calculus, homological techniques, and applications.Acta Numerica152006, 1–155DOIback to text
- 37 articleConsistency of Trace Norm Minimization.J. Mach. Learn. Res.9June 2008, 1019--1048URL: http://dl.acm.org/citation.cfm?id=1390681.1390716back to text
- 38 articleConsistency of the Group Lasso and Multiple Kernel Learning.J. Mach. Learn. Res.9June 2008, 1179--1225URL: http://dl.acm.org/citation.cfm?id=1390681.1390721back to text
- 39 miscApplications of weak transport theory.2020back to text
- 40 unpublishedA spatial Pareto exchange economy problem.December 2021, working paper or preprintHALback to text
- 41 articleA mean field game model for the evolution of cities.Journal of Dynamics and Games2021HALback to text
- 42 articleComputing Large Deformation Metric Mappings via Geodesic Flows of Diffeomorphisms.International Journal of Computer Vision612February 2005, 139--157URL: http://dx.doi.org/10.1023/B:VISI.0000043755.93987.aaback to textback to text
- 43 articleModel-independent bounds for option prices mass transport approach.Finance and Stochastics1732013, 477--501URL: http://dx.doi.org/10.1007/s00780-013-0205-8DOIback to text
- 44 articleA computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem.Numer. Math.8432000, 375--393URL: http://dx.doi.org/10.1007/s002110050002DOIback to textback to text
- 45 articleWeak existence for the semigeostrophic equations formulated as a coupled Monge-Ampère/transport problem.SIAM J. Appl. Math.5851998, 1450--1461back to textback to text
- 46 articleAugmented Lagrangian algorithms for variational problems with divergence constraints.JOTA2015back to text
- 47 articleIterative Bregman Projections for Regularized Transportation Problems.SIAM J. Sci. Comp.to appear2015back to textback to textback to text
- 48 articlePoint Source Regularization of the Finite Source Reflector Problem.Journal of Computational PhysicsMay 2022HALback to text
- 49 articleNumerical solution of the optimal transportation problem using the Monge--Ampere equation.Journal of Computational Physics2602014, 107--126back to text
- 50 articleTwo numerical methods for the elliptic Monge-Ampère equation.M2AN Math. Model. Numer. Anal.4442010, 737--758back to text
- 51 articleConsistent estimation of a population barycenter in the Wasserstein space.Preprint arXiv:1212.25622012back to text
- 52 articleAdditive manufacturing scanning paths optimization using shape optimization tools.Struct. Multidiscip. Optim.6162020, 2437--2466URL: https://doi.org/10.1007/s00158-020-02614-3DOIback to text
- 53 articleSliced and Radon Wasserstein Barycenters of Measures.Journal of Mathematical Imaging and Vision5112015, 22--45URL: http://hal.archives-ouvertes.fr/hal-00881872/back to text
- 54 articleCharacterization of optimal shapes and masses through Monge-Kantorovich equation.J. Eur. Math. Soc. (JEMS)322001, 139--168URL: http://dx.doi.org/10.1007/s100970000027DOIback to text
- 55 articleDeterministic guarantees for Burer-Monteiro factorizations of smooth semidefinite programs.preprinthttps://arxiv.org/abs/1804.020082018back to text
- 56 articleA generalized conditional gradient method for dynamic inverse problems with optimal transport regularization.arXiv preprint arXiv:2012.117062020back to textback to text
- 57 articleDécomposition polaire et réarrangement monotone des champs de vecteurs.C. R. Acad. Sci. Paris Sér. I Math.305191987, 805--808back to text
- 58 articleReconstruction of the early universe as a convex optimization problem.Mon. Not. Roy. Astron. Soc.3462003, 501--524URL: http://arxiv.org/pdf/astro-ph/0304214.pdfback to text
- 59 articleGeneralized solutions and hydrostatic approximation of the Euler equations.Phys. D23714-172008, 1982--1988URL: http://dx.doi.org/10.1016/j.physd.2008.02.026DOIback to text
- 60 articlePolar factorization and monotone rearrangement of vector-valued functions.Comm. Pure Appl. Math.4441991, 375--417URL: http://dx.doi.org/10.1002/cpa.3160440402DOIback to text
- 61 articleA guide to the TV zoo.Level-Set and PDE-based Reconstruction Methods, Springer2013back to text
- 62 bookOptimal Urban Networks via Mass Transportation.1961Lecture Notes in MathematicsBerlin, HeidelbergSpringer Berlin Heidelberg2009, URL: http://link.springer.com/10.1007/978-3-540-85799-0DOIback to text
- 63 incollectionOn the numerical solution of the problem of reflector design with given far-field scattering data.Monge Ampère equation: applications to geometry and optimization (Deerfield Beach, FL, 1997)226Contemp. Math.Providence, RIAmer. Math. Soc.1999, 13--32URL: http://dx.doi.org/10.1090/conm/226/03233DOIback to text
- 64 articleThe regularity of mappings with a convex potential.J. Amer. Math. Soc.511992, 99--104URL: http://dx.doi.org/10.2307/2152752DOIback to text
- 65 articleComputational Analysis of LDDMM for Brain Mapping.Frontiers in Neuroscience72013back to text
- 66 articleSuper-Resolution from Noisy Data.Journal of Fourier Analysis and Applications1962013, 1229--1254back to text
- 67 articleTowards a Mathematical Theory of Super-Resolution.Communications on Pure and Applied Mathematics6762014, 906--956back to text
- 68 articleAn Introduction to Compressive Sensing.IEEE Signal Processing Magazine2522008, 21--30back to text
- 69 unpublishedSISTA: Learning Optimal Transport Costs under Sparsity Constraints.October 2020, working paper or preprintHALback to text
- 70 articleConvergence rate of general entropic optimal transport costs.Calculus of Variations and Partial Differential Equations624May 2023, 116HALDOIback to text
- 71 articleOn the total variation Wasserstein gradient flow and the TV-JKO scheme.ESAIM: Control, Optimisation and Calculus of Variations2019HALback to text
- 72 articleTowards Off-the-grid Algorithms for Total Variation Regularized Inverse Problems.Journal of Mathematical Imaging and VisionJuly 2022HALDOIback to text
- 73 articleKinetic models for chemotaxis and their drift-diffusion limits.Monatsh. Math.1421-22004, 123--141URL: http://dx.doi.org/10.1007/s00605-004-0234-7DOIback to text
-
74
articleVariational approximation of size-mass energies for
-dimensional currents.ESAIM Control Optim. Calc. Var.252019, Paper No. 43, 39URL: https://doi.org/10.1051/cocv/2018027DOIback to text - 75 articleCrouzeix-Raviart approximation of the total variation on simplicial meshes.J. Math. Imaging Vision626-72020, 872--899URL: https://doi.org/10.1007/s10851-019-00939-3DOIback to text
- 76 articleLearning consistent discretizations of the total variation.SIAM J. Imaging Sci.1422021, 778--813URL: https://doi.org/10.1137/20M1377199DOIback to text
-
77
articleThe
-Wasserstein Distance: Local Solutions and Existence of Optimal Transport Maps.401January 2008, 1--20URL: https://epubs.siam.org/doi/10.1137/07069938XDOIback to text - 78 articleAtomic decomposition by basis pursuit.SIAM journal on scientific computing2011999, 33--61back to text
- 79 inproceedingsFaster Wasserstein Distance Estimation with the Sinkhorn Divergence.Neural Information Processing SystemsAdvances in Neural Information Processing SystemsVancouver, CanadaDecember 2020HALback to text
- 80 articleDiscrete total variation: new definition and minimization.SIAM J. Imaging Sci.1032017, 1258--1290URL: https://doi.org/10.1137/16M1075247back to text
- 81 articleDensity Functional Theory and Optimal Transportation with Coulomb Cost.Communications on Pure and Applied Mathematics6642013, 548--599URL: http://dx.doi.org/10.1002/cpa.21437DOIback to text
- 82 bookA Mathematical Theory of Large-Scale Atmosphere/Ocean Flow.Imperial College Press2006, URL: https://books.google.fr/books?id=JxBqDQAAQBAJback to text
- 83 articleThe semigeostrophic equations discretized in reference and dual variables.Arch. Ration. Mech. Anal.18522007, 341--363URL: http://dx.doi.org/10.1007/s00205-006-0040-6DOIback to text
- 84 articleGeneralised Lagrangian solutions for atmospheric and oceanic flows.SIAM J. Appl. Math.5111991, 20--31back to text
- 85 inproceedingsSinkhorn Distances: Lightspeed Computation of Optimal Transport.Proc. NIPS2013, 2292--2300back to textback to text
-
86
articleA Study of the Dual Problem of the One-Dimensional
-Optimal Transport Problem with Applications.27611June 2019, 3304--3324URL: https://www.sciencedirect.com/science/article/pii/S0022123619300643DOIback to text - 87 articleNumerical methods for fully nonlinear elliptic equations of the Monge-Ampère type.Comput. Methods Appl. Mech. Engrg.19513-162006, 1344--1386back to text
- 88 articleExact Support Recovery for Sparse Spikes Deconvolution.Foundations of Computational Mathematics2014, 1--41URL: http://dx.doi.org/10.1007/s10208-014-9228-6DOIback to text
- 89 articleSupport detection in super-resolution.Proc. Proceedings of the 10th International Conference on Sampling Theory and Applications2013, 145--148back to text
- 90 unpublishedInterpolating between Optimal Transport and MMD using Sinkhorn Divergences.October 2018, working paper or preprintHALback to text
- 91 articleBreaking the Curse of Dimension in Multi-Marginal Kantorovich Optimal Transport on Finite State Spaces.SIAM Journal on Mathematical Analysis5042018, 3996--4019URL: https://doi.org/10.1137/17M1150025DOIback to text
- 92 articleMonge-Ampère-Kantorovitch (MAK) reconstruction of the eary universe.Nature4172602002back to text
- 93 articleA stochastic control approach to No-Arbitrage bounds given marginals, with an application to Loopback options.submitted to Annals of Applied Probability2011back to text
- 94 articleThe geometry of optimal transportation.Acta Math.17721996, 113--161URL: http://dx.doi.org/10.1007/BF02392620DOIback to text
- 95 articleGaspard Monge, Le mémoire sur les déblais et les remblais.Image des mathématiques, CNRS2012, URL: http://images.math.cnrs.fr/Gaspard-Monge,1094.htmlback to text
-
96
articlePath Dependent Optimal Transport and Model Calibration on Exotic Derivatives.SSRN Electron.~J.Available at
doi:10.2139/ssrn.3302384 01 2018DOIback to text - 97 articleNotes on a PDE System for Biological Network Formation.138June 2016, 127--155URL: https://www.sciencedirect.com/science/article/pii/S0362546X15004344DOIback to text
- 98 articleSoliton dynamics in computational anatomy.NeuroImage232004, S170--S178back to text
- 99 incollectionThe mathematical theory of frontogenesis.Annual review of fluid mechanics, Vol. 14Palo Alto, CAAnnual Reviews1982, 131--151back to text
- 100 articleThe back-and-forth method for Wasserstein gradient flows.ESAIM: Control, Optimisation and Calculus of Variations272021, 28back to text
- 101 articleA fast approach to optimal transport: The back-and-forth method.Numer. Math.1462020, 513--544URL: https://doi.org/10.1007/s00211-020-01154-8DOIback to text
- 102 articleOn explosions of solutions to a system of partial differential equations modelling chemotaxis.Trans. Amer. Math. Soc.32921992, 819--824URL: http://dx.doi.org/10.2307/2153966DOIback to text
- 103 articleThe variational formulation of the Fokker-Planck equation.SIAM J. Math. Anal.2911998, 1--17back to text
- 104 articleOn the translocation of masses.C. R. (Doklady) Acad. Sci. URSS (N.S.)371942, 199--201back to text
- 105 articleMinimizing within convex bodies using a convex hull method.SIAM Journal on Optimization162January 2005, 368--379HALback to text
- 106 articleMean field games.Jpn. J. Math.212007, 229--260URL: http://dx.doi.org/10.1007/s11537-007-0657-8DOIback to text
- 107 articleOptimal Transport Approximation of 2-Dimensional Measures.SIAM Journal on Imaging Sciences122January 2019, 762--787URL: https://epubs.siam.org/doi/10.1137/18M1193736DOIback to text
- 108 articleA survey of the Schrödinger problem and some of its connections with optimal transport.Discrete Contin. Dyn. Syst.3442014, 1533--1574URL: http://dx.doi.org/10.3934/dcds.2014.34.1533DOIback to text
- 109 articleActive sets, nonsmoothness, and sensitivity.SIAM Journal on Optimization1332003, 702--725back to text
- 110 articleOptimal transportation meshfree approximation schemes for Fluid and plastic Flows.Int. J. Numer. Meth. Engng 83:1541--579832010, 1541--1579back to text
- 111 articleA fully nonlinear version of the incompressible Euler equations: the semigeostrophic system.SIAM J. Math. Anal.3832006, 795--823 (electronic)back to text
- 112 articleNumerical solution of the Monge-Ampére equation by a Newton's algorithm.C. R. Math. Acad. Sci. Paris34042005, 319--324back to text
- 113 bookA wavelet tour of signal processing.Elsevier/Academic Press, Amsterdam2009back to text
- 114 articleA macroscopic crowd motion model of gradient flow type.Math. Models Methods Appl. Sci.20102010, 1787--1821URL: http://dx.doi.org/10.1142/S0218202510004799DOIback to text
- 115 articleA multiscale approach to optimal transport.Computer Graphics Forum3052011, 1583--1592back to text
- 116 articleGeodesic Shooting for Computational Anatomy.Journal of Mathematical Imaging and Vision242March 2006, 209--228URL: http://dx.doi.org/10.1007/s10851-005-3624-0back to text
- 117 articleAdaptive, Anisotropic and Hierarchical cones of Discrete Convex functions.Numerische Mathematik132435 pages, 11 figures. (Second version fixes a small bug in Lemma 3.2. Modifications are anecdotic.)2016, 807--853HALback to text
- 118 articleA Modica-Mortola Approximation for Branched Transport and Applications.Archive for Rational Mechanics and Analysis2011July 2011, 115--142URL: http://link.springer.com/10.1007/s00205-011-0402-6DOIback to text
- 119 articleOptimal transport: From moving soil to same-sex marriage.CMS Notes452013, 14--15back to text
- 120 articleUniqueness and Monge Solutions in the Multimarginal Optimal Transportation Problem.SIAM Journal on Mathematical Analysis4362011, 2758--2775back to text
- 121 miscRegularized Optimal Transport is Ground Cost Adversarial.2020back to text
- 122 articleA Generalized Forward-Backward Splitting.SIAM Journal on Imaging Sciences632013, 1199--1226URL: http://hal.archives-ouvertes.fr/hal-00613637/DOIback to text
- 123 articleNonlinear total variation based noise removal algorithms.Physica D: Nonlinear Phenomena6011992, 259--268URL: http://dx.doi.org/10.1016/0167-2789(92)90242-Fback to text
- 124 articleConvolutional Wasserstein Distances: Efficient Optimal Transportation on Geometric Domains.ACM Transaction on Graphics, Proc. SIGGRAPH'15to appear2015back to textback to text
- 125 articleRegression shrinkage and selection via the Lasso.Journal of the Royal Statistical Society. Series B. Methodological5811996, 267--288back to textback to text
- 126 unpublishedDynamical Programming for off-the-grid dynamic Inverse Problems.December 2022, working paper or preprintHALback to text
- 127 inproceedingsParameter tuning and model selection in optimal transport with semi-dual Brenier formulation.NeurIPSNew Orleans, France2022HALback to text
- 128 articleModel Selection with Piecewise Regular Gauges.Information and Inferenceto appear2015, URL: http://hal.archives-ouvertes.fr/hal-00842603/back to text
- 129 bookOptimal transport.338Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]Old and newBerlinSpringer-Verlag2009, xxii+973URL: http://dx.doi.org/10.1007/978-3-540-71050-9DOIback to text
- 130 bookTopics in optimal transportation.58Graduate Studies in MathematicsAmerican Mathematical Society, Providence, RI2003, xvi+370back to text
- 131 articleOn the design of a reflector antenna. II.Calc. Var. Partial Differential Equations2032004, 329--341URL: http://dx.doi.org/10.1007/s00526-003-0239-4DOIback to text
- 132 articleA continuum mechanical approach to geodesics in shape space.International Journal of Computer Vision9332011, 293--318back to text
- 133 articleSparse representation for computer vision and pattern recognition.Proceedings of the IEEE9862010, 1031--1044back to text
- 134 phdthesisEntropic Unbalanced Optimal Transport: Application to Full-Waveform Inversion and Numerical Illustration.Université de ParisDecember 2021HALback to text