Convex discretization of functionals involving the Monge-Ampère operator
|
|
- Katherine Wells
- 5 years ago
- Views:
Transcription
1 Convex discretization of functionals involving the Monge-Ampère operator Quentin Mérigot CNRS / Université Paris-Dauphine Joint work with J.D. Benamou, G. Carlier and É. Oudet GeMeCod Conference October 26-30, 2015, Marne-la-Vallée 1
2 2 1. Discretization of the space of convex functions
3 Motivation: variational problems over convex functions Principal-agent problem in economy min φ K,φ 0 [1,2] ( φ(x) x 2 φ(x)) d x K := {φ convex } [Rochet-Choné 98] φ # Leb [1,2] 2 3
4 Motivation: variational problems over convex functions Principal-agent problem in economy min φ K,φ 0 [1,2] ( φ(x) x 2 φ(x)) d x K := {φ convex } [Rochet-Choné 98] φ # Leb [1,2] 2 Gradient flows in Wasserstein space [Otto] [Jordan-Kinderlehrer-Otto 98] min φ KX 1 2τ W2 2(ρ, φ # ρ) + E( φ # ρ) K X := {φ convex ; φ X} 3
5 Motivation: variational problems over convex functions Principal-agent problem in economy min φ K,φ 0 [1,2] ( φ(x) x 2 φ(x)) d x K := {φ convex } [Rochet-Choné 98] φ # Leb [1,2] 2 Gradient flows in Wasserstein space [Otto] [Jordan-Kinderlehrer-Otto 98] min φ KX 1 2τ W2 2(ρ, φ # ρ) + E( φ # ρ) K X := {φ convex ; φ X} Numerical exploration of conjectures in convex geometry, e.g. Meissner. K := { support function of convex bodies } 3
6 Motivation: variational problems over convex functions Principal-agent problem in economy min φ K,φ 0 [1,2] ( φ(x) x 2 φ(x)) d x K := {φ convex } [Rochet-Choné 98] φ # Leb [1,2] 2 Gradient flows in Wasserstein space [Otto] [Jordan-Kinderlehrer-Otto 98] min φ KX 1 2τ W2 2(ρ, φ # ρ) + E( φ # ρ) K X := {φ convex ; φ X} Numerical exploration of conjectures in convex geometry, e.g. Meissner. K := { support function of convex bodies } 3
7 Prior work: Discretizing the Space of Convex Functions Finite elements: piecewise-linear convex functions over a fixed mesh non-density result [Choné-Le Meur 99] h standard grid 4
8 Prior work: Discretizing the Space of Convex Functions Finite elements: piecewise-linear convex functions over a fixed mesh non-density result [Choné-Le Meur 99] h standard grid Convex interpolates: [Carlier-Lachand-Robert-Maury 01] [Ekeland Moreno-Bromberg 10] Def: Given a finite subset P R d, K(P ) := { restriction of convex functions to P } P R d 4
9 Prior work: Discretizing the Space of Convex Functions Finite elements: piecewise-linear convex functions over a fixed mesh non-density result [Choné-Le Meur 99] h standard grid Convex interpolates: [Carlier-Lachand-Robert-Maury 01] [Ekeland Moreno-Bromberg 10] Def: Given a finite subset P R d, K(P ) := { restriction of convex functions to P } good convergence properties under convergence of P P R d 4
10 Prior work: Discretizing the Space of Convex Functions Finite elements: piecewise-linear convex functions over a fixed mesh non-density result [Choné-Le Meur 99] h standard grid Convex interpolates: [Carlier-Lachand-Robert-Maury 01] [Ekeland Moreno-Bromberg 10] Def: Given a finite subset P R d, K(P ) := { restriction of convex functions to P } good convergence properties under convergence of P K(P ) is a finite-dimensional convex polytope defined by Card(P ) 2 linear inequalities P R d adaptive method exterior parameterization [Oberman 14] [Mirebeau 14] [Oudet-M. 14] 4
11 Prior work: Discretizing the Space of Convex Functions Finite elements: piecewise-linear convex functions over a fixed mesh non-density result [Choné-Le Meur 99] h standard grid Convex interpolates: [Carlier-Lachand-Robert-Maury 01] [Ekeland Moreno-Bromberg 10] Def: Given a finite subset P R d, K(P ) := { restriction of convex functions to P } good convergence properties under convergence of P K(P ) is a finite-dimensional convex polytope defined by Card(P ) 2 linear inequalities P R d adaptive method exterior parameterization [Oberman 14] [Mirebeau 14] [Oudet-M. 14] This talk: Using the Monge-Ampère operator to describe K(P ) more conveniently. 4
12 5 2. Space-discretization of Wasserstein gradient flows
13 Background: Optimal transport P 2 (R d ) = prob. measures with finite second moment P 2 ac(r d ) = P 2 (R d ) L 1 (R d ) ν R d p 2 π Wasserstein distance between µ, ν P 2 (R d ), Γ(µ, ν) := {π P(R d R d ); p 1# π = µ, p 2# π = ν} Definition: W 2 2(µ, ν) := min π Γ(µ,ν) x y 2 d π(x, y). µ p 1 R d 6
14 Background: Optimal transport P 2 (R d ) = prob. measures with finite second moment P 2 ac(r d ) = P 2 (R d ) L 1 (R d ) ν R d p 2 π Wasserstein distance between µ, ν P 2 (R d ), Γ(µ, ν) := {π P(R d R d ); p 1# π = µ, p 2# π = ν} Definition: W 2 2(µ, ν) := min π Γ(µ,ν) x y 2 d π(x, y). µ p 1 R d Relation to convex functions: Def: K := convex functions on R d 6
15 Background: Optimal transport P 2 (R d ) = prob. measures with finite second moment P 2 ac(r d ) = P 2 (R d ) L 1 (R d ) ν R d p 2 π Wasserstein distance between µ, ν P 2 (R d ), Γ(µ, ν) := {π P(R d R d ); p 1# π = µ, p 2# π = ν} Definition: W 2 2(µ, ν) := min π Γ(µ,ν) x y 2 d π(x, y). µ p 1 R d Relation to convex functions: Def: K := convex functions on R d Theorem (Brenier): Given µ P 2 ac(r d ) the map φ K φ # µ P 2 (R d ) is surjective and moreover, W 2 2(µ, φ # µ) = R d x φ(x) 2 d µ(x) [Brenier 91] 6
16 Background: Optimal transport P 2 (R d ) = prob. measures with finite second moment P 2 ac(r d ) = P 2 (R d ) L 1 (R d ) ν R d p 2 π Wasserstein distance between µ, ν P 2 (R d ), Γ(µ, ν) := {π P(R d R d ); p 1# π = µ, p 2# π = ν} Definition: W 2 2(µ, ν) := min π Γ(µ,ν) x y 2 d π(x, y). µ p 1 R d Relation to convex functions: Def: K := convex functions on R d Theorem (Brenier): Given µ P 2 ac(r d ) the map φ K φ # µ P 2 (R d ) is surjective and moreover, W 2 2(µ, φ # µ) = R d x φ(x) 2 d µ(x) [Brenier 91] Lagrangian parameterization of P 2 (R d ), as seen from µ P ac 2 (R d ). 6
17 Motivation 1: Crowd Motion Under Congestion JKO scheme for crowd motion with hard congestion: [Maury-Roudneff-Chupin-Santambrogio 10] ρ τ k+1 = min σ P 2 (X) 1 2τ W2 2(ρ τ k, σ) + E(σ) + U(σ) ( ) where X R d is convex and bounded, and U(ν) := 0 if ν P ac and d ν d x <= 1 E(ν) := V (x) d ν(x) R d + if not congestion potential energy 7
18 Motivation 1: Crowd Motion Under Congestion JKO scheme for crowd motion with hard congestion: [Maury-Roudneff-Chupin-Santambrogio 10] ρ τ k+1 = min σ P 2 (X) 1 2τ W2 2(ρ τ k, σ) + E(σ) + U(σ) ( ) where X R d is convex and bounded, and U(ν) := 0 if ν P ac and d ν d x <= 1 E(ν) := V (x) d ν(x) R d + if not congestion potential energy Assuming σ = φ # ρ τ k with φ convex, the Wasserstein term becomes explicit: ( ) min φ 1 2τ R d x φ(x) 2 ρ τ k (x) d x + E( φ #ρ τ k ) + U( φ #ρ τ k ) On the other hand, the constraint becomes strongly nonlinear: U( φ # ρ τ k ) < + det D2 φ(x) ρ k (x) 7
19 Motivation 2: Nonlinear Diffusion ( ) ρ t = div [ρ (U (ρ) + V + W ρ)] ρ(0,.) = ρ 0 ρ(t,.) P ac (R d ) 8
20 Motivation 2: Nonlinear Diffusion ( ) ρ t = div [ρ (U (ρ) + V + W ρ)] ρ(0,.) = ρ 0 ρ(t,.) P ac (R d ) Formally, ( ) can be interpreted as the W 2 -gradient flow of U + E, with U(ν) := U( d ν R d x ) d x if ν P2 ac(r d ) d + if not internal energy, ex: U(r) = r log r = entropy E(ν) := R d V (x) d ν(x) + R d W (x y) d[ν ν](x, y) potential energy interaction energy 8
21 Motivation 2: Nonlinear Diffusion ( ) ρ t = div [ρ (U (ρ) + V + W ρ)] ρ(0,.) = ρ 0 ρ(t,.) P ac (R d ) Formally, ( ) can be interpreted as the W 2 -gradient flow of U + E, with U(ν) := U( d ν R d x ) d x if ν P2 ac(r d ) d + if not internal energy, ex: U(r) = r log r = entropy E(ν) := R d V (x) d ν(x) + R d W (x y) d[ν ν](x, y) potential energy interaction energy JKO time discrete scheme: for τ > 0, [Jordan-Kinderlehrer-Otto 98] ρ τ k+1 = arg min σ P(R d ) 1 2τ W 2(ρ τ k, σ)2 + U(σ) + E(σ) 8
22 Motivation 2: Nonlinear Diffusion ( ) ρ t = div [ρ (U (ρ) + V + W ρ)] ρ(0,.) = ρ 0 ρ(t,.) P ac (R d ) Formally, ( ) can be interpreted as the W 2 -gradient flow of U + E, with U(ν) := U( d ν R d x ) d x if ν P2 ac(r d ) d + if not internal energy, ex: U(r) = r log r = entropy E(ν) := R d V (x) d ν(x) + R d W (x y) d[ν ν](x, y) potential energy interaction energy JKO time discrete scheme: for τ > 0, [Jordan-Kinderlehrer-Otto 98] ρ τ k+1 = arg min σ P(R d ) 1 2τ W 2(ρ τ k, σ)2 + U(σ) + E(σ) 8 Many applications: porous medium equation, cell movement via chemotaxis, generalization to higher order PDEs, etc.
23 Displacement Convex Setting min σ P(R d ) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) min φ K 1 2τ W2 2(µ, φ # µ) + U( φ # µ) + E( φ # µ) 9
24 Displacement Convex Setting min σ P(R d ) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) min φ K 1 2τ W2 2(µ, φ # µ) + U( φ # µ) + E( φ # µ) When is this minimization problem convex in φ? 9
25 Displacement Convex Setting min σ P(R d ) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) min φ K 1 2τ W2 2(µ, φ # µ) + U( φ # µ) + E( φ # µ) When is this minimization problem convex in φ? Formally, the pushforward σ = φ # µ satisfies the Monge-Ampère equation σ( φ(x))ma[φ](x) = µ(x) MA[φ](x) := det(d 2 φ(x)) 9
26 Displacement Convex Setting min σ P(R d ) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) min φ K 1 2τ W2 2(µ, φ # µ) + U( φ # µ) + E( φ # µ) When is this minimization problem convex in φ? Formally, the pushforward σ = φ # µ satisfies the Monge-Ampère equation σ( φ(x))ma[φ](x) = µ(x) MA[φ](x) := det(d 2 φ(x)) With a change-of-variable formula, U( φ # µ) = R d U(σ(x)) d x = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x 9
27 Displacement Convex Setting min σ P(R d ) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) min φ K 1 2τ W2 2(µ, φ # µ) + U( φ # µ) + E( φ # µ) When is this minimization problem convex in φ? Formally, the pushforward σ = φ # µ satisfies the Monge-Ampère equation σ( φ(x))ma[φ](x) = µ(x) MA[φ](x) := det(d 2 φ(x)) With a change-of-variable formula, U( φ # µ) = R d U(σ(x)) d x = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x Using concavity of det 1/d over SDP, this gives: Theorem: Given µ P 2 ac(r 2 ), for the functional φ K U( φ #µ ) to be convex, it suffices that r > 0 r d U(r d ) is convex non-increasing and U(0) = 0. [McCann 94] 9
28 Discretizing the Space of Convex Functions We replace K with K X := {φ convex ; φ X} Definition: Given a finite subset P X, let K X (P ) := {ψ P ; ψ K X } R P 10
29 Discretizing the Space of Convex Functions We replace K with K X := {φ convex ; φ X} Definition: Given a finite subset P X, let K X (P ) := {ψ P ; ψ K X } R P = φ : P R P X 10
30 Discretizing the Space of Convex Functions We replace K with K X := {φ convex ; φ X} Definition: Given a finite subset P X, let K X (P ) := {ψ P ; ψ K X } R P K X envelope of a function φ on P : φ : P R ˆφ : R d R ˆφ = φ : P R ˆφ := max{ψ; ψ K X and ψ P φ} K X P X 10
31 Discretizing the Space of Convex Functions We replace K with K X := {φ convex ; φ X} Definition: Given a finite subset P X, let K X (P ) := {ψ P ; ψ K X } R P K X envelope of a function φ on P : φ : P R ˆφ : R d R ˆφ = φ : P R ˆφ := max{ψ; ψ K X and ψ P φ} K X P X Given µ = p P µ pδ p on P and φ K X (P ) we need to define φ # µ. 10
32 Discretizing the Space of Convex Functions We replace K with K X := {φ convex ; φ X} Definition: Given a finite subset P X, let K X (P ) := {ψ P ; ψ K X } R P K X envelope of a function φ on P : φ : P R ˆφ : R d R ˆφ = φ : P R ˆφ := max{ψ; ψ K X and ψ P φ} K X V p P X Given µ = p P µ pδ p on P and φ K X (P ) we need to define φ # µ. Definition: For φ K X (P ), let V p := ˆφ(p) and V p p X 10
33 Discretizing the Space of Convex Functions We replace K with K X := {φ convex ; φ X} Definition: Given a finite subset P X, let K X (P ) := {ψ P ; ψ K X } R P K X envelope of a function φ on P : φ : P R ˆφ : R d R ˆφ = φ : P R ˆφ := max{ψ; ψ K X and ψ P φ} K X V p P X Given µ = p P µ pδ p on P and φ K X (P ) we need to define φ # µ. Definition: For φ K X (P ), let V p := ˆφ(p) and V p p X Gφ # µ := p P µ p Leb(V p ) Leb V p P ac (X) 10
34 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p 11
35 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p Internal energy of the push-forward Gφ # µ = p P (µ p/ma[φ](p)) Leb ˆφ(p) : U(Gφ # µ) = U(Gφ # µ(x)) d x, = p P ˆφ(p) U (µ p/ma[φ](p))) d x 11
36 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p Internal energy of the push-forward Gφ # µ = p P (µ p/ma[φ](p)) Leb ˆφ(p) : U(Gφ # µ) = U(Gφ # µ(x)) d x, = p P ˆφ(p) U (µ p/ma[φ](p))) d x = p P U (µ p/ma[φ](p)) MA[φ](p) (Compare to the continuous formula: U( φ # µ) = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x) 11
37 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p Internal energy of the push-forward Gφ # µ = p P (µ p/ma[φ](p)) Leb ˆφ(p) : U(Gφ # µ) = U(Gφ # µ(x)) d x, = p P ˆφ(p) U (µ p/ma[φ](p))) d x = p P U (µ p/ma[φ](p)) MA[φ](p) (Compare to the continuous formula: U( φ # µ) = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x) Concavity of discrete MA operator: φ K(P ) MA[φ](p) 1/d is concave 11
38 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p Internal energy of the push-forward Gφ # µ = p P (µ p/ma[φ](p)) Leb ˆφ(p) : U(Gφ # µ) = U(Gφ # µ(x)) d x, = p P ˆφ(p) U (µ p/ma[φ](p))) d x = p P U (µ p/ma[φ](p)) MA[φ](p) (Compare to the continuous formula: U( φ # µ) = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x) Concavity of discrete MA operator: φ K(P ) MA[φ](p) 1/d is concave Leb( ˆφ t ) 1/d Leb((1 t) ˆφ 0 (p) t ˆφ 1 (p)) 1/d. 11
39 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p Internal energy of the push-forward Gφ # µ = p P (µ p/ma[φ](p)) Leb ˆφ(p) : U(Gφ # µ) = U(Gφ # µ(x)) d x, = p P ˆφ(p) U (µ p/ma[φ](p))) d x = p P U (µ p/ma[φ](p)) MA[φ](p) (Compare to the continuous formula: U( φ # µ) = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x) Concavity of discrete MA operator: φ K(P ) MA[φ](p) 1/d is concave Brunn-Minkowski Leb( ˆφ t ) 1/d Leb((1 t) ˆφ 0 (p) t ˆφ 1 (p)) 1/d. (1 t)leb( ˆφ 0 (p)) 1/d + tleb( ˆφ 1 (p)) 1/d. 11
40 Internal Energy and Discrete Monge-Ampère Operator Discrete Monge-Ampère operator: MA[φ](p) := Leb( ˆφ(p)). 11 Related to Alexandrov s Monge-Ampère measure of ˆφ, i.e. p P MA[φ](p)δ p Internal energy of the push-forward Gφ # µ = p P (µ p/ma[φ](p)) Leb ˆφ(p) : U(Gφ # µ) = U(Gφ # µ(x)) d x, = p P ˆφ(p) U (µ p/ma[φ](p))) d x = p P U (µ p/ma[φ](p)) MA[φ](p) (Compare to the continuous formula: U( φ # µ) = R d U ( ) µ(x) MA[φ](x) MA[φ](x) d x) Concavity of discrete MA operator: φ K(P ) MA[φ](p) 1/d is concave Brunn-Minkowski Leb( ˆφ t ) 1/d Leb((1 t) ˆφ 0 (p) t ˆφ 1 (p)) 1/d. (1 t)leb( ˆφ 0 (p)) 1/d + tleb( ˆφ 1 (p)) 1/d. Proposition: Under McCann s assumption, φ K X (P ) U(Gφ # µ P ) is convex.
41 Convex Space-Discretization of One Gradient Step Theorem: Under McCann s hypotheses the discrete problem min φ KX (P ) 1 2τ W2 2(µ, Hφ # µ) + E(Hφ # µ) + U(Gφ # µ) is convex, and the minimum is unique if r d U(r d ) is strictly convex. [Benamou, Carlier, M., Oudet 14] 12
42 Convex Space-Discretization of One Gradient Step Theorem: Under McCann s hypotheses the discrete problem min φ KX (P ) 1 2τ W2 2(µ, Hφ # µ) + E(Hφ # µ) + U(Gφ # µ) is convex, and the minimum is unique if r d U(r d ) is strictly convex. [Benamou, Carlier, M., Oudet 14] Unfortunately, two different definitions of push-forward seem necessary to get a convex discretisation (i.e. φ E(Gφ # µ) is non-convex). 12
43 Convex Space-Discretization of One Gradient Step Theorem: Under McCann s hypotheses the discrete problem min φ KX (P ) 1 2τ W2 2(µ, Hφ # µ) + E(Hφ # µ) + U(Gφ # µ) is convex, and the minimum is unique if r d U(r d ) is strictly convex. [Benamou, Carlier, M., Oudet 14] Unfortunately, two different definitions of push-forward seem necessary to get a convex discretisation (i.e. φ E(Gφ # µ) is non-convex). If lim r U(r)/r = +, the internal energy U is a barrier for convexity, i.e. U(Gφ # µ P ) < + = p P, MA[φ](p) = Leb( ˆφ(p)) > 0 = φ is in the interior of K X (P ). NB: P non-linear constraints vs P 2 linear constraints 12
44 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) 13
45 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) (C1) E continuous, U l.s.c on (P(X), W 2 ) 13
46 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) (C1) E continuous, U l.s.c on (P(X), W 2 ) (C2) U(ρ) = U(ρ(x)) d x, where U is convex with superlinear growth. 13
47 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) (C1) E continuous, U l.s.c on (P(X), W 2 ) (C2) U(ρ) = U(ρ(x)) d x, where U is convex with superlinear growth. McCann s condition for displacement-convexity 13
48 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) (C1) E continuous, U l.s.c on (P(X), W 2 ) (C2) U(ρ) = U(ρ(x)) d x, where U is convex with superlinear growth. McCann s condition for displacement-convexity Theorem: Let P n X finite, µ n P(P n ) with lim W 2 (µ n, µ) = 0, and: min φ KX (P n ) 1 2τ W 2(µ n, Hφ # µ n ) + E(Hφ # µ n ) + U(Gφ # µ n ) ( ) n If φ n minimizes ( ) n, then (Gφ n #µ n ) is a minimizing sequence for ( ). [Benamou, Carlier, M., Oudet 14] 13
49 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) (C1) E continuous, U l.s.c on (P(X), W 2 ) (C2) U(ρ) = U(ρ(x)) d x, where U is convex with superlinear growth. McCann s condition for displacement-convexity Theorem: Let P n X finite, µ n P(P n ) with lim W 2 (µ n, µ) = 0, and: min φ KX (P n ) 1 2τ W 2(µ n, Hφ # µ n ) + E(Hφ # µ n ) + U(Gφ # µ n ) ( ) n If φ n minimizes ( ) n, then (Gφ n #µ n ) is a minimizing sequence for ( ). Proof relies on Caffarelli s regularity theorem. [Benamou, Carlier, M., Oudet 14] 13
50 Convergence of the Space-Discretization Setting: X bounded and convex, µ P ac (X) with density c 1 ρ c min σ P(X) 1 2τ W2 2(µ, σ) + E(σ) + U(σ) ( ) (C1) E continuous, U l.s.c on (P(X), W 2 ) (C2) U(ρ) = U(ρ(x)) d x, where U is convex with superlinear growth. McCann s condition for displacement-convexity Theorem: Let P n X finite, µ n P(P n ) with lim W 2 (µ n, µ) = 0, and: min φ KX (P n ) 1 2τ W 2(µ n, Hφ # µ n ) + E(Hφ # µ n ) + U(Gφ # µ n ) ( ) n If φ n minimizes ( ) n, then (Gφ n #µ n ) is a minimizing sequence for ( ). Proof relies on Caffarelli s regularity theorem. [Benamou, Carlier, M., Oudet 14] When P n is a regular grid, there is an alternative (and quantitative) argument 13
51 14 3. Numerical results
52 Monge-Ampère operator and Computational Geometry U(Gφ # µ P ) = p P U (µ p/ma[φ](p)) MA[φ](p) with MA[φ](p) := Leb( ˆφ(p)). Goal: fast computation of MA[φ] and its 1st/2nd derivatives w.r.t φ 15
53 Monge-Ampère operator and Computational Geometry U(Gφ # µ P ) = p P U (µ p/ma[φ](p)) MA[φ](p) with MA[φ](p) := Leb( ˆφ(p)). Goal: fast computation of MA[φ] and its 1st/2nd derivatives w.r.t φ We rely on the notion of Laguerre diagram in computational geometry Definition: Given a function φ : P R, Lag φ P (p) := {y Rd ; q P, φ(q) φ(p) + q p y } Lag φ P (p) X p 15
54 Monge-Ampère operator and Computational Geometry U(Gφ # µ P ) = p P U (µ p/ma[φ](p)) MA[φ](p) with MA[φ](p) := Leb( ˆφ(p)). Goal: fast computation of MA[φ] and its 1st/2nd derivatives w.r.t φ We rely on the notion of Laguerre diagram in computational geometry Definition: Given a function φ : P R, Lag φ P (p) := {y Rd ; q P, φ(q) φ(p) + q p y } Lag φ P (p) X For φ K X (P ), ˆφ(p) = Lag φ P (p) X. p 15
55 Monge-Ampère operator and Computational Geometry U(Gφ # µ P ) = p P U (µ p/ma[φ](p)) MA[φ](p) with MA[φ](p) := Leb( ˆφ(p)). Goal: fast computation of MA[φ] and its 1st/2nd derivatives w.r.t φ We rely on the notion of Laguerre diagram in computational geometry Definition: Given a function φ : P R, Lag φ P (p) := {y Rd ; q P, φ(q) φ(p) + q p y } Lag φ P (p) X For φ K X (P ), ˆφ(p) = Lag φ P (p) X. p For φ(p) = p 2 /2, one gets the Voronoi cell: Lag φ P (p) := {y; q P, q y 2 p y 2 } 15
56 Monge-Ampère operator and Computational Geometry U(Gφ # µ P ) = p P U (µ p/ma[φ](p)) MA[φ](p) with MA[φ](p) := Leb( ˆφ(p)). Goal: fast computation of MA[φ] and its 1st/2nd derivatives w.r.t φ We rely on the notion of Laguerre diagram in computational geometry Definition: Given a function φ : P R, Lag φ P (p) := {y Rd ; q P, φ(q) φ(p) + q p y } Lag φ P (p) X For φ K X (P ), ˆφ(p) = Lag φ P (p) X. p For φ(p) = p 2 /2, one gets the Voronoi cell: Lag φ P (p) := {y; q P, q y 2 p y 2 } Computation in time O( P log P ) in 2D 15
57 Example: semi-lagrangian scheme for crowd motion Gradient flow model of crowd motion with congestion, with a JKO scheme: [Maury-Roudneff-Chupin-Santambrogio 10] µ k+1 = min ν P(X) 1 2τ W2 2(µ k, ν) + E(ν) + U(ν) E(ν) := V (x) d ν(x) X { 0 if d ν/ d x 1, U(ν) := + if not 16
58 Example: semi-lagrangian scheme for crowd motion Gradient flow model of crowd motion with congestion, with a JKO scheme: [Maury-Roudneff-Chupin-Santambrogio 10] µ k+1 = min ν P(X) 1 2τ W2 2(µ k, ν) + E(ν) + U(ν) E(ν) := V (x) d ν(x) X { 0 if d ν/ d x 1, U(ν) := + if not Convex optimization problem if V is λ-convex (V + λ. 2 convex) and τ λ/2. We solve this problem with a relaxed hard congestion term: 6 r r α log(1 r 1/d ) U α (ρ) := ρ(x) α log(1 ρ(x) 1/d ) d x 4 2 α=1 α=5 α=10 16 r = 0 r = 1
59 Example: semi-lagrangian scheme for crowd motion Initial density on X = [ 2, 2] 2 Potential 0 1 (2, 2) ( 2, 2) P = regular grid. V (x) = x (2, 0) exp( 5 x 2 /2) Algorithm: Input: µ 0 P(P ), τ > 0, α > 0, β 1. For k {0,..., T } φ arg min φ K G X (P k ) 1 2τ W2 2(µ k, Hφ # µ k ) + E(Hφ # µ k ) + αu β (Gφ # µ k ) ν G φ# µ k ; µ k+1 projection of ν [ 2,2) [ 2,2] on P. 17
60 4. Application to Geometry? Ongoing work with J.M. Mirebeau 18
61 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. 19
62 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. interior point method for shape optimization problem: Minkowski, Meissner, etc. 19
63 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. Definition: Given a convex body K, h K : u S d 1 max p K u p. 19
64 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. Definition: Given a convex body K, h K : u S d 1 max p K u p. K (P ) := {h K P ; K bounded convex body } n 1 n 2 n 5 n 3 n 4 P = {n 1,..., n N } S d 1 R d 19
65 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. Definition: Given a convex body K, h K : u S d 1 max p K u p. K (P ) := {h K P ; K bounded convex body } K(h) n 1 n 2 h 1 h 2 n 5 n 3 h 5 h 3 h 4 n 4 P = {n 1,..., n N } S d 1 R d K(h) := N i=1 {x; x n i h i } 19
66 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. Definition: Given a convex body K, h K : u S d 1 max p K u p. K (P ) := {h K P ; K bounded convex body } K(h) n 1 n 2 h 1 h 2 n 5 n 3 h 5 h 3 h 4 n 4 P = {n 1,..., n N } S d 1 R d K(h) := N i=1 {x; x n i h i } Prop: Φ(h) := N i=1 log(hd 1 (ith face of K(h))) is a convex barrier for K (P ). 19
67 Support Functions of Convex Bodies Objective: Logarithmic barrier for the space of support functions of convex bodies. Definition: Given a convex body K, h K : u S d 1 max p K u p. K (P ) := {h K P ; K bounded convex body } K(h) n 1 n 2 h 1 h 2 n 5 n 3 h 5 h 3 h 4 n 4 P = {n 1,..., n N } S d 1 R d K(h) := N i=1 {x; x n i h i } Prop: Φ(h) := N i=1 log(hd 1 (ith face of K(h))) is a convex barrier for K (P ). Extension to other spaces, such as c-convex functions in optimal transport. 19
68 Application: convex bodies with constant-width Definition: The width of K in direction u is w K (u) = h K (u) + h K ( u). K has constant width α if w K (u) = α. 20
69 Application: convex bodies with constant-width α Definition: The width of K in direction u is w K (u) = h K (u) + h K ( u). K has constant width α if w K (u) = α. Example: Reuleaux triangle in R 2 20
70 Application: convex bodies with constant-width α Definition: The width of K in direction u is w K (u) = h K (u) + h K ( u). K has constant width α if w K (u) = α. Example: Reuleaux triangle in R 2 Theorem: Reuleaux triangles minimize the area over convex sets of the plane with constant width α. [Blaschke-Lebesgue] 20
71 Application: convex bodies with constant-width Definition: The width of K in direction u is α w K (u) = h K (u) + h K ( u). K has constant width α if w K (u) = α. Example: Reuleaux triangle in R 2 Theorem: Reuleaux triangles minimize the area over convex sets of the plane with constant width α. [Blaschke-Lebesgue] Bonnensen-Fenchel conjecture in R 3 : Meissner s body minimize the volume among convex sets with fixed constant width. 20
72 Application: convex bodies with constant-width Definition: The width of K in direction u is α w K (u) = h K (u) + h K ( u). K has constant width α if w K (u) = α. Example: Reuleaux triangle in R 2 Theorem: Reuleaux triangles minimize the area over convex sets of the plane with constant width α. [Blaschke-Lebesgue] Bonnensen-Fenchel conjecture in R 3 : Meissner s body minimize the volume among convex sets with fixed constant width. few theoretical results, but no proof of first-order optimality. a role to be played by numerical experimentation? [Lachand-Robert, Oudet] 20
73 Application: convex bodies with constant-width Definition: The width of K in direction u is α w K (u) = h K (u) + h K ( u). K has constant width α if w K (u) = α. Example: Reuleaux triangle in R 2 Theorem: Reuleaux triangles minimize the area over convex sets of the plane with constant width α. [Blaschke-Lebesgue] Bonnensen-Fenchel conjecture in R 3 : Meissner s body minimize the volume among convex sets with fixed constant width. few theoretical results, but no proof of first-order optimality. a role to be played by numerical experimentation? [Lachand-Robert, Oudet] 20
74 21 Appendix
75 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) Cost function: c : X Y R 22
76 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) K c (P ) := {φ P ; φ : X R is c-concave } Cost function: c : X Y R 22
77 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) K c (P ) := {φ P ; φ : X R is c-concave } Cost function: c : X Y R Characterization of the costs functions such that the space K c is convex Application: generalized principal-agent problem [Figalli, Kim, McCann 10] 22
78 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) K c (P ) := {φ P ; φ : X R is c-concave } Cost function: c : X Y R Characterization of the costs functions such that the space K c is convex Application: generalized principal-agent problem [Figalli, Kim, McCann 10] Under the same hypotheses, there exists a convex logarithmic barrier for K c (P ): 22
79 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) K c (P ) := {φ P ; φ : X R is c-concave } Cost function: c : X Y R Characterization of the costs functions such that the space K c is convex Application: generalized principal-agent problem [Figalli, Kim, McCann 10] Under the same hypotheses, there exists a convex logarithmic barrier for K c (P ): c-laguerre diagram: Lag ψ c (p) = {x X; q P, c(x, p) + ψ(p) c(x, q) + ψ(q)} Lag ψ c (p) X 22
80 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) K c (P ) := {φ P ; φ : X R is c-concave } Cost function: c : X Y R Characterization of the costs functions such that the space K c is convex Application: generalized principal-agent problem [Figalli, Kim, McCann 10] Under the same hypotheses, there exists a convex logarithmic barrier for K c (P ): c-laguerre diagram: Lag ψ c (p) = {x X; q P, c(x, p) + ψ(p) c(x, q) + ψ(q)} exp 1 p (Lag ψ c (p)) R d is convex exp p := [ c(., p)] 1 Lag ψ c (p) X 22
81 c-concave Functions Definition: φ : X R is c-concave if ψ : Y R s.t. φ = min y c(x, y) + ψ(y) K c (P ) := {φ P ; φ : X R is c-concave } Cost function: c : X Y R Characterization of the costs functions such that the space K c is convex Application: generalized principal-agent problem [Figalli, Kim, McCann 10] Under the same hypotheses, there exists a convex logarithmic barrier for K c (P ): c-laguerre diagram: Lag ψ c (p) = {x X; q P, c(x, p) + ψ(p) c(x, q) + ψ(q)} exp 1 p (Lag ψ c (p)) R d is convex exp p := [ c(., p)] 1 Lag ψ c (p) X Prop: Φ(h) := N p P log(hd (exp 1 p (Lag ψ c (p))) is a convex barrier for K c (P ). 22
Convex discretization of functionals involving the Monge-Ampère operator
Convex discretization of functionals involving the Monge-Ampère operator Quentin Mérigot CNRS / Université Paris-Dauphine Joint work with J.D. Benamou, G. Carlier and É. Oudet Workshop on Optimal Transport
More informationLagrangian discretization of incompressibility, congestion and linear diffusion
Lagrangian discretization of incompressibility, congestion and linear diffusion Quentin Mérigot Joint works with (1) Thomas Gallouët (2) Filippo Santambrogio Federico Stra (3) Hugo Leclerc Seminaire du
More informationDiscrete Ricci curvature via convexity of the entropy
Discrete Ricci curvature via convexity of the entropy Jan Maas University of Bonn Joint work with Matthias Erbar Simons Institute for the Theory of Computing UC Berkeley 2 October 2013 Starting point McCann
More informationNumerical Optimal Transport and Applications
Numerical Optimal Transport and Applications Gabriel Peyré Joint works with: Jean-David Benamou, Guillaume Carlier, Marco Cuturi, Luca Nenna, Justin Solomon www.numerical-tours.com Histograms in Imaging
More informationThe semi-geostrophic equations - a model for large-scale atmospheric flows
The semi-geostrophic equations - a model for large-scale atmospheric flows Beatrice Pelloni, University of Reading with M. Cullen (Met Office), D. Gilbert, T. Kuna INI - MFE Dec 2013 Introduction - Motivation
More informationGradient Flows: Qualitative Properties & Numerical Schemes
Gradient Flows: Qualitative Properties & Numerical Schemes J. A. Carrillo Imperial College London RICAM, December 2014 Outline 1 Gradient Flows Models Gradient flows Evolving diffeomorphisms 2 Numerical
More informationNumerical Methods for Optimal Transportation 13frg167
Numerical Methods for Optimal Transportation 13frg167 Jean-David Benamou (INRIA) Adam Oberman (McGill University) Guillaume Carlier (Paris Dauphine University) July 2013 1 Overview of the Field Optimal
More informationA few words about the MTW tensor
A few words about the Ma-Trudinger-Wang tensor Université Nice - Sophia Antipolis & Institut Universitaire de France Salah Baouendi Memorial Conference (Tunis, March 2014) The Ma-Trudinger-Wang tensor
More informationOn the regularity of solutions of optimal transportation problems
On the regularity of solutions of optimal transportation problems Grégoire Loeper April 25, 2008 Abstract We give a necessary and sufficient condition on the cost function so that the map solution of Monge
More informationA Modest Proposal for MFG with Density Constraints
A Modest Proposal for MFG with Density Constraints Mean Field Games and related topics, Rome, May 011 Filippo Santambrogio October 31, 011 Abstract We consider a typical problem in Mean Field Games: the
More informationMinimal time mean field games
based on joint works with Samer Dweik and Filippo Santambrogio PGMO Days 2017 Session Mean Field Games and applications EDF Lab Paris-Saclay November 14th, 2017 LMO, Université Paris-Sud Université Paris-Saclay
More informationMoment Measures. Bo az Klartag. Tel Aviv University. Talk at the asymptotic geometric analysis seminar. Tel Aviv, May 2013
Tel Aviv University Talk at the asymptotic geometric analysis seminar Tel Aviv, May 2013 Joint work with Dario Cordero-Erausquin. A bijection We present a correspondence between convex functions and Borel
More informationRepresentation of the polar cone of convex functions and applications
Representation of the polar cone of convex functions and applications G. Carlier, T. Lachand-Robert October 23, 2006 version 2.1 Abstract Using a result of Y. Brenier [1], we give a representation of the
More informationCONGESTED AGGREGATION VIA NEWTONIAN INTERACTION
CONGESTED AGGREGATION VIA NEWTONIAN INTERACTION KATY CRAIG, INWON KIM, AND YAO YAO Abstract. We consider a congested aggregation model that describes the evolution of a density through the competing effects
More informationSobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations
Sobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations Alessio Figalli Abstract In this note we review some recent results on the Sobolev regularity of solutions
More informationLogarithmic Sobolev Inequalities
Logarithmic Sobolev Inequalities M. Ledoux Institut de Mathématiques de Toulouse, France logarithmic Sobolev inequalities what they are, some history analytic, geometric, optimal transportation proofs
More informationA VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION 1. INTRODUCTION
A VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION GERARD AWANOU AND LEOPOLD MATAMBA MESSI ABSTRACT. We give a proof of existence of a solution to the discrete problem
More informationC 1 regularity of solutions of the Monge-Ampère equation for optimal transport in dimension two
C 1 regularity of solutions of the Monge-Ampère equation for optimal transport in dimension two Alessio Figalli, Grégoire Loeper Abstract We prove C 1 regularity of c-convex weak Alexandrov solutions of
More informationPARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION
PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION ALESSIO FIGALLI AND YOUNG-HEON KIM Abstract. Given Ω, Λ R n two bounded open sets, and f and g two probability densities concentrated
More informationFULL CHARACTERIZATION OF OPTIMAL TRANSPORT PLANS FOR CONCAVE COSTS
FULL CHARACTERIZATION OF OPTIMAL TRANSPORT PLANS FOR CONCAVE COSTS PAUL PEGON, DAVIDE PIAZZOLI, FILIPPO SANTAMBROGIO Abstract. This paper slightly improves a classical result by Gangbo and McCann (1996)
More informationIntroduction to Optimal Transport Theory
Introduction to Optimal Transport Theory Filippo Santambrogio Grenoble, June 15th 2009 These very short lecture notes do not want to be an exhaustive presentation of the topic, but only a short list of
More informationFree Energy, Fokker-Planck Equations, and Random walks on a Graph with Finite Vertices
Free Energy, Fokker-Planck Equations, and Random walks on a Graph with Finite Vertices Haomin Zhou Georgia Institute of Technology Jointly with S.-N. Chow (Georgia Tech) Wen Huang (USTC) Yao Li (NYU) Research
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More informationA DIFFUSIVE MODEL FOR MACROSCOPIC CROWD MOTION WITH DENSITY CONSTRAINTS
A DIFFUSIVE MODEL FOR MACROSCOPIC CROWD MOTION WITH DENSITY CONSTRAINTS Abstract. In the spirit of the macroscopic crowd motion models with hard congestion (i.e. a strong density constraint ρ ) introduced
More informationIntroduction to optimal transport
Introduction to optimal transport Nicola Gigli May 20, 2011 Content Formulation of the transport problem The notions of c-convexity and c-cyclical monotonicity The dual problem Optimal maps: Brenier s
More informationBranched transport limit of the Ginzburg-Landau functional
Branched transport limit of the Ginzburg-Landau functional Michael Goldman CNRS, LJLL, Paris 7 Joint work with S. Conti, F. Otto and S. Serfaty Introduction Superconductivity was first observed by Onnes
More informationA description of transport cost for signed measures
A description of transport cost for signed measures Edoardo Mainini Abstract In this paper we develop the analysis of [AMS] about the extension of the optimal transport framework to the space of real measures.
More informationNumerical methods for Mean Field Games: additional material
Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical
More informationREGULARITY OF MONOTONE TRANSPORT MAPS BETWEEN UNBOUNDED DOMAINS
REGULARITY OF MONOTONE TRANSPORT MAPS BETWEEN UNBOUNDED DOMAINS DARIO CORDERO-ERAUSQUIN AND ALESSIO FIGALLI A Luis A. Caffarelli en su 70 años, con amistad y admiración Abstract. The regularity of monotone
More informationMoment Measures. D. Cordero-Erausquin 1 and B. Klartag 2
Moment Measures D. Cordero-Erausquin 1 and B. Klartag 2 Abstract With any convex function ψ on a finite-dimensional linear space X such that ψ goes to + at infinity, we associate a Borel measure µ on X.
More informationOptimal Transport for Applied Mathematicians
Optimal Transport for Applied Mathematicians Calculus of Variations, PDEs and Modelling Filippo Santambrogio 1 1 Laboratoire de Mathématiques d Orsay, Université Paris Sud, 91405 Orsay cedex, France filippo.santambrogio@math.u-psud.fr,
More informationGradient Flow. Chang Liu. April 24, Tsinghua University. Chang Liu (THU) Gradient Flow April 24, / 91
Gradient Flow Chang Liu Tsinghua University April 24, 2017 Chang Liu (THU) Gradient Flow April 24, 2017 1 / 91 Contents 1 Introduction 2 Gradient flow in the Euclidean space Variants of Gradient Flow in
More informationDisplacement convexity of the relative entropy in the discrete h
Displacement convexity of the relative entropy in the discrete hypercube LAMA Université Paris Est Marne-la-Vallée Phenomena in high dimensions in geometric analysis, random matrices, and computational
More informationTutorial on Optimal Mass Transprot for Computer Vision
Tutorial on Optimal Mass Transprot for Computer Vision David Gu State University of New York at Stony Brook Graduate Summer School: Computer Vision, IPAM, UCLA July 31, 2013 Joint work with Shing-Tung
More informationUNIQUENESS ISSUES FOR EVOLUTIVE EQUATIONS WITH DENSITY CONSTRAINTS
UNIQUENESS ISSUES FOR EVOLUTIVE EQUATIONS WITH DENSITY CONSTRAINTS SIMONE DI MARINO AND ALPÁR RICHÁRD MÉSZÁROS Abstract. In this paper we present some basic uniqueness results for evolutive equations under
More informationarxiv: v1 [math.ap] 10 Apr 2013
QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT DAMON ALEXANDER, INWON KIM, AND YAO YAO arxiv:1304.3072v1 [math.ap] 10 Apr 2013 Abstract. We consider the relationship between Hele-Shaw evolution with
More informationOptimal Transport: A Crash Course
Optimal Transport: A Crash Course Soheil Kolouri and Gustavo K. Rohde HRL Laboratories, University of Virginia Introduction What is Optimal Transport? The optimal transport problem seeks the most efficient
More informationCurvature and the continuity of optimal transportation maps
Curvature and the continuity of optimal transportation maps Young-Heon Kim and Robert J. McCann Department of Mathematics, University of Toronto June 23, 2007 Monge-Kantorovitch Problem Mass Transportation
More informationFast convergent finite difference solvers for the elliptic Monge-Ampère equation
Fast convergent finite difference solvers for the elliptic Monge-Ampère equation Adam Oberman Simon Fraser University BIRS February 17, 211 Joint work [O.] 28. Convergent scheme in two dim. Explicit solver.
More informationDynamic and Stochastic Brenier Transport via Hopf-Lax formulae on Was
Dynamic and Stochastic Brenier Transport via Hopf-Lax formulae on Wasserstein Space With many discussions with Yann Brenier and Wilfrid Gangbo Brenierfest, IHP, January 9-13, 2017 ain points of the
More informationVariational approach to mean field games with density constraints
1 / 18 Variational approach to mean field games with density constraints Alpár Richárd Mészáros LMO, Université Paris-Sud (based on ongoing joint works with F. Santambrogio, P. Cardaliaguet and F. J. Silva)
More informationAPPLICATIONS DU TRANSPORT OPTIMAL DANS LES JEUX À POTENTIEL
APPLICATIONS DU TRANSPORT OPTIMAL DANS LES JEUX À POTENTIEL Adrien Blanchet IAST/TSE/Université Toulouse 1 Capitole Journées SMAI-MODE 216 En collaboration avec G. Carlier & P. Mossay & F. Santambrogio
More informationA variational proof of partial regularity for optimal transportation maps between sets
A variational proof of partial regularity for optimal transportation maps between sets M. Goldman F. Otto October 24, 27 Abstract We provide a new proof of the known partial regularity result for the optimal
More informationarxiv:math/ v3 [math.dg] 30 Jul 2007
OPTIMAL TRANSPORT AND RICCI CURVATURE FOR METRIC-MEASURE SPACES ariv:math/0610154v3 [math.dg] 30 Jul 2007 JOHN LOTT Abstract. We survey work of Lott-Villani and Sturm on lower Ricci curvature bounds for
More informationESTIMATES FOR THE MONGE-AMPERE EQUATION
GLOBAL W 2,p ESTIMATES FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We use a localization property of boundary sections for solutions to the Monge-Ampere equation obtain global W 2,p estimates under
More informationStochastic representation of random positive-definite tensor-valued properties: application to 3D anisotropic permeability random fields
Sous la co-tutelle de : LABORATOIRE DE MODÉLISATION ET SIMULATION MULTI ÉCHELLE CNRS UPEC UNIVERSITÉ PARIS-EST CRÉTEIL UPEM UNIVERSITÉ PARIS-EST MARNE-LA-VALLÉE Stochastic representation of random positive-definite
More informationPATH FUNCTIONALS OVER WASSERSTEIN SPACES. Giuseppe Buttazzo. Dipartimento di Matematica Università di Pisa.
PATH FUNCTIONALS OVER WASSERSTEIN SPACES Giuseppe Buttazzo Dipartimento di Matematica Università di Pisa buttazzo@dm.unipi.it http://cvgmt.sns.it ENS Ker-Lann October 21-23, 2004 Several natural structures
More informationarxiv: v2 [math.na] 7 Mar 2017
CONVERGENCE OF A NEWTON ALGORITHM FOR SEMI-DISCRETE OPTIMAL TRANSPORT JUN KITAGAWA, QUENTIN MÉRIGOT, AND BORIS THIBERT arxiv:1603.05579v2 [math.na] 7 Mar 2017 Abstract. A popular way to solve optimal transport
More informationOn a Class of Multidimensional Optimal Transportation Problems
Journal of Convex Analysis Volume 10 (2003), No. 2, 517 529 On a Class of Multidimensional Optimal Transportation Problems G. Carlier Université Bordeaux 1, MAB, UMR CNRS 5466, France and Université Bordeaux
More informationOptimal Transport for Data Analysis
Optimal Transport for Data Analysis Bernhard Schmitzer 2017-05-16 1 Introduction 1.1 Reminders on Measure Theory Reference: Ambrosio, Fusco, Pallara: Functions of Bounded Variation and Free Discontinuity
More information1. Introduction Boundary estimates for the second derivatives of the solution to the Dirichlet problem for the Monge-Ampere equation
POINTWISE C 2,α ESTIMATES AT THE BOUNDARY FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We prove a localization property of boundary sections for solutions to the Monge-Ampere equation. As a consequence
More informationRicci curvature for metric-measure spaces via optimal transport
Annals of athematics, 169 (2009), 903 991 Ricci curvature for metric-measure spaces via optimal transport By John Lott and Cédric Villani* Abstract We define a notion of a measured length space having
More informationSeismic imaging and optimal transport
Seismic imaging and optimal transport Bjorn Engquist In collaboration with Brittany Froese, Sergey Fomel and Yunan Yang Brenier60, Calculus of Variations and Optimal Transportation, Paris, January 10-13,
More informationarxiv: v3 [stat.co] 22 Apr 2016
A fixed-point approach to barycenters in Wasserstein space Pedro C. Álvarez-Esteban1, E. del Barrio 1, J.A. Cuesta-Albertos 2 and C. Matrán 1 1 Departamento de Estadística e Investigación Operativa and
More informationA LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION O. SAVIN. Introduction In this paper we study the geometry of the sections for solutions to the Monge- Ampere equation det D 2 u = f, u
More information1 An Elementary Introduction to Monotone Transportation
An Elementary Introduction to Monotone Transportation after K. Ball [3] A summary written by Paata Ivanisvili Abstract We outline existence of the Brenier map. As an application we present simple proofs
More informationOptimal Transport for Data Analysis
Optimal Transport for Data Analysis Bernhard Schmitzer 2017-05-30 1 Introduction 1.1 Reminders on Measure Theory Reference: Ambrosio, Fusco, Pallara: Functions of Bounded Variation and Free Discontinuity
More informationUniqueness of the solution to the Vlasov-Poisson system with bounded density
Uniqueness of the solution to the Vlasov-Poisson system with bounded density Grégoire Loeper December 16, 2005 Abstract In this note, we show uniqueness of weak solutions to the Vlasov- Poisson system
More informationOn the approximation of the principal eigenvalue for a class of nonlinear elliptic operators
On the approximation of the principal eigenvalue for a class of nonlinear elliptic operators Fabio Camilli ("Sapienza" Università di Roma) joint work with I.Birindelli ("Sapienza") I.Capuzzo Dolcetta ("Sapienza")
More informationTransport Continuity Property
On Riemannian manifolds satisfying the Transport Continuity Property Université de Nice - Sophia Antipolis (Joint work with A. Figalli and C. Villani) I. Statement of the problem Optimal transport on Riemannian
More informationThe optimal partial transport problem
The optimal partial transport problem Alessio Figalli Abstract Given two densities f and g, we consider the problem of transporting a fraction m [0, min{ f L 1, g L 1}] of the mass of f onto g minimizing
More informationTopology-preserving diffusion equations for divergence-free vector fields
Topology-preserving diffusion equations for divergence-free vector fields Yann BRENIER CNRS, CMLS-Ecole Polytechnique, Palaiseau, France Variational models and methods for evolution, Levico Terme 2012
More informationsoftened probabilistic
Justin Solomon MIT Understand geometry from a softened probabilistic standpoint. Somewhere over here. Exactly here. One of these two places. Query 1 2 Which is closer, 1 or 2? Query 1 2 Which is closer,
More informationOptimal Transportation and Economic Applications
Lecture Notes Optimal Transportation and Economic Applications Guillaume CARLIER IMA, New Mathematical Models in Economics and Finance, June 2010 2 3 Contents 1 Introduction to optimal transportation 6
More informationGenerative Models and Optimal Transport
Generative Models and Optimal Transport Marco Cuturi Joint work / work in progress with G. Peyré, A. Genevay (ENS), F. Bach (INRIA), G. Montavon, K-R Müller (TU Berlin) Statistics 0.1 : Density Fitting
More informationNear-Potential Games: Geometry and Dynamics
Near-Potential Games: Geometry and Dynamics Ozan Candogan, Asuman Ozdaglar and Pablo A. Parrilo September 6, 2011 Abstract Potential games are a special class of games for which many adaptive user dynamics
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationarxiv: v2 [math.oc] 25 Sep 2012
From Knothe s Rearrangement to Brenier s Optimal Transport Map arxiv:1205.1099v2 [math.oc] 25 Sep 2012 Nicolas Bonnotte Université Paris-Sud, Orsay Scuola Normale Superiore, Pisa May 1, 2014 (revised version)
More informationLinear Programming Methods
Chapter 11 Linear Programming Methods 1 In this chapter we consider the linear programming approach to dynamic programming. First, Bellman s equation can be reformulated as a linear program whose solution
More informationGradient Flow in the Wasserstein Metric
3 2 1 Time Position -0.5 0.0 0.5 0 Graient Flow in the Wasserstein Metric Katy Craig University of California, Santa Barbara JMM, AMS Metric Geometry & Topology January 11th, 2018 graient flow an PDE Examples:
More informationJean-David Benamou 1, Guillaume Carlier 2 and Maxime Laborde 3
ESAIM: PROCEEDINGS AND SURVEYS, June 6, Vol. 54, p. -7 B. Düring, C.-B. Schönlieb and M.-T. Wolfram, Editors AN AUGMENTED LAGRANGIAN APPROACH TO WASSERSTEIN GRADIENT FLOWS AND APPLICATIONS Jean-David Benamou,
More informationMean-Field optimization problems and non-anticipative optimal transport. Beatrice Acciaio
Mean-Field optimization problems and non-anticipative optimal transport Beatrice Acciaio London School of Economics based on ongoing projects with J. Backhoff, R. Carmona and P. Wang Thera Stochastics
More informationMartingale optimal transport with Monge s cost function
Martingale optimal transport with Monge s cost function Martin Klimmek, joint work with David Hobson (Warwick) klimmek@maths.ox.ac.uk Crete July, 2013 ca. 1780 c(x, y) = x y No splitting, transport along
More informationMean field games via probability manifold II
Mean field games via probability manifold II Wuchen Li IPAM summer school, 2018. Introduction Consider a nonlinear Schrödinger equation hi t Ψ = h2 1 2 Ψ + ΨV (x) + Ψ R d W (x, y) Ψ(y) 2 dy. The unknown
More informationNumerical Approximation of L 1 Monge-Kantorovich Problems
Rome March 30, 2007 p. 1 Numerical Approximation of L 1 Monge-Kantorovich Problems Leonid Prigozhin Ben-Gurion University, Israel joint work with J.W. Barrett Imperial College, UK Rome March 30, 2007 p.
More informationFrom Euler to Monge and vice versa. Y.B. Lecture 2: From Euler to Vlasov through Monge-Ampère and Kantorovich.
From Euler to Monge and vice versa. Y.B. Lecture 2: From Euler to Vlasov through Monge-Ampère and Kantorovich. Yann Brenier, Mikaela Iacobelli, Filippo Santambrogio, Paris, Durham, Lyon. MFO SEMINAR 1842,
More informationBodies of constant width in arbitrary dimension
Bodies of constant width in arbitrary dimension Thomas Lachand-Robert, Edouard Oudet To cite this version: Thomas Lachand-Robert, Edouard Oudet. Bodies of constant width in arbitrary dimension. Mathematische
More informationMinimizing within convex bodies using a convex hull method
Minimizing within convex bodies using a convex hull method Thomas Lachand-Robert, Edouard Oudet To cite this version: Thomas Lachand-Robert, Edouard Oudet. Minimizing within convex bodies using a convex
More informationNumerical Schemes for Calculating the Discrete Wasserstein Distance
Numerical Schemes for Calculating the Discrete Wasserstein Distance Benjamin Cabrera Born 9th April 99 in Groß-Gerau, Germany February 6, 6 Master s Thesis Mathematics Advisor: Prof. Dr. Martin Rumpf Mathematical
More informationFrom Knothe s transport to Brenier s map and a continuation method for optimal transport
From Knothe s transport to Brenier s map and a continuation method for optimal transport G. Carlier, A. Galichon, F. Santambrogio November 4, 29 Abstract A simple procedure to map two probability measures
More informationarxiv: v1 [math.cv] 30 Jun 2008
EQUIDISTRIBUTION OF FEKETE POINTS ON COMPLEX MANIFOLDS arxiv:0807.0035v1 [math.cv] 30 Jun 2008 ROBERT BERMAN, SÉBASTIEN BOUCKSOM Abstract. We prove the several variable version of a classical equidistribution
More information12. Interior-point methods
12. Interior-point methods Convex Optimization Boyd & Vandenberghe inequality constrained minimization logarithmic barrier function and central path barrier method feasibility and phase I methods complexity
More informationSpaces with Ricci curvature bounded from below
Spaces with Ricci curvature bounded from below Nicola Gigli February 23, 2015 Topics 1) On the definition of spaces with Ricci curvature bounded from below 2) Analytic properties of RCD(K, N) spaces 3)
More informationA Dirichlet problem in the strip
Electronic Journal of Differential Equations, Vol. 1996(1996), No. 10, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp (login: ftp) 147.26.103.110 or 129.120.3.113
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More informationComputing regularization paths for learning multiple kernels
Computing regularization paths for learning multiple kernels Francis Bach Romain Thibaux Michael Jordan Computer Science, UC Berkeley December, 24 Code available at www.cs.berkeley.edu/~fbach Computing
More informationarxiv: v1 [math.ca] 20 Sep 2010
Introduction to Optimal Transport Theory arxiv:1009.3856v1 [math.ca] 20 Sep 2010 Filippo Santambrogio Grenoble, June 15th, 2009 Revised version : March 23rd, 2010 Abstract These notes constitute a sort
More informationc 2017 Michael J. Lindsey and Yanir A. Rubinstein
SIAM J. MATH. ANAL. Vol. 49, No. 4, pp. 3073 3124 c 2017 Michael J. Lindsey and Yanir A. Rubinstein OPTIMAL TRANSPORT VIA A MONGE AMPÈRE OPTIMIZATION PROBLEM MICHAEL LINDSEY AND YANIR A. RUBINSTEIN Abstract.
More informationThe structure of cities
The structure of cities G Université Paris-Dauphine et University of British Columbia Paris, July 6, 2009 (Université Paris-Dauphine The structure et University of cities of British Columbia) Paris, July
More informationGeneralized Ricci Bounds and Convergence of Metric Measure Spaces
Generalized Ricci Bounds and Convergence of Metric Measure Spaces Bornes Généralisées de la Courbure Ricci et Convergence des Espaces Métriques Mesurés Karl-Theodor Sturm a a Institut für Angewandte Mathematik,
More informationA STRATEGY FOR NON-STRICTLY CONVEX TRANSPORT COSTS AND THE EXAMPLE OF x y p IN R 2
A STRATEGY FOR NON-STRICTLY CONVEX TRANSPORT COSTS AND THE EXAMPLE OF x y p IN R 2 GUILLAUME CARLIER, LUIGI DE PASCALE, AND FILIPPO SANTAMBROGIO Abstract. This paper deals with the existence of optimal
More informationA Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth
A Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth E. DiBenedetto 1 U. Gianazza 2 C. Klaus 1 1 Vanderbilt University, USA 2 Università
More informationREGULARITY OF POTENTIAL FUNCTIONS IN OPTIMAL TRANSPORTATION. Centre for Mathematics and Its Applications The Australian National University
ON STRICT CONVEXITY AND C 1 REGULARITY OF POTENTIAL FUNCTIONS IN OPTIMAL TRANSPORTATION Neil Trudinger Xu-Jia Wang Centre for Mathematics and Its Applications The Australian National University Abstract.
More informationConvergence to equilibrium of Markov processes (eventually piecewise deterministic)
Convergence to equilibrium of Markov processes (eventually piecewise deterministic) A. Guillin Université Blaise Pascal and IUF Rennes joint works with D. Bakry, F. Barthe, F. Bolley, P. Cattiaux, R. Douc,
More informationOptimal transport for Seismic Imaging
Optimal transport for Seismic Imaging Bjorn Engquist In collaboration with Brittany Froese and Yunan Yang ICERM Workshop - Recent Advances in Seismic Modeling and Inversion: From Analysis to Applications,
More informationGRADIENT FLOWS FOR NON-SMOOTH INTERACTION POTENTIALS
GRADIENT FLOWS FOR NON-SMOOTH INTERACTION POTENTIALS J. A. CARRILLO, S. LISINI, E. MAININI Abstract. We deal with a nonlocal interaction equation describing the evolution of a particle density under the
More informationCharacterization of Self-Polar Convex Functions
Characterization of Self-Polar Convex Functions Liran Rotem School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel Abstract In a work by Artstein-Avidan and Milman the concept of
More informationSMOOTH OPTIMAL TRANSPORTATION ON HYPERBOLIC SPACE
SMOOTH OPTIMAL TRANSPORTATION ON HYPERBOLIC SPACE JIAYONG LI A thesis completed in partial fulfilment of the requirement of Master of Science in Mathematics at University of Toronto. Copyright 2009 by
More informationApplications of the time derivative of the L 2 -Wasserstein distance and the free entropy dissipation
Applications of the time derivative of the L 2 -Wasserstein distance and the free entropy dissipation Hiroaki YOSHIDA Ochanomizu University Tokyo, Japan at Fields Institute 23 July, 2013 Plan of talk 1.
More informationThe quasi-periodic Frenkel-Kontorova model
The quasi-periodic Frenkel-Kontorova model Philippe Thieullen (Bordeaux) joint work with E. Garibaldi (Campinas) et S. Petite (Amiens) Korean-French Conference in Mathematics Pohang, August 24-28,2015
More information