The automatic solution of PDEs using a global spectral method
|
|
- Aldous Chapman
- 5 years ago
- Views:
Transcription
1 The automatic solution of PDEs using a global spectral method 2014 CBMS-NSF Conference, 29th June 2014 Alex Townsend PhD student University of Oxford (with Sheehan Olver) Supervised by Nick Trefethen, University of Oxford. Supported by EPSRC grant EP/P505666/1.
2 Introduction Chebfun Piecewise smooth (2010) Endpoint singularities (2010*) Reducing regularity Blow up functions (2011) Chebfun (2004) Linear ODEs (2008) Nonlinear ODEs (2012) Chebgui (2011*) Ordinary differential equations Chebfun2 (2013) Piecewise smooth Linear PDEs Two dimensions Arbitrary domains Nonlinear PDEs Alex Oxford 1/21
3 Introduction Chebop: Spectral collocation for ODEs In 2008: Overload the MATLAB backslash command \ for operators [Driscoll, Bornemann, & Trefethen 2008]. L = chebop(@(x,u) diff(u,2)-x.*u,[-30 30]); % Airy equation L.lbc = 1; L.rbc = 0; % Set boundary conditions u = L \ 0; plot(u) % Solve and plot Alex Oxford 2/21
4 Introduction Spectral collocation basics Given values on a grid, what are the values of the derivative on that same grid?: u u 1 2 u 5 u 9 u 1 u 2 u 5 u 9 x 1 x 2 x 5 x 9 D n u 1 u n u 1 u n. =., D n = diffmat(n). For example, u (x) + cos(x)u(x) is represented as L n = D n + diag (cos(x 1 ),..., cos(x n )) R n n. Alex Oxford 3/21
5 Introduction Why do spectral methods get a bad press? 1. Dense matrices. 2. Ill-conditioned matrices. 3. When has it converged? Tricky. See, for example: [Canuto et al. 07],, [Fornberg 98], [Trefethen 00]. Alex Oxford 4/21
6 Introduction Why do spectral methods get a bad press? 1. Dense matrices. 2. Ill-conditioned matrices. 3. When has it converged? Tricky. See, for example: [Canuto et al. 07],, [Fornberg 98], [Trefethen 00]. Alex Oxford 4/21
7 Introduction Why do spectral methods get a bad press? 1. Dense matrices. 2. Ill-conditioned matrices. 3. When has it converged? Tricky Condition number See, for example: [Canuto et al. 07],, [Fornberg 98], [Trefethen 00]. Alex Oxford 4/21
8 Introduction Why do spectral methods get a bad press? 1. Dense matrices. 2. Ill-conditioned matrices. 3. When has it converged? Tricky Condition number 10 0 Error in solution See, for example: [Canuto et al. 07],, [Fornberg 98], [Trefethen 00]. Alex Oxford 4/21
9 A fast and well-conditioned spectral method Differentiation operator Work with coefficients: Spectral methods do not have to result in dense, illconditioned matrices. (Just don t discretize the differentiation operator faithfully.) The idea is to use simple relations between Chebyshev polynomials: 1 2 dt k dx = ku k 1, k 1, (U k U k 2 ), k 2, 1 T k = 2 0, k = 0, U 1, k = 1, U 0, k = 0. D = , S = Olver & T., A fast and well-conditioned spectral method, SIAM Review, Alex Oxford 5/21
10 A fast and well-conditioned spectral method Multiplication operator T j T k = 1 2 T j k + 2a 0 a 1 a 2 a 3... a 1 2a 0 a 1 a T j+k a 1 a 2 a 3 a 4... M[a] = 1 a 2 a 1 2a 0 a a 2 a 3 a 4 a a 3 a 2 a 1 2a 0... a 3 a 4 a 5 a } {{ } Toeplitz } {{ } Hankel + rank-1 Multiplication is not a dense operator in finite precision. It is m-banded: m a(x) = a k T k (x) = ã k T k (x) + O(ɛ), k=0 Alex Oxford 6/21 k=0
11 A fast and well-conditioned spectral method What about this new spectral method? 1. Almost banded matrices. 2. Well-conditioned matrices. 3. When has it converged? Trivial. Other approaches: [Clenshaw 57], [Greengard 91], [Shen 03]. Alex Oxford 7/21
12 A fast and well-conditioned spectral method What about this new spectral method? 1. Almost banded matrices. 2. Well-conditioned matrices. 3. When has it converged? Trivial. Other approaches: [Clenshaw 57], [Greengard 91], [Shen 03]. Alex Oxford 7/21
13 A fast and well-conditioned spectral method What about this new spectral method? 1. Almost banded matrices. 2. Well-conditioned matrices. 3. When has it converged? Trivial Condition number Other approaches: [Clenshaw 57], [Greengard 91], [Shen 03]. Alex Oxford 7/21
14 A fast and well-conditioned spectral method What about this new spectral method? 1. Almost banded matrices. 2. Well-conditioned matrices. 3. When has it converged? Trivial Condition number 10 0 Error in solution Other approaches: [Clenshaw 57], [Greengard 91], [Shen 03]. Alex Oxford 7/21
15 A fast and well-conditioned spectral method First example u (x) + x 3 u(x) = 100 sin(20,000x 2 ), u( 1) = 0. The exact solution is ( x ) u(x) = e x e t4 4 sin(20,000t 2 )dt. 1 u(x) degree(u) = time = 15.5s x N = chebop(@(x,u) ); N.lbc = 0; u = N \ f; Adaptively selects the discretisation size. Forms a chebfun object [Chebfun V4.2]. ũ u = Alex Oxford 8/21
16 A fast and well-conditioned spectral method Another example u 1 (x) + u(x) = 0, u( 1) = ,000x2 The exact solution with a = 50,000 is ( u(x) = exp tan 1 ( ax) + tan 1 ( ) a). a u(x) degree(u) = x Absolute error Old chebop Preconditioned collocation US method x Alex Oxford 9/21
17 A fast and well-conditioned spectral method A high-order example u (10) (x) + cosh(x)u (8) (x) + cos(x)u (2) (x) + x 2 u(x) = 0 u(±1) = 0, u (±1) = 1, u (k) (±1) = 0, k = 2, 3, u(x) u n+1 u n ( 1 1 (ũ(x) + ũ( x)) 2 ) 1 2 = x n Alex Oxford 10/21
18 Chebop and Chebop2 Convenience for the user L = chebop(@(x,u) diff(u,2)-x.*u,[-30 30]); % Airy equation L.lbc = 1; L.rbc = 0; % Set boundary conditions u = L \ 0; % u = chebfun Convert handle into discretisation instructions Construct discretisation A R n n Impose bcs and solve Ãx = b Converged to the solution? yes Construct a chebfun no, increase n L = chebop2(@(x,y,u) laplacian(u)+(1000+y)*u);% Helmholtz with gravity L.lbc = 1; L.rbc = 1; L.ubc = 1; L.dbc = 1;% Set boundary conditions u = L \ 0; % u = chebfun2 Convert handle into discretisation instructions Construct a matrix equation with an n y n x solution matrix Impose bcs and solve matrix equation Converged to the solution? yes Construct a chebfun2 no, increase n x or n y or both Alex Oxford 11/21
19 Chebop and Chebop2 Convenience for the user L = chebop(@(x,u) diff(u,2)-x.*u,[-30 30]); % Airy equation L.lbc = 1; L.rbc = 0; % Set boundary conditions u = L \ 0; % u = chebfun Convert handle into discretisation instructions Construct discretisation A R n n Impose bcs and solve Ãx = b Converged to the solution? yes Construct a chebfun no, increase n L = chebop2(@(x,y,u) laplacian(u)+(1000+y)*u);% Helmholtz with gravity L.lbc = 1; L.rbc = 1; L.ubc = 1; L.dbc = 1;% Set boundary conditions u = L \ 0; % u = chebfun2 Convert handle into discretisation instructions Construct a matrix equation with an n y n x solution matrix Impose bcs and solve matrix equation Converged to the solution? yes Construct a chebfun2 no, increase n x or n y or both Alex Oxford 11/21
20 Interpreting user-defined input Automatic differentiation Implemented by forward-mode operator overloading u xx + u yy + 50u + yu + Interpret anonymous function as a sequence of elementary operations + + * y.* u Can also calculate Fréchet derivatives u xx u yy u Key people: Ásgeir Birkisson and Toby Driscoll u u Alex Oxford 12/21
21 Low rank approximation Numerical rank For A C m n, SVD gives best rank k wrt 2-norm [Eckart & Young 1936] A = min(m,n) σ j u j v j k σ j u j v j, σ k+1 < tol. For Lipschitz smooth bivariate functions [Schmidt 1909, Smithies 1937] f(x, y) = σ j u j (y)v j (x) k σ j u j (y)v j (x). For compact linear operators acting on functions of two variables, L! = σ j L y L x j j k σ j L y L x j j. Alex Oxford 13/21
22 Low rank approximation Numerical rank For A C m n, SVD gives best rank k wrt 2-norm [Eckart & Young 1936] A = min(m,n) σ j u j v j k σ j u j v j, σ k+1 < tol. For Lipschitz smooth bivariate functions [Schmidt 1909, Smithies 1937] f(x, y) = σ j u j (y)v j (x) k σ j u j (y)v j (x). For compact linear operators acting on functions of two variables, L! = σ j L y L x j j k σ j L y L x j j. Alex Oxford 13/21
23 Low rank approximation Numerical rank For A C m n, SVD gives best rank k wrt 2-norm [Eckart & Young 1936] A = min(m,n) σ j u j v j k σ j u j v j, σ k+1 < tol. For Lipschitz smooth bivariate functions [Schmidt 1909, Smithies 1937] f(x, y) = σ j u j (y)v j (x) k σ j u j (y)v j (x). For compact linear operators acting on functions of two variables, L! = σ j L y L x j j k σ j L y L x j j. Alex Oxford 13/21
24 Low rank approximation Do the low rank stuff before discretization Low rank-then-discretize: Instead of low rank techniques after discretization, do them before. For example, Helmholtz is of rank 2 2 u + K 2 u = (u xx + K 2 2 u) + (u yy + K 2 2 u) = (D2 + K 2 2 I) I + I (D2 + K 2 Let A be your favourite ODE discretization of D 2 + K 2 2 I, then (typically) In general, if L is of rank k we have AXI + IXA T. 2 I). k A j XB T j = F Alex Oxford 14/21
25 Low rank approximation Do the low rank stuff before discretization Low rank-then-discretize: Instead of low rank techniques after discretization, do them before. For example, Helmholtz is of rank 2 2 u + K 2 u = (u xx + K 2 2 u) + (u yy + K 2 2 u) = (D2 + K 2 2 I) I + I (D2 + K 2 Let A be your favourite ODE discretization of D 2 + K 2 2 I, then (typically) In general, if L is of rank k we have AXI + IXA T. 2 I). k A j XB T j = F Alex Oxford 14/21
26 Low rank approximation Do the low rank stuff before discretization Low rank-then-discretize: Instead of low rank techniques after discretization, do them before. For example, Helmholtz is of rank 2 2 u + K 2 u = (u xx + K 2 2 u) + (u yy + K 2 2 u) = (D2 + K 2 2 I) I + I (D2 + K 2 Let A be your favourite ODE discretization of D 2 + K 2 2 I, then (typically) In general, if L is of rank k we have AXI + IXA T. 2 I). k A j XB T j = F Alex Oxford 14/21
27 Low rank approximation Computing the rank of a partial differential operator Recast differential operators as polynomials: Once you have polynomials computing the rank is easy. The rank of N y N x L = i=0 j=0 a ij (x, y) i y i j x j equals a TT-rank [Oseledets 2011] (between {x, s} and {y, t}) of N y N x k h(x, s, y, t) = a ij (s, t)y i x j = c j (t, y)r j (s, x). i=0 j=0 Rank 1: ODEs Trivial PDEs Rank 2: Laplace, Helmholtz Transport, Heat, Wave Black-Scholes Rank 3: Biharmonic Lots here. Alex Oxford 15/21
28 Low rank approximation Computing the rank of a partial differential operator Recast differential operators as polynomials: Once you have polynomials computing the rank is easy. The rank of N y N x L = i=0 j=0 a ij (x, y) i y i j x j equals a TT-rank [Oseledets 2011] (between {x, s} and {y, t}) of N y N x k h(x, s, y, t) = a ij (s, t)y i x j = c j (t, y)r j (s, x). i=0 j=0 Rank 1: ODEs Trivial PDEs Rank 2: Laplace, Helmholtz Transport, Heat, Wave Black-Scholes Rank 3: Biharmonic Lots here. Alex Oxford 15/21
29 Low rank approximation Computing the rank of a partial differential operator Recast differential operators as polynomials: Once you have polynomials computing the rank is easy. The rank of N y N x L = i=0 j=0 a ij (x, y) i y i j x j equals a TT-rank [Oseledets 2011] (between {x, s} and {y, t}) of N y N x k h(x, s, y, t) = a ij (s, t)y i x j = c j (t, y)r j (s, x). i=0 j=0 Rank 1: ODEs Trivial PDEs Rank 2: Laplace, Helmholtz Transport, Heat, Wave Black-Scholes Rank 3: Biharmonic Lots here. Alex Oxford 15/21
30 Low rank approximation Construct a n x by n y generalised Sylvester matrix equation If the PDE is Lu = f, where L is of rank-k then we solve for X C n y n x in, k X = solution s coefficients σ j A j XB T j = F, A j C n y n y, B j C n x n x. A j = A j, B j = 1D spectral discretization of L y j, Lx j Alex Oxford 16/21
31 Low rank approximation Matrix equation solvers Rank 1: A 1 XB T 1 = F. Solve A 1Y = F, then B 1 X T = Y T. Rank 2: A 1 XB T 1 + A 2XB T = F. Generalised Sylvester solver (RECSY) 2 [Jonsson & Kågström, 2002]. Rank k, k 3: Solve N N system using almost banded structure Execution time O(N 2 ) O(N 3/2 ) O(N) blue = rank 1 green = rank 2 red = rank N Alex Oxford 17/21
32 Examples Helmholtz equation 2 u + 2ω2 u = 0, u(x, ±1) = f (x, ±1), u(±1, y) = f (±1, y), where f (x, y) = cos(ωx) cos(ωy). Wave numbers vs. polynomial degree ω = n y Alex Oxford x ω 18/21
33 Examples Variable helmholtz equation N = chebop2(@(x,y,u) laplacian(u) (1/2+sin(x)ˆ2).*cos(y)ˆ2.*u); N.lbc = 1; N.rbc = 1; N.ubc = 1; N.dbc = 1; u = N \ chebfun2(@(x,y) cos(x.*y)); 10 0 Absolute maximum error N N = 1,050,625, error , time = 44.2s. Alex Oxford 19/21
34 Examples Wave and Klein Gordon equation N = chebop2(@(u) diff(u,2,1) - diff(u,2,2) + 5*u); % u_tt - u_xx + 5u N.dbc [u-exp(-10*x) diff(u)]; N.lbc = 0; N.rbc = 0; u = N \ 0; Alex Oxford 20/21
35 Conclusion Spectral methods do not have to be ill-conditioned. (Don t discretize differentiation faithfully.) Spectral methods are extremely convenient and flexible. As of 2014, global spectral methods are heavily restricted to a few geometries. Thank you for listening Alex Oxford 21/21
Chebfun2: Exploring Constant Coefficient PDEs on Rectangles
Chebfun2: Exploring Constant Coefficient PDEs on Rectangles Alex Townsend Oxford University with Sheehan Olver, Sydney University FUN13 workshop, 11th of April 2013 Alex Townsend Constant coefficient PDEs
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationN y. N x. y i x j, (1.1)
THE AUTOMATIC SOLUTION OF PARTIAL DIFFERENTIAL EQUATIONS USING A GLOBAL SPECTRAL METHOD ALEX TOWNSEND AND SHEEHAN OLVER Abstract. A spectral method for solving linear partial differential equations (PDEs)
More informationAlex Townsend Cornell University
Continuous analogues of Krylov methods for differential operators Alex Townsend Cornell University townsend@cornell.edu Joint work with: Marc Aurèle Gilles Based on Continuous analogues of Krylov methods
More informationarxiv: v1 [math.na] 4 Nov 2015
ON WELL-CONDITIONED SPECTRAL COLLOCATION AND SPECTRAL METHODS BY THE INTEGRAL REFORMULATION KUI DU arxiv:529v [mathna 4 Nov 25 Abstract Well-conditioned spectral collocation and spectral methods have recently
More informationComputing the common zeros of two bivariate functions via Bézout resultants Colorado State University, 26th September 2013
Work supported by supported by EPSRC grant EP/P505666/1. Computing the common zeros of two bivariate functions via Bézout resultants Colorado State University, 26th September 2013 Alex Townsend PhD student
More informationContinuous analogues of matrix factorizations NASC seminar, 9th May 2014
Continuous analogues of matrix factorizations NSC seminar, 9th May 2014 lex Townsend DPhil student Mathematical Institute University of Oxford (joint work with Nick Trefethen) Many thanks to Gil Strang,
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More informationTHE CHEBOP SYSTEM FOR AUTOMATIC SOLUTION OF DIFFERENTIAL EQUATIONS
BIT Numerical Mathematics (28)46:- c Springer 28. DOI:1.17/s1543---x THE CHEBOP SYSTEM FOR AUTOMATIC SOLUTION OF DIFFERENTIAL EQUATIONS TOBIN A. DRISCOLL 1, FOLKMAR BORNEMANN 2 and LLOYD N. TREFETHEN 3
More informationFast inverse transform sampling in one and two dimensions
Fast inverse transform sampling in one and two dimensions Sheehan Olver Alex Townsend Abstract We develop a computationally efficient and robust algorithm for generating pseudo-random samples from a broad
More informationChebfun: A New Kind of Numerical Computing
Chebfun: A New Kind of Numerical Computing R.B. Platte 1 and L.N. Trefethen 2 1 Oxford University Computing Laboratory, Oxford, OX1 3QD, UK (rodp@comlab.ox.ac.uk). After December 29: Department of Mathematics
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationChebfun: A New Kind of Numerical Computing
Report no. 08/13 Chebfun: A New Kind of Numerical Computing by Rodrigo B. Platte and Lloyd N. Trefethen The functionalities of the chebfun and chebop systems are surveyed. The chebfun system is a collection
More informationNumerical tensor methods and their applications
Numerical tensor methods and their applications 8 May 2013 All lectures 4 lectures, 2 May, 08:00-10:00: Introduction: ideas, matrix results, history. 7 May, 08:00-10:00: Novel tensor formats (TT, HT, QTT).
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 43: RBF-PS Methods in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 43 1 Outline
More informationSolving Boundary Value Problems (with Gaussians)
What is a boundary value problem? Solving Boundary Value Problems (with Gaussians) Definition A differential equation with constraints on the boundary Michael McCourt Division Argonne National Laboratory
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationRectangular spectral collocation
IMA Journal of Numerical Analysis Advance Access published February 6, 2015 IMA Journal of Numerical Analysis (2015) Page 1 of 25 doi:10.1093/imanum/dru062 Rectangular spectral collocation Tobin A. Driscoll
More informationTHE CHEBOP SYSTEM FOR AUTOMATIC SOLUTION OF DIFFERENTIAL EQUATIONS
BIT Numerical Mathematics (2008) 48: 701 723 DOI: 10.1007/s10543-008-0198-4 Published online: 22 November 2008 c Springer 2008 THE CHEBOP SYSTEM FOR AUTOMATIC SOLUTION OF DIFFERENTIAL EQUATIONS TOBIN A.
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationMath 671: Tensor Train decomposition methods
Math 671: Eduardo Corona 1 1 University of Michigan at Ann Arbor December 8, 2016 Table of Contents 1 Preliminaries and goal 2 Unfolding matrices for tensorized arrays The Tensor Train decomposition 3
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationBoundary Layer Problems and Applications of Spectral Methods
Boundary Layer Problems and Applications of Spectral Methods Sumi Oldman University of Massachusetts Dartmouth CSUMS Summer 2011 August 4, 2011 Abstract Boundary layer problems arise when thin boundary
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 43: RBF-PS Methods in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 43 1 Outline
More informationRBF Collocation Methods and Pseudospectral Methods
RBF Collocation Methods and Pseudospectral Methods G. E. Fasshauer Draft: November 19, 24 Abstract We show how the collocation framework that is prevalent in the radial basis function literature can be
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationMATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:
MATH 5 Final Examination May 3, 07 FORM A Name: Student Number: Section: This exam has 6 questions for a total of 50 points. In order to obtain full credit for partial credit problems, all work must be
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationProblem Set 1. This week. Please read all of Chapter 1 in the Strauss text.
Math 425, Spring 2015 Jerry L. Kazdan Problem Set 1 Due: Thurs. Jan. 22 in class. [Late papers will be accepted until 1:00 PM Friday.] This is rust remover. It is essentially Homework Set 0 with a few
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationAn Introduction to Partial Differential Equations
An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)
More informationSyllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department
Syllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department Alex Barnett, Scott Pauls, Dan Rockmore August 12, 2011 We aim to touch upon many topics that a professional
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationMATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:
MATH 5 Final Examination December 6, 5 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 5 points. In order to obtain full credit for partial credit problems, all work must
More informationA Automatic reformulation of ODEs to systems of first order equations
A Automatic reformulation of ODEs to systems of first order equations Á. BIRKISSON, University of Oxford Most numerical ODE solvers require problems to be written as systems of first order erential equations.
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationName: INSERT YOUR NAME HERE. Due to dropbox by 6pm PDT, Wednesday, December 14, 2011
AMath 584 Name: INSERT YOUR NAME HERE Take-home Final UWNetID: INSERT YOUR NETID Due to dropbox by 6pm PDT, Wednesday, December 14, 2011 The main part of the assignment (Problems 1 3) is worth 80 points.
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationScientific Computing: Dense Linear Systems
Scientific Computing: Dense Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 February 9th, 2012 A. Donev (Courant Institute)
More informationNumerical Computing without Discretization Woes
Numerical Computing without Discretization Woes Alex Townsend Cornell University townsend@cornell.edu give a visionary computational science talk Ben Zhang & Ricardo Baptista Discretization alleviating
More informationLeplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math 16.
Leplace s Analyzing the Analyticity of Analytic Analysis Engineering Math 16.364 1 The Laplace equations are built on the Cauchy- Riemann equations. They are used in many branches of physics such as heat
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationLinear Algebra. Session 12
Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)
More informationNumerical solution of linear time delay systems using Chebyshev-tau spectral method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Applications and Applied Mathematics: An International Journal (AAM) Vol., Issue (June 07), pp. 445 469 Numerical solution of linear time
More informationProper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,
More informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationQuestion 9: PDEs Given the function f(x, y), consider the problem: = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1. x 2 u. u(x, 0) = u(x, 1) = 0 for 0 x 1
Question 9: PDEs Given the function f(x, y), consider the problem: 2 u x 2 u = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1 u(x, 0) = u(x, 1) = 0 for 0 x 1 u(0, y) = u(1, y) = 0 for 0 y 1. a. Discuss how you
More informationarxiv: v1 [math.na] 2 Jul 2015
A FAST AND WELL-CONDITIONED SPECTRAL METHOD FOR SINGULAR INTEGRAL EQUATIONS RICHARD MIKAEL SLEVINSKY AND SHEEHAN OLVER arxiv:507.00596v [math.na] Jul 05 Abstract. We develop a spectral method for solving
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More information11.3 MATLAB for Partial Differential Equations
276 3. Generate the shape functions N (i) j = a (i) j where i =1, 2, 3,..., m and j =1, 2, 3,..., m. + b (i) j x + c(i) j y (11.2.17) 4. Compute the integrals for matrix elements α ij and vector elements
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationA HIERARCHICAL 3-D DIRECT HELMHOLTZ SOLVER BY DOMAIN DECOMPOSITION AND MODIFIED FOURIER METHOD
SIAM J. SCI. COMPUT. Vol. 26, No. 5, pp. 1504 1524 c 2005 Society for Industrial and Applied Mathematics A HIERARCHICAL 3-D DIRECT HELMHOLTZ SOLVER BY DOMAIN DECOMPOSITION AND MODIFIED FOURIER METHOD E.
More informationChapter 3. Second Order Linear PDEs
Chapter 3. Second Order Linear PDEs 3.1 Introduction The general class of second order linear PDEs are of the form: ax, y)u xx + bx, y)u xy + cx, y)u yy + dx, y)u x + ex, y)u y + f x, y)u = gx, y). 3.1)
More informationRadial basis function partition of unity methods for PDEs
Radial basis function partition of unity methods for PDEs Elisabeth Larsson, Scientific Computing, Uppsala University Credit goes to a number of collaborators Alfa Ali Alison Lina Victor Igor Heryudono
More informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
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 informationMath 671: Tensor Train decomposition methods II
Math 671: Tensor Train decomposition methods II Eduardo Corona 1 1 University of Michigan at Ann Arbor December 13, 2016 Table of Contents 1 What we ve talked about so far: 2 The Tensor Train decomposition
More informationCHAPTER 4. Introduction to the. Heat Conduction Model
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 4 A COLLECTION OF HANDOUTS ON PARTIAL DIFFERENTIAL EQUATIONS
More informationSingular Value Decomposition of the frequency response operator
Measures of input-output amplification (across frequency) Largest singular value Hilbert-Schmidt norm (power spectral density) Systems with one spatial variable Two point boundary value problems EE 8235:
More informationPROJECTION METHODS FOR DYNAMIC MODELS
PROJECTION METHODS FOR DYNAMIC MODELS Kenneth L. Judd Hoover Institution and NBER June 28, 2006 Functional Problems Many problems involve solving for some unknown function Dynamic programming Consumption
More informationLecture 1: Introduction
CBMS Conference on Fast Direct Solvers Dartmouth College June June 7, 4 Lecture : Introduction Gunnar Martinsson The University of Colorado at Boulder Research support by: Many thanks to everyone who made
More informationSolution of Linear Equations
Solution of Linear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 7, 07 We have discussed general methods for solving arbitrary equations, and looked at the special class of polynomial equations A subclass
More informationExplicit construction of rectangular differentiation matrices
IMA Journal of Numerical Analysis (206) 36, 68 632 doi:0.093/imanum/drv03 Advance Access publication on April 26, 205 Explicit construction of rectangular differentiation matrices Kuan Xu Mathematical
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationApplied Mathematics 505b January 22, Today Denitions, survey of applications. Denition A PDE is an equation of the form F x 1 ;x 2 ::::;x n ;~u
Applied Mathematics 505b January 22, 1998 1 Applied Mathematics 505b Partial Dierential Equations January 22, 1998 Text: Sobolev, Partial Dierentail Equations of Mathematical Physics available at bookstore
More informationMATH 251 Final Examination December 16, 2014 FORM A. Name: Student Number: Section:
MATH 2 Final Examination December 6, 204 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 0 points. In order to obtain full credit for partial credit problems, all work must
More informationA WELL-CONDITIONED COLLOCATION METHOD USING A PSEUDOSPECTRAL INTEGRATION MATRIX
A WELL-CONDITIONED COLLOCATION METHOD USING A PSEUDOSPECTRAL INTEGRATION MATRIX LI-LIAN WANG, MICHAEL DANIEL SAMSON, AND XIAODAN ZHAO Abstract. In this paper, a well-conditioned collocation method is constructed
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationSpeeding up numerical computations via conformal maps
Speeding up numerical computations via conformal maps Nick Trefethen, Oxford University Thanks to Nick Hale, Nick Higham and Wynn Tee 1/35 SIAM 1997 SIAM 2000 Cambridge 2003 2/35 Princeton 2005 Bornemann
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationMATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination May 4, 2015 FORM A Name: Student Number: Section: This exam has 16 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work must
More informationLeplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math EECE
Leplace s Analyzing the Analyticity of Analytic Analysis Engineering Math EECE 3640 1 The Laplace equations are built on the Cauchy- Riemann equations. They are used in many branches of physics such as
More information2.3. VECTOR SPACES 25
2.3. VECTOR SPACES 25 2.3 Vector Spaces MATH 294 FALL 982 PRELIM # 3a 2.3. Let C[, ] denote the space of continuous functions defined on the interval [,] (i.e. f(x) is a member of C[, ] if f(x) is continuous
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More information