Numerical Minimization of Potential Energies on Specific Manifolds
|
|
- Annice Andrews
- 5 years ago
- Views:
Transcription
1 Numerical Minimization of Potential Energies on Specific Manifolds SIAM Conference on Applied Linear Algebra June 2012, Valencia Manuel Gra f 1 1 Chemnitz University of Technology, Germany, supported by DFG
2 Outline 1 The General Problem 2 Optimization on Riemannian Manifolds 3 Efficient Evaluation Approaches 4 Discrepancies and Potential Energies 5 Numerical Examples
3 The General Problem General Task For a given compact set X R n, and functions K : X X R, h : X R minimize the potential energy E(P, w) := M M w i w j K(p i, p j ) 2 w i h(p i ) + C i,j=1 for points P := (p i ) M i=1 X M and weights w := (w i ) M i=1 Ω RM Applications and Relations optimal quadrature points in reproducing kernel Hilbert spaces (spherical t-designs, low-discrepancy points, Gauß-quadrature) halftoning/dithering of images i=1 optimal point configurations (Thomson problem on S 2 )...
4 Optimization on Riemannian Manifolds minimization of E : X M Ω R for manifolds X R n, Ω R M by a nonlinear CG method on Riemannian manifolds Advantages no boundary issues for compact manifolds the method is invariant under change of coordinate system no need of local coordinates for specific manifolds and optimization problems group symmetries may be easily incorporated superlinear convergence rate Disadvantages calculations of geodesics may be too expansive or impractically
5 Optimization on Riemannian Manifolds Tangent Spaces and Geodesic Curves let M R n be a compact d-dimensional manifold M is characterized by smooth curves γ : R M tangent vector v R n at point x M v := γ(0) := d dt γ(t) t=0, γ(0) = x tangent space of M at x M defined by T x M := {v = γ(0) : γ(0) = x} R n is a d-dimensional subspace of R n geodesic curves defined by the relation γ(t) γ(t) = 0, t R for any point x M and vector v T x M there exists a unique geodesic γ x,v : R M with γ x,v (0) = x, γ x,v (0) = v
6 Optimization on Riemannian Manifolds Example of a geodesic on M = S 2
7 Optimization on Riemannian Manifolds Example of a geodesic on M = S 2
8 Optimization on Riemannian Manifolds Example of a geodesic on M = S 2
9 Optimization on Riemannian Manifolds Gradient and Hessian smooth function f : M R gradient of f at x is the unique tangent vector M f (x) T x M satisfying M f (x) v = d dt f γ x,v(t), t=0 v T x M Hessian of f at x is the unique symmetric bilinear form H M f (x) : T x M T x M R satisfying H M f (x)(v, v) = d2 dt 2 f γ x,v(t), t=0 v T x M
10 Optimization on Riemannian Manifolds Representations of the Gradient and Hessian smooth function f := f M induced by function f : R n R gradient and Hessian of f at x := (x 1,..., x n ) R n ( ) n ( f f (x) := (x) R n 2 f, H f (x) := (x) x i x i x j i=1 gradient of f at x M M f (x) = P TxM f (x) T x M ) n i,j=1 where P TxM : R n T x M is the ortho. projection onto T x M matrix representation of the Hessian of f at x M H M f (x) = H f (x) + N f (x) R n n where N f (x) R n n is a symmetric matrix satisfying R n n v N f (x)v = f (x) γ x,v (0), v T x M
11 Optimization on Riemannian Manifolds Remarks optimization methods on Riemannian manifolds M rely on closed form expressions of geodesics γ x,v, x M, v T x M for M = S d := {x R d+1 : x 2 = 1} the geodesics are v γ x,v (t) := cos(t v 2 )x+sin(t v 2 ), x M, v T x M, v 2 closed form expressions are known for specific manifolds, e.g., matrix Lie groups M R n n = R n2 have the geodesics γ X,V (t) := X exp X 1 Vt, X M, V T X M, where the matrix exponential is defined by exp A := k=0 A k k!, A Rn n
12 Optimization on Riemannian Manifolds nonlinear CG method 1 : local minimization of f : M R n R Input: initial guess x 0 M, accuracy ε > 0 1: g 0 := M f (x 0 ) T x0 M, d 0 = g 0, k = 0 2: while ε < g k 2 do 3: determine α k > 0 by a line search in direction d k : M f (γ xk,d k (α k )) γ xk,d k (α k ) = 0 4: x k+1 := γ xk,d k (α k ) M 5: dk := γ xk,d k (α k ) T xk+1 M 6: g k+1 := M f (x k+1 ) T xk+1 M H M f (x k+1 )(g k+1, d k ) 7: β k :=, k + 1 mod d 0 H M f (x k+1 )( d k, d k ) 0 k + 1 mod d 0 8: d k+1 := g k+1 + β k dk 9: k := k : end while Output: approximate stationary point x k 1 S. T. Smith (1994)
13 Optimization on Riemannian Manifolds Example of a CG-iteration on M = S 2
14 Optimization on Riemannian Manifolds Example of a CG-iteration on M = S 2
15 Optimization on Riemannian Manifolds Example of a CG-iteration on M = S 2
16 Optimization on Riemannian Manifolds Example of a CG-iteration on M = S 2
17 Efficient Evaluation Approaches Minimization of the energy E(P, w), P := (p 1,..., p M ) X M, w := (w 1,..., w M ) Ω needs optimization on the product manifold M := X M Ω R M(n+1) x := (x 1,..., x M, w) M, x j X, j = 1,..., M v := (v 1,..., v M, v 0 ) T x M, v j T xj X, v 0 T w Ω geodesics are given by γ x,v (t) = (γ x1,v 1 (t),..., γ xm,v M (t), γ w,v0 (t)) M, t R gradient M E(P, w) and Hessian H M E(P, w) can be computed by the geodesic γ x,v
18 Efficient Evaluation Approaches Naive Evaluation fixed compact manifold X R n energy E(P, w) needs O(M 2 ) function evaluations of K : X X R, h : X R, cf. E(P, w) = w Kw 2h w + C, K(p 1, p 1 )... K(p 1, p M ) h(p K := )..., h :=... K(p M, p 1 )... K(p M, p M ) h(p M ) the gradient M E(P, w) and the Hessian H M E(P, w) can be computed using O(M 2 ) function evaluation of the derivatives of K and h
19 Efficient Evaluation Approaches Nonequispaced Fast Fourier Trasnsforms 1 X = S 2, SO(3), or T d, d N polynomial functions of the form K(x, y) = λ l ψ l (x)ψ l (y), h(y) = ˆν l ψ X,l (y) x, y X, N d l=0 for certain polynomial bases ψ l L 2 (X ), l N, on X evaluation of the energy E(P, w), its gradient M E(P, w), and the matrix-vector multiplication H M E(P, w)v, v T x M, need 2 S 2 : O(M log 2 M), (M N 2 ) N d l=0 SO(3) : O(M log 2 M), (M N 3 ) T d : O(M log M), (M N d ) arithmetic operations 1 NFFT Library (J. Keiner, S. Kunis, D. Potts): http// potts/nfft 2 MG, D. Potts (2011); MG, D. Potts, G. Steidl (2011)
20 Efficient Evaluation Approaches Local Kernels fixed compact d-dimensional manifold X R n kernel satisfies for R > 0 K R (x, y) = 0, x y 2 R, x, y X for point set P := {p 1,..., p M } X M define N R (x) := {p P : p x 2 R} points P X M are well distributed if, for fixed c > 0, M R d N R (p i ) < c, i = 1,..., M for well distributed points, evaluation of the energy E(P, w), M E(P, w), and H M E(P, w)v, v T x M : the determination of the sets N R (p i ), i = 1,..., M, in O(M log M) operations O(M) function evaluations of K, h, and their derivatives
21 Discrepancies and Potential Energies compact set X R n with standard Borel σ-algebra σ(x ) Borel measure ν : σ(x ) R with ν(x ) = 1 normalized point measure at P X M { ν P := 1 M 1, x Ω, δ pi, δ x (Ω) := Ω X, x X M 0, else, i=1 discrepancy system 1 D σ(x ) of measurable test sets L 2 -discrepancy between ν P and ν (formally) ( ( ) 1 2dµD 2 D D (P, ν) := ν P (B) ν(b)) (B) D for certain choices of K : X X R, h : X R, C one calculates E(P) := E (P, 1M ) e = DD(P, 2 ν), P X M 1 M. Drmota, R. F. Tichy Sequences, Discrepancies and Applications
22 Discrepancies and Potential Energies L 2 -Discrepancies over Halfspaces discrepancy system D consists of halfspaces h + (n, d) := {x R n : x n d}, (n, d) D := S n 1 R such that the correpsonding halfplanes have non-empty intersection with X, i.e., D := H X := {h + (n, d) X : (n, d) D X }, D X := {(n, d) D : h + (n, d) h ( n, d) X } canonical measure on D induced by µ S n 1 µ R the L 2 -discrepancy over halfspaces reads as 1 D 2 H X (P, ν) = 1 M 2 M K H (p i, p j ) 2 M i,j=1 with the Euclidean distance kernel 1 R. Alexander (1975) M i=1 K H (x, y) := C X c n x y 2, X K H (p i, x)d ν (x)+c K,ν x, y X
23 Discrepancies and Potential Energies L 2 -Discrepancies over Euclidean Balls for fixed radius R > 0, the discrepancy system D consists of Euclidean balls B n (c, r) := {x R n : c x 2 r}, c D := R n, r > 0, which have non-empty intersection with X, i.e., D := B R := {B n (c, R) X : c D} canonical measure on D induced by µ R n the L 2 -discrepancy over Euclidean balls reads as DB 2 R (P, ν) = 1 M M 2 K BR (p i, p j ) 2 M K BR (p M i, x)d ν (x)+c K,ν with kernel i,j=1 i=1 K BR (x, y) := µ R n( Bn (x, R) B n (y, R) ), x, y X X
24 Discrepancies and Potential Energies Remarks the kernels K BR, R > 0, are local kernels with K BR (x, y) = 0, x y 2 2R, x, y X, positive definite, and radial, i.e, K BR (x, y) = k n R ( x y 2), x, y X R n, where kr n : [0, ) R is n 2-times differentiable for m n the kernels defined by K R (x, y) := k m R ( x y 2), x, y X R n, are positive definite and correspond to some L 2 -discrepancy the kernels K R, R > 0, are used for computing low-discrepancy points of D H (P, ν)
25 Numerical Examples Low-Discrepancy Points on the Sphere S d X = S d R d+1 with measure ν := µ S d ω d, ω d := µ S d (S d ) E(P) = DH 2 (P, ν) for S d K(x, y) := ω d 1 d x y 2, h(x) := C := 2ω 2dω 2 d 1 dω 2d 1 ω d optimal points P M X M of the L 2 -Discrepancy over halfspaces satisfy 1 c d M d+1 2d D HS d (P M, ν) C dm d+1 2d, M N, cd, C d > 0 low-discrepancy points P M X M have optimal order M d+1 2d random points with respect to ν have expected order M W. M. Schmidt (1969); J. Beck (1984)
26 Numerical Examples Computation of Low-Discrepancy Points Fourier approximation for X = S 2 : for given M N minimize the energy E(P) with h C = 0 for polynomial kernels K N (x, y) := N n n=1 k= n 16π 2 Y k n (x)y k n (y) (2n 1)(2n + 1)(2n + 3), where N := cm 1 2 for some fixed constant c 1 Heuristic with local kernels for X = S d, d N: for given M N minimize the energy E(P) with h C = 0 for local kernels K R (x, y) := k m R ( x y 2), m := max{d + 1, 4}, where R := cm 1 d for some fixed constant c ( ) d ω 1 d d ω d 1
27 Numerical Examples Low-Discrepancy Points by Polynomial Kernels K N on S 2 L 2 -discrepancy DH S (P, ν) Number of points M Average order: 4 3 π M 1 2 Random points Optimal order: c M 3 4, c := Local min.: N = 2M J. S. Brauchart, D. P. Hardin, E. B. Saff (2012)
28 Numerical Examples Low-Discrepancy Points by Local Kernels K R on S d, d = 2, 3, 4, (P, ν) L 2 -discrepancy DH S d Number of points M S 5 : 1.78 M 3 5 Local min.: R = 1.6/M 1 5 S 4 : 1.79 M 5 8 Local min.: R = 1.7/M 1 4 S 3 : 1.74 M 2 3 Local min.: R = 1.8/M 1 3 S 2 : 1.58 M 3 4 Local min.: R = 2.4/M 1 2
29 Numerical Examples Different Geometric Characteristics of Local Minimizers on S 2 polynomial kernels K N for M = 400 points P N = 2M 1 2 N = 3M 1 2 D (P, ν) = D HS 2 H S 2 (P, ν) =
30 Numerical Examples Different Geometric Characteristics of Local Minimizers on S 2 local kernels K R for M = 400 points P R = 2.9/M 1 2 R = 2.4/M 1 2 D (P, ν) = D HS 2 H S 2 (P, ν) =
31 Conclusion general framework for optimizing potential energies on Riemannian manifolds efficient optimization by the nonlinear CG method efficient evaluation of the energy, its gradient, and Hessian times vector multiplication: for polynomial functions on S 2, SO(3), T d, d N for local kernels on S d, d N computation of low-discrepancy points on the sphere S d, d = 2, 3, 4, 5 publications: grman/
32 Conclusion general framework for optimizing potential energies on Riemannian manifolds efficient optimization by the nonlinear CG method efficient evaluation of the energy, its gradient, and Hessian times vector multiplication: for polynomial functions on S 2, SO(3), T d, d N for local kernels on S d, d N computation of low-discrepancy points on the sphere S d, d = 2, 3, 4, 5 publications: grman/ Thank you for your attention!
Efficient Algorithms for the Computation of Optimal Quadrature Points on Riemannian Manifolds
Efficient Algorithms for the Computation of Optimal Quadrature Points on Riemannian Manifolds Inaugural Dissertation zur Erlangung des Akademischen Grades Doktor eines Wissenschaftszweiges doctor rerum
More informationOn the computation of spherical designs by a new optimization approach based on fast spherical Fourier transforms
On the computation of spherical designs by a new optimization approach based on fast spherical Fourier transforms anuel Gräf Daniel Potts In this paper we consider the problem of finding numerical spherical
More informationGradient Estimates and Sobolev Inequality
Gradient Estimates and Sobolev Inequality Jiaping Wang University of Minnesota ( Joint work with Linfeng Zhou) Conference on Geometric Analysis in honor of Peter Li University of California, Irvine January
More informationThe Karcher Mean of Points on SO n
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 Contents Introduction CESAME
More informationAdvances in Manifold Learning Presented by: Naku Nak l Verm r a June 10, 2008
Advances in Manifold Learning Presented by: Nakul Verma June 10, 008 Outline Motivation Manifolds Manifold Learning Random projection of manifolds for dimension reduction Introduction to random projections
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 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 informationHyperbolic Geometry on Geometric Surfaces
Mathematics Seminar, 15 September 2010 Outline Introduction Hyperbolic geometry Abstract surfaces The hemisphere model as a geometric surface The Poincaré disk model as a geometric surface Conclusion Introduction
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 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 informationOptimisation 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 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 informationLecture 2: Some basic principles of the b-calculus
Lecture 2: Some basic principles of the b-calculus Daniel Grieser (Carl von Ossietzky Universität Oldenburg) September 20, 2012 Summer School Singular Analysis Daniel Grieser (Oldenburg) Lecture 2: Some
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 informationHandlebody Decomposition of a Manifold
Handlebody Decomposition of a Manifold Mahuya Datta Statistics and Mathematics Unit Indian Statistical Institute, Kolkata mahuya@isical.ac.in January 12, 2012 contents Introduction What is a handlebody
More informationFourier extension and sampling on the sphere
Fourier extension and sampling on the sphere Daniel Potts Technische Universität Chemnitz Faculty of Mathematics 97 Chemnitz, Germany Email: potts@mathematiktu-chemnitzde iel Van Buggenhout KU Leuven Faculty
More informationQUADRATURE ERRORS, DISCREPANCIES AND THEIR RELATIONS TO HALFTONING ON THE TORUS AND THE SPHERE
QUADRATURE ERRORS, DISCREPANCIES AND THEIR RELATIONS TO HALFTONING ON THE TORUS AND THE SPHERE MANUEL GRÄF, DANIEL POTTS, AND GABRIELE STEIDL Abstract. This paper deals with a special halftoning process,
More informationQUADRATURE ERRORS, DISCREPANCIES AND THEIR RELATIONS TO HALFTONING ON THE TORUS AND THE SPHERE
QUADRATURE ERRORS, DISCREPANCIES AND THEIR RELATIONS TO HALFTONING ON THE TORUS AND THE SPHERE MANUEL GRÄF, DANIEL POTTS, AND GABRIELE STEIDL Abstract. This paper deals with continuous-domain quantization,
More informationQuasi-invariant Measures on Path Space. Denis Bell University of North Florida
Quasi-invariant Measures on Path Space Denis Bell University of North Florida Transformation of measure under the flow of a vector field Let E be a vector space (or a manifold), equipped with a finite
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 informationNotation. General. Notation Description See. Sets, Functions, and Spaces. a b & a b The minimum and the maximum of a and b
Notation General Notation Description See a b & a b The minimum and the maximum of a and b a + & a f S u The non-negative part, a 0, and non-positive part, (a 0) of a R The restriction of the function
More informationHeat Kernel and Analysis on Manifolds Excerpt with Exercises. Alexander Grigor yan
Heat Kernel and Analysis on Manifolds Excerpt with Exercises Alexander Grigor yan Department of Mathematics, University of Bielefeld, 33501 Bielefeld, Germany 2000 Mathematics Subject Classification. Primary
More informationCOBORDISM AND FORMAL POWER SERIES
COBORDISM AND FORMAL POWER SERIES NEIL STRICKLAND Thom s cobordism theorem The graded ring of cobordism classes of manifolds is Z/2[x 2, x 4, x 5, x 6, x 8, x 9, x 10, x 11, x 12, x 13, x 14, x 16, x 17,...
More informationEfficient Spherical Designs with Good Geometric Properties
Efficient Spherical Designs with Good Geometric Properties Rob Womersley, R.Womersley@unsw.edu.au School of Mathematics and Statistics, University of New South Wales 45-design with N = 1059 and symmetrtic
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
More informationPICARD S THEOREM STEFAN FRIEDL
PICARD S THEOREM STEFAN FRIEDL Abstract. We give a summary for the proof of Picard s Theorem. The proof is for the most part an excerpt of [F]. 1. Introduction Definition. Let U C be an open subset. A
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 informationWARPED PRODUCTS PETER PETERSEN
WARPED PRODUCTS PETER PETERSEN. Definitions We shall define as few concepts as possible. A tangent vector always has the local coordinate expansion v dx i (v) and a function the differential df f dxi We
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 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 informationSparse and Fast Fourier Transforms in Biomedical Imaging
Sparse and Fast Fourier Transforms in Biomedical Imaging Stefan Kunis 1 1 joint work with Ines Melzer, TU Chemnitz and Helmholtz Center Munich, supported by DFG and Helmholtz association Outline 1 Trigonometric
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 informationCurvature homogeneity of type (1, 3) in pseudo-riemannian manifolds
Curvature homogeneity of type (1, 3) in pseudo-riemannian manifolds Cullen McDonald August, 013 Abstract We construct two new families of pseudo-riemannian manifolds which are curvature homegeneous of
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 1 Part 3: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu
More informationDistance to curves and surfaces in the Heisenberg group
Distance to curves and surfaces in the Heisenberg group Università di Bologna 3 giugno 2012 Survey based on... Fausto Ferrari, N.A. Metric normal and distance function in the Heisenberg group, Math.Z.
More informationSYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS. 1. Introduction
SYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS CRAIG JACKSON 1. Introduction Generally speaking, geometric quantization is a scheme for associating Hilbert spaces
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 informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationThe Kernel Trick. Robert M. Haralick. Computer Science, Graduate Center City University of New York
The Kernel Trick Robert M. Haralick Computer Science, Graduate Center City University of New York Outline SVM Classification < (x 1, c 1 ),..., (x Z, c Z ) > is the training data c 1,..., c Z { 1, 1} specifies
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 informationMath 215B: Solutions 3
Math 215B: Solutions 3 (1) For this problem you may assume the classification of smooth one-dimensional manifolds: Any compact smooth one-dimensional manifold is diffeomorphic to a finite disjoint union
More informationWarped Products. by Peter Petersen. We shall de ne as few concepts as possible. A tangent vector always has the local coordinate expansion
Warped Products by Peter Petersen De nitions We shall de ne as few concepts as possible. A tangent vector always has the local coordinate expansion a function the di erential v = dx i (v) df = f dxi We
More information1 First and second variational formulas for area
1 First and second variational formulas for area In this chapter, we will derive the first and second variational formulas for the area of a submanifold. This will be useful in our later discussion on
More informationGaussian Fields and Percolation
Gaussian Fields and Percolation Dmitry Beliaev Mathematical Institute University of Oxford RANDOM WAVES IN OXFORD 18 June 2018 Berry s conjecture In 1977 M. Berry conjectured that high energy eigenfunctions
More informationCHAPTER 3. Gauss map. In this chapter we will study the Gauss map of surfaces in R 3.
CHAPTER 3 Gauss map In this chapter we will study the Gauss map of surfaces in R 3. 3.1. Surfaces in R 3 Let S R 3 be a submanifold of dimension 2. Let {U i, ϕ i } be a DS on S. For any p U i we have a
More informationOn local normal forms of completely integrable systems
On local normal forms of completely integrable systems ENCUENTRO DE Răzvan M. Tudoran West University of Timişoara, România Supported by CNCS-UEFISCDI project PN-II-RU-TE-2011-3-0103 GEOMETRÍA DIFERENCIAL,
More informationSome topics in sub-riemannian geometry
Some topics in sub-riemannian geometry Luca Rizzi CNRS, Institut Fourier Mathematical Colloquium Universität Bern - December 19 2016 Sub-Riemannian geometry Known under many names: Carnot-Carathéodory
More informationA local time scaling exponent for compact metric spaces
A local time scaling exponent for compact metric spaces John Dever School of Mathematics Georgia Institute of Technology Fractals 6 @ Cornell, June 15, 2017 Dever (GaTech) Exit time exponent Fractals 6
More informationContribution from: Springer Verlag Berlin Heidelberg 2005 ISBN
Contribution from: Mathematical Physics Studies Vol. 7 Perspectives in Analysis Essays in Honor of Lennart Carleson s 75th Birthday Michael Benedicks, Peter W. Jones, Stanislav Smirnov (Eds.) Springer
More information1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0
4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationSymmetric Spaces Toolkit
Symmetric Spaces Toolkit SFB/TR12 Langeoog, Nov. 1st 7th 2007 H. Sebert, S. Mandt Contents 1 Lie Groups and Lie Algebras 2 1.1 Matrix Lie Groups........................ 2 1.2 Lie Group Homomorphisms...................
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 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 informationControllability, Observability & Local Decompositions
ontrollability, Observability & Local Decompositions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Lie Bracket Distributions ontrollability ontrollability Distributions
More informationAlgebras of singular integral operators with kernels controlled by multiple norms
Algebras of singular integral operators with kernels controlled by multiple norms Alexander Nagel Conference in Harmonic Analysis in Honor of Michael Christ This is a report on joint work with Fulvio Ricci,
More informationEigenvalues and eigenfunctions of the Laplacian. Andrew Hassell
Eigenvalues and eigenfunctions of the Laplacian Andrew Hassell 1 2 The setting In this talk I will consider the Laplace operator,, on various geometric spaces M. Here, M will be either a bounded Euclidean
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
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 informationCalderón-Zygmund inequality on noncompact Riem. manifolds
The Calderón-Zygmund inequality on noncompact Riemannian manifolds Institut für Mathematik Humboldt-Universität zu Berlin Geometric Structures and Spectral Invariants Berlin, May 16, 2014 This talk is
More informationEuler Characteristic of Two-Dimensional Manifolds
Euler Characteristic of Two-Dimensional Manifolds M. Hafiz Khusyairi August 2008 In this work we will discuss an important notion from topology, namely Euler Characteristic and we will discuss several
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 informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationPACKING DIMENSIONS, TRANSVERSAL MAPPINGS AND GEODESIC FLOWS
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 489 500 PACKING DIMENSIONS, TRANSVERSAL MAPPINGS AND GEODESIC FLOWS Mika Leikas University of Jyväskylä, Department of Mathematics and
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
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 informationVolume comparison theorems without Jacobi fields
Volume comparison theorems without Jacobi fields Dominique Bakry Laboratoire de Statistique et Probabilités Université Paul Sabatier 118 route de Narbonne 31062 Toulouse, FRANCE Zhongmin Qian Mathematical
More informationKernel-Based Contrast Functions for Sufficient Dimension Reduction
Kernel-Based Contrast Functions for Sufficient Dimension Reduction Michael I. Jordan Departments of Statistics and EECS University of California, Berkeley Joint work with Kenji Fukumizu and Francis Bach
More informationRiemannian Stein Variational Gradient Descent for Bayesian Inference
Riemannian Stein Variational Gradient Descent for Bayesian Inference Chang Liu, Jun Zhu 1 Dept. of Comp. Sci. & Tech., TNList Lab; Center for Bio-Inspired Computing Research State Key Lab for Intell. Tech.
More informationCMS winter meeting 2008, Ottawa. The heat kernel on connected sums
CMS winter meeting 2008, Ottawa The heat kernel on connected sums Laurent Saloff-Coste (Cornell University) December 8 2008 p(t, x, y) = The heat kernel ) 1 x y 2 exp ( (4πt) n/2 4t The heat kernel is:
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 informationNOTES ON THE REGULARITY OF QUASICONFORMAL HOMEOMORPHISMS
NOTES ON THE REGULARITY OF QUASICONFORMAL HOMEOMORPHISMS CLARK BUTLER. Introduction The purpose of these notes is to give a self-contained proof of the following theorem, Theorem.. Let f : S n S n be a
More informationHI CAMBRIDGE n S P UNIVERSITY PRESS
Infinite-Dimensional Dynamical Systems An Introduction to Dissipative Parabolic PDEs and the Theory of Global Attractors JAMES C. ROBINSON University of Warwick HI CAMBRIDGE n S P UNIVERSITY PRESS Preface
More informationHomework for Math , Spring 2012
Homework for Math 6170 1, Spring 2012 Andres Treibergs, Instructor April 24, 2012 Our main text this semester is Isaac Chavel, Riemannian Geometry: A Modern Introduction, 2nd. ed., Cambridge, 2006. Please
More informationRandom measures, intersections, and applications
Random measures, intersections, and applications Ville Suomala joint work with Pablo Shmerkin University of Oulu, Finland Workshop on fractals, The Hebrew University of Jerusalem June 12th 2014 Motivation
More informationInformation geometry for bivariate distribution control
Information geometry for bivariate distribution control C.T.J.Dodson + Hong Wang Mathematics + Control Systems Centre, University of Manchester Institute of Science and Technology Optimal control of stochastic
More informationPoincaré Duality Angles on Riemannian Manifolds with Boundary
Poincaré Duality Angles on Riemannian Manifolds with Boundary Clayton Shonkwiler Department of Mathematics University of Pennsylvania June 5, 2009 Realizing cohomology groups as spaces of differential
More informationGenerated sets as sampling schemes for hyperbolic cross trigonometric polynomials
Generated sets as sampling schemes for hyperbolic cross trigonometric polynomials Lutz Kämmerer Chemnitz University of Technology joint work with S. Kunis and D. Potts supported by DFG-SPP 1324 1 / 16
More informationGeometry/Topology 868
Geometry/Topology 868 Homework Samuel Otten CLAIM: A bijective map between manifolds is not necessarily a diffeomorphism. Consider the example f : R R, x x 3. We know that R is a manifold by example from
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 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 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 informationKernels A Machine Learning Overview
Kernels A Machine Learning Overview S.V.N. Vishy Vishwanathan vishy@axiom.anu.edu.au National ICT of Australia and Australian National University Thanks to Alex Smola, Stéphane Canu, Mike Jordan and Peter
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
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 informationHomework in Topology, Spring 2009.
Homework in Topology, Spring 2009. Björn Gustafsson April 29, 2009 1 Generalities To pass the course you should hand in correct and well-written solutions of approximately 10-15 of the problems. For higher
More informationPositive Definite Functions on Spheres
Positive Definite Functions on Spheres Tilmann Gneiting Institute for Applied Mathematics Heidelberg University, Germany t.gneiting@uni-heidelberg.de www.math.uni-heidelberg.de/spatial/tilmann/ International
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
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 informationMicrolocal Methods in X-ray Tomography
Microlocal Methods in X-ray Tomography Plamen Stefanov Purdue University Lecture I: Euclidean X-ray tomography Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Methods
More informationL 2 Geometry of the Symplectomorphism Group
University of Notre Dame Workshop on Innite Dimensional Geometry, Vienna 2015 Outline 1 The Exponential Map on D s ω(m) 2 Existence of Multiplicity of Outline 1 The Exponential Map on D s ω(m) 2 Existence
More informationIntroduction to Practical FFT and NFFT
Introduction to Practical FFT and NFFT Michael Pippig and Daniel Potts Faculty of Mathematics Chemnitz University of Technology 07.09.2010 supported by BMBF grant 01IH08001B Table of Contents 1 Serial
More informationBredon, Introduction to compact transformation groups, Academic Press
1 Introduction Outline Section 3: Topology of 2-orbifolds: Compact group actions Compact group actions Orbit spaces. Tubes and slices. Path-lifting, covering homotopy Locally smooth actions Smooth actions
More informationIntroduction to Practical FFT and NFFT
Introduction to Practical FFT and NFFT Michael Pippig and Daniel Potts Department of Mathematics Chemnitz University of Technology September 14, 211 supported by BMBF grant 1IH81B Table of Contents 1 Serial
More informationLocal radial basis function approximation on the sphere
Local radial basis function approximation on the sphere Kerstin Hesse and Q. T. Le Gia School of Mathematics and Statistics, The University of New South Wales, Sydney NSW 05, Australia April 30, 007 Abstract
More informationWeek 6: Differential geometry I
Week 6: Differential geometry I Tensor algebra Covariant and contravariant tensors Consider two n dimensional coordinate systems x and x and assume that we can express the x i as functions of the x i,
More informationQualifying Exams I, 2014 Spring
Qualifying Exams I, 2014 Spring 1. (Algebra) Let k = F q be a finite field with q elements. Count the number of monic irreducible polynomials of degree 12 over k. 2. (Algebraic Geometry) (a) Show that
More informationHadamard s Theorem. Rich Schwartz. September 10, The purpose of these notes is to prove the following theorem.
Hadamard s Theorem Rich Schwartz September 10, 013 1 The Result and Proof Outline The purpose of these notes is to prove the following theorem. Theorem 1.1 (Hadamard) Let M 1 and M be simply connected,
More informationHilbert Space Methods for Reduced-Rank Gaussian Process Regression
Hilbert Space Methods for Reduced-Rank Gaussian Process Regression Arno Solin and Simo Särkkä Aalto University, Finland Workshop on Gaussian Process Approximation Copenhagen, Denmark, May 2015 Solin &
More information