The Karcher Mean of Points on SO n
|
|
- Anissa Hart
- 5 years ago
- Views:
Transcription
1 The Karcher Mean of Points on SO n Knut Hüper joint work with Jonathan Manton (Univ. Melbourne) Knut.Hueper@nicta.com.au National ICT Australia Ltd. CESAME LLN, 15/7/04 p.1/25
2 Contents Introduction CESAME LLN, 15/7/04 p.2/25
3 Contents Introduction Centroids, Karcher mean, Fréchet mean CESAME LLN, 15/7/04 p.2/25
4 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case CESAME LLN, 15/7/04 p.2/25
5 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n CESAME LLN, 15/7/04 p.2/25
6 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n Radii of convexity and injectivity CESAME LLN, 15/7/04 p.2/25
7 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n Radii of convexity and injectivity Karcher mean on SO n CESAME LLN, 15/7/04 p.2/25
8 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n Radii of convexity and injectivity Karcher mean on SO n Cost function, gradient and Hessian CESAME LLN, 15/7/04 p.2/25
9 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n Radii of convexity and injectivity Karcher mean on SO n Cost function, gradient and Hessian Newton-type algorithm CESAME LLN, 15/7/04 p.2/25
10 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n Radii of convexity and injectivity Karcher mean on SO n Cost function, gradient and Hessian Newton-type algorithm Convergence results CESAME LLN, 15/7/04 p.2/25
11 Contents Introduction Centroids, Karcher mean, Fréchet mean The Euclidean case Motivation, why SO n Radii of convexity and injectivity Karcher mean on SO n Cost function, gradient and Hessian Newton-type algorithm Convergence results Discussion, outlook CESAME LLN, 15/7/04 p.2/25
12 Given x 1,..., x k R n. Several ways to define a centroid x C CESAME LLN, 15/7/04 p.3/25
13 Given x 1,..., x k R n. 1) As the sum Several ways to define a centroid x C x c := 1 n k i=1 x i. CESAME LLN, 15/7/04 p.3/25
14 Given x 1,..., x k R n. 1) As the sum Several ways to define a centroid x C x c := 1 n k i=1 x i. 2) Equivalently, to ask the vector sum to vanish. xx xx k CESAME LLN, 15/7/04 p.3/25
15 Several ways to define a centroid x C 3) (Appolonius of Perga) As unique minimum of x c := argmin x R k i=1 x x i 2. CESAME LLN, 15/7/04 p.4/25
16 Several ways to define a centroid x C 3) (Appolonius of Perga) As unique minimum of x c := argmin x R k i=1 x x i 2. 4) More generally, assign to each x i a mass m i, mi = 1. By induction x c1,2 = m 1x 1 + m 2 x 2 m 1 + m 2 x c1,2,3 = (m 1 + m 2 )x c1,2 + m 3 x 3 m 1 + m 2 + m 3,... Also works on spheres. CESAME LLN, 15/7/04 p.4/25
17 Several ways to define a centroid x C 5) Axiomatically: Let Φ : R n R n Ξ R n be a rule mapping points to its centroid. CESAME LLN, 15/7/04 p.5/25
18 Several ways to define a centroid x C 5) Axiomatically: Let Φ : R n R n Ξ R n be a rule mapping points to its centroid. Axioms: (A1) Φ is symmetric in its arguments. (A2) Φ is smooth. (A3) Φ commutes with the induced action of SE n on R n R n. (A4) If Ω R n is an open convex ball then Φ maps Ω Ω into Ω. CESAME LLN, 15/7/04 p.5/25
19 Axioms (A1) Centroid is independent of the ordering of the points. (A2) Small changes in the location of the points causes only small changes in x c. (A3) Invariance w.r.t. translation and rotation. (A4) Centroid lies in the "same region" as the points themselves. Especially, Φ(x,., x) = x. CESAME LLN, 15/7/04 p.6/25
20 Why centroids on manifolds? Engineering, Mathematics, Physics CESAME LLN, 15/7/04 p.7/25
21 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds CESAME LLN, 15/7/04 p.7/25
22 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics CESAME LLN, 15/7/04 p.7/25
23 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking CESAME LLN, 15/7/04 p.7/25
24 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking fuzzy control on manifolds (defuzzification) CESAME LLN, 15/7/04 p.7/25
25 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking fuzzy control on manifolds (defuzzification) smoothing data CESAME LLN, 15/7/04 p.7/25
26 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking fuzzy control on manifolds (defuzzification) smoothing data plate tectonics CESAME LLN, 15/7/04 p.7/25
27 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking fuzzy control on manifolds (defuzzification) smoothing data plate tectonics sequence dep. continuum modeling of DNA CESAME LLN, 15/7/04 p.7/25
28 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking fuzzy control on manifolds (defuzzification) smoothing data plate tectonics sequence dep. continuum modeling of DNA comparison theorems (diff. geometry) CESAME LLN, 15/7/04 p.7/25
29 Why centroids on manifolds? Engineering, Mathematics, Physics statistical inferences on manifolds pose estimation in vision and robotics shape analysis and shape tracking fuzzy control on manifolds (defuzzification) smoothing data plate tectonics sequence dep. continuum modeling of DNA comparison theorems (diff. geometry) stochastic flows of mass distributions on manifolds (jets in gravitational field) CESAME LLN, 15/7/04 p.7/25
30 The special orthogonal group SO n SO n := {X R n n X X = I, det X = 1}. CESAME LLN, 15/7/04 p.8/25
31 The special orthogonal group SO n Facts: SO n := {X R n n X X = I, det X = 1}. CESAME LLN, 15/7/04 p.8/25
32 The special orthogonal group SO n SO n := {X R n n X X = I, det X = 1}. Facts: a) SO n is a Lie group, CESAME LLN, 15/7/04 p.8/25
33 The special orthogonal group SO n SO n := {X R n n X X = I, det X = 1}. Facts: a) SO n is a Lie group, b) is in general not diffeomorphic to a sphere, CESAME LLN, 15/7/04 p.8/25
34 The special orthogonal group SO n SO n := {X R n n X X = I, det X = 1}. Facts: a) SO n is a Lie group, b) is in general not diffeomorphic to a sphere, c) can be equipped with a Riemannian metric, therefore notion of distance is available, CESAME LLN, 15/7/04 p.8/25
35 The special orthogonal group SO n SO n := {X R n n X X = I, det X = 1}. Facts: a) SO n is a Lie group, b) is in general not diffeomorphic to a sphere, c) can be equipped with a Riemannian metric, therefore notion of distance is available, d) is compact and connected, but in general not simply connected. CESAME LLN, 15/7/04 p.8/25
36 Geometry of SO n a) We think of SO n as a submanifold of R n n. CESAME LLN, 15/7/04 p.9/25
37 Geometry of SO n a) We think of SO n as a submanifold of R n n. b) Tangent space T X SO n = {XA A R n n, A = A}. CESAME LLN, 15/7/04 p.9/25
38 Geometry of SO n a) We think of SO n as a submanifold of R n n. b) Tangent space T X SO n = {XA A R n n, A = A}. c) (Scaled) Frobenius inner product on R n n restricts to U, V = 1 2 tr(v U) XU, XV = 1 2 tr(v U), U, V T X SO n. Gives Riemannian metric on SO n. CESAME LLN, 15/7/04 p.9/25
39 Geometry of SO n cont d d) Let X SO n, Ω = Ω R n n. CESAME LLN, 15/7/04 p.10/25
40 Geometry of SO n cont d d) Let X SO n, Ω = Ω R n n. γ : R SO n, t X e t Ω is a geodesic through X = γ(0). CESAME LLN, 15/7/04 p.10/25
41 Geometry of SO n cont d d) Let X SO n, Ω = Ω R n n. γ : R SO n, t X e t Ω is a geodesic through X = γ(0). T 0 γ(t), γ(t) 1 2 d t is minimal (for T not too large..) CESAME LLN, 15/7/04 p.10/25
42 Geometry of SO n cont d e) Squared distance between any two points X, Y SO n CESAME LLN, 15/7/04 p.11/25
43 Geometry of SO n cont d e) Squared distance between any two points X, Y SO n d 2 (X, Y ) = 1 2 min A = A exp(a)=x Y tr(aa ) = 1 2 tr(log(x Y )) 2 CESAME LLN, 15/7/04 p.11/25
44 Centroid of SO n by axioms Let Ξ SO n SO n be open and consider Φ : Ξ SO n. CESAME LLN, 15/7/04 p.12/25
45 Centroid of SO n by axioms Let Ξ SO n SO n be open and consider Φ : Ξ SO n. (A1) Φ is symmetric in its arguments. CESAME LLN, 15/7/04 p.12/25
46 Centroid of SO n by axioms Let Ξ SO n SO n be open and consider Φ : Ξ SO n. (A1) Φ is symmetric in its arguments. (A2) Φ is smooth. CESAME LLN, 15/7/04 p.12/25
47 Centroid of SO n by axioms Let Ξ SO n SO n be open and consider Φ : Ξ SO n. (A1) Φ is symmetric in its arguments. (A2) Φ is smooth. (A3) Φ commutes with left and right translation. CESAME LLN, 15/7/04 p.12/25
48 Centroid of SO n by axioms Let Ξ SO n SO n be open and consider Φ : Ξ SO n. (A1) Φ is symmetric in its arguments. (A2) Φ is smooth. (A3) Φ commutes with left and right translation. (A4) If Ω SO n is an open convex ball then Φ maps Ω Ω into Ω. CESAME LLN, 15/7/04 p.12/25
49 Notion of convexity Ω SO n is defined to be convex if for any X, Y SO n there is a unique geodesic wholly contained in Ω connecting X to Y and such that it is also the unique minimising geodesic in SO n connecting X to Y. CESAME LLN, 15/7/04 p.13/25
50 Notion of convexity Ω SO n is defined to be convex if for any X, Y SO n there is a unique geodesic wholly contained in Ω connecting X to Y and such that it is also the unique minimising geodesic in SO n connecting X to Y. A function f : Ω R is convex if for any geodesic γ : [0, 1] Ω, the function f γ : [0, 1] R is convex in the usual sense, that is, f(γ(t)) (1 t)f(γ(0)) + tf(γ(1)), t [0, 1]. CESAME LLN, 15/7/04 p.13/25
51 Notion of convexity cont d Maximal convex ball (centered at the identity I n ) CESAME LLN, 15/7/04 p.14/25
52 Notion of convexity cont d Maximal convex ball (centered at the identity I n ) B(I, r) = {X SO n d(i, X) < r}. r conv is the largest r s.t. B(I, r) is convex and d(i, X) is convex on B(I, r). CESAME LLN, 15/7/04 p.14/25
53 Notion of convexity cont d Maximal convex ball (centered at the identity I n ) B(I, r) = {X SO n d(i, X) < r}. r conv is the largest r s.t. B(I, r) is convex and d(i, X) is convex on B(I, r). Theorem: For SO n it holds r conv = π 2. CESAME LLN, 15/7/04 p.14/25
54 Injectivity radius For so n := {A R n n A = A} let exp : so n SO n, Ψ exp(ψ), and B(0, ρ) = {A so n 1 2 tr A A < ρ 2 }. CESAME LLN, 15/7/04 p.15/25
55 Injectivity radius For so n := {A R n n A = A} let and exp : so n SO n, Ψ exp(ψ), B(0, ρ) = {A so n 1 2 tr A A < ρ 2 }. The injectivity radius r inj of so n is the largest ρ s.t. exp B(0,ρ) is a diffeomorphism onto its image. CESAME LLN, 15/7/04 p.15/25
56 Injectivity radius For so n := {A R n n A = A} let and exp : so n SO n, Ψ exp(ψ), B(0, ρ) = {A so n 1 2 tr A A < ρ 2 }. The injectivity radius r inj of so n is the largest ρ s.t. exp B(0,ρ) is a diffeomorphism onto its image. Theorem: For so n it holds r inj = π. CESAME LLN, 15/7/04 p.15/25
57 Let Ω SO n be open. Karcher mean on SO n CESAME LLN, 15/7/04 p.16/25
58 Karcher mean on SO n Let Ω SO n be open. A Karcher mean of Q 1,..., Q k SO n is defined to be a minimiser of f : Ω R, f(x) = k i=1 d 2 (Q i, X). CESAME LLN, 15/7/04 p.16/25
59 Karcher mean on SO n Let Ω SO n be open. A Karcher mean of Q 1,..., Q k SO n is defined to be a minimiser of f : Ω R, f(x) = k i=1 d 2 (Q i, X). Existence, uniqueness? CESAME LLN, 15/7/04 p.16/25
60 Results Theorem (MH 04): The critical points of f : Ω R, f(x) = k i=1 d 2 (Q i, X) are precisely the solutions of k i=1 log(q i X) = 0. CESAME LLN, 15/7/04 p.17/25
61 Results Theorem (MH 04): The Karcher mean is well defined and satisfies axioms (A1)-(A4) of a centroid on the open set Ξ = Y SO n B(Y, π/2) B(Y, π/2). CESAME LLN, 15/7/04 p.18/25
62 Results Theorem (MH 04): The Hessian of f represented along geodesics d 2 d t2(f γ)(t) is always positive definite. t=0 CESAME LLN, 15/7/04 p.19/25
63 f, grad f and Hessian explicitly f : Ω R, f(x) = k d 2 (Q i, X) = k i=1 i=1 1 2 tr(log(x Q i )) 2. CESAME LLN, 15/7/04 p.20/25
64 f, grad f and Hessian explicitly f : Ω R, f(x) = k d 2 (Q i, X) = k i=1 i=1 1 2 tr(log(x Q i )) 2. D f(x)xa = k = i=1 2X k tr(log(q i X)A) log(q i X) i=1 }{{} =grad f(x), XA. CESAME LLN, 15/7/04 p.20/25
65 f, grad f and Hessian explicitly d 2 ( d ε 2f X e εa ) = ε=0 vec A H(X) vec A CESAME LLN, 15/7/04 p.21/25
66 f, grad f and Hessian explicitly d 2 ( d ε 2f X e εa ) = ε=0 vec A H(X) vec A with (n 2 n 2 ) matrix H(X) := k i=1 Z i (X) coth(z i (X)) CESAME LLN, 15/7/04 p.21/25
67 f, grad f and Hessian explicitly d 2 ( d ε 2f X e εa ) = ε=0 vec A H(X) vec A with (n 2 n 2 ) matrix H(X) := k i=1 Z i (X) coth(z i (X)) and Z i (X) := I n log(q i X) + log(q i X) I n 2 CESAME LLN, 15/7/04 p.21/25
68 Algorithm Given Q 1,., Q k SO n, compute a local minimum of f. Step 1: Set X SO n to an initial estimate. Step 2: Compute k i=1 log(q i X). Step 3: Stop if k i=1 log(q i X) is suff. small. Step 4: Compute the update direction vec A opt = (H(X)) 1 k i=1 vec(log(q i X)) Step 5: Set X := X e A opt. Step 6: Go to Step 2. CESAME LLN, 15/7/04 p.22/25
69 Results Theorem (MH 04): The algorithm is an intrinsic Newton method. CESAME LLN, 15/7/04 p.23/25
70 Results Theorem (MH 04): The algorithm is an intrinsic Newton method. Theorem: If the algorithm converges, then it converges locally quadratically fast. CESAME LLN, 15/7/04 p.23/25
71 Discussion, outlook Need simple test to ensure that update step in algorithm remains in open convex ball global convergence. CESAME LLN, 15/7/04 p.24/25
72 Discussion, outlook Need simple test to ensure that update step in algorithm remains in open convex ball global convergence. Different RM, e.g. Cayley-like, gives different function, geodesics, etc.., but typically KM cay KM exp << 1. CESAME LLN, 15/7/04 p.24/25
73 Discussion, outlook Need simple test to ensure that update step in algorithm remains in open convex ball global convergence. Different RM, e.g. Cayley-like, gives different function, geodesics, etc.., but typically KM cay KM exp << 1. (H(X)) 1 via EVD. CESAME LLN, 15/7/04 p.24/25
74 Discussion, outlook Need simple test to ensure that update step in algorithm remains in open convex ball global convergence. Different RM, e.g. Cayley-like, gives different function, geodesics, etc.., but typically KM cay KM exp << 1. (H(X)) 1 via EVD. Quasi-Newton (rank-one updates). CESAME LLN, 15/7/04 p.24/25
75 Discussion, outlook Linear convergent algorithm (joint work with Robert Orsi, ANU) X i+1 = X i e 1 k k j=1 log(x i Q j) CESAME LLN, 15/7/04 p.25/25
76 Discussion, outlook Linear convergent algorithm (joint work with Robert Orsi, ANU) X i+1 = X i e 1 k k j=1 log(x i Q j) Centroids on homogeneous (symmetric) spaces. CESAME LLN, 15/7/04 p.25/25
77 Discussion, outlook Linear convergent algorithm (joint work with Robert Orsi, ANU) X i+1 = X i e 1 k k j=1 log(x i Q j) Centroids on homogeneous (symmetric) spaces. Project with NICTA vision/robotic program (Richard Hartley) to treat SE 3 case. CESAME LLN, 15/7/04 p.25/25
Optimisation on Manifolds
Optimisation on Manifolds K. Hüper MPI Tübingen & Univ. Würzburg K. Hüper (MPI Tübingen & Univ. Würzburg) Applications in Computer Vision Grenoble 18/9/08 1 / 29 Contents 2 Examples Essential matrix estimation
More informationEssential Matrix Estimation via Newton-type Methods
Essential Matrix Estimation via Newton-type Methods Uwe Helmke 1, Knut Hüper, Pei Yean Lee 3, John Moore 3 1 Abstract In this paper camera parameters are assumed to be known and a novel approach for essential
More informationCALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =
CALCULUS ON MANIFOLDS 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = a M T am, called the tangent bundle, is itself a smooth manifold, dim T M = 2n. Example 1.
More informationII. DIFFERENTIABLE MANIFOLDS. Washington Mio CENTER FOR APPLIED VISION AND IMAGING SCIENCES
II. DIFFERENTIABLE MANIFOLDS Washington Mio Anuj Srivastava and Xiuwen Liu (Illustrations by D. Badlyans) CENTER FOR APPLIED VISION AND IMAGING SCIENCES Florida State University WHY MANIFOLDS? Non-linearity
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationINTRINSIC MEAN ON MANIFOLDS. Abhishek Bhattacharya Project Advisor: Dr.Rabi Bhattacharya
INTRINSIC MEAN ON MANIFOLDS Abhishek Bhattacharya Project Advisor: Dr.Rabi Bhattacharya 1 Overview Properties of Intrinsic mean on Riemannian manifolds have been presented. The results have been applied
More informationA few words about the MTW tensor
A few words about the Ma-Trudinger-Wang tensor Université Nice - Sophia Antipolis & Institut Universitaire de France Salah Baouendi Memorial Conference (Tunis, March 2014) The Ma-Trudinger-Wang tensor
More informationGeodesic Equivalence in sub-riemannian Geometry
03/27/14 Equivalence in sub-riemannian Geometry Supervisor: Dr.Igor Zelenko Texas A&M University, Mathematics Some Preliminaries: Riemannian Metrics Let M be a n-dimensional surface in R N Some Preliminaries:
More informationMicrolocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries
Microlocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries David Dos Santos Ferreira LAGA Université de Paris 13 Wednesday May 18 Instituto de Ciencias Matemáticas,
More informationPractice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009.
Practice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009. Solutions (1) Let Γ be a discrete group acting on a manifold M. (a) Define what it means for Γ to act freely. Solution: Γ acts
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 informationfy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))
1. Basic algebra of vector fields Let V be a finite dimensional vector space over R. Recall that V = {L : V R} is defined to be the set of all linear maps to R. V is isomorphic to V, but there is no canonical
More informationDIFFERENTIAL GEOMETRY HW 12
DIFFERENTIAL GEOMETRY HW 1 CLAY SHONKWILER 3 Find the Lie algebra so(n) of the special orthogonal group SO(n), and the explicit formula for the Lie bracket there. Proof. Since SO(n) is a subgroup of GL(n),
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More informationMany of the exercises are taken from the books referred at the end of the document.
Exercises in Geometry I University of Bonn, Winter semester 2014/15 Prof. Christian Blohmann Assistant: Néstor León Delgado The collection of exercises here presented corresponds to the exercises for the
More informationThe nonsmooth Newton method on Riemannian manifolds
The nonsmooth Newton method on Riemannian manifolds C. Lageman, U. Helmke, J.H. Manton 1 Introduction Solving nonlinear equations in Euclidean space is a frequently occurring problem in optimization and
More informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationLECTURE 16: CONJUGATE AND CUT POINTS
LECTURE 16: CONJUGATE AND CUT POINTS 1. Conjugate Points Let (M, g) be Riemannian and γ : [a, b] M a geodesic. Then by definition, exp p ((t a) γ(a)) = γ(t). We know that exp p is a diffeomorphism near
More informationMath Topology II: Smooth Manifolds. Spring Homework 2 Solution Submit solutions to the following problems:
Math 132 - Topology II: Smooth Manifolds. Spring 2017. Homework 2 Solution Submit solutions to the following problems: 1. Let H = {a + bi + cj + dk (a, b, c, d) R 4 }, where i 2 = j 2 = k 2 = 1, ij = k,
More informationSYMPLECTIC GEOMETRY: LECTURE 5
SYMPLECTIC GEOMETRY: LECTURE 5 LIAT KESSLER Let (M, ω) be a connected compact symplectic manifold, T a torus, T M M a Hamiltonian action of T on M, and Φ: M t the assoaciated moment map. Theorem 0.1 (The
More informationEinstein-Hilbert action on Connes-Landi noncommutative manifolds
Einstein-Hilbert action on Connes-Landi noncommutative manifolds Yang Liu MPIM, Bonn Analysis, Noncommutative Geometry, Operator Algebras Workshop June 2017 Motivations and History Motivation: Explore
More informationRepresentation of Markov chains. Stochastic stability AIMS 2014
Markov chain model Joint with J. Jost and M. Kell Max Planck Institute for Mathematics in the Sciences Leipzig 0th July 204 AIMS 204 We consider f : M! M to be C r for r 0 and a small perturbation parameter
More informationTransport Continuity Property
On Riemannian manifolds satisfying the Transport Continuity Property Université de Nice - Sophia Antipolis (Joint work with A. Figalli and C. Villani) I. Statement of the problem Optimal transport on Riemannian
More informationDifferential Topology Final Exam With Solutions
Differential Topology Final Exam With Solutions Instructor: W. D. Gillam Date: Friday, May 20, 2016, 13:00 (1) Let X be a subset of R n, Y a subset of R m. Give the definitions of... (a) smooth function
More informationIntersection of stable and unstable manifolds for invariant Morse functions
Intersection of stable and unstable manifolds for invariant Morse functions Hitoshi Yamanaka (Osaka City University) March 14, 2011 Hitoshi Yamanaka (Osaka City University) ()Intersection of stable and
More informationDEVELOPMENT OF MORSE THEORY
DEVELOPMENT OF MORSE THEORY MATTHEW STEED Abstract. In this paper, we develop Morse theory, which allows us to determine topological information about manifolds using certain real-valued functions defined
More informationLecture 11. Geodesics and completeness
Lecture 11. Geodesics and completeness In this lecture we will investigate further the metric properties of geodesics of the Levi-Civita connection, and use this to characterise completeness of a Riemannian
More informationTheorem 2. Let n 0 3 be a given integer. is rigid in the sense of Guillemin, so are all the spaces ḠR n,n, with n n 0.
This monograph is motivated by a fundamental rigidity problem in Riemannian geometry: determine whether the metric of a given Riemannian symmetric space of compact type can be characterized by means of
More informationManifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds
MA 755 Fall 05. Notes #1. I. Kogan. Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds Definition 1 An n-dimensional C k -differentiable manifold
More informationLECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups
LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS 1. Lie groups A Lie group is a special smooth manifold on which there is a group structure, and moreover, the two structures are compatible. Lie groups are
More informationOnline Learning on Spheres
Online Learning on Spheres Krzysztof A. Krakowski 1, Robert E. Mahony 1, Robert C. Williamson 2, and Manfred K. Warmuth 3 1 Department of Engineering, Australian National University, Canberra, ACT 0200,
More informationTravel Time Tomography and Tensor Tomography, I
Travel Time Tomography and Tensor Tomography, I Plamen Stefanov Purdue University Mini Course, MSRI 2009 Plamen Stefanov (Purdue University ) Travel Time Tomography and Tensor Tomography, I 1 / 17 Alternative
More informationNumerical Minimization of Potential Energies on Specific Manifolds
Numerical Minimization of Potential Energies on Specific Manifolds SIAM Conference on Applied Linear Algebra 2012 22 June 2012, Valencia Manuel Gra f 1 1 Chemnitz University of Technology, Germany, supported
More informationLecture 8. Connections
Lecture 8. Connections This lecture introduces connections, which are the machinery required to allow differentiation of vector fields. 8.1 Differentiating vector fields. The idea of differentiating vector
More informationKilling fields of constant length on homogeneous Riemannian manifolds
Killing fields of constant length on homogeneous Riemannian manifolds Southern Mathematical Institute VSC RAS Poland, Bedlewo, 21 October 2015 1 Introduction 2 3 4 Introduction Killing vector fields (simply
More informationPogorelov Klingenberg theorem for manifolds homeomorphic to R n (translated from [4])
Pogorelov Klingenberg theorem for manifolds homeomorphic to R n (translated from [4]) Vladimir Sharafutdinov January 2006, Seattle 1 Introduction The paper is devoted to the proof of the following Theorem
More informationPseudo-Poincaré Inequalities and Applications to Sobolev Inequalities
Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities Laurent Saloff-Coste Abstract Most smoothing procedures are via averaging. Pseudo-Poincaré inequalities give a basic L p -norm control
More informationH-convex Riemannian submanifolds
H-convex Riemannian submanifolds Constantin Udrişte and Teodor Oprea Abstract. Having in mind the well known model of Euclidean convex hypersurfaces [4], [5] and the ideas in [1], many authors defined
More informationMostow Rigidity. W. Dison June 17, (a) semi-simple Lie groups with trivial centre and no compact factors and
Mostow Rigidity W. Dison June 17, 2005 0 Introduction Lie Groups and Symmetric Spaces We will be concerned with (a) semi-simple Lie groups with trivial centre and no compact factors and (b) simply connected,
More informationGeodesic shooting on shape spaces
Geodesic shooting on shape spaces Alain Trouvé CMLA, Ecole Normale Supérieure de Cachan GDR MIA Paris, November 21 2014 Outline Riemannian manifolds (finite dimensional) Spaces spaces (intrinsic metrics)
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
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 informationWelcome to Copenhagen!
Welcome to Copenhagen! Schedule: Monday Tuesday Wednesday Thursday Friday 8 Registration and welcome 9 Crash course on Crash course on Introduction to Differential and Differential and Information Geometry
More informationChoice of Riemannian Metrics for Rigid Body Kinematics
Choice of Riemannian Metrics for Rigid Body Kinematics Miloš Žefran1, Vijay Kumar 1 and Christopher Croke 2 1 General Robotics and Active Sensory Perception (GRASP) Laboratory 2 Department of Mathematics
More informationLECTURE 10: THE ATIYAH-GUILLEMIN-STERNBERG CONVEXITY THEOREM
LECTURE 10: THE ATIYAH-GUILLEMIN-STERNBERG CONVEXITY THEOREM Contents 1. The Atiyah-Guillemin-Sternberg Convexity Theorem 1 2. Proof of the Atiyah-Guillemin-Sternberg Convexity theorem 3 3. Morse theory
More informationThe theory of manifolds Lecture 3. Tf : TR n TR m
The theory of manifolds Lecture 3 Definition 1. The tangent space of an open set U R n, TU is the set of pairs (x, v) U R n. This should be thought of as a vector v based at the point x U. Denote by T
More informationHOMEWORK 2 - RIEMANNIAN GEOMETRY. 1. Problems In what follows (M, g) will always denote a Riemannian manifold with a Levi-Civita connection.
HOMEWORK 2 - RIEMANNIAN GEOMETRY ANDRÉ NEVES 1. Problems In what follows (M, g will always denote a Riemannian manifold with a Levi-Civita connection. 1 Let X, Y, Z be vector fields on M so that X(p Z(p
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 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 information4.7 The Levi-Civita connection and parallel transport
Classnotes: Geometry & Control of Dynamical Systems, M. Kawski. April 21, 2009 138 4.7 The Levi-Civita connection and parallel transport In the earlier investigation, characterizing the shortest curves
More informationThe Square Root Velocity Framework for Curves in a Homogeneous Space
The Square Root Velocity Framework for Curves in a Homogeneous Space Zhe Su 1 Eric Klassen 1 Martin Bauer 1 1 Florida State University October 8, 2017 Zhe Su (FSU) SRVF for Curves in a Homogeneous Space
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More information1 Differentiable manifolds and smooth maps
1 Differentiable manifolds and smooth maps Last updated: April 14, 2011. 1.1 Examples and definitions Roughly, manifolds are sets where one can introduce coordinates. An n-dimensional manifold is a set
More informationDetecting submanifolds of minimum volume with calibrations
Detecting submanifolds of minimum volume with calibrations Marcos Salvai, CIEM - FaMAF, Córdoba, Argentina http://www.famaf.unc.edu.ar/ salvai/ EGEO 2016, Córdoba, August 2016 It is not easy to choose
More information1. Geometry of the unit tangent bundle
1 1. Geometry of the unit tangent bundle The main reference for this section is [8]. In the following, we consider (M, g) an n-dimensional smooth manifold endowed with a Riemannian metric g. 1.1. Notations
More informationOn isometries of Finsler manifolds
On isometries of Finsler manifolds International Conference on Finsler Geometry February 15-16, 2008, Indianapolis, IN, USA László Kozma (University of Debrecen, Hungary) Finsler metrics, examples isometries
More informationHow curvature shapes space
How curvature shapes space Richard Schoen University of California, Irvine - Hopf Lecture, ETH, Zürich - October 30, 2017 The lecture will have three parts: Part 1: Heinz Hopf and Riemannian geometry Part
More information2 Lie Groups. Contents
2 Lie Groups Contents 2.1 Algebraic Properties 25 2.2 Topological Properties 27 2.3 Unification of Algebra and Topology 29 2.4 Unexpected Simplification 31 2.5 Conclusion 31 2.6 Problems 32 Lie groups
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 informationGeneralized Weiszfeld Algorithms for Lq Optimization
1 Generalized Weiszfeld Algorithms for Lq Optimization Khurrum Aftab, Richard Hartley, Jochen Trumpf Abstract In many computer vision applications, a desired model of some type is computed by minimizing
More informationA LOWER BOUND ON THE SUBRIEMANNIAN DISTANCE FOR HÖLDER DISTRIBUTIONS
A LOWER BOUND ON THE SUBRIEMANNIAN DISTANCE FOR HÖLDER DISTRIBUTIONS SLOBODAN N. SIMIĆ Abstract. Whereas subriemannian geometry usually deals with smooth horizontal distributions, partially hyperbolic
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationB.7 Lie Groups and Differential Equations
96 B.7. LIE GROUPS AND DIFFERENTIAL EQUATIONS B.7 Lie Groups and Differential Equations Peter J. Olver in Minneapolis, MN (U.S.A.) mailto:olver@ima.umn.edu The applications of Lie groups to solve differential
More informationRecent progress on the explicit inversion of geodesic X-ray transforms
Recent progress on the explicit inversion of geodesic X-ray transforms François Monard Department of Mathematics, University of Washington. Geometric Analysis and PDE seminar University of Cambridge, May
More informationLECTURE 9: MOVING FRAMES IN THE NONHOMOGENOUS CASE: FRAME BUNDLES. 1. Introduction
LECTURE 9: MOVING FRAMES IN THE NONHOMOGENOUS CASE: FRAME BUNDLES 1. Introduction Until now we have been considering homogenous spaces G/H where G is a Lie group and H is a closed subgroup. The natural
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 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 informationRiemannian geometry of surfaces
Riemannian geometry of surfaces In this note, we will learn how to make sense of the concepts of differential geometry on a surface M, which is not necessarily situated in R 3. This intrinsic approach
More informationη = (e 1 (e 2 φ)) # = e 3
Research Statement My research interests lie in differential geometry and geometric analysis. My work has concentrated according to two themes. The first is the study of submanifolds of spaces with riemannian
More informationChanging sign solutions for the CR-Yamabe equation
Changing sign solutions for the CR-Yamabe equation Ali Maalaoui (1) & Vittorio Martino (2) Abstract In this paper we prove that the CR-Yamabe equation on the Heisenberg group has infinitely many changing
More informationTRANSITIVE HOLONOMY GROUP AND RIGIDITY IN NONNEGATIVE CURVATURE. Luis Guijarro and Gerard Walschap
TRANSITIVE HOLONOMY GROUP AND RIGIDITY IN NONNEGATIVE CURVATURE Luis Guijarro and Gerard Walschap Abstract. In this note, we examine the relationship between the twisting of a vector bundle ξ over a manifold
More informationTerse Notes on Riemannian Geometry
Terse Notes on Riemannian Geometry Tom Fletcher January 26, 2010 These notes cover the basics of Riemannian geometry, Lie groups, and symmetric spaces. This is just a listing of the basic definitions and
More informationThe theory of manifolds Lecture 2
The theory of manifolds Lecture 2 Let X be a subset of R N, Y a subset of R n and f : X Y a continuous map. We recall Definition 1. f is a C map if for every p X, there exists a neighborhood, U p, of p
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More information10. The subgroup subalgebra correspondence. Homogeneous spaces.
10. The subgroup subalgebra correspondence. Homogeneous spaces. 10.1. The concept of a Lie subgroup of a Lie group. We have seen that if G is a Lie group and H G a subgroup which is at the same time a
More informationDistances, volumes, and integration
Distances, volumes, and integration Scribe: Aric Bartle 1 Local Shape of a Surface A question that we may ask ourselves is what significance does the second fundamental form play in the geometric characteristics
More informationQuasi-conformal minimal Lagrangian diffeomorphisms of the
Quasi-conformal minimal Lagrangian diffeomorphisms of the hyperbolic plane (joint work with J.M. Schlenker) January 21, 2010 Quasi-symmetric homeomorphism of a circle A homeomorphism φ : S 1 S 1 is quasi-symmetric
More informationCoordinate Systems and Canonical Forms
Appendix D Coordinate Systems and Canonical Forms D.1. Local Coordinates Let O be an open set in R n. We say that an n-tuple of smooth realvalued functions defined in O, (φ 1,...,φ n ), forms a local coordinate
More informationMILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES
MILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES NILAY KUMAR In these lectures I want to introduce the Chern-Weil approach to characteristic classes on manifolds, and in particular, the Chern classes.
More informationAdapted complex structures and Riemannian homogeneous spaces
ANNALES POLONICI MATHEMATICI LXX (1998) Adapted complex structures and Riemannian homogeneous spaces by Róbert Szőke (Budapest) Abstract. We prove that every compact, normal Riemannian homogeneous manifold
More informationDifferential Geometry Exercises
Differential Geometry Exercises Isaac Chavel Spring 2006 Jordan curve theorem We think of a regular C 2 simply closed path in the plane as a C 2 imbedding of the circle ω : S 1 R 2. Theorem. Given the
More informationclass # MATH 7711, AUTUMN 2017 M-W-F 3:00 p.m., BE 128 A DAY-BY-DAY LIST OF TOPICS
class # 34477 MATH 7711, AUTUMN 2017 M-W-F 3:00 p.m., BE 128 A DAY-BY-DAY LIST OF TOPICS [DG] stands for Differential Geometry at https://people.math.osu.edu/derdzinski.1/courses/851-852-notes.pdf [DFT]
More informationAFFINE MAXIMAL HYPERSURFACES. Xu-Jia Wang. Centre for Mathematics and Its Applications The Australian National University
AFFINE MAXIMAL HYPERSURFACES Xu-Jia Wang Centre for Mathematics and Its Applications The Australian National University Abstract. This is a brief survey of recent works by Neil Trudinger and myself on
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 informationThe Inverse Function Theorem via Newton s Method. Michael Taylor
The Inverse Function Theorem via Newton s Method Michael Taylor We aim to prove the following result, known as the inverse function theorem. Theorem 1. Let F be a C k map (k 1) from a neighborhood of p
More informationarxiv: v4 [math.dg] 10 Feb 2018
CYLINDERS AS LEFT INVARIANT CMC SURFACES IN Sol 3 AND E(κ, τ)-spaces DIFFEOMORPHIC TO R 3 MIROSLAV VRŽINA arxiv:1412.6820v4 [math.dg] 10 Feb 2018 ABSTRACT. In the present paper we give a geometric proof
More informationDifferential Geometry qualifying exam 562 January 2019 Show all your work for full credit
Differential Geometry qualifying exam 562 January 2019 Show all your work for full credit 1. (a) Show that the set M R 3 defined by the equation (1 z 2 )(x 2 + y 2 ) = 1 is a smooth submanifold of R 3.
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 informationTHE INVERSE FUNCTION THEOREM FOR LIPSCHITZ MAPS
THE INVERSE FUNCTION THEOREM FOR LIPSCHITZ MAPS RALPH HOWARD DEPARTMENT OF MATHEMATICS UNIVERSITY OF SOUTH CAROLINA COLUMBIA, S.C. 29208, USA HOWARD@MATH.SC.EDU Abstract. This is an edited version of a
More informationTHE EULER CHARACTERISTIC OF A LIE GROUP
THE EULER CHARACTERISTIC OF A LIE GROUP JAY TAYLOR 1 Examples of Lie Groups The following is adapted from [2] We begin with the basic definition and some core examples Definition A Lie group is a smooth
More informationREGULAR TRIPLETS IN COMPACT SYMMETRIC SPACES
REGULAR TRIPLETS IN COMPACT SYMMETRIC SPACES MAKIKO SUMI TANAKA 1. Introduction This article is based on the collaboration with Tadashi Nagano. In the first part of this article we briefly review basic
More informationThe Laplace-Beltrami-Operator on Riemannian Manifolds. 1 Why do we need the Laplace-Beltrami-Operator?
Frank Schmidt Computer Vision Group - Technische Universität ünchen Abstract This report mainly illustrates a way to compute the Laplace-Beltrami-Operator on a Riemannian anifold and gives information
More informationCHAPTER 1 PRELIMINARIES
CHAPTER 1 PRELIMINARIES 1.1 Introduction The aim of this chapter is to give basic concepts, preliminary notions and some results which we shall use in the subsequent chapters of the thesis. 1.2 Differentiable
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 17
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 17 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory May 29, 2012 Andre Tkacenko
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationAPPLICATIONS OF DIFFERENTIABILITY IN R n.
APPLICATIONS OF DIFFERENTIABILITY IN R n. MATANIA BEN-ARTZI April 2015 Functions here are defined on a subset T R n and take values in R m, where m can be smaller, equal or greater than n. The (open) ball
More informationNOTES ON DIFFERENTIAL FORMS. PART 5: DE RHAM COHOMOLOGY
NOTES ON DIFFERENTIAL FORMS. PART 5: DE RHAM COHOMOLOGY 1. Closed and exact forms Let X be a n-manifold (not necessarily oriented), and let α be a k-form on X. We say that α is closed if dα = 0 and say
More informationEigenvalue (mis)behavior on manifolds
Bucknell University Lehigh University October 20, 2010 Outline 1 Isoperimetric inequalities 2 3 4 A little history Rayleigh quotients The Original Isoperimetric Inequality The Problem of Queen Dido: maximize
More informationLecture 4 - The Basic Examples of Collapse
Lecture 4 - The Basic Examples of Collapse July 29, 2009 1 Berger Spheres Let X, Y, and Z be the left-invariant vector fields on S 3 that restrict to i, j, and k at the identity. This is a global frame
More informationComparison for infinitesimal automorphisms. of parabolic geometries
Comparison techniques for infinitesimal automorphisms of parabolic geometries University of Vienna Faculty of Mathematics Srni, January 2012 This talk reports on joint work in progress with Karin Melnick
More information