Splitting Enables Overcoming the Curse of Dimensionality
|
|
- Amberly Bryan
- 5 years ago
- Views:
Transcription
1 Splitting Enables Overcoming the Curse of Dimensionality Jerome Darbon Stanley Osher In this short note we briefly outline a new and remarkably fast algorithm for solving a a large class of high dimensional Hamilton-Jacobi (H-J) initial value problems arising in optimal control and elsewhere [1]. This is done without the use of grids or numerial approximations. Moreover, by using the level set method [8] we can rapidly compute projections of a point in R n, n large to a fairly arbitrary compact set [2]. The method seems to generalize widely beyond what will we present here to some nonconvex Hamiltonians, state dependent Hamiltonians, differential games and perhaps new linear programming algorithms. We begin with the HJ initial value problem ϕ t (x, t) + H( xϕ(x, t)) = in R n (, + ) (1) ϕ(x, ) = J(x) x R n. We assume J : R n R is convex and one coercive, i.e., lim x 2 + J(x) x 2 = +, H : R n R is convex and positively one homogeneous (we sometimes relax all these assumptions). A good example of this is H(v) = v 2. Here v p = (Σ n i=1 v i p ) 1 p for p 1 and x, v = Σ n i=1 x iv i. If we take for J a convex Lipschitz function having the property that, for Ω a convex compact set of R n J(x) < for any x int Ω J(x) = for any x (Ω \ int Ω) J(x) > for any x (R n \ Ω). CNRS/CMLA-Ecole Normale Supérieure de Cachan Mathematics, UCLA, Los Angeles, CA 99-1, Research supported by ONR grants N111683, N and DOE grant DE-SC
2 We call this level set initial data. Then the set of points for which ϕ(x, t) =, t > are exactly those at a distance t, from the boundary of Ω. In fact given x / Ω, then the closest point x opt from x to Ω is exactly To solve (1) we use the Hopf formula [] ϕ( x, t) x opt = x t. (2) ϕ( x, t) 2 ϕ(x, t) = (J + th) (x) = min v R nj (v) + th(v) x, v }, where the Fenchel-Legendre transform f : R n R (+ ) of the convex function f is defined by f (v) = sup x R n v, x f(x)}. Moreover, for free we get that the minimizer satisfies arg min v R nj (v) + th(v) x, v } = x ϕ(x, t). (3) whenever ϕ(, t) is differentiable at x. Let us note here that our algorithm computes ϕ(x, t) but also x ϕ(x, t). Also, we can use the Hopf-Lax formula [, 6] to solve (1). ϕ(x, t) = min z R n J(z) + th ( x z t for convex H. From (4) it is easy to show that if we have k different initial value problems i = 1,... k ϕ i t (x, t) + H( xϕ i (x, t)) =, in R n (, + ) ϕ i (x, ) = J i (x) x R n with the usual hypotheses, then (4) implies, for any x R n, t > ( )} x z ϕ i (x, t) = min J i (z) + th. z R n t )} (4) So min ϕ i(x, t) = i=1,k min z R n min i=1,...,k ( )}} x z J i (z) + th t 2
3 solves the initial value problem ϕ t (x, t) + H( xϕ(x, t)) =, in R n (, + ) ϕ(x, ) = min J i(x) x R n. i=1,...,k This means that if Ω = i=1,...,k Ω i, each Ω i is compact and convex and may overlap, then we can easily compute the set of all points at distance t from Ω which is exactly the solution to () where each J i is a level set function for Ω i. Moreover, at every point x outside of Ω for which there is one i such that ϕ i ( x, t) < ϕ i ( x, t) for any i i, then the closest point x opt to x and Ω is again x opt = x t xϕ i ( x, t) x ϕ i ( x, t). If there are several i for which ϕ i ( x, t) is the minimum among all k of them, then x ϕ will be multivalued, i.e. it will have jumps, but any of the x opt defined above will be a closest point on Ω to x. We solve the optimization problem (3) by using the split Bregman algorithm [4, 3, 9] as follows v k+1 = arg min v R nj (v) x, v + λ 2 dk v b k 2 2}, (6) d k+1 = arg min th(d) + λ } d R n 2 d vk+1 b k 2 2 (7) b k+1 = b k + v k+1 d k+1. (8) Here the sequences (v k ) k N, (d k ) k N both converge to x ϕ(x, t). Let us emphasize again that our numerical algorithm not only computes the solution ϕ(x, t) but also computes x ϕ(x, t) when ϕ(, t) is differentiable. Both (6) and (7), up to change of variables, can be reformulated as finding the unique minimizer of arg min w αf(w) w z 2 2 which is the proximal map of f. Equation (6) can be solved if either J or J have easily computable proximal maps, which often occurs, especially if one of them is smooth. Equation (7) can be easily solved if H(d) = d 2 via the shrink 2 operator defined by z z shrink 2 (z, α) = 2 max( z 2 α, ) if z if z = 3 } ()
4 and we have arg min α w } w 2 w z 2 2 = shrink 2 (z, α) If H(d) = d 1 we use shrink 1 operator defined as follows for any i = 1,..., n z i α if z i > α (shrink 1 (z, α)) i = if z i α z i + α if z i < α and we have arg min α w } w 2 w z 2 2 = shrink 1 (z, α). To solve (7) for more general H(d) convex one homogeneous or to find the proximal map for f of that type we use the fact that H is the characteristic function of a closed convex set C R n H = I c. By using the Moreau identity [7] we realize that the proximal map of H can be obtained by projecting onto C. To do this projection, we merely solve the eikonal equation with level set initial data for C via split Bregman as above in (6), (7), (8) with H(d) = d 2. This is easy using the shrink 2 operator. We then use (2) to obtain the projection and repeat the entire iteration. Numerical experiments on an Intel Laptop Core i-3u running at 2.3 GHz are now presented. We consider diagonal matrices D defined by D ii = i n for i = 1,..., n. We also consider matrices A defined by A ii = 2 for i = 1,..., n and A ij = 1 for i, j = 1,..., n. Table 1 presents the average time (in seconds) to evaluate the solution over 1,, samples (x, t) uniformly drawn in [ 1, 1] n [, 1]. The convergence is remarkably rapid: 1 6 to 1 8 seconds on a standard laptop, per function evaluation. Figure 1 depicts 2-dimensional slices at different times for the (H-J) equation with a weighted l 1 Hamiltonian H = D 1, initial data J = and n = 8. References [1] J. Darbon and S. Osher, Algorithms for Overcoming the Curse of Dimensionality for Certain Hamilton-Jacobi Equations Arising in Control Theory and Elsewhere UCLA CAM report 1-. 4
5 n y 1 y 2 y y D y A e-8 1.2e e-7 7.e e e e e-7 1.7e-6 1.7e e-8 1.6e-7 7.9e-7 1.9e e e-8 1.e e-7 2.4e-6 2.9e-6 Table 1: Time results in seconds for the average time per call for evaluting the solution of the HJ-PDE with the initial data J = , several Hamiltonians and various dimensions n. [2] J. Darbon and S. Osher, Fast projections onto compact sets in high dimensions using the level set method, Hopf formulas and optimization. To appear. [3] R. Glowinski and A. Marrocco, Sur l approximation par éléments finis d ordre, et la resolution par pénalzation-dualité, d une classe de problemémes de Dirichlet non-linéares, C.R. hebd. Séanc. Acad. Sci., Paris 278, série A, (1974), pp [4] T. Goldstein and S. Osher, The Split Bregman Method for L 1 Regularized Problems, SIAM J. on Imaging Science, v.2, (2), 29, pp [] E. Hopf, Generalized Solutions of Nonlinear Equations of the First Order, J. Math. Mech., v.14, (196), pp [6] P. D. Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Conservation Laws, SIAM, 1973 [7] J.-J. Moreau. Proximité et dualité dans un espace hilbertien. Bulletin Spc. Math. France, 93, pp , 196. [8] S. Osher and J.A. Sethian, Fronts Propagating with Curvature Dependent Speech: Algorithms Based on Hamilton-Jacobi Formulations, J. Comput. Phys., v.79, (1), (1988), pp [9] S. Osher and W. Yin, Error Forgetting of Bregman Iteration, J. Sci. Comput., v.4, (2), (213), pp
6 (a) (b) (c) (d) Figure 1: Evaluation of the solution φ((x 1, x 2,,,,,, ), t) of the HJ- PDE with initial data J = and Hamiltonian H = D 1 for (x 1, x 2 ) [ 2, 2] 2 for different times t. Plots for t =, 3,, 8 and respectively depicted in (a), (b), (c) and (d). The level lines multiple of 1 are superimposed on the plots. 6
Algorithm for Overcoming the Curse of Dimensionality for Certain Non-convex Hamilton-Jacobi Equations, Projections and Differential Games
Algorithm for Overcoming the Curse of Dimensionality for Certain Non-convex Hamilton-Jacobi Equations, Projections and Differential Games Yat Tin Chow Jérôme Darbon Stanley Osher Wotao Yin Abstract In
More informationA STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS *
Journal of Computational Mathematics Vol.8, No.4, 010, 55 568. http://www.global-sci.org/jcm doi:10.408/jcm.1003-m0016 A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI
More information/00 $ $.25 per page
Contemporary Mathematics Volume 00, 0000 Domain Decomposition For Linear And Nonlinear Elliptic Problems Via Function Or Space Decomposition UE-CHENG TAI Abstract. In this article, we use a function decomposition
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationBrøndsted-Rockafellar property of subdifferentials of prox-bounded functions. Marc Lassonde Université des Antilles et de la Guyane
Conference ADGO 2013 October 16, 2013 Brøndsted-Rockafellar property of subdifferentials of prox-bounded functions Marc Lassonde Université des Antilles et de la Guyane Playa Blanca, Tongoy, Chile SUBDIFFERENTIAL
More informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
More informationc 1998 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 9, No. 1, pp. 179 189 c 1998 Society for Industrial and Applied Mathematics WEAK SHARP SOLUTIONS OF VARIATIONAL INEQUALITIES PATRICE MARCOTTE AND DAOLI ZHU Abstract. In this work we
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationMerit functions and error bounds for generalized variational inequalities
J. Math. Anal. Appl. 287 2003) 405 414 www.elsevier.com/locate/jmaa Merit functions and error bounds for generalized variational inequalities M.V. Solodov 1 Instituto de Matemática Pura e Aplicada, Estrada
More informationJournal of Convex Analysis (accepted for publication) A HYBRID PROJECTION PROXIMAL POINT ALGORITHM. M. V. Solodov and B. F.
Journal of Convex Analysis (accepted for publication) A HYBRID PROJECTION PROXIMAL POINT ALGORITHM M. V. Solodov and B. F. Svaiter January 27, 1997 (Revised August 24, 1998) ABSTRACT We propose a modification
More informationDuality and dynamics in Hamilton-Jacobi theory for fully convex problems of control
Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control RTyrrell Rockafellar and Peter R Wolenski Abstract This paper describes some recent results in Hamilton- Jacobi theory
More informationSelf-dual Smooth Approximations of Convex Functions via the Proximal Average
Chapter Self-dual Smooth Approximations of Convex Functions via the Proximal Average Heinz H. Bauschke, Sarah M. Moffat, and Xianfu Wang Abstract The proximal average of two convex functions has proven
More informationIterative Convex Optimization Algorithms; Part One: Using the Baillon Haddad Theorem
Iterative Convex Optimization Algorithms; Part One: Using the Baillon Haddad Theorem Charles Byrne (Charles Byrne@uml.edu) http://faculty.uml.edu/cbyrne/cbyrne.html Department of Mathematical Sciences
More informationAuxiliary-Function Methods in Iterative Optimization
Auxiliary-Function Methods in Iterative Optimization Charles L. Byrne April 6, 2015 Abstract Let C X be a nonempty subset of an arbitrary set X and f : X R. The problem is to minimize f over C. In auxiliary-function
More informationGENERAL NONCONVEX SPLIT VARIATIONAL INEQUALITY PROBLEMS. Jong Kyu Kim, Salahuddin, and Won Hee Lim
Korean J. Math. 25 (2017), No. 4, pp. 469 481 https://doi.org/10.11568/kjm.2017.25.4.469 GENERAL NONCONVEX SPLIT VARIATIONAL INEQUALITY PROBLEMS Jong Kyu Kim, Salahuddin, and Won Hee Lim Abstract. In this
More informationA Dykstra-like algorithm for two monotone operators
A Dykstra-like algorithm for two monotone operators Heinz H. Bauschke and Patrick L. Combettes Abstract Dykstra s algorithm employs the projectors onto two closed convex sets in a Hilbert space to construct
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 informationNUMERICAL STUDIES ON HAMILTON-JACOBI EQUATION IN NON-ORTHOGONAL COORDINATES
J. Appl. Math. & Computing Vol. 21(2006), No. 1-2, pp. 603-622 NUMERICAL STUDIES ON HAMILTON-JACOBI EQUATION IN NON-ORTHOGONAL COORDINATES MINYOUNG YUN Abstract. The level-set method is quite efficient
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationarxiv: v1 [math.oc] 8 Apr 2017
Algorithm for Overcoming the Curse of Dimensionality for State-dependent Hamilton-Jacobi equations Yat Tin Chow Jérôme Darbon Stanley Osher Wotao Yin arxiv:74.54v [math.oc] 8 Apr 7 Abstract In this paper,
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationRecent Developments of Alternating Direction Method of Multipliers with Multi-Block Variables
Recent Developments of Alternating Direction Method of Multipliers with Multi-Block Variables Department of Systems Engineering and Engineering Management The Chinese University of Hong Kong 2014 Workshop
More informationConvergence of a first order scheme for a non local eikonal equation
Convergence of a first order scheme for a non local eikonal equation O. Alvarez, E. Carlini, R. Monneau, E. Rouy March 22, 2005 Abstract We prove the convergence of a first order finite difference scheme
More informationA General Framework for a Class of Primal-Dual Algorithms for TV Minimization
A General Framework for a Class of Primal-Dual Algorithms for TV Minimization Ernie Esser UCLA 1 Outline A Model Convex Minimization Problem Main Idea Behind the Primal Dual Hybrid Gradient (PDHG) Method
More informationThe Hopf - Lax solution for state dependent Hamilton - Jacobi equations
The Hopf - Lax solution for state dependent Hamilton - Jacobi equations Italo Capuzzo Dolcetta Dipartimento di Matematica, Università di Roma - La Sapienza 1 Introduction Consider the Hamilton - Jacobi
More informationOn convergence rate of the Douglas-Rachford operator splitting method
On convergence rate of the Douglas-Rachford operator splitting method Bingsheng He and Xiaoming Yuan 2 Abstract. This note provides a simple proof on a O(/k) convergence rate for the Douglas- Rachford
More informationarxiv: v3 [cs.sy] 13 Jun 2018
A Primal-Dual Method for Optimal Control and Trajectory Generation in High-Dimensional Systems* Matthew R Kirchner 12 Gary Hewer 1 Jérôme Darbon 3 and Stanley Osher 4 arxiv:17128226v3 [cssy] 13 Jun 218
More informationCoordinate Update Algorithm Short Course Proximal Operators and Algorithms
Coordinate Update Algorithm Short Course Proximal Operators and Algorithms Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 36 Why proximal? Newton s method: for C 2 -smooth, unconstrained problems allow
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationVISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS, AND ASYMPTOTICS FOR HAMILTONIAN SYSTEMS
VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS, AND ASYMPTOTICS FOR HAMILTONIAN SYSTEMS DIOGO AGUIAR GOMES U.C. Berkeley - CA, US and I.S.T. - Lisbon, Portugal email:dgomes@math.ist.utl.pt Abstract.
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationContraction Methods for Convex Optimization and Monotone Variational Inequalities No.16
XVI - 1 Contraction Methods for Convex Optimization and Monotone Variational Inequalities No.16 A slightly changed ADMM for convex optimization with three separable operators Bingsheng He Department of
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationA characterization of essentially strictly convex functions on reflexive Banach spaces
A characterization of essentially strictly convex functions on reflexive Banach spaces Michel Volle Département de Mathématiques Université d Avignon et des Pays de Vaucluse 74, rue Louis Pasteur 84029
More informationAUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION
New Horizon of Mathematical Sciences RIKEN ithes - Tohoku AIMR - Tokyo IIS Joint Symposium April 28 2016 AUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION T. Takeuchi Aihara Lab. Laboratories for Mathematics,
More informationAn Algorithmic Framework of Generalized Primal-Dual Hybrid Gradient Methods for Saddle Point Problems
An Algorithmic Framework of Generalized Primal-Dual Hybrid Gradient Methods for Saddle Point Problems Bingsheng He Feng Ma 2 Xiaoming Yuan 3 January 30, 206 Abstract. The primal-dual hybrid gradient method
More informationThuong Nguyen. SADCO Internal Review Metting
Asymptotic behavior of singularly perturbed control system: non-periodic setting Thuong Nguyen (Joint work with A. Siconolfi) SADCO Internal Review Metting Rome, Nov 10-12, 2014 Thuong Nguyen (Roma Sapienza)
More informationScaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations
Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.
More informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationSome Inexact Hybrid Proximal Augmented Lagrangian Algorithms
Some Inexact Hybrid Proximal Augmented Lagrangian Algorithms Carlos Humes Jr. a, Benar F. Svaiter b, Paulo J. S. Silva a, a Dept. of Computer Science, University of São Paulo, Brazil Email: {humes,rsilva}@ime.usp.br
More informationIn collaboration with J.-C. Pesquet A. Repetti EC (UPE) IFPEN 16 Dec / 29
A Random block-coordinate primal-dual proximal algorithm with application to 3D mesh denoising Emilie CHOUZENOUX Laboratoire d Informatique Gaspard Monge - CNRS Univ. Paris-Est, France Horizon Maths 2014
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
More informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
More informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationANALYSIS IN SOBOLEV AND BV SPACES
ж ш Л Fi I /V T I Jf"\ Ik 1 Л 1 ANALYSIS IN SOBOLEV AND BV SPACES APPLICATIONS TO PDES AND OPTIMIZATION Hedy Attouch Universite Montpellier I! Montpellier, France Giuseppe Buttazzo Universitä di Pisa Pisa,
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More informationNew Identities for Weak KAM Theory
New Identities for Weak KAM Theory Lawrence C. Evans Department of Mathematics University of California, Berkeley Abstract This paper records for the Hamiltonian H = p + W (x) some old and new identities
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationE. YABLONOVITCH photonic crystals by using level set methods
Appl. Phys. B 81, 235 244 (2005) DOI: 10.1007/s00340-005-1877-3 Applied Physics B Lasers and Optics C.Y. KAO 1, Maximizing band gaps in two-dimensional S. OSHER 2 E. YABLONOVITCH photonic crystals by using
More informationHyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems. Joanna Zwierzchowska
Decision Making in Manufacturing and Services Vol. 9 05 No. pp. 89 00 Hyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems Joanna Zwierzchowska Abstract.
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 informationCONVERGENCE OF APPROXIMATING FIXED POINTS FOR MULTIVALUED NONSELF-MAPPINGS IN BANACH SPACES. Jong Soo Jung. 1. Introduction
Korean J. Math. 16 (2008), No. 2, pp. 215 231 CONVERGENCE OF APPROXIMATING FIXED POINTS FOR MULTIVALUED NONSELF-MAPPINGS IN BANACH SPACES Jong Soo Jung Abstract. Let E be a uniformly convex Banach space
More informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationRobust control and applications in economic theory
Robust control and applications in economic theory In honour of Professor Emeritus Grigoris Kalogeropoulos on the occasion of his retirement A. N. Yannacopoulos Department of Statistics AUEB 24 May 2013
More informationFlux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks
Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks C. Imbert and R. Monneau June 24, 2014 Abstract We study Hamilton-Jacobi equations on networks in the case where Hamiltonians
More informationDevelopments on Variational Inclusions
Advance Physics Letter Developments on Variational Inclusions Poonam Mishra Assistant Professor (Mathematics), ASET, AMITY University, Raipur, Chhattisgarh Abstract - The purpose of this paper is to study
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 4, No. 6, pp. 61 63 c 5 Society for Industrial and Applied Mathematics FAST SWEEPING METHODS FOR STATIC HAMILTON JACOBI EQUATIONS CHIU-YEN KAO, STANLEY OSHER, AND YEN-HSI TSAI
More informationExistence of viscosity solutions for a nonlocal equation modelling polymer
Existence of viscosity solutions for a nonlocal modelling Institut de Recherche Mathématique de Rennes CNRS UMR 6625 INSA de Rennes, France Joint work with P. Cardaliaguet (Univ. Paris-Dauphine) and A.
More informationA splitting minimization method on geodesic spaces
A splitting minimization method on geodesic spaces J.X. Cruz Neto DM, Universidade Federal do Piauí, Teresina, PI 64049-500, BR B.P. Lima DM, Universidade Federal do Piauí, Teresina, PI 64049-500, BR P.A.
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationVISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS
VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS LUIS SILVESTRE These are the notes from the summer course given in the Second Chicago Summer School In Analysis, in June 2015. We introduce the notion of viscosity
More informationA LEVEL SET FRAMEWORK FOR TRACKING MULTI-VALUED SOLUTIONS OF NONLINEAR FIRST-ORDER EQUATIONS
A LEVEL SET FRAMEWORK FOR TRACKING MULTI-VALUED SOLUTIONS OF NONLINEAR FIRST-ORDER EQUATIONS HAILIANG LIU, LI-TIEN CHENG, AND STANLEY OSHER Abstract. We introduce a level set method for the computation
More informationConvergence Theorems for Bregman Strongly Nonexpansive Mappings in Reflexive Banach Spaces
Filomat 28:7 (2014), 1525 1536 DOI 10.2298/FIL1407525Z Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Convergence Theorems for
More informationConvexity of level sets for solutions to nonlinear elliptic problems in convex rings. Paola Cuoghi and Paolo Salani
Convexity of level sets for solutions to nonlinear elliptic problems in convex rings Paola Cuoghi and Paolo Salani Dip.to di Matematica U. Dini - Firenze - Italy 1 Let u be a solution of a Dirichlet problem
More informationMOSCO STABILITY OF PROXIMAL MAPPINGS IN REFLEXIVE BANACH SPACES
MOSCO STABILITY OF PROXIMAL MAPPINGS IN REFLEXIVE BANACH SPACES Dan Butnariu and Elena Resmerita Abstract. In this paper we establish criteria for the stability of the proximal mapping Prox f ϕ =( ϕ+ f)
More informationShape from Shading with discontinuous image brightness
Shape from Shading with discontinuous image brightness Fabio Camilli Dip. di Matematica Pura e Applicata, Università dell Aquila (Italy), camilli@ing.univaq.it Emmanuel Prados UCLA Vision Lab. (USA), eprados@cs.ucla.edu
More informationOn the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (06), 5536 5543 Research Article On the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
More informationCoordinate Update Algorithm Short Course Subgradients and Subgradient Methods
Coordinate Update Algorithm Short Course Subgradients and Subgradient Methods Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 30 Notation f : H R { } is a closed proper convex function domf := {x R n
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationExistence and Approximation of Fixed Points of. Bregman Nonexpansive Operators. Banach Spaces
Existence and Approximation of Fixed Points of in Reflexive Banach Spaces Department of Mathematics The Technion Israel Institute of Technology Haifa 22.07.2010 Joint work with Prof. Simeon Reich General
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 informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationON GAP FUNCTIONS OF VARIATIONAL INEQUALITY IN A BANACH SPACE. Sangho Kum and Gue Myung Lee. 1. Introduction
J. Korean Math. Soc. 38 (2001), No. 3, pp. 683 695 ON GAP FUNCTIONS OF VARIATIONAL INEQUALITY IN A BANACH SPACE Sangho Kum and Gue Myung Lee Abstract. In this paper we are concerned with theoretical properties
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationDeep Relaxation: PDEs for optimizing Deep Neural Networks
Deep Relaxation: PDEs for optimizing Deep Neural Networks IPAM Mean Field Games August 30, 2017 Adam Oberman (McGill) Thanks to the Simons Foundation (grant 395980) and the hospitality of the UCLA math
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 informationON THE ASYMPTOTIC BEHAVIOR OF ELLIPTIC PROBLEMS IN PERIODICALLY PERFORATED DOMAINS WITH MIXED-TYPE BOUNDARY CONDITIONS
Bulletin of the Transilvania University of Braşov Series III: Mathematics, Informatics, Physics, Vol 5(54) 01, Special Issue: Proceedings of the Seventh Congress of Romanian Mathematicians, 73-8, published
More informationOn the Bellman equation for control problems with exit times and unbounded cost functionals 1
On the Bellman equation for control problems with exit times and unbounded cost functionals 1 Michael Malisoff Department of Mathematics, Hill Center-Busch Campus Rutgers University, 11 Frelinghuysen Road
More informationA proximal-like algorithm for a class of nonconvex programming
Pacific Journal of Optimization, vol. 4, pp. 319-333, 2008 A proximal-like algorithm for a class of nonconvex programming Jein-Shan Chen 1 Department of Mathematics National Taiwan Normal University Taipei,
More informationVisco-penalization of the sum of two monotone operators
Visco-penalization of the sum of two monotone operators Patrick L. Combettes a and Sever A. Hirstoaga b a Laboratoire Jacques-Louis Lions, Faculté de Mathématiques, Université Pierre et Marie Curie Paris
More informationAPPROXIMATING SOLUTIONS FOR THE SYSTEM OF REFLEXIVE BANACH SPACE
Bulletin of Mathematical Analysis and Applications ISSN: 1821-1291, URL: http://www.bmathaa.org Volume 2 Issue 3(2010), Pages 32-39. APPROXIMATING SOLUTIONS FOR THE SYSTEM OF φ-strongly ACCRETIVE OPERATOR
More informationSingular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity
Singular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity Sylvie Méléard, Sepideh Mirrahimi September 2, 214 Abstract We perform an asymptotic analysis
More informationThe Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control
The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control A.Younes 1 A. Jarray 2 1 Faculté des Sciences de Tunis, Tunisie. e-mail :younesanis@yahoo.fr 2 Faculté des Sciences
More informationLevel-set convex Hamilton-Jacobi equations on networks
Level-set convex Hamilton-Jacobi equations on networks C. Imbert and R. Monneau January 17, 2014 Abstract The paper deals with Hamilton-Jacobi equations on networks with level-set convex (in the gradient
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT -3 Stability analysis of the numerical density-enthalpy model Ibrahim, F J Vermolen, C Vui ISSN 389-65 Reports of the Department of Applied Mathematical Analysis Delft
More informationThe Fast Sweeping Method. Hongkai Zhao. Department of Mathematics University of California, Irvine
The Fast Sweeping Method Hongkai Zhao Department of Mathematics University of California, Irvine Highlights of fast sweeping mehtod Simple (A Gauss-Seidel type of iterative method). Optimal complexity.
More informationExistence and convergence theorems for the split quasi variational inequality problems on proximally smooth sets
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (206), 2364 2375 Research Article Existence and convergence theorems for the split quasi variational inequality problems on proximally smooth
More informationOn the Convergence and O(1/N) Complexity of a Class of Nonlinear Proximal Point Algorithms for Monotonic Variational Inequalities
STATISTICS,OPTIMIZATION AND INFORMATION COMPUTING Stat., Optim. Inf. Comput., Vol. 2, June 204, pp 05 3. Published online in International Academic Press (www.iapress.org) On the Convergence and O(/N)
More informationContraction Methods for Convex Optimization and monotone variational inequalities No.12
XII - 1 Contraction Methods for Convex Optimization and monotone variational inequalities No.12 Linearized alternating direction methods of multipliers for separable convex programming Bingsheng He Department
More informationDual and primal-dual methods
ELE 538B: Large-Scale Optimization for Data Science Dual and primal-dual methods Yuxin Chen Princeton University, Spring 2018 Outline Dual proximal gradient method Primal-dual proximal gradient method
More informationPrimal Dual Methods in Imaging Science
Primal Dual Methods in Imaging Science Ernie Esser UC Irvine Joint work with Xiaoqun Zhang (Shanghai Jiaotong University) and Tony Chan (HKUST) ICIAM 2011 1 Outline Motivation to Study Primal Dual Methods
More informationARock: an algorithmic framework for asynchronous parallel coordinate updates
ARock: an algorithmic framework for asynchronous parallel coordinate updates Zhimin Peng, Yangyang Xu, Ming Yan, Wotao Yin ( UCLA Math, U.Waterloo DCO) UCLA CAM Report 15-37 ShanghaiTech SSDS 15 June 25,
More informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationA double projection method for solving variational inequalities without monotonicity
A double projection method for solving variational inequalities without monotonicity Minglu Ye Yiran He Accepted by Computational Optimization and Applications, DOI: 10.1007/s10589-014-9659-7,Apr 05, 2014
More informationScaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations
Journal of Statistical Physics, Vol. 122, No. 2, January 2006 ( C 2006 ) DOI: 10.1007/s10955-005-8006-x Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations Jan
More informationAccelerated primal-dual methods for linearly constrained convex problems
Accelerated primal-dual methods for linearly constrained convex problems Yangyang Xu SIAM Conference on Optimization May 24, 2017 1 / 23 Accelerated proximal gradient For convex composite problem: minimize
More informationExtremal Solutions of Differential Inclusions via Baire Category: a Dual Approach
Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Alberto Bressan Department of Mathematics, Penn State University University Park, Pa 1682, USA e-mail: bressan@mathpsuedu
More information