We investigate the convergence rate of multi-marginal optimal transport costs that are regularized with the Boltzmann–Shannon entropy, as the noise parameter $\varepsilon $ tends to $0$. We establish lower and upper bounds on the difference with the unregularized cost of the form $C\varepsilon \log (1/\varepsilon )+O(\varepsilon )$ for some explicit dimensional constants C depending on the marginals and on the ground cost, but not on the optimal transport plans themselves. Upper bounds are obtained for Lipschitz costs or locally semiconcave costs for a finer estimate, and lower bounds for $\mathscr {C}^2$ costs satisfying some signature condition on the mixed second derivatives that may include degenerate costs, thus generalizing results previously in the two marginals case and for nondegenerate costs. We obtain in particular matching bounds in some typical situations where the optimal plan is deterministic.

In modelling joint probability distributions it is often desirable to incorporate standard marginal distributions and match a set of key observed mixed moments. At the same time it may also be prudent to avoid additional unwarranted assumptions. The problem is to find the least ordered distribution that respects the prescribed constraints. In this paper we will construct a suitable joint probability distribution by finding the checkerboard copula of maximum entropy that allows us to incorporate the appropriate marginal distributions and match the nominated set of observed moments.

In this paper we study the compact and convex setsK ⊆ Ω ⊆ ℝ2 that minimize

\begin{equation*} \int_{\Omega} \dist(\x,K) \,{\rm d}\x +\lambda_1 {\rm Vol}(K)+\lambda_2 {\rm Per}(K) \end{equation*}${\mathrm{\int }}_{\mathit{\Omega }}\mathrm{dist}\mathrm{\left(}{x}\mathit{,K}\mathrm{\right)}\mathrm{d}{x}\mathrm{+}{\mathit{\lambda }}_{\mathrm{1}}\mathrm{Vol}\mathrm{\left(}\mathit{K}\mathrm{\right)}\mathrm{+}{\mathit{\lambda }}_{\mathrm{2}}\mathrm{Per}\mathrm{\left(}\mathit{K}\mathrm{\right)}$

The aim of this paper is to develop a crowd motion model designed to handle highly packed situations. The model we propose rests on two principles: we first define a spontaneous velocity which corresponds to the velocity each individual would like to have in the absence of other people. The actual velocity is then computed as the projection of the spontaneous velocity onto the set of admissible velocities (i.e. velocities which do not violate the non-overlapping constraint). We describe here the underlying mathematical framework, and we explain how recent results by J.F. Edmond and L. Thibault on the sweeping process by uniformly prox-regular sets can be adapted to handle this situation in terms of well-posedness. We propose a numerical scheme for this contact dynamics model, based on a prediction-correction algorithm. Numerical illustrations are finally presented and discussed.

We apply conjugate duality to establish the existence of optimal portfolios in an asset-allocation problem, with the goal of minimizing the variance of the final wealth which results from trading over a fixed, finite horizon in a continuous-time, complete market, subject to the constraints that the expected final wealth equal a specified target value and the portfolio of the investor (defined by the dollar amount invested in each stock) take values in a given closed, convex set. The asset prices are modelled by Itô processes, for which the market parameters are random processes adapted to the information filtration available to the investor. We synthesize a dual optimization problem and establish a set of optimality relations, similar to the Euler-Lagrange and transversality relations of calculus of variations, giving necessary and sufficient conditions for the given optimization problem and its dual to each have a solution, with zero duality gap. We then solve these relations, to establish the existence of an optimal portfolio.

Approximation theory in the context of probability density function turns out to go beyond the classical idea of orthogonal projection. Special tools have to be designed so as to respect the nonnegativity of the approximate function. We develop here and justify from the theoretical point of view an approximation procedure introduced by Levermore [Levermore, J. Stat. Phys.83 (1996) 1021–1065] and based on an entropy minimization principle under moment constraints. We prove in particular a global existence theorem for such an approximation and derive as a by-product a necessary and sufficient condition for the so-called problem of moment realizability. Applications of the above result are given in kinetic theory: first in the context of Levermore's approach and second to design generalized BGK models for Maxwellian molecules.

A homogeneous real polynomial $p$ is hyperbolic with respect to a given vector $d$ if the univariate polynomial $t\,\mapsto \,p(x\,-\,td)$ has all real roots for all vectors $x$. Motivated by partial differential equations, Gårding proved in 1951 that the largest such root is a convex function of $x$, and showed various ways of constructing new hyperbolic polynomials. We present a powerful new such construction, and use it to generalize Gårding’s result to arbitrary symmetric functions of the roots. Many classical and recent inequalities follow easily. We develop various convex-analytic tools for such symmetric functions, of interest in interior-point methods for optimization problems over related cones.

The computation of glacier movements leads to a system of nonlinear partial differentialequations. The existence and uniqueness of a weak solution is established by using the calculus ofvariations. A discretization by the finite element method is done. Thesolution of the discrete problem is proved to be convergent to the exactsolution. A first simple numerical algorithm is proposed and its convergence numericallystudied.

An important question in the study of moment problems is to determine when a fixed point in ℝn lies in the moment cone of vectors , with μ a nonnegative measure. In associated optimization problems it is also important to be able to distinguish between the interior and boundary of the moment cone. Recent work of Dachuna-Castelle, Gamboa and Gassiat derived elegant computational characterizations for these problems, and for related questions with an upper bound on μ. Their technique involves a probabilistic interpretation and large deviations theory. In this paper a purely convex analytic approach is used, giving a more direct understanding of the underlying duality, and allowing the relaxation of their assumptions.