Proximal point algorithm in Hadamard spaces
|
|
- Barnaby Flynn
- 5 years ago
- Views:
Transcription
1 Proximal point algorithm in Hadamard spaces Miroslav Bacak Télécom ParisTech Optimisation Géométrique sur les Variétés - Paris, 21 novembre 2014
2 Contents of the talk 1 Basic facts on Hadamard spaces 2 Proximal point algorithm 3 Applications to computational phylogenetics
3 Proximal point algorithm in Hadamard spaces Why? Well... it is used in: Phylogenetics: computing medians and means of phylogenetic trees. diffusion tensor imaging: the space P (n, R) of symmetric positive definite matrices n n with real entries is a Hadamard space if it is equipped with the Riemannian metric X, Y A := Tr ( A 1 XA 1 Y ), X, Y T A (P (n, R)), for every A P (n, R). Computational biology: shape analyses of tree-like structures:
4 Tree-like structures in organisms Figure: Bronchial tubes in lungs Figure: Transport system in plants Figure: Human circulatory system
5 Definition of Hadamard space Let (H, d) be a complete metric space where: 1 any two points x 0 and x 1 are connected by a geodesic x: [0, 1] H: t x t, 2 and, d (y, x t ) 2 (1 t)d (y, x 0 ) 2 + td (y, x 1 ) 2 t(1 t)d (x 0, x 1 ) 2, for every y H. Then (H, d) is called a Hadamard space. For today: assume that local compactness.
6 Geodesic space x 2 x 1
7 Geodesic space x 2 x 1
8 Geodesic space (1 λ)x 1 + λx 2 x 2 x 1
9 Definition of nonpositive curvature A geodesic triangle in a geodesic space: x 3 x 2 x 1
10 Terminology remark CAT(κ) spaces, for κ R, were introduced in 1987 by Michail Gromov C = Cartan A = Alexandrov T = Toponogov We are particularly interested in CAT(0) spaces.
11 Examples of Hadamard spaces 1 Hilbert spaces, the Hilbert ball 2 complete simply connected Riemannian manifolds with Sec 0 3 R-trees: a metric space T is an R-tree if for x, y T there is a unique geodesic [x, y] if [x, y] [y, z] = {y}, then [x, z] = [x, y] [y, z] 4 Euclidean buildings 5 the BHV tree space (space of phylogenetic trees) 6 L 2 (M, H), where (M, µ) is a probability space: ( d 2 (u, v) := M ) 1 d (u(x), v(x)) 2 2 dµ(x), u, v L 2 (M, H)
12 Convexity in Hadamard spaces Let (H, d) be a Hadamard space. These spaces allow for a natural definition of convexity: Definition A set C H is convex if, given x, y C, we have [x, y] C. Definition A function f : H (, ] is convex if f γ is a convex function for each geodesic γ : [0, 1] H.
13 Convexity in Hadamard spaces Let (H, d) be a Hadamard space. These spaces allow for a natural definition of convexity: Definition A set C H is convex if, given x, y C, we have [x, y] C. Definition A function f : H (, ] is convex if f γ is a convex function for each geodesic γ : [0, 1] H.
14 Convexity in Hadamard spaces Let (H, d) be a Hadamard space. These spaces allow for a natural definition of convexity: Definition A set C H is convex if, given x, y C, we have [x, y] C. Definition A function f : H (, ] is convex if f γ is a convex function for each geodesic γ : [0, 1] H.
15 Examples of convex functions 1 The indicator function of a convex closed set C H : ι C (x) := 0, if x C, and ι C (x) :=, if x / C. 2 The distance function to a closed convex subset C H : d C (x) := inf d(x, c), x H. c C 3 The displacement function of an isometry T : H H : δ T (x) := d(x, T x), x H.
16 Examples of convex functions 1 The indicator function of a convex closed set C H : ι C (x) := 0, if x C, and ι C (x) :=, if x / C. 2 The distance function to a closed convex subset C H : d C (x) := inf d(x, c), x H. c C 3 The displacement function of an isometry T : H H : δ T (x) := d(x, T x), x H.
17 Examples of convex functions 1 The indicator function of a convex closed set C H : ι C (x) := 0, if x C, and ι C (x) :=, if x / C. 2 The distance function to a closed convex subset C H : d C (x) := inf d(x, c), x H. c C 3 The displacement function of an isometry T : H H : δ T (x) := d(x, T x), x H.
18 Examples of convex functions 4 Let c: [0, ) H be a geodesic ray. The function b c : H R defined by b c (x) := lim t [d (x, c(t)) t], x H, is called the Busemann function associated to the ray c. 5 The energy of a mapping u: M H given by E(u) := d (u(x), u(y)) 2 p(x, dy)dµ(x), M M where (M, µ) is a measure space with a Markov kernel p(x, dy). E is convex continuous on L 2 (M, H).
19 Examples of convex functions 4 Let c: [0, ) H be a geodesic ray. The function b c : H R defined by b c (x) := lim t [d (x, c(t)) t], x H, is called the Busemann function associated to the ray c. 5 The energy of a mapping u: M H given by E(u) := d (u(x), u(y)) 2 p(x, dy)dµ(x), M M where (M, µ) is a measure space with a Markov kernel p(x, dy). E is convex continuous on L 2 (M, H).
20 Examples of convex functions 6 Given a 1,..., a N H and w 1,..., w N > 0, set f(x) := N w n d (x, a n ) p, x H, n=1 where p [1, ). If p = 1, we get Fermat-Weber problem for optimal facility location. A minimizer of f is called a median. If p = 2, then a minimizer of f is the barycenter of µ := or the mean of a 1,..., a N. N w n δ an, n=1
21 Examples of convex functions 6 Given a 1,..., a N H and w 1,..., w N > 0, set f(x) := N w n d (x, a n ) p, x H, n=1 where p [1, ). If p = 1, we get Fermat-Weber problem for optimal facility location. A minimizer of f is called a median. If p = 2, then a minimizer of f is the barycenter of µ := or the mean of a 1,..., a N. N w n δ an, n=1
22 Strongly convex functions A function f : H (, ] is strongly convex with parameter β > 0 if f ((1 t)x + ty) (1 t)f(x) + tf(y) βt(1 t)d(x, y) 2, for any x, y H and t [0, 1]. Each strongly has a unique minimizer. Example Given y H, the function f := d(y, ) 2 is strongly convex. Indeed, d (y, x t ) 2 (1 t)d (y, x 0 ) 2 + td (y, x 1 ) 2 t(1 t)d (x 0, x 1 ) 2, for each geodesic x: [0, 1] H.
23 1 Basic facts on Hadamard spaces 2 Proximal point algorithm 3 Applications to computational phylogenetics
24 Proximal point algorithm Let f : H (, ] be convex lsc. Optimization problem: min f(x). x H Recall: no (sub)differential, no shooting (singularities). Implicit methods are appropriate. The PPA generates a sequence [ x i := J λi (x i 1 ) := arg min f(y) + 1 ] d(y, x i 1 ) 2, y H 2λ i where x 0 H is a given starting point and λ i > 0, for each i N.
25 Proximal point algorithm Let f : H (, ] be convex lsc. Optimization problem: min f(x). x H Recall: no (sub)differential, no shooting (singularities). Implicit methods are appropriate. The PPA generates a sequence [ x i := J λi (x i 1 ) := arg min f(y) + 1 ] d(y, x i 1 ) 2, y H 2λ i where x 0 H is a given starting point and λ i > 0, for each i N.
26 Convergence of proximal point algorithm Theorem (M.B., 2011) Let f : H (, ] be a convex lsc function attaining its minimum. Given x 0 H and (λ i ) such that 1 λ i =, the PPA sequence (x i ) converges to a minimizer of f. (Resolvents are firmly nonexpansive - cheap version for λ i = λ.) Disadvantage: The resolvents x i := J λi (x i 1 ) := arg min y H are often difficult to compute. [ f(y) + 1 ] d(y, x i 1 ) 2, 2λ i
27 Convergence of proximal point algorithm Theorem (M.B., 2011) Let f : H (, ] be a convex lsc function attaining its minimum. Given x 0 H and (λ i ) such that 1 λ i =, the PPA sequence (x i ) converges to a minimizer of f. (Resolvents are firmly nonexpansive - cheap version for λ i = λ.) Disadvantage: The resolvents x i := J λi (x i 1 ) := arg min y H are often difficult to compute. [ f(y) + 1 ] d(y, x i 1 ) 2, 2λ i
28 Splitting proximal point algorithm Let f 1,..., f N be convex lsc and consider N f(x) := f n (x), x H. n=1 Example (Median and mean) f n := d (, a n ), f n := d (, a n ) 2. Key idea: apply resolvents J n λ s of f n s in a cyclic or random order.
29 Splitting proximal point algorithm Let f 1,..., f N be convex lsc and consider N f(x) := f n (x), x H. n=1 Example (Median and mean) f n := d (, a n ), f n := d (, a n ) 2. Key idea: apply resolvents J n λ s of f n s in a cyclic or random order.
30 Splitting proximal point algorithm Let f 1,..., f N be convex lsc and consider N f(x) := f n (x), x H. n=1 Example (Median and mean) f n := d (, a n ), f n := d (, a n ) 2. Key idea: apply resolvents J n λ s of f n s in a cyclic or random order.
31 Splitting proximal point algorithm Let x 0 H be a starting point. For each k N 0 we apply resolvents in cyclic order: x kn+1 := J 1 λ k (x kn ), x kn+2 := J 2 λ k (x kn+1 ),. x kn+n := J N λ k (x kn+n 1 ), or in random order: x i+1 := J r i λ i (x i ), where (r i ) are random variables with values in {1,..., N}.
32 Convergence of splitting proximal point algorithm Theorem (Cyclic order version + Random order version) Assume that f n are Lipschitz (or locally Lipschitz and the minimizing sequence is bounded). Then 1 the cyclic PPA sequence converges to a minimizer of f 2 the random PPA sequence converges to a minimizer of f almost surely. Assumptions are satisfied for f(x) := where p [1, ). N w n d (x, a n ) p, x H, n=1
33 Convergence of splitting proximal point algorithm Theorem (Cyclic order version + Random order version) Assume that f n are Lipschitz (or locally Lipschitz and the minimizing sequence is bounded). Then 1 the cyclic PPA sequence converges to a minimizer of f 2 the random PPA sequence converges to a minimizer of f almost surely. Assumptions are satisfied for f(x) := where p [1, ). N w n d (x, a n ) p, x H, n=1
34 Splitting proximal point algorithm (for mean) Hence instead of computing (the usual PPA) x i+1 := arg min z H [ N ] d (z, a n ) d (z, x i ) 2, 2λ i n=1 we are to minimize the function x i+1 := arg min z H [ d (z, a n ) λ i d (z, x i ) 2 where a n are chosen in a cyclic or random order. This is one-dimensional problem! = x i+1 is a convex combination of a n and x i. ],
35 Splitting proximal point algorithm (for mean) Hence instead of computing (the usual PPA) x i+1 := arg min z H [ N ] d (z, a n ) d (z, x i ) 2, 2λ i n=1 we are to minimize the function x i+1 := arg min z H [ d (z, a n ) λ i d (z, x i ) 2 where a n are chosen in a cyclic or random order. This is one-dimensional problem! = x i+1 is a convex combination of a n and x i. ],
36 1 Basic facts on Hadamard spaces 2 Proximal point algorithm 3 Applications to computational phylogenetics
37 Left: One of many evolutionary trees Right: A picture of an evolutionary tree by Charles Darwin (1837)
38 Billera-Holmes-Vogtmann tree space T d Figure: 5 out of 15 orthants of T 4
39 Billera-Holmes-Vogtmann tree space T d Figure: A finite set of trees in T 4
40 Billera-Holmes-Vogtmann tree space T d Figure: Consider the most frequent tree topology only
41 Computing the mean: Random order version Algorithm (SPPA with f n := d(, T n ) 2 ) Input: T 1,..., T N T d Step 1: S 1 := T 1 and i := 1 Step 2: choose r {1,..., N} at random Step 3: S i+1 := 1 i+1 T r + i i+1 S i Step 4: i := i + 1 Step 5: go to Step 2 Geodesics can be computed in polynomial time: The Owen-Provan algorithm (2011)
42 Computing the mean: Random order version Algorithm (SPPA with f n := d(, T n ) 2 ) Input: T 1,..., T N T d Step 1: S 1 := T 1 and i := 1 Step 2: choose r {1,..., N} at random Step 3: S i+1 := 1 i+1 T r + i i+1 S i Step 4: i := i + 1 Step 5: go to Step 2 Geodesics can be computed in polynomial time: The Owen-Provan algorithm (2011)
43 Computing the mean: Random order version Algorithm (SPPA with f n := d(, T n ) 2 ) Input: T 1,..., T N T d Step 1: S 1 := T 1 and i := 1 Step 2: choose r {1,..., N} at random Step 3: S i+1 := 1 i+1 T r + i i+1 S i Step 4: i := i + 1 Step 5: go to Step 2 Geodesics can be computed in polynomial time: The Owen-Provan algorithm (2011)
44 Computing the mean: Random order version T 5 T 4 T 6 T 3 T 7 T 1 T 2
45 Computing the mean: Random order version T 5 T 4 T 6 T 7 S 1 T 3 T 1 T 2
46 Computing the mean: Random order version T 5 T 4 T 6 T 7 S 1 S 2 T 3 T 1 T 2
47 Computing the mean: Random order version T 5 T 4 T 6 T 7 S 1 S 3 S 2 T 3 T 1 T 2
48 Computing the mean: Random order version T 5 T 4 T 6 T 7 S 1 S 4 S 3 S 2 T 3 T 1 T 2
49 Computing the mean: Random order version T 5 T 4 T 6 T 3 T 7 T 1 T 2
50 Computing the mean: Random order version T 5 T 4 T 6 T 3 T 7 T 1 T 2
51 Computing the mean: Random order version T 5 T 4 T 6 T 3 T 7 T 1 T 2
52 Computing the mean: Random order version T 5 T 4 T 6 T 3 T 7 T 1 T 2
53 Computing the mean: Random order version T 5 T 4 T 6 T 3 T 7 T 1 T 2
54 Le Big Data Space T d : orthant dimension = d 2, # of orthants = (2d 3)!! The actual dimension of T d is d (d 2)(2d 3)!! # of trees = N (e.g. coming from an MCMC simulation) Our computation: d = 12 hence dim and N = 10 5
55 Le Big Data Space T d : orthant dimension = d 2, # of orthants = (2d 3)!! The actual dimension of T d is d (d 2)(2d 3)!! # of trees = N (e.g. coming from an MCMC simulation) Our computation: d = 12 hence dim and N = 10 5
56 Resuls (courtesy of Philipp Benner) Pika Rabbit Orangutan Baboon Marmoset Gorilla Human Chimpanzee Kangaroo rat Squirrel Guinea pig Rat Mouse Pika Rabbit Squirrel Orangutan Gorilla Human Chimpanzee Baboon Guinea pig Marmoset Kangaroo rat Rat Mouse Pika Rabbit Kangaroo rat Rat Mouse Squirrel Guinea pig Orangutan Baboon Marmoset Gorilla Human Chimpanzee Pika Rabbit Squirrel Kangaroo rat Guinea pig Rat Mouse Orangutan Gorilla Human Chimpanzee Baboon Marmoset
57 Resuls (courtesy of Philipp Benner)... continued. Pika Rabbit Baboon Marmoset Orangutan Gorilla Chimpanzee Rat Human Mouse Kangaroo rat Guinea pig Squirrel Figure: Approximation of the mean of the 100,000 trees.
58 References 1 M.B.: Computing medians and means in Hadamard spaces, SIAM J. Optim. 24 (2014), no. 3, Benner, Bacak, Bourguignon: Point estimates in phylogenetic reconstructions, Bioinformatics, Vol 30 (2014), Issue M.B.: Convex analysis and optimization in Hadamard spaces, De Gruyter, Berlin, 2014.
A NOTION OF NONPOSITIVE CURVATURE FOR GENERAL METRIC SPACES
A NOTION OF NONPOSITIVE CURVATURE FOR GENERAL METRIC SPACES MIROSLAV BAČÁK, BOBO HUA, JÜRGEN JOST, MARTIN KELL, AND ARMIN SCHIKORRA Abstract. We introduce a new definition of nonpositive curvature in metric
More informationThe Metric Geometry of the Multivariable Matrix Geometric Mean
Trieste, 2013 p. 1/26 The Metric Geometry of the Multivariable Matrix Geometric Mean Jimmie Lawson Joint Work with Yongdo Lim Department of Mathematics Louisiana State University Baton Rouge, LA 70803,
More informationResults from MathSciNet: Mathematical Reviews on the Web c Copyright American Mathematical Society 2000
2000k:53038 53C23 20F65 53C70 57M07 Bridson, Martin R. (4-OX); Haefliger, André (CH-GENV-SM) Metric spaces of non-positive curvature. Grundlehren der Mathematischen Wissenschaften [Fundamental Principles
More informationAvailable online at J. Nonlinear Sci. Appl., 10 (2017), Research Article
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 (2017), 2719 2726 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa An affirmative answer to
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 informationProjection methods in geodesic metric spaces
Projection methods in geodesic metric spaces Computer Assisted Research Mathematics & Applications University of Newcastle http://carma.newcastle.edu.au/ 7th Asian Conference on Fixed Point Theory and
More informationPROXIMAL POINT ALGORITHMS INVOLVING FIXED POINT OF NONSPREADING-TYPE MULTIVALUED MAPPINGS IN HILBERT SPACES
PROXIMAL POINT ALGORITHMS INVOLVING FIXED POINT OF NONSPREADING-TYPE MULTIVALUED MAPPINGS IN HILBERT SPACES Shih-sen Chang 1, Ding Ping Wu 2, Lin Wang 3,, Gang Wang 3 1 Center for General Educatin, China
More informationViscosity approximation methods for the implicit midpoint rule of nonexpansive mappings in CAT(0) Spaces
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 017, 386 394 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Viscosity approximation methods
More informationGeneralized monotonicity and convexity for locally Lipschitz functions on Hadamard manifolds
Generalized monotonicity and convexity for locally Lipschitz functions on Hadamard manifolds A. Barani 1 2 3 4 5 6 7 Abstract. In this paper we establish connections between some concepts of generalized
More informationBasic Properties of Metric and Normed Spaces
Basic Properties of Metric and Normed Spaces Computational and Metric Geometry Instructor: Yury Makarychev The second part of this course is about metric geometry. We will study metric spaces, low distortion
More informationMONOTONE OPERATORS ON BUSEMANN SPACES
MONOTONE OPERATORS ON BUSEMANN SPACES David Ariza-Ruiz Genaro López-Acedo Universidad de Sevilla Departamento de Análisis Matemático V Workshop in Metric Fixed Point Theory and Applications 15-17 Noviembre,
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 informationSpaces with Ricci curvature bounded from below
Spaces with Ricci curvature bounded from below Nicola Gigli March 10, 2014 Lessons Basics of optimal transport Definition of spaces with Ricci curvature bounded from below Analysis on spaces with Ricci
More informationA new proof of Gromov s theorem on groups of polynomial growth
A new proof of Gromov s theorem on groups of polynomial growth Bruce Kleiner Courant Institute NYU Groups as geometric objects Let G a group with a finite generating set S G. Assume that S is symmetric:
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
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 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 informationA Second Order Non-smooth Variational Model for Restoring Manifold-valued Images
A Second Order Non-smooth Variational Model for Restoring Manifold-valued Images Ronny Bergmann Mathematical Image Processing and Data Analysis Department of Mathematics University of Kaiserslautern July
More informationNegative sectional curvature and the product complex structure. Harish Sheshadri. Department of Mathematics Indian Institute of Science Bangalore
Negative sectional curvature and the product complex structure Harish Sheshadri Department of Mathematics Indian Institute of Science Bangalore Technical Report No. 2006/4 March 24, 2006 M ath. Res. Lett.
More informationLecture 1: Background on Convex Analysis
Lecture 1: Background on Convex Analysis John Duchi PCMI 2016 Outline I Convex sets 1.1 Definitions and examples 2.2 Basic properties 3.3 Projections onto convex sets 4.4 Separating and supporting hyperplanes
More informationScientific report for the period December November 2014
Scientific report for the period December 203 - November 204 Publications. Published papers D. Ariza-Ruiz, L. Leuştean, G. Lopez-Acedo, Firmly nonexpansive mappings in classes of geodesic spaces, Transactions
More informationICM 2014: The Structure and Meaning. of Ricci Curvature. Aaron Naber ICM 2014: Aaron Naber
Outline of Talk Background and Limit Spaces Structure of Spaces with Lower Ricci Regularity of Spaces with Bounded Ricci Characterizing Ricci Background: s (M n, g, x) n-dimensional pointed Riemannian
More informationg 2 (x) (1/3)M 1 = (1/3)(2/3)M.
COMPACTNESS If C R n is closed and bounded, then by B-W it is sequentially compact: any sequence of points in C has a subsequence converging to a point in C Conversely, any sequentially compact C R n is
More informationarxiv: v3 [q-bio.pe] 23 Apr 2018
CONVERGENCE OF RANDOM WALKS TO BROWNIAN MOTION ON CUBICAL COMPLEXES arxiv:1508.02906v3 [q-bio.pe] 23 Apr 2018 TOM M. W. NYE Abstract. Cubical complexes are metric spaces constructed by gluing together
More informationRANDOM FIELDS AND GEOMETRY. Robert Adler and Jonathan Taylor
RANDOM FIELDS AND GEOMETRY from the book of the same name by Robert Adler and Jonathan Taylor IE&M, Technion, Israel, Statistics, Stanford, US. ie.technion.ac.il/adler.phtml www-stat.stanford.edu/ jtaylor
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 informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 64 (2012) 643 650 Contents lists available at SciVerse ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa
More informationALEXANDER LYTCHAK & VIKTOR SCHROEDER
AFFINE FUNCTIONS ON CAT (κ)-spaces ALEXANDER LYTCHAK & VIKTOR SCHROEDER Abstract. We describe affine functions on spaces with an upper curvature bound. 1. introduction A map f : X Y between geodesic metric
More informationZair Ibragimov CSUF. Talk at Fullerton College July 14, Geometry of p-adic numbers. Zair Ibragimov CSUF. Valuations on Rational Numbers
integers Talk at Fullerton College July 14, 2011 integers Let Z = {..., 2, 1, 0, 1, 2,... } denote the integers. Let Q = {a/b : a, b Z and b > 0} denote the rational. We can add and multiply rational :
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationSection 6. Laplacian, volume and Hessian comparison theorems
Section 6. Laplacian, volume and Hessian comparison theorems Weimin Sheng December 27, 2009 Two fundamental results in Riemannian geometry are the Laplacian and Hessian comparison theorems for the distance
More informationFixed Point Theory in Reflexive Metric Spaces
Babeş-Bolyai University Cluj-Napoca Department of Mathematics and University of Seville Department of Mathematical Analysis Fixed Point Theory in Reflexive Metric Spaces Ph.D. Thesis Summary Adriana-Maria
More informationStochastic calculus and martingales on trees
Stochastic calculus and martingales on trees (Calcul stochastique et martingales sur les arbres) Jean Picard Laboratoire de Mathématiques Appliquées (CNRS UMR 6620) Université Blaise Pascal 63177 Aubière
More informationarxiv:math/ v1 [math.dg] 24 May 2005
arxiv:math/0505496v1 [math.dg] 24 May 2005 INF-CONVOLUTION AND REGULARIZATION OF CONVEX FUNCTIONS ON RIEMANNIAN MANIFOLDS OF NONPOSITIVE CURVATURE DANIEL AZAGRA AND JUAN FERRERA Abstract. We show how an
More informationKrasnoselskii-Mann Method for Multi-valued Non-Self Mappings in CAT(0) Spaces
Filomat 31:14 (2017), 4629 4640 https://doi.org/10.2298/fil1714629t Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Krasnoselskii-Mann
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 informationTree Space: Algorithms & Applications Part I. Megan Owen University of Waterloo
Tree Space: lgorithms & pplications Part I Megan Owen University of Waterloo Phylogenetic Trees a phylogenetic tree: GTTTGT GTTGTT GTGT GGTTGTT? 7 5 1 6 questions: how do we infer a tree from data? how
More informationVariational inequalities for set-valued vector fields on Riemannian manifolds
Variational inequalities for set-valued vector fields on Riemannian manifolds Chong LI Department of Mathematics Zhejiang University Joint with Jen-Chih YAO Chong LI (Zhejiang University) VI on RM 1 /
More informationOn the convergence of SP-iterative scheme for three multivalued nonexpansive mappings in CAT (κ) spaces
Palestine Journal of Mathematics Vol. 4(1) (2015), 73 83 Palestine Polytechnic University-PPU 2015 On the convergence of SP-iterative scheme for three multivalued nonexpansive mappings in CAT () spaces
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationProximal Point Method for a Class of Bregman Distances on Riemannian Manifolds
Proximal Point Method for a Class of Bregman Distances on Riemannian Manifolds E. A. Papa Quiroz and P. Roberto Oliveira PESC-COPPE Federal University of Rio de Janeiro Rio de Janeiro, Brazil erik@cos.ufrj.br,
More informationStrictly convex functions on complete Finsler manifolds
Proc. Indian Acad. Sci. (Math. Sci.) Vol. 126, No. 4, November 2016, pp. 623 627. DOI 10.1007/s12044-016-0307-2 Strictly convex functions on complete Finsler manifolds YOE ITOKAWA 1, KATSUHIRO SHIOHAMA
More informationarxiv: v4 [math.dg] 18 Jun 2015
SMOOTHING 3-DIMENSIONAL POLYHEDRAL SPACES NINA LEBEDEVA, VLADIMIR MATVEEV, ANTON PETRUNIN, AND VSEVOLOD SHEVCHISHIN arxiv:1411.0307v4 [math.dg] 18 Jun 2015 Abstract. We show that 3-dimensional polyhedral
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationLoos Symmetric Cones. Jimmie Lawson Louisiana State University. July, 2018
Louisiana State University July, 2018 Dedication I would like to dedicate this talk to Joachim Hilgert, whose 60th birthday we celebrate at this conference and with whom I researched and wrote a big blue
More informationParaconvex functions and paraconvex sets
STUDIA MATHEMATICA 184 (1) (2008) Paraconvex functions and paraconvex sets by Huynh Van Ngai (Qui Nhon) and Jean-Paul Penot (Pau) Abstract. We study a class of functions which contains both convex functions
More informationVariational methods for restoration of phase or orientation data
Variational methods for restoration of phase or orientation data Martin Storath joint works with Laurent Demaret, Michael Unser, Andreas Weinmann Image Analysis and Learning Group Universität Heidelberg
More informationGoebel and Kirk fixed point theorem for multivalued asymptotically nonexpansive mappings
CARPATHIAN J. MATH. 33 (2017), No. 3, 335-342 Online version at http://carpathian.ubm.ro Print Edition: ISSN 1584-2851 Online Edition: ISSN 1843-4401 Goebel and Kirk fixed point theorem for multivalued
More informationMATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:
MATH 51H Section 4 October 16, 2015 1 Continuity Recall what it means for a function between metric spaces to be continuous: Definition. Let (X, d X ), (Y, d Y ) be metric spaces. A function f : X Y is
More informationContinuity. Matt Rosenzweig
Continuity Matt Rosenzweig Contents 1 Continuity 1 1.1 Rudin Chapter 4 Exercises........................................ 1 1.1.1 Exercise 1............................................. 1 1.1.2 Exercise
More informationSOME ELEMENTARY GENERAL PRINCIPLES OF CONVEX ANALYSIS. A. Granas M. Lassonde. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 5, 1995, 23 37 SOME ELEMENTARY GENERAL PRINCIPLES OF CONVEX ANALYSIS A. Granas M. Lassonde Dedicated, with admiration,
More informationTopological properties
CHAPTER 4 Topological properties 1. Connectedness Definitions and examples Basic properties Connected components Connected versus path connected, again 2. Compactness Definition and first examples Topological
More informationFirst-order Methods for Geodesically Convex Optimization
JMLR: Workshop and Conference Proceedings vol 49:, 06 First-order Methods for Geodesically Convex Optimization Hongyi Zhang Suvrit Sra Laboratory for Information & Decision Systems, Massachusetts Institute
More information8.8. Codimension one isoperimetric inequalities Distortion of a subgroup in a group 283
Contents Preface xiii Chapter 1. Geometry and topology 1 1.1. Set-theoretic preliminaries 1 1.1.1. General notation 1 1.1.2. Growth rates of functions 2 1.1.3. Jensen s inequality 3 1.2. Measure and integral
More informationMetric Thickenings of Euclidean Submanifolds
Metric Thickenings of Euclidean Submanifolds Advisor: Dr. Henry Adams Committe: Dr. Chris Peterson, Dr. Daniel Cooley Joshua Mirth Masters Thesis Defense October 3, 2017 Introduction Motivation Can we
More informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationReflexive Metric Spaces and The Fixed Point Property
Reflexive Metric Spaces and The Fixed Point Property M.A. Khamsi Department of Mathematical Sciences University of Texas at El Paso El Paso, TX 79968-0514 mohamed@math.utep.edu 1 Introduction As for the
More informationMATH 4200 HW: PROBLEM SET FOUR: METRIC SPACES
MATH 4200 HW: PROBLEM SET FOUR: METRIC SPACES PETE L. CLARK 4. Metric Spaces (no more lulz) Directions: This week, please solve any seven problems. Next week, please solve seven more. Starred parts of
More informationHopf-Rinow and Hadamard Theorems
Summersemester 2015 University of Heidelberg Riemannian geometry seminar Hopf-Rinow and Hadamard Theorems by Sven Grützmacher supervised by: Dr. Gye-Seon Lee Prof. Dr. Anna Wienhard Contents Introduction..........................................
More informationSOME REMARKS ON THE SPACE OF DIFFERENCES OF SUBLINEAR FUNCTIONS
APPLICATIONES MATHEMATICAE 22,3 (1994), pp. 419 426 S. G. BARTELS and D. PALLASCHKE (Karlsruhe) SOME REMARKS ON THE SPACE OF DIFFERENCES OF SUBLINEAR FUNCTIONS Abstract. Two properties concerning the space
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 informationStrongly convex functions, Moreau envelopes and the generic nature of convex functions with strong minimizers
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part B Faculty of Engineering and Information Sciences 206 Strongly convex functions, Moreau envelopes
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 informationConvex Optimization Notes
Convex Optimization Notes Jonathan Siegel January 2017 1 Convex Analysis This section is devoted to the study of convex functions f : B R {+ } and convex sets U B, for B a Banach space. The case of B =
More informationFixed point sets of parabolic isometries of CAT(0)-spaces
Comment. Math. Helv. 81 (2006), 305 335 Commentarii Mathematici Helvetici Swiss Mathematical Society Fixed point sets of parabolic isometries of CAT(0)-spaces Koji Fujiwara, Koichi Nagano and Takashi Shioya
More informationTHE FUNDAMENTAL GROUP OF NON-NEGATIVELY CURVED MANIFOLDS David Wraith The aim of this article is to oer a brief survey of an interesting, yet accessib
THE FUNDAMENTAL GROUP OF NON-NEGATIVELY CURVED MANIFOLDS David Wraith The aim of this article is to oer a brief survey of an interesting, yet accessible line of research in Dierential Geometry. A fundamental
More informationOn a result of Pazy concerning the asymptotic behaviour of nonexpansive mappings
On a result of Pazy concerning the asymptotic behaviour of nonexpansive mappings Heinz H. Bauschke, Graeme R. Douglas, and Walaa M. Moursi September 8, 2015 (final version) Abstract In 1971, Pazy presented
More informationLaplace Operator and Heat Kernel for Shape Analysis
Laplace Operator and Heat Kernel for Shape Analysis Jian Sun Mathematical Sciences Center, Tsinghua University R kf := 2 f x 2 1 Laplace Operator on R k, the standard Laplace operator: R kf := div f +
More informationDIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17
DIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17 6. Geodesics A parametrized line γ : [a, b] R n in R n is straight (and the parametrization is uniform) if the vector γ (t) does not depend on t. Thus,
More informationNON POSITIVELY CURVED METRIC IN THE SPACE OF POSITIVE DEFINITE INFINITE MATRICES
REVISTA DE LA UNIÓN MATEMÁTICA ARGENTINA Volumen 48, Número 1, 2007, Páginas 7 15 NON POSITIVELY CURVED METRIC IN THE SPACE OF POSITIVE DEFINITE INFINITE MATRICES ESTEBAN ANDRUCHOW AND ALEJANDRO VARELA
More informationWARPED PRODUCT METRICS ON R 2. = f(x) 1. The Levi-Civita Connection We did this calculation in class. To quickly recap, by metric compatibility,
WARPED PRODUCT METRICS ON R 2 KEVIN WHYTE Let M be R 2 equipped with the metric : x = 1 = f(x) y and < x, y >= 0 1. The Levi-Civita Connection We did this calculation in class. To quickly recap, by metric
More informationMATH 54 - TOPOLOGY SUMMER 2015 FINAL EXAMINATION. Problem 1
MATH 54 - TOPOLOGY SUMMER 2015 FINAL EXAMINATION ELEMENTS OF SOLUTION Problem 1 1. Let X be a Hausdorff space and K 1, K 2 disjoint compact subsets of X. Prove that there exist disjoint open sets U 1 and
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 informationPLURISUBHARMONIC FUNCTIONS AND THE STRUCTURE OF COMPLETE KÄHLER MANIFOLDS WITH NONNEGATIVE CURVATURE
j. differential geometry 64 (2003) 457-524 PLURISUBHARMONIC FUNCTIONS AND THE STRUCTURE OF COMPLETE KÄHLER MANIFOLDS WITH NONNEGATIVE CURVATURE LEI NI & LUEN-FAI TAM Abstract In this paper, we study global
More informationRicci Curvature and Bochner Formula on Alexandrov Spaces
Ricci Curvature and Bochner Formula on Alexandrov Spaces Sun Yat-sen University March 18, 2013 (work with Prof. Xi-Ping Zhu) Contents Alexandrov Spaces Generalized Ricci Curvature Geometric and Analytic
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationBeyond Scalar Affinities for Network Analysis or Vector Diffusion Maps and the Connection Laplacian
Beyond Scalar Affinities for Network Analysis or Vector Diffusion Maps and the Connection Laplacian Amit Singer Princeton University Department of Mathematics and Program in Applied and Computational Mathematics
More informationConvergence of three-step iterations for Ciric-quasi contractive operator in CAT(0) spaces
Acta Univ. Sapientiae, Mathematica, 7, 1 (15) 89 15 DOI: 1.1515/ausm-15-6 Convergence of three-step iterations for Ciric-quasi contractive operator in CAT() spaces Gurucharan S. Saluja Department of Mathematics,
More informationExtending Lipschitz and Hölder maps between metric spaces
Extending Lipschitz and Hölder maps between metric spaces Shin-ichi OHTA Department of Mathematics, Faculty of Science, Kyoto University, Kyoto 606-8502, JAPAN e-mail: sohta@math.kyoto-u.ac.jp Abstract
More informationSmooth Dynamics 2. Problem Set Nr. 1. Instructor: Submitted by: Prof. Wilkinson Clark Butler. University of Chicago Winter 2013
Smooth Dynamics 2 Problem Set Nr. 1 University of Chicago Winter 2013 Instructor: Submitted by: Prof. Wilkinson Clark Butler Problem 1 Let M be a Riemannian manifold with metric, and Levi-Civita connection.
More informationOPTIMALITY CONDITIONS FOR GLOBAL MINIMA OF NONCONVEX FUNCTIONS ON RIEMANNIAN MANIFOLDS
OPTIMALITY CONDITIONS FOR GLOBAL MINIMA OF NONCONVEX FUNCTIONS ON RIEMANNIAN MANIFOLDS S. HOSSEINI Abstract. A version of Lagrange multipliers rule for locally Lipschitz functions is presented. Using Lagrange
More informationDifferential Geometry II Lecture 1: Introduction and Motivation
Differential Geometry II Lecture 1: Introduction and Motivation Robert Haslhofer 1 Content of this lecture This course is on Riemannian geometry and the geometry of submanifol. The goal of this first lecture
More informationPythagorean Property and Best Proximity Pair Theorems
isibang/ms/2013/32 November 25th, 2013 http://www.isibang.ac.in/ statmath/eprints Pythagorean Property and Best Proximity Pair Theorems Rafa Espínola, G. Sankara Raju Kosuru and P. Veeramani Indian Statistical
More informationGroups up to quasi-isometry
OSU November 29, 2007 1 Introduction 2 3 Topological methods in group theory via the fundamental group. group theory topology group Γ, a topological space X with π 1 (X) = Γ. Γ acts on the universal cover
More informationASYMPTOTIC INVARIANTS OF HADAMARD MANIFOLDS. Mohamad A. Hindawi. A Dissertation. Mathematics
ASYMPTOTIC INVARIANTS OF HADAMARD MANIFOLDS Mohamad A. Hindawi A Dissertation in Mathematics Presented to the Faculties of the University of Pennsylvania in Partial Fulfillment of the Requirements for
More informationGeometric Analysis, Karlovassi, Samos
Departments of Mathematics Geometric Analysis, Karlovassi, Samos 31 May - 4 June 2010 University of the Aegean Program Monday 31/05 Tuesday 1/06 Wed. 2 Thursday 3 Friday 4 9:00-10:00 9:30 Welcome Knieper
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationThe problems I left behind
The problems I left behind Kazimierz Goebel Maria Curie-Sk lodowska University, Lublin, Poland email: goebel@hektor.umcs.lublin.pl During over forty years of studying and working on problems of metric
More informationConvex Analysis and Economic Theory Winter 2018
Division of the Humanities and Social Sciences Ec 181 KC Border Convex Analysis and Economic Theory Winter 2018 Supplement A: Mathematical background A.1 Extended real numbers The extended real number
More informationConvergence Theorems of Approximate Proximal Point Algorithm for Zeroes of Maximal Monotone Operators in Hilbert Spaces 1
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 4, 175-186 Convergence Theorems of Approximate Proximal Point Algorithm for Zeroes of Maximal Monotone Operators in Hilbert Spaces 1 Haiyun Zhou Institute
More informationFunctional Analysis I
Functional Analysis I Course Notes by Stefan Richter Transcribed and Annotated by Gregory Zitelli Polar Decomposition Definition. An operator W B(H) is called a partial isometry if W x = X for all x (ker
More informationMA3051: Mathematical Analysis II
MA3051: Mathematical Analysis II Course Notes Stephen Wills 2014/15 Contents 0 Assumed Knowledge 2 Sets and maps................................ 2 Sequences and series............................. 2 Metric
More informationTEST CODE: PMB SYLLABUS
TEST CODE: PMB SYLLABUS Convergence and divergence of sequence and series; Cauchy sequence and completeness; Bolzano-Weierstrass theorem; continuity, uniform continuity, differentiability; directional
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 informationChapter 2 Convex Analysis
Chapter 2 Convex Analysis The theory of nonsmooth analysis is based on convex analysis. Thus, we start this chapter by giving basic concepts and results of convexity (for further readings see also [202,
More informationThe Hopf argument. Yves Coudene. IRMAR, Université Rennes 1, campus beaulieu, bat Rennes cedex, France
The Hopf argument Yves Coudene IRMAR, Université Rennes, campus beaulieu, bat.23 35042 Rennes cedex, France yves.coudene@univ-rennes.fr slightly updated from the published version in Journal of Modern
More informationThe Structure of Singular Spaces of Dimension 2
The Structure of Singular Spaces of Dimension 2 Chikako Mese Department of Mathematics Box 5657 Connecticut College 270 Mohegan Avenue New London, CT 06320 cmes@conncoll.edu Abstract: In this paper, we
More informationChapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves
Chapter 3 Riemannian Manifolds - I The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves embedded in Riemannian manifolds. A Riemannian manifold is an abstraction
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationBoundaries in non positive curvature
Boundaries in non positive curvature Maciej Czarnecki Uniwersytet Lódzki, Katedra Geometrii Lódź, Poland maczar@math.uni.lodz.pl January 10, 2017 1 Contents 0 Introduction 3 1 Preliminaries 4 1.1 Riemannian
More information