Numerical Approximation of L 1 Monge-Kantorovich Problems
|
|
- Marcus Osborne
- 5 years ago
- Views:
Transcription
1 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
2 Rome March 30, 2007 p. 2 Monge, 1781 Given two distributions of mass with the same finite total, embankment f + and excavation f, find a transport map T which carries f + into f, f + = f B R n Borel, T 1 (B) and minimizes the total transportation cost c(x,t(x))f + (x) B where c(x,y) is the cost of transporting one unit of mass from x to y and f ± M + (R n ), the set of nonnegative measures.
3 Rome March 30, 2007 p. 3 Kantorovich, 1942 Find an optimal transport plan, a measure γ M + (R n R n ) satisfying γ = f +, γ = B R n B R n B and minimizing the transportation cost R n R n c(x,y)γ(x,y) B f The difference is that splitting of mass is allowed: the mass from a point x is distributed as γ(x, ). This is a linear programming problem. The infimum is known as Kantorovich-Wasserstein distance.
4 Rome March 30, 2007 p. 4 Recent much progress is due to Ambrosio, Brenier, Caffarelli, Evans, Feldman, Gangbo, McCann,... and was obtained for very general cost functions. The most studied: c(x,y) = x y p, with p = 1, 2. p = 2: The optimal plan is unique, it is an optimal map, and this map is a gradient of a convex function. Numerical algorithms: Benamou and Brenier 00, Angenent, Haker, Tannenbaum 03. p = 1: The cost is not strictly convex, not every optimal plan is a map, an optimal map exists but is not unique. Numerical algorithms: the issue of this talk.
5 Rome March 30, 2007 p. 5 transport density and potential Dual formulation of L 1 MK problem, Kantorovich 42, { } max uf : u Lip 1 (R n ), R n where f = f + f. This problem is equivalent to the following PDE (MK equation, Evans 97:) (a u) = f u 1, a 0, u < 1 a = 0 Here u is Kantorovich potential, a - transport density. A weak solution {u,a} Lip 1 (R n ) M + (R n ) exists for any f with finite total variation and zero average, Bouchitté, Buttazzo, and Seppecher 97.
6 Rome March 30, 2007 p. 6 transport flux Although optimal transport plan is not unique, the density a is the same for all such plans (under mild regularity conditions on f, Feldman and McCann 02). The transport flux q = a u is also the same for all optimal plans and contains all information on the local transport. The Kantorovich-Wasserstein distance between f + and f (the minimal cost of transportation) can be found as C = a = q Our goal: computing the optimal transport flux.
7 Rome March 30, 2007 p. 7 dual formulation for flux Let q = a u, where {a,u} solve the MK equation, (a u) = f u 1, a 0, u < 1 a = 0. For any test vector function v, u (v q) u v + u q v q. Hence we arrive at q = f and v q u (v q) v, a mixed formulation of the equivalent dual to MK problem (Bouchitté, Buttazzo, and Seppecher 97): min { q : q [M(R n )] n and q = f }
8 Rome March 30, 2007 p. 8 generalization the transport should be confined to the closure of a Lipschitz domain Ω R n with spt(f) Ω. the cost of transporting one unit of mass to the distance dx near a point x is k(x)dx, where k : Ω R >0 is a given function. In this case the equivalent flux formulation of MK problem, Pratelli 05, reads: { } min k(x) q q = f in Ω, q n = 0 on Ω Ω where the constraint is understood in the sense of distributions, (q, φ) = (f,φ), φ C 1.
9 Rome March 30, 2007 p. 9 mixed method Let V M 0 = v [M(Ω)]n v M(Ω), v,φ = v, φ φ C 1 (Ω) Problem (Q): Find q V M 0 and u C(Ω) such that q,φ = f,φ φ C(Ω), v,k q,k (v q),u v V M 0
10 Rome March 30, 2007 p. 10 regularization Replace by 1 r r with r > 1 and consider the regularized Problem (Q r ): Find q r V0 r and u r L p I (Ω) such that ( q r,φ) = (f r,φ) φ L p (Ω), (k q r r 2 q r,v) = ( q r,u r ) v V0 r Here V0 r = {v [Lr (Ω)] n : v L r (Ω), v n Ω = 0} L p I (Ω) = {ψ Lp (Ω) : Ω ψ = 0}, is a regularization of f. 1 r + 1 p = 1, and f r L r I
11 Rome March 30, 2007 p. 11 Theorem 1. The problem (Q r ), r > 1, has a unique solution. As r 1 there exists a subsequence of {q r,u r } such that q r q -weakly in [M(Ω)] n ; q r q -weakly in M(Ω); u r u strongly in C(Ω). Moreover, {q,u} V M 0 C I (Ω) and solves (Q).
12 Rome March 30, 2007 p. 12 approximation of (Q r ) q r approximated using divergence-conforming Raviart-Thomas finite elements of the lowest order, RT h 0, and vertex sampling on nonlinear term; u r by piecewise constant functions, their set is denoted S h, and S h I = {φ Sh : Ω φ = 0}. Here h is the maximal element diameter. Problem (Q h r): Find q h r V h 0 and uh r Sh I such that ( q h r,φh ) = (f r,φ h ) φ h S h, (k q h r r 2 q h r,vh ) h = ( q h r,uh r ) vh RT h 0.
13 Theorem 2. The unique solution {q h r,u h r } of (Q h r) converges, for r > 1 fixed and h 0, to {q r,u r }, the unique solution of (Q r ). Implementation: iterations, ( Q j+1 2,φ h ) =(f r,φ h ) φ h S h, (k{ Q j r 2 Q j + Q j r 2 ε (Q j+1 2 Q j )},v h ) h = (U j+1, v h ) v h RT h 0 where the regularization Q ε = ( Q 2 + ε 2 ) 1 2 makes all terms well defined, successive over-relaxation Q j+1 = γq j (1 γ)q j, and adaptive mesh generation. Rome March 30, 2007 p. 13
14 Rome March 30, 2007 p. 14 inhomogeneous domain f + k 1 f k = 1 Ω is a square, f + and f are uniform in their supports, left and right ellipses, respectively. The cost parameter: k = 1 outside the inclined ellipse, inside the ellipse k = k 1.
15 Rome March 30, 2007 p. 15 k 1 = 0.01 and k 1 = 1 mesh flux q q -levels
16 Rome March 30, 2007 p. 16 k 1 = 3 and k 1 = mesh flux q q -levels
17 f+ = Pl 1 i=1 l δ(x yi ), f = Pn 1 i=1 n δ(x xi ); yi, xi, Upper row: exact optimal transport rays (linear programming); lower row: approximate solution (mixed scheme). Transport cost errors: 0.05% left (l = 9, n = 6); 0.3% right (l = 11, n = 7). Rome March 30, 2007 p. 17
18 Rome March 30, 2007 p. 18 opt. couplings: f + 1/ Ω, f = n i=1 α iδ(x x i ), αi = 1 Solution characterization, Rüschendorf and Uckelmann 00: Let F(x) = max i { x x i + a i } and A i = {x Ω : F(x) = x x i + a i }. The coupling {A i } is optimal iff a i are such that A i = α i, i. Left: all a i are zero, points are coupled to the closest sink, exact boundaries can be found. Right: all sinks are equal, the boundaries are found as small-level contours of q.
19 Rome March 30, 2007 p. 19 partial MK problem Let only a given amount, m min( f +, f ), should be optimally transported from f + to f ; the distributions can be unbalanced. It is needed to find both the optimal transportation domains and the optimal transport plan for them. Free boundary problem in optimal transportation: Caffarelli and McCann 07, cost = (distance) 2. We now derive a variational formulation for the partial L 1 MK problems, cost = distance, and solve them numerically.
20 Rome March 30, 2007 p. 20 partial L 1 MK problem We assume again transport is allowed only in Ω, the local cost of transporting one unit of mass is k(x)dx. Let us define Γ m M+ (Ω Ω) as the set of nonnegative measures satisfying Ω Ω γ = m and dominated by f ±, i.e. for every Borel set B Ω B Ω γ B f+, Ω B γ B f. The optimal transport plan minimizes the total cost of transportation in Γ m. The plan can be presented as a part of an optimal plan for a balanced problem.
21 Rome March 30, 2007 p. 21 auxiliary balanced problem imbed R n into R n+1 as a subspace x n+1 = 0; introduce two auxiliary sources, w and w +, supported in hyperplanes x n+1 = l and x n+1 = l, respectively; postulate free transportation, k = 0, in auxiliary planes; everywhere else we assume k = 1. The MK problem obtained is balanced. It contains Dirichlet regions, Bouchitté and Buttazzo 01. Flux formulation holds true, Pratelli 05.
22 Rome March 30, 2007 p. 22 auxiliary balanced problem Transport in auxiliary planes is free; If l > l 0 = 1 2 max{d Ω(x,y) : x spt(f + ), y spt(f )} transport from f + to f via an auxiliary plane is not paying; The cost of interplane transport is the same for all reasonable plans. If q is an optimal flux, Q is a solution to the partial { problem. Flux formulation } min k q : q = f Ω [ l,l] can be simplified.
23 Rome March 30, 2007 p. 23 new variational formulation Problem (PMK): min K P Ω { Q + l( q + + q ) } where K P V M 0 [M(Ω)] 2 is the set of triples {Q,q +,q } satisfying three mass balance conditions: Q + q + q = f + f, q + = f + m, Ω Ω q = f m. Ω Ω
24 Rome March 30, 2007 p. 24 special case Let m = f + < f Then w = 0 and q + = 0, so q = Q f + + f and q = f m, Q V M 0 Problem (PMK 0 ): { Q +l Q f + +f } min Q V M 0 Same if m = f < f +.
25 Rome March 30, 2007 p. 25 back to balance In the balanced case, m = f + = f, the fluxes q +, q are zero. Hence Q = f + f and we return to the known flux formulation, { min Q V M 0 Q } Q = f+ f Alternatively, for any l > l 0 we have also { Q + l Q f + + f } min Q V M 0 as an equivalent flux formulation of L 1 MK problem.
26 Rome March 30, 2007 p. 26 existence for partial MK Barrett and P., in preparation: Regularization: 1 r r, where 0 < r 1 1; For any r > 1 the Euler-Lagrange equations for regularized problem have a unique solution: the fluxes Q r V r 0, q ± r L r (Ω ± ) and Lagrange multipliers u r L p I (Ω), λ± r R. Here spt(f ± ) Ω ± Ω, p = r r 1. As r 1 the solutions converge (for a subsequence) to a solution of Euler-Lagrange equations for the flux formulation of the partial L 1 MK problem.
27 Rome March 30, 2007 p. 27 solution of partial problems: regularization: 1 r r where 0 < r 1 1; constraint optimization: augmented Lagrangian method; approximation: Raviart-Thomas elements of the least order and vertex sampling on nonlinear terms for Q r, piecewise constant elements for q r ± and the Lagrange multiplier u r of the pointwise divergence constraint; nonlinearity: iterations, successive overrelaxation; singularities: adaptive mesh generation.
28 Rome March 30, 2007 p. 28 example Left: f + = 1 in the left rectangle, f = 1 in the right one; m = 1 3 f + = 1 2 An optimal plan is easy to find, transportation cost error is less than 0.05%; red line " " is the free boundary. Right: elliptic obstacle is placed between the rectangles. f.
29 Rome March 30, 2007 p. 29 example: partial coupling f + = 1 in the whole unit square Ω f is a sum of five equal point sinks of total mass m < 1. " " is the free boundary. Left: m = 0.5; right: m = 0.75;
30 Rome March 30, 2007 p. 30 partial coupling: mesh Results of three mesh adaptions, elements.
Optimal 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 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 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 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 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 informationExample 1. Hamilton-Jacobi equation. In particular, the eikonal equation. for some n( x) > 0 in Ω. Here 1 / 2
Oct. 1 0 Viscosity S olutions In this lecture we take a glimpse of the viscosity solution theory for linear and nonlinear PDEs. From our experience we know that even for linear equations, the existence
More informationContinuum Modeling of Transportation Networks with Differential Equations
with Differential Equations King Abdullah University of Science and Technology Thuwal, KSA Examples of transportation networks The Silk Road Examples of transportation networks Painting by Latifa Echakhch
More informationHomework Set #6 - Solutions
EE 15 - Applications of Convex Optimization in Signal Processing and Communications Dr Andre Tkacenko JPL Third Term 11-1 Homework Set #6 - Solutions 1 a The feasible set is the interval [ 4] The unique
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 informationTHREE SINGULAR VARIATIONAL PROBLEMS. Lawrence C. Evans Department of Mathematics University of California Berkeley, CA 94720
THREE SINGLAR VARIATIONAL PROBLEMS By Lawrence C. Evans Department of Mathematics niversity of California Berkeley, CA 9470 Some of the means I use are trivial and some are quadrivial. J. Joyce Abstract.
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 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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
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 informationExam in TMA4180 Optimization Theory
Norwegian University of Science and Technology Department of Mathematical Sciences Page 1 of 11 Contact during exam: Anne Kværnø: 966384 Exam in TMA418 Optimization Theory Wednesday May 9, 13 Tid: 9. 13.
More informationFree boundaries in optimal transport and Monge-Ampère obstacle problems
Annals of Mathematics, 0 (2006), 1 57 Free boundaries in optimal transport and Monge-Ampère obstacle problems By Luis A. Caffarelli and Robert J. McCann* E-mail address: caffarel@math.utexas.edu and mccann@math.toronto.edu
More informationBoundary value problems for the infinity Laplacian. regularity and geometric results
: regularity and geometric results based on joint works with Graziano Crasta, Roma La Sapienza Calculus of Variations and Its Applications - Lisboa, December 2015 on the occasion of Luísa Mascarenhas 65th
More informationThin Film Problems in Superconductivity
Thin Film Problems in Superconductivity L. Prigozhin 1 J. W. Barrett 2 V. Sokolovsky 1 1 Ben-Gurion University of the Negev, Israel 2 Imperial College London, UK Sept 2016 Outlook: Thin film magnetization
More informationMeasure-valued - strong uniqueness for hyperbolic systems
Measure-valued - strong uniqueness for hyperbolic systems joint work with Tomasz Debiec, Eduard Feireisl, Ondřej Kreml, Agnieszka Świerczewska-Gwiazda and Emil Wiedemann Institute of Mathematics Polish
More informationConvex Optimization and Modeling
Convex Optimization and Modeling Duality Theory and Optimality Conditions 5th lecture, 12.05.2010 Jun.-Prof. Matthias Hein Program of today/next lecture Lagrangian and duality: the Lagrangian the dual
More informationSur le tas de Sable. Noureddine Igbida (LAMFA CNRS-UMR 6140, Université de Picardie Jules Sur le Verne, tas de SableAmiens)
Sur le tas de Sable Noureddine Igbida LAMFA CNRS-UMR 6140, Université de Picardie Jules Verne, 80000 Amiens Nancy, 4 Mars 2008 Nancy, 4 Mars 2008 1 / 26 Une source de sable qui tourne : résultat d une
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 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 informationOptimal transportation on non-compact manifolds
Optimal transportation on non-compact manifolds Albert Fathi, Alessio Figalli 07 November 2007 Abstract In this work, we show how to obtain for non-compact manifolds the results that have already been
More informationGeodesic Density Tracking with Applications to Data Driven Modeling
Geodesic Density Tracking with Applications to Data Driven Modeling Abhishek Halder and Raktim Bhattacharya Department of Aerospace Engineering Texas A&M University College Station, TX 77843-3141 American
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 informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationLagrange duality. The Lagrangian. We consider an optimization program of the form
Lagrange duality Another way to arrive at the KKT conditions, and one which gives us some insight on solving constrained optimization problems, is through the Lagrange dual. The dual is a maximization
More informationLagrangian Duality Theory
Lagrangian Duality Theory Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye Chapter 14.1-4 1 Recall Primal and Dual
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 informationRecent developments in elliptic partial differential equations of Monge Ampère type
Recent developments in elliptic partial differential equations of Monge Ampère type Neil S. Trudinger Abstract. In conjunction with applications to optimal transportation and conformal geometry, there
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 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 informationNeedle decompositions and Ricci curvature
Tel Aviv University CMC conference: Analysis, Geometry, and Optimal Transport KIAS, Seoul, June 2016. A trailer (like in the movies) In this lecture we will not discuss the following: Let K 1, K 2 R n
More informationAnalyse d un Modèle de tas de Sable
Analyse d un Modèle de tas de Sable Noureddine Igbida LAMFA CNRS-UMR 6140, Université de Picardie Jules Verne, 80000 Amiens Premier congrès de la SM2A, Rabat 5-8 Février 2008 de Sable Premier congrés de
More informationRegularity and compactness for the DiPerna Lions flow
Regularity and compactness for the DiPerna Lions flow Gianluca Crippa 1 and Camillo De Lellis 2 1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. g.crippa@sns.it 2 Institut für Mathematik,
More informationOptimal Martingale Transport problem in higher dimensions
University of British Columbia Based on joint work with Young-Heon Kim & Tongseok Lim (UBC) September, 2015 Varese Main Points of the talk: Comparison between optimal transport problem and optimal martingale
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 informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
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 informationMath 742: Geometric Analysis
Math 742: Geometric Analysis Lecture 5 and 6 Notes Jacky Chong jwchong@math.umd.edu The following notes are based upon Professor Yanir ubenstein s lectures with reference to Variational Methods 4 th edition
More informationarxiv: v1 [math.oc] 19 Dec 2013
THE MONGE-KANTOROVICH PROBLEM FOR DISTRIBUTIONS AND APPLICATIONS GUY BOUCHITTÉ, GIUSEPPE BUTTAZZO, LUIGI DE PASCALE arxiv:1312.5453v1 [math.oc] 19 Dec 2013 Abstract. We study the Kantorovich-Rubinstein
More informationOutline. Relaxation. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Lagrangian Relaxation. Lecture 12 Single Machine Models, Column Generation
Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING 1. Lagrangian Relaxation Lecture 12 Single Machine Models, Column Generation 2. Dantzig-Wolfe Decomposition Dantzig-Wolfe Decomposition Delayed Column
More informationOptimal transportation and action-minimizing measures
Scuola Normale Superiore of Pisa and École Normale Supérieure of Lyon Phd thesis 24th October 27 Optimal transportation and action-minimizing measures Alessio Figalli a.figalli@sns.it Advisor Prof. Luigi
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationOn Lagrangian solutions for the semi-geostrophic system with singular initial data
On Lagrangian solutions for the semi-geostrophic system with singular initial data Mikhail Feldman Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA feldman@math.wisc.edu
More informationOptimal transportation for a quadratic cost with convex constraints and applications
Optimal transportation for a quadratic cost with convex constraints and applications C. Jimenez, F. Santambrogio October 13, 2011 Abstract We prove existence of an optimal transport map in the Monge- Kantorovich
More informationautomating the analysis and design of large-scale optimization algorithms
automating the analysis and design of large-scale optimization algorithms Laurent Lessard University of Wisconsin Madison Joint work with Ben Recht, Andy Packard, Bin Hu, Bryan Van Scoy, Saman Cyrus October,
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationA primal-dual approach for a total variation Wasserstein flow
A primal-dual approach for a total variation Wasserstein flow Martin Benning 1, Luca Calatroni 2, Bertram Düring 3, Carola-Bibiane Schönlieb 4 1 Magnetic Resonance Research Centre, University of Cambridge,
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 informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationConvex Optimization. Ofer Meshi. Lecture 6: Lower Bounds Constrained Optimization
Convex Optimization Ofer Meshi Lecture 6: Lower Bounds Constrained Optimization Lower Bounds Some upper bounds: #iter μ 2 M #iter 2 M #iter L L μ 2 Oracle/ops GD κ log 1/ε M x # ε L # x # L # ε # με f
More informationOn the Optimal Insulation of Conductors 1
JOURNAL OF OPTIMIZATION THEORY AND APPLICATIONS: Vol. 100, No. 2, pp. 253-263, FEBRUARY 1999 On the Optimal Insulation of Conductors 1 S. J. COX, 2 B. KAWOHL, 3 AND P. X. UHLIG 4 Communicated by K. A.
More informationOutline. Roadmap for the NPP segment: 1 Preliminaries: role of convexity. 2 Existence of a solution
Outline Roadmap for the NPP segment: 1 Preliminaries: role of convexity 2 Existence of a solution 3 Necessary conditions for a solution: inequality constraints 4 The constraint qualification 5 The Lagrangian
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
More informationNon Linear Elliptic Theory and the Monge-Ampere Equation
ICM 2002 Vol. I 179^187 Non Linear Elliptic Theory and the Monge-Ampere Equation Luis A. Caffarelli* Abstract The Monge-Ampere equation, plays a central role in the theory of fully non linear equations.
More informationLecture 3: Lagrangian duality and algorithms for the Lagrangian dual problem
Lecture 3: Lagrangian duality and algorithms for the Lagrangian dual problem Michael Patriksson 0-0 The Relaxation Theorem 1 Problem: find f := infimum f(x), x subject to x S, (1a) (1b) where f : R n R
More informationMultiphysics optimal transportation and image analysis
Multiphysics optimal transportation and image analysis Afaf Bouharguane, Emmanuel Maitre, Edouard Oudet, Nicolas Papadakis To cite this version: Afaf Bouharguane, Emmanuel Maitre, Edouard Oudet, Nicolas
More informationAlgorithms for Nonsmooth Optimization
Algorithms for Nonsmooth Optimization Frank E. Curtis, Lehigh University presented at Center for Optimization and Statistical Learning, Northwestern University 2 March 2018 Algorithms for Nonsmooth Optimization
More informationNumerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation
Numerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation Sören Bartels University of Freiburg, Germany Joint work with G. Buttazzo (U Pisa), L. Diening (U Bielefeld),
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 informationCONSTRUCTING OPTIMAL MAPS FOR MONGE'S TRANSPORT PROBLEM AS A LIMIT OF STRICTLY CONVEX COSTS
JOURNAL OF THE AMERICAN MATHEMATICAL SOCIETY Volume 15, Number 1, Pages 1{6 S 0894-0347(01)00376-9 Article electronically published on July 31, 001 CONSTRUCTING OPTIMAL MAPS FOR MONGE'S TRANSPORT PROBLEM
More informationThe Fundamental Theorem of Linear Inequalities
The Fundamental Theorem of Linear Inequalities Lecture 8, Continuous Optimisation Oxford University Computing Laboratory, HT 2006 Notes by Dr Raphael Hauser (hauser@comlab.ox.ac.uk) Constrained Optimisation
More informationPartial regularity for fully nonlinear PDE
Partial regularity for fully nonlinear PDE Luis Silvestre University of Chicago Joint work with Scott Armstrong and Charles Smart Outline Introduction Intro Review of fully nonlinear elliptic PDE Our result
More informationBeyond Heuristics: Applying Alternating Direction Method of Multipliers in Nonconvex Territory
Beyond Heuristics: Applying Alternating Direction Method of Multipliers in Nonconvex Territory Xin Liu(4Ð) State Key Laboratory of Scientific and Engineering Computing Institute of Computational Mathematics
More informationUpon successful completion of MATH 220, the student will be able to:
MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient
More informationAnisotropic congested transport
Anisotropic congested transport Lorenzo Brasco LATP Aix-Marseille Université lbrasco@gmail.com http://www.latp.univ-mrs.fr/~brasco/ Sankt Peterburg, 04/06/2012 References Part of the results here presented
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
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 informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
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 informationA primal-dual approach for a total variation Wasserstein flow
A primal-dual approach for a total variation Wasserstein flow Martin Benning 1, Luca Calatroni 2, Bertram Düring 3, Carola-Bibiane Schönlieb 4 1 Magnetic Resonance Research Centre, University of Cambridge,
More informationSymmetry breaking for a problem in optimal insulation
Symmetry breaking for a problem in optimal insulation Giuseppe Buttazzo Dipartimento di Matematica Università di Pisa buttazzo@dm.unipi.it http://cvgmt.sns.it Geometric and Analytic Inequalities Banff,
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 informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationJuan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationMAT 578 FUNCTIONAL ANALYSIS EXERCISES
MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.
More informationConvex 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 GeMeCod Conference October
More informationAlgorithms for constrained local optimization
Algorithms for constrained local optimization Fabio Schoen 2008 http://gol.dsi.unifi.it/users/schoen Algorithms for constrained local optimization p. Feasible direction methods Algorithms for constrained
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
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 informationCalculus of Variations
Calc. Var. 20, 329 341 (2004) DOI: 10.1007/s00526-003-0239-4 Calculus of Variations Xu-Jia Wang On the design of a reflector antenna II Received: 7 July 2003 / Accepted: 26 August 2003 Published online:
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
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 informationCOMP 652: Machine Learning. Lecture 12. COMP Lecture 12 1 / 37
COMP 652: Machine Learning Lecture 12 COMP 652 Lecture 12 1 / 37 Today Perceptrons Definition Perceptron learning rule Convergence (Linear) support vector machines Margin & max margin classifier Formulation
More informationConvex Optimization & Lagrange Duality
Convex Optimization & Lagrange Duality Chee Wei Tan CS 8292 : Advanced Topics in Convex Optimization and its Applications Fall 2010 Outline Convex optimization Optimality condition Lagrange duality KKT
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationINTEGRAL ESTIMATES FOR TRANSPORT DENSITIES
INTEGRAL ESTIMATES FOR TRANSPORT DENSITIES L. DE PASCALE, L. C. EVANS, AND A. PRATELLI Abstract. We introduce some integration-by-parts methods that improve upon the L p estimates on transport densitites
More informationCalculus of Variations Summer Term 2014
Calculus of Variations Summer Term 2014 Lecture 12 26. Juni 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture 12 26. Juni 2014 1 / 25 Purpose of Lesson Purpose of Lesson: To discuss numerical
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 informationNumerische Mathematik
Numer. Math. (2) 84: 375 393 Digital Object Identifier (DOI) 1.17/s21199117 Numerische Mathematik c Springer-Verlag 2 A computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem
More informationMEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS
MEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS KARL KUNISCH, KONSTANTIN PIEPER, AND BORIS VEXLER Abstract. A directional sparsity framework allowing for measure valued controls
More informationXingye Yue. Soochow University, Suzhou, China.
Relations between the multiscale methods p. 1/30 Relationships beteween multiscale methods for elliptic homogenization problems Xingye Yue Soochow University, Suzhou, China xyyue@suda.edu.cn Relations
More information