Towards a Faster Spherical Harmonic transform
|
|
- Miranda McKinney
- 5 years ago
- Views:
Transcription
1 Towards a Faster Spherical Harmonic transform Nathanaël Schaeffer LGIT / CNRS / Université Joseph Fourier IUGG CMG, Pisa, 7 May 2010 N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
2 1 Introduction 2 Taking advantage of the DCT 3 The SHTns library N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
3 Spherical Harmonic Transform (SHT) Definition Spherical Harmonic of degree l and order m, defined on the sphere : Y m l (θ, φ) Eigenfunction of the Laplace operator on the sphere : Y m l = l(l + 1)Y m l N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
4 Spherical Harmonic Transform (SHT) Definition Spherical Harmonic of degree l and order m, defined on the sphere : Y m l (θ, φ) Eigenfunction of the Laplace operator on the sphere : Y m l They form an orthonormal basis f (θ, φ) = l,m = l(l + 1)Y m l Q m l Y m l (θ, φ) Q m l = f (θ, φ) Y m l (θ, φ) sin(θ)dθdφ N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
5 Application for numerical simulations Advantages Spectral convergence Exact derivatives Boundary conditions for a magnetic field are straightforward. N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
6 Application for numerical simulations Advantages Spectral convergence Exact derivatives Boundary conditions for a magnetic field are straightforward. Drawback As of today, Gauss-Legendre algorithm is the best choice (complexity N θ (N lm + N φ log N φ )), which means that for high resolution you ll spend most of your time performing SHT! N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
7 Application for numerical simulations Advantages Spectral convergence Exact derivatives Boundary conditions for a magnetic field are straightforward. Drawback As of today, Gauss-Legendre algorithm is the best choice (complexity N θ (N lm + N φ log N φ )), which means that for high resolution you ll spend most of your time performing SHT! A fast algorithm would not harm! N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
8 Fast algorithms There are a few fast algorithm : Healy (1996,2003) : FPM (source available) Mohlenkamp (1997,1999) : Wavelet approach (approximate) Reiji Suda & Masayasu Takami (2001) : FMM (approximate) Pots (2003, 2006) : DCT, FPM (source available) Tygert (2008) : FMM (approximate) N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
9 Fast algorithms There are a few fast algorithm : Healy (1996,2003) : FPM (source available) Mohlenkamp (1997,1999) : Wavelet approach (approximate) Reiji Suda & Masayasu Takami (2001) : FMM (approximate) Pots (2003, 2006) : DCT, FPM (source available) Tygert (2008) : FMM (approximate) They all claim to be faster for N 512. But... N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
10 Fast algorithms There are a few fast algorithm : Healy (1996,2003) : FPM (source available) Mohlenkamp (1997,1999) : Wavelet approach (approximate) Reiji Suda & Masayasu Takami (2001) : FMM (approximate) Pots (2003, 2006) : DCT, FPM (source available) Tygert (2008) : FMM (approximate) They all claim to be faster for N 512. But... Timing for 1 scalar SHT l max Pots Healy Gauss Healy/Gauss s 5.5 ms 0.69 ms s 27 ms 6.5 ms 4.2 N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
11 Fast algorithms There are a few fast algorithm : Healy (1996,2003) : FPM (source available) Mohlenkamp (1997,1999) : Wavelet approach (approximate) Reiji Suda & Masayasu Takami (2001) : FMM (approximate) Pots (2003, 2006) : DCT, FPM (source available) Tygert (2008) : FMM (approximate) They all claim to be faster for N 512. But... Timing for 1 scalar SHT l max Pots Healy Gauss Healy/Gauss s 5.5 ms 0.69 ms s 27 ms 6.5 ms s 143 ms 49 ms 2.9 N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
12 Fast algorithms There are a few fast algorithm : Healy (1996,2003) : FPM (source available) Mohlenkamp (1997,1999) : Wavelet approach (approximate) Reiji Suda & Masayasu Takami (2001) : FMM (approximate) Pots (2003, 2006) : DCT, FPM (source available) Tygert (2008) : FMM (approximate) They all claim to be faster for N 512. But... Timing for 1 scalar SHT l max Pots Healy Gauss Healy/Gauss s 5.5 ms 0.69 ms s 27 ms 6.5 ms s 143 ms 49 ms s 860 ms 297 ms 2.9 N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
13 Gauss beats them all! Gauss-Legendre Quadrature P k (x)dx N P k (x j )w j For polynomials of degree l, the quadrature is exact if 2N > k. j=1 N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
14 Gauss beats them all! Gauss-Legendre Quadrature P k (x)dx N P k (x j )w j For polynomials of degree l, the quadrature is exact if 2N > k. j=1 This means that our integral f m (θ) l (θ) sin θdθ }{{}}{{} l l Y m is exactly computed by the quadrature if we use 2N θ > 2l max Gauss points. N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
15 Gauss beats them all! Gauss-Legendre Quadrature P k (x)dx N P k (x j )w j For polynomials of degree l, the quadrature is exact if 2N > k. j=1 This means that our integral f m (θ) l (θ) sin θdθ }{{}}{{} l l Y m is exactly computed by the quadrature if we use 2N θ > 2l max Gauss points. The Gauss-Legendre quadrature rule keeps the number of colocation points low N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
16 Gauss on steroids θ Parity separation : all Y m l have defined parity (x2 speedup) N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
17 Gauss on steroids θ Parity separation : all Yl m have defined parity (x2 speedup) Polar optimization : for large m, Yl m (θ) is almost zero near the poles ( 15% speed-up) N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
18 Go to DCT space For m even with Y m l (θ, φ) = P m l (cos θ)e imφ P m l (cos θ) = l k=0 ak lm cos(kθ) N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
19 Go to DCT space For m even with Y m l (θ, φ) = P m l (cos θ)e imφ P m l (cos θ) = l k=0 ak lm cos(kθ) + Less terms in the sum (independant of N θ ) : Efficient for small m only. l max fk m = Ql m ak lm l=k + Fast DCT : f m (θ i ) = DCT [f m k ] + Fast DCT requires N θ > l max equally spaced nodes. We can no longer use the Gauss quadrature. N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
20 Fejér (or Clenshaw-Curtis) quadrature Quadrature in DCT space l max Ql m = f (θ)yl m (θ) = f k k=0 l max T k Yl m = k=0 f k W lm k With T k (cos θ) = cos(kθ) the Tchebychev polynomial. This is te Fejér quadrature with weight Y m l, which is exact for N θ > l max. As the Wk lm matrix is dense, there is no advantage in computing the quadrature in the cosine domain, so we go back to real domain : Q m l = k f m (θ k ) IDCT [W lm k ] which is still exact for N θ > l max as the Gauss quadrature We can of course precompute the W lm k coefficients. N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
21 For non-linear simulations? Aliasing For an equation with quadratic non-linear terms : f m (θ) Yl m (θ) sin θdθ }{{}}{{} 2l l is exactly computed by a Gauss quadrature if we use 2N θ > 3l max Gauss points. a Fejér quadrature if we use N θ > 2l max regular grid points. Aliasing for non-linear problems penalizes the Fejér quadrature. N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
22 a library for Spherical Harmonic Transform Aimed at Numerical Simulations, it focuses on speed (but with high double precision accuracy). both scalar and vector transforms flexible truncation latitude-major or longitude-major spatial data ordering various conventions (normalizations, Condon-Shortley phase) can be called from fortran or c/c++ programs Highly optimized Gauss-Legendre algorithm DCT algorithm for regular nodes (generalized Fejér quadrature) support for SSE2 & SSE3 vectorization with gcc N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
23 Timing comparison The inverse SHT Ql m f (θ, φ) is much faster with the DCT : l max m max perf. gain /Gauss x4 x1.9 x1.15 N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
24 Timing comparison The inverse SHT Ql m f (θ, φ) is much faster with the DCT : l max m max perf. gain /Gauss x4 x1.9 x1.15 However, the analysis step involving quadrature is just the same. N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
25 Timing comparison The inverse SHT Ql m f (θ, φ) is much faster with the DCT : l max m max perf. gain /Gauss x4 x1.9 x1.15 However, the analysis step involving quadrature is just the same. Hence the real-life results are not as good : l max = 1021, m max = 10, NL order=2, 2x isht, 1x SHT DCT is % faster N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
26 Conclusions & future developments Conclusions Asymptotically fast algorithms do not appear fast yet! A DCT based SHT is more efficient for anisotropic truncations. A highly optimized Gauss SHT is still the overall winner! Anyone producing a fast algorithm should compare with it. Freely available SHTns library that picks the fastest for you N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
27 Conclusions & future developments Conclusions Asymptotically fast algorithms do not appear fast yet! A DCT based SHT is more efficient for anisotropic truncations. A highly optimized Gauss SHT is still the overall winner! Anyone producing a fast algorithm should compare with it. Freely available SHTns library that picks the fastest for you Future work SHTns will use the tremendous computing power of GPUs. Maybe parallelization. Perhaps somebody needs to writes a good implementation of a fast algorithm... N. Schaeffer (LGIT/CNRS/UJF) Towards a faster SHT IUGG CMG, Pisa, 7 May / 13
Angular Momentum. Classically the orbital angular momentum with respect to a fixed origin is. L = r p. = yp z. L x. zp y L y. = zp x. xpz L z.
Angular momentum is an important concept in quantum theory, necessary for analyzing motion in 3D as well as intrinsic properties such as spin Classically the orbital angular momentum with respect to a
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationThe spectral transform method
The spectral transform method by Nils Wedi European Centre for Medium-Range Weather Forecasts wedi@ecmwf.int Advanced Numerical Methods for Earth-System Modelling Slide 1 Advanced Numerical Methods for
More informationFast Algorithms for Spherical Harmonic Expansions
Fast Algorithms for Spherical Harmonic Expansions Vladimir Rokhlin and Mark Tygert Research Report YALEU/DCS/RR-1309 December 17, 004 Abstract An algorithm is introduced for the rapid evaluation at appropriately
More informationA FAST SPHERICAL HARMONICS TRANSFORM ALGORITHM
MATHEMATICS OF COMPUTATION Volume 71, Number 238, Pages 703 715 S 0025-5718(01)01386-2 Article electronically published on November 28, 2001 A FAST SPHERICAL HARMONICS TRANSFORM ALGORITHM REIJI SUDA AND
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationTurbulence in geodynamo simulations
Turbulence in geodynamo simulations Nathanaël Schaeffer, A. Fournier, D. Jault, J. Aubert, H-C. Nataf,... ISTerre / CNRS / Université Grenoble Alpes Journée des utilisateurs CIMENT, Grenoble, 23 June 2016
More informationFast algorithms based on asymptotic expansions of special functions
Fast algorithms based on asymptotic expansions of special functions Alex Townsend MIT Cornell University, 23rd January 2015, CAM Colloquium Introduction Two trends in the computing era Tabulations Software
More informationEigenfunctions on the surface of a sphere. In spherical coordinates, the Laplacian is. u = u rr + 2 r u r + 1 r 2. sin 2 (θ) + 1
Eigenfunctions on the surface of a sphere In spherical coordinates, the Laplacian is u = u rr + 2 r u r + 1 r 2 [ uφφ sin 2 (θ) + 1 sin θ (sin θ u θ) θ ]. Eigenfunctions on the surface of a sphere In spherical
More informationA Fast Spherical Filter with Uniform Resolution
JOURNAL OF COMPUTATIONAL PHYSICS 136, 580 584 (1997) ARTICLE NO. CP975782 A Fast Spherical Filter with Uniform Resolution Rüdiger Jakob-Chien*, and Bradley K. Alpert *Department of Computer Science & Engineering,
More informationComparison of Clenshaw-Curtis and Gauss Quadrature
WDS'11 Proceedings of Contributed Papers, Part I, 67 71, 211. ISB 978-8-7378-184-2 MATFYZPRESS Comparison of Clenshaw-Curtis and Gauss Quadrature M. ovelinková Charles University, Faculty of Mathematics
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationA fast Chebyshev Legendre transform using an asymptotic formula
A fast Chebyshev Legendre transform using an asymptotic formula Alex Townsend MIT (Joint work with ick Hale) SIAM OPSFA, 5th June 2015 Introduction What is the Chebyshev Legendre transform? Suppose we
More informationA Fast Algorithm for Spherical Grid Rotations and its Application to Singular Quadrature
A Fast Algorithm for Spherical Grid Rotations and its Application to Singular Quadrature Zydrunas Gimbutas Shravan Veerapaneni September 11, 2013 Abstract We present a fast and accurate algorithm for evaluating
More informationSpectral transforms. Contents. ARPEGE-Climat Version 5.1. September Introduction 2
Spectral transforms ARPEGE-Climat Version 5.1 September 2008 Contents 1 Introduction 2 2 Spectral representation 2 2.1 Spherical harmonics....................... 2 2.2 Collocation grid.........................
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 7, No. 6, pp. 1903 198 c 006 Society for Industrial and Applied Mathematics FAST ALGORITHMS FOR SPHERICAL HARMONIC EXPANSIONS VLADIMIR ROKHLIN AND MARK TYGERT Abstract. An algorithm
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More informationThe Spherical Harmonics
Physics 6C Fall 202 The Spherical Harmonics. Solution to Laplace s equation in spherical coordinates In spherical coordinates, the Laplacian is given by 2 = ( r 2 ) ( + r 2 r r r 2 sin 2 sinθ ) + θ θ θ
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More information( )e ikx. ( ) = ˆq k
! Revised December 15, 2015 6:57 PM! 1 Chapter 13: Spectral Methods Copyright 2015, David A. Randall Introduction Assume that q( x,t) is real and integrable. If the domain is periodic, with period L, we
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More information(ii) Observational Geomagnetism. Lecture 5: Spherical harmonic field models
(ii) Observational Geomagnetism Lecture 5: Spherical harmonic field models Lecture 5: Spherical harmonic field models 5.1 Introduction 5.2 How to represent functions on a spherical surface 5.3 Spherical
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationFFTs for the 2-Sphere Improvements and Variations
FFTs for the -Sphere Improvements and Variations D.M. Healy Jr., D. Rockmore and P. J. Kostelec Department of Mathematics Dartmouth College Hanover, NH 03755 Sean S. B. Moore BBN Systems and Technologies
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationSummary: angular momentum derivation
Summary: angular momentum derivation L = r p L x = yp z zp y, etc. [x, p y ] = 0, etc. (-) (-) (-3) Angular momentum commutation relations [L x, L y ] = i hl z (-4) [L i, L j ] = i hɛ ijk L k (-5) Levi-Civita
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationThe Rigid Rotor Problem: A quantum par7cle confined to a sphere Spherical Harmonics. Reading: McIntyre 7.6
The Rigid Rotor Problem: A quantum par7cle confined to a sphere Spherical Harmonics Reading: McIntyre 7.6 Legendre s equa7on (m = 0) The series must be finite! If the series is not finite, the polynomial
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationSlepian functions on the sphere, generalized Gaussian quadrature rule
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems (4 877 89 INVERSE PROBLEMS PII: S66-56(4769-7 Slepian functions on the sphere, generalized Gaussian quadrature rule L Miranian Department of Mathematics,
More informationOrdinary Least Squares and its applications
Ordinary Least Squares and its applications Dr. Mauro Zucchelli University Of Verona December 5, 2016 Dr. Mauro Zucchelli Ordinary Least Squares and its applications December 5, 2016 1 / 48 Contents 1
More informationHigh precision pixel-based simulations of CMB lensing
Laboratoire APC AstroParticule et Cosmologie University of Paris V Denis Diderot 22nd April 2011 Berkeley CMB lensing workshop Outline 1 Pixel Based methods 2 LenS 2 HAT 3 Application 4 Conclusions How
More informationQuantum Mechanics in 3-Dimensions
Quantum Mechanics in 3-Dimensions Pavithran S Iyer, 2nd yr BSc Physics, Chennai Mathematical Institute Email: pavithra@cmi.ac.in August 28 th, 2009 1 Schrodinger equation in Spherical Coordinates 1.1 Transforming
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
More informationH atom solution. 1 Introduction 2. 2 Coordinate system 2. 3 Variable separation 4
H atom solution Contents 1 Introduction 2 2 Coordinate system 2 3 Variable separation 4 4 Wavefunction solutions 6 4.1 Solution for Φ........................... 6 4.2 Solution for Θ...........................
More informationFun With Carbon Monoxide. p. 1/2
Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationPhysics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics
Physics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics We now consider the spatial degrees of freedom of a particle moving in 3-dimensional space, which of course is an important
More informationAPPLICATION OF MVDR BEAMFORMING TO SPHERICAL ARRAYS
AMBISONICS SYMPOSIUM 29 June 2-27, Graz APPLICATION OF MVDR BEAMFORMING TO SPHERICAL ARRAYS Anton Schlesinger 1, Marinus M. Boone 2 1 University of Technology Delft, The Netherlands (a.schlesinger@tudelft.nl)
More informationarxiv: v1 [math.ca] 21 Mar 2014
SLEPIAN SPATIAL-SPECTRAL CONCENTRATION ON THE BALL ZUBAIR KHALID, RODNEY A. KENNEDY, AND JASON D. MCEWEN arxiv:1403.5553v1 [math.ca] 21 Mar 2014 Abstract. We formulate and solve the Slepian spatial-spectral
More informationFFTs in Graphics and Vision. Fast Alignment of Spherical Functions
FFTs in Graphics and Vision Fast Alignment of Spherical Functions Outline Math Review Fast Rotational Alignment Review Recall 1: We can represent any rotation R in terms of the triplet of Euler angles
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 information(a) Determine the general solution for φ(ρ) near ρ = 0 for arbitary values E. (b) Show that the regular solution at ρ = 0 has the series expansion
Problem 1. Curious Wave Functions The eigenfunctions of a D7 brane in a curved geometry lead to the following eigenvalue equation of the Sturm Liouville type ρ ρ 3 ρ φ n (ρ) = E n w(ρ)φ n (ρ) w(ρ) = where
More informationIV. Electronic Spectroscopy, Angular Momentum, and Magnetic Resonance
IV. Electronic Spectroscopy, Angular Momentum, and Magnetic Resonance The foundation of electronic spectroscopy is the exact solution of the time-independent Schrodinger equation for the hydrogen atom.
More informationChapter 4 Multipole model of the Earth's magnetic field
Chapter 4 Multipole model of the Earth's magnetic field 1 Previously A measurement of the geomagnetic field at any given point and time consists of a superposition of fields from different sources: Internal
More informationAtoms 2010 update -- start with single electron: H-atom
VII 33 Atoms 2010 update -- start with single electron: H-atom x z φ θ e -1 y 3-D problem - free move in x, y, z - easier if change coord. systems: Cartesian Spherical Coordinate (x, y, z) (r, θ, φ) Reason:
More informationThe Hydrogen atom. Chapter The Schrödinger Equation. 2.2 Angular momentum
Chapter 2 The Hydrogen atom In the previous chapter we gave a quick overview of the Bohr model, which is only really valid in the semiclassical limit. cf. section 1.7.) We now begin our task in earnest
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherical Bessel Functions We would like to solve the free Schrödinger equation [ h d l(l + ) r R(r) = h k R(r). () m r dr r m R(r) is the radial wave function ψ( x)
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 informationS2kit: A Lite Version of SpharmonicKit
February 23, 2004 S2kit: A Lite Version of SpharmonicKit Peter J. Kostelec Daniel N. Rockmore Department of Mathematics Dartmouth College Hanover, NH 03755 S2kit, version 1.0, is a collection of C routines
More informationAtoms 2012 update -- start with single electron: H-atom
Atoms 2012 update -- start with single electron: H-atom x z φ θ e -1 y 3-D problem - free move in x, y, z - easier if change coord. systems: Cartesian Spherical Coordinate (x, y, z) (r, θ, φ) Reason: V(r)
More informationSpherical Harmonics and Related Topics. David Randall 2 S = 0, r 2 r r S 2. S = r n Y n
! Revised April 10, 2017 1:32 PM! 1 Spherical Harmonics and Related Topics David Randall The spherical surface harmonics are convenient functions for representing the distribution of geophysical quantities
More informationTotal Angular Momentum for Hydrogen
Physics 4 Lecture 7 Total Angular Momentum for Hydrogen Lecture 7 Physics 4 Quantum Mechanics I Friday, April th, 008 We have the Hydrogen Hamiltonian for central potential φ(r), we can write: H r = p
More informationarxiv: v3 [math.na] 30 Nov 2016
Fast Fourier Transforms for Spherical Gauss-Laguerre Basis Functions Jürgen Prestin and Christian Wüler arxiv:164.514v3 [math.na] 3 Nov 216 Abstract Spherical Gauss-Laguerre (SGL) basis functions, i.e.,
More informationSeparation of Variables in Polar and Spherical Coordinates
Separation of Variables in Polar and Spherical Coordinates Polar Coordinates Suppose we are given the potential on the inside surface of an infinitely long cylindrical cavity, and we want to find the potential
More informationRecent Developments in Numerical Methods for 4d-Var
Recent Developments in Numerical Methods for 4d-Var Mike Fisher Slide 1 Recent Developments Numerical Methods 4d-Var Slide 2 Outline Non-orthogonal wavelets on the sphere: - Motivation: Covariance Modelling
More informationSIGNALS are inherently defined on the sphere in a variety
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 62, NO. 17, SEPTEMBER 1, 2014 4597 An Optimal-Dimensionality Sampling Scheme on the Sphere With Fast Spherical Harmonic Transforms Zubair Khalid, Member,IEEE,RodneyA.Kennedy,Fellow,IEEE,and
More informationElectrodynamics I Midterm - Part A - Closed Book KSU 2005/10/17 Electro Dynamic
Electrodynamics I Midterm - Part A - Closed Book KSU 5//7 Name Electro Dynamic. () Write Gauss Law in differential form. E( r) =ρ( r)/ɛ, or D = ρ, E= electricfield,ρ=volume charge density, ɛ =permittivity
More informationExploring the physics of the Earth s core with numerical simulations
Exploring the physics of the Earth s core with numerical simulations Nathanaël Schaeffer To cite this version: Nathanaël Schaeffer. Exploring the physics of the Earth s core with numerical simulations.
More informationFast Evaluation and Irregular Sampling of Band-limited Functions on the Sphere
CONSTRUCTIVE THEORY OF FUNCTIONS, Sozopol 2013 (K. Ivanov, G. Nikolov and R. Uluchev, Eds.), pp. 115-140 Prof. Marin Drinov Academic Publishing House, Sofia, 2014 Fast Evaluation and Irregular Sampling
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More information5.61 Lecture #17 Rigid Rotor I
561 Fall 2017 Lecture #17 Page 1 561 Lecture #17 Rigid Rotor I Read McQuarrie: Chapters E and 6 Rigid Rotors molecular rotation and the universal angular part of all central force problems another exactly
More informationSpherical Coordinates
Spherical Coordinates Bo E. Sernelius 4:6 SPHERICAL COORDINATES Φ = 1 Φ + 1 Φ + 1 Φ = 0 r r r θ sin r r sinθ θ θ r sin θ ϕ ( ) = ( ) ( ) ( ) Φ r, θϕ, R r P θq ϕ 1 d 1 1 0 r R dr r dr d dp d Q dr + sinθ
More informationOn an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University
On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials
More informationFinal Examination Solutions
Math. 42 Fulling) 4 December 25 Final Examination Solutions Calculators may be used for simple arithmetic operations only! Laplacian operator in polar coordinates: Some possibly useful information 2 u
More informationIndex. for Ɣ(a, z), 39. convergent asymptotic representation, 46 converging factor, 40 exponentially improved, 39
Index Abramowitz function computed by Clenshaw s method, 74 absolute error, 356 Airy function contour integral for, 166 Airy functions algorithm, 359 asymptotic estimate of, 18 asymptotic expansions, 81,
More informationSlow Growth for Gauss Legendre Sparse Grids
Slow Growth for Gauss Legendre Sparse Grids John Burkardt, Clayton Webster April 4, 2014 Abstract A sparse grid for multidimensional quadrature can be constructed from products of 1D rules. For multidimensional
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter ifor a one-dimenional problem to a specific and very important
More informationDifferentiation and Integration
Differentiation and Integration (Lectures on Numerical Analysis for Economists II) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 12, 2018 1 University of Pennsylvania 2 Boston College Motivation
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
More informationWavelet Frames on the Sphere for Sparse Representations in High Angular Resolution Diusion Imaging. Chen Weiqiang
Wavelet Frames on the Sphere for Sparse Representations in High Angular Resolution Diusion Imaging Chen Weiqiang Overview 1. Introduction to High Angular Resolution Diusion Imaging (HARDI). 2. Wavelets
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationLegendre Polynomials and Angular Momentum
University of Connecticut DigitalCommons@UConn Chemistry Education Materials Department of Chemistry August 006 Legendre Polynomials and Angular Momentum Carl W. David University of Connecticut, Carl.David@uconn.edu
More informationPHYS 502 Lecture 8: Legendre Functions. Dr. Vasileios Lempesis
PHYS 502 Lecture 8: Legendre Functions Dr. Vasileios Lempesis Introduction Legendre functions or Legendre polynomials are the solutions of Legendre s differential equation that appear when we separate
More informationSOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS
SOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS NATASHA FLYER 1. Introduction Mathematical modeling of space and climate phenomena generally requires
More informationSeismic compression. François G. Meyer Department of Electrical Engineering University of Colorado at Boulder
Seismic compression François G. Meyer Department of Electrical Engineering University of Colorado at Boulder francois.meyer@colorado.edu http://ece-www.colorado.edu/ fmeyer IPAM, MGA 2004 Seismic Compression
More informationBOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann,
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More informationFastMMLib: a generic Fast Multipole Method library. Eric DARRIGRAND. IRMAR Université de Rennes 1
FastMMLib: a generic Fast Multipole Method library Eric DARRIGRAND joint work with Yvon LAFRANCHE IRMAR Université de Rennes 1 Outline Introduction to FMM Historically SVD principle A first example of
More informationLecture 9. Angular momentum - 2
Lecture 9 Angular momentum - 2 83 84 LECTURE 9. ANGULAR MOMENTUM - 2 9.1 Spectroscopic Notation We saw in the last section that the only allowed values of the angular momentum quantum number,, are =0,
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More informationTowards safe and effective high-order Legendre transforms with applications to FFTs for the 2-sphere
Advances in Computational Mathematics 21: 59 105, 2004. 2004 Kluwer Academic Publishers. Printed in the Netherlands. Towards safe and effective high-order Legendre transforms with applications to FFTs
More informationLecture 11 Spin, orbital, and total angular momentum Mechanics. 1 Very brief background. 2 General properties of angular momentum operators
Lecture Spin, orbital, and total angular momentum 70.00 Mechanics Very brief background MATH-GA In 9, a famous experiment conducted by Otto Stern and Walther Gerlach, involving particles subject to a nonuniform
More informationQuantum Mechanics in Multidimensions
Quantum Mechanics in Multidimensions In this chapter we discuss bound state solutions of the Schrödinger equation in more than one dimension. In this chapter we will discuss some particularly straightforward
More informationPose Detection of 3-D Objects Using Images Sampled on SO(3), Spherical Harmonics, and Wigner-D Matrices
4th IEEE Conference on Automation Science and Engineering Key Bridge Marriott, Washington DC, USA August 3-6, 008 Pose Detection of 3-D Objects Using Images Sampled on SO(3), Spherical Harmonics, and Wigner-D
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More informationAtoms 09 update-- start with single electron: H-atom
Atoms 09 update-- start with single electron: H-atom VII 33 x z φ θ e -1 y 3-D problem - free move in x, y, z - handy to change systems: Cartesian Spherical Coordinate (x, y, z) (r, θ, φ) Reason: V(r)
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationFinal Examination Solutions
Math. 42 Fulling) December 22 Final Examination Solutions Calculators may be used for simple arithmetic operations only! Laplacian operator in polar coordinates: Some possibly useful information 2 u =
More information