Chebfun2: Exploring Constant Coefficient PDEs on Rectangles
|
|
- Myra Boone
- 5 years ago
- Views:
Transcription
1 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 0/11
2 A simplified history of the Chebfun project Piecewise smooth (2010) Endpoint singularities (2010*) Reducing regularity Blow up functions (2011) Chebfun (2004) Linear ODEs (2008) Ṅonlinear ODEs (2012) Chebgui (2011*) Ordinary differential equations Chebfun2 (2013) Piecewise smooth Linear PDEs Arbitrary domains Nonlinear PDEs Two dimensions Alex Townsend Constant coefficient PDEs 1/11
3 Chebop and Chebop2: Overall workflow L = chebop(@(x,u) diff(u,2)-x*u, [-30 30]); Llbc = 1; Lrbc = 0; u = L \ 0; % Airy eqn % bcs % 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 [Driscoll, Bornemann, & Trefethen, 2008], [Birkisson & Driscoll, 2012] Alex Townsend Constant coefficient PDEs 2/11
4 Chebop and Chebop2: Overall workflow L = chebop(@(x,u) diff(u,2)-x*u, [-30 30]); Llbc = 1; Lrbc = 0; u = L \ 0; % Airy eqn % bcs % 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 [Driscoll, Bornemann, & Trefethen, 2008], [Birkisson & Driscoll, 2012] L = chebop2(@(u) diffx(u,2)+diffy(u,2)+100*u); Llbc = 1; Lrbc = 1; Lubc = 1; Ldbc = 1; u = L \ 0; % Helmholtz % bcs % u = chebfun2 Convert handle into discretisation instructions Construct a m by n generalised Sylvester matrix equation Impose bcs and solve matrix equation Converged to the solution? yes Construct a chebfun2 no, increase m or n or both Alex Townsend Constant coefficient PDEs 2/11
5 The numerical rank of a matrix, function, and operator For an m n matrix, Numerical Rank A = min(m,n) j=1 σ j u j v T j k σ j u j vj T, σ k+1 < tol For a function of two variables (this is what Chebfun2 uses), j=1 f (x, y)! = k σ j u j (y)v j (x) σ j u j (y)v j (x) j=1 j=1 For operators acting on functions of two variables, k L =! σ j L y j L x j σ j L y j L x j j=1 j=1 Discretise by a 1D spectral method [Olver & T, 2013] Alex Townsend Constant coefficient PDEs 3/11
6 Determining the rank of a PDE operator Given a PDE with constant coefficients of the form, L = N N j j=0 i=0 a ij i+j i x j y, A = (a ij), We can rewrite this as Find by AD, thanks A Birkisson 0 / x 0 L = D(y) T AD(x), D(x) = N / x N Hence, if A = UΣV T is the truncated SVD of A then ( ) ( ) ( ) L = D(y) T UΣV T D(x) = D(y) T U Σ V T D(x) This low rank representation for L tells us how to discretize and solve the PDE Alex Townsend Constant coefficient PDEs 4/11
7 Examples Helmholtz is of rank 2, u xx + u yy + 100u = (u xx + 50u) + (u yy + 50u) The SVD obtains the rank 2 expression, u xx + u yy + 100u = ( 0009u yy 0999u ) ( 0009u xx 0999u) ( 0999u yy 0009u ) (0009u xx 0999u) Rank 1: ODEs Trivial PDEs Rank 3: Almost everything else Rank 2: Laplace Helmholtz Transport Heat Wave The method of separation of variables is useful Alex Townsend Constant coefficient PDEs 5/11
8 Construct a m by n generalised Sylvester matrix equation If the PDE is Lu = f, where L is of rank-k then we solve for X R m n in, k σ j A j XBj T = F, A j R m m, B j R n n j=1 solution s ċoefficients 1D spectral ḋiscretization of Ly j, Lx j A j = ( ) Alex Townsend Constant coefficient PDEs 6/11
9 Matrix equation solvers Rank-1: A 1 XB1 T = F Solve A 1Y = F, then B 1 X T = Y T Rank-2: A 1 XB1 T + A 2XB2 T = F Generalised Bartels Stewart algorithm [Gardener, Laub, Amato, & Moler, 1992] Rank-k, k 3: Solve mn mn system with \ in Matlab 10 2 Time to solve matrix equation O(m 2 n 2 ) O(m 3 + n 3 ) Solve time O(mn) k = 1 k = 2 k m & n Alex Townsend Constant coefficient PDEs 7/11
10 Have we converged? (simplified version) Increase m and n until the tail of the solution s Chebyshev coefficients is near/below machine precision, and then truncate Always 2 j + 1 Truncated coefficients These coefficients can then be used by Chebfun2 to construct a low rank function approximation [T & Trefethen, 2013] Alex Townsend Constant coefficient PDEs 8/11
11 Overview and recap L = chebop2(@(u) diffx(u,2)+diffy(u,2)+100*u); Llbc = 1; Lrbc = 1; Lubc = 1; Ldbc = 1; u = L \ 0; % Helmholtz % bcs % u = chebfun2 Convert handle into discretisation instructions Construct a m by n generalised Sylvester matrix equation Impose bcs and solve matrix equation Converged to the solution? yes Construct a chebfun2 no, increase m or n or both Alex Townsend Constant coefficient PDEs 9/11
12 Demo DEMO Alex Townsend Constant coefficient PDEs 10/11
13 Future Work I would like to be able to do vector-valued PDEs such as Stokes equation: For Ω = [ 1, 1] [ 1, 1], u + p = f in Ω, u = 0 in Ω, u = 0 on Ω, p(x, y)dxdy = 0 in Ω Aiming for the following Chebfun2 syntax: N = chebop2v(@(u,p) [-lap(u) + grad(p); div(u)]); Nlbc = 0; Nrbc = 0; Nubc = 0; Ndbc = 0; Nbc sum2(p); uu = N \ [f;0]; Alex Townsend Constant coefficient PDEs 11/11
14 References [Birkisson & Driscoll, 2012] Automatic Fréchet differentiation for the numerical solution of boundary-value problems, (ACM Trans Math Softw, 2012) [Driscoll, Bornemann, & Trefethen, 2008] The chebop system for automatic solution of differential equations (BIT Numer Math, 2008) [Gardener, Laub, Amato, & Moler, 1992] Solution of the Sylvester matrix equation AXB T + CXD T = E, (ACM Trans Math Softw, 1992) [Olver & T, 2013] A fast and well-conditioned spectral method, (to appear in SIAM Review, 2013) [T & Trefethen, 2013] An extension of Chebfun to two dimensions, (submitted, 2013) Alex Townsend Constant coefficient PDEs 12/11
The automatic solution of PDEs using a global spectral method
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,
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 informationA Continuation Approach to a Quadratic Matrix Equation
A Continuation Approach to a Quadratic Matrix Equation Nils Wagner nwagner@mecha.uni-stuttgart.de Institut A für Mechanik, Universität Stuttgart GAMM Workshop Applied and Numerical Linear Algebra September
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 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 informationPositive Solution to a Generalized Lyapunov Equation via a Coupled Fixed Point Theorem in a Metric Space Endowed With a Partial Order
Filomat 29:8 2015, 1831 1837 DOI 10.2298/FIL1508831B Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Positive Solution to a Generalized
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 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 informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
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 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 informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Numerical Linear Algebra Background Cho-Jui Hsieh UC Davis May 15, 2018 Linear Algebra Background Vectors A vector has a direction and a magnitude
More informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Lecture 5: Numerical Linear Algebra Cho-Jui Hsieh UC Davis April 20, 2017 Linear Algebra Background Vectors A vector has a direction and a magnitude
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
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 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 informationComplex Analysis Homework 4
Complex Analysis Homework 4 Isaac Defrain Steve Clanton David Holz April 9, 2009 Problem 1 Problem. Show that if f = u + iv is analytic then ( v) ( v) = 0. Explain this geometrically. Since f is analytic,
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
More information1 Sylvester equations
1 Sylvester equations Notes for 2016-11-02 The Sylvester equation (or the special case of the Lyapunov equation) is a matrix equation of the form AX + XB = C where A R m m, B R n n, B R m n, are known,
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
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 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 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 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 informationMath 322. Spring 2015 Review Problems for Midterm 2
Linear Algebra: Topic: Linear Independence of vectors. Question. Math 3. Spring Review Problems for Midterm Explain why if A is not square, then either the row vectors or the column vectors of A are linearly
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 informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
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 informationMath 126 Winter 2012: Rough lecture notes
Math 126 Winter 2012: Rough lecture notes Alex Barnett, Dartmouth College March 1, 2012 1 Introduction Numerical mathematics is at the intersection of analysis (devising and proving theorems), computation
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 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 informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
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 informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More information194 B. TSANG, S.W. TAYLOR AND G.C. WAKE and Forray [4]. The question as to which boundary value problems have a variational formulation in which all o
cfl Journal of Applied Mathematics & Decision Sciences, 4(2), 193 24 (2) Reprints Available directly from the Editor. Printed in New Zealand. Variational methods for boundary value problems B. TSANG b.tsang@math.auckland.ac.nz
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 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 informationLinear Algebra Math 221
Linear Algebra Math 221 Open Book Exam 1 Open Notes 3 Sept, 24 Calculators Permitted Show all work (except #4) 1 2 3 4 2 1. (25 pts) Given A 1 2 1, b 2 and c 4. 1 a) (7 pts) Bring matrix A to echelon form.
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
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 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 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 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 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 informationPropagation of discontinuities in solutions of First Order Partial Differential Equations
Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in
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 informationChapter 2: First Order DE 2.6 Exact DE and Integrating Fa
Chapter 2: First Order DE 2.6 Exact DE and Integrating Factor First Order DE Recall the general form of the First Order DEs (FODE): dy dx = f(x, y) (1) (In this section x is the independent variable; not
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 informationLecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018
Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent
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 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 informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationTHE SINGULAR VALUE DECOMPOSITION MARKUS GRASMAIR
THE SINGULAR VALUE DECOMPOSITION MARKUS GRASMAIR 1. Definition Existence Theorem 1. Assume that A R m n. Then there exist orthogonal matrices U R m m V R n n, values σ 1 σ 2... σ p 0 with p = min{m, n},
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 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 informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationLinear maps. Matthew Macauley. Department of Mathematical Sciences Clemson University Math 8530, Spring 2017
Linear maps Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 8530, Spring 2017 M. Macauley (Clemson) Linear maps Math 8530, Spring 2017
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationLecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems
Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now
More informationLeast Squares. Tom Lyche. October 26, Centre of Mathematics for Applications, Department of Informatics, University of Oslo
Least Squares Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo October 26, 2010 Linear system Linear system Ax = b, A C m,n, b C m, x C n. under-determined
More informationMATH 304 Linear Algebra Lecture 15: Linear transformations (continued). Range and kernel. Matrix transformations.
MATH 304 Linear Algebra Lecture 15: Linear transformations (continued). Range and kernel. Matrix transformations. Linear mapping = linear transformation = linear function Definition. Given vector spaces
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
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 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 informationFactorized Solution of Sylvester Equations with Applications in Control
Factorized Solution of Sylvester Equations with Applications in Control Peter Benner Abstract Sylvester equations play a central role in many areas of applied mathematics and in particular in systems and
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
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 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 informationMathematical Computing
IMT2b2β Department of Mathematics University of Ruhuna A.W.L. Pubudu Thilan Differential Equations Types of Differential Equations Differential equations can basically be classified as ordinary differential
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 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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
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 informationMaths III - Numerical Methods
Maths III - Numerical Methods Matt Probert matt.probert@york.ac.uk 4 Solving Differential Equations 4.1 Introduction Many physical problems can be expressed as differential s, but solving them is not always
More informationFirst and Second Order ODEs
Civil Engineering 2 Mathematics Autumn 211 M. Ottobre First and Second Order ODEs Warning: all the handouts that I will provide during the course are in no way exhaustive, they are just short recaps. Notation
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
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 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 informationMATH 612 Computational methods for equation solving and function minimization Week # 2
MATH 612 Computational methods for equation solving and function minimization Week # 2 Instructor: Francisco-Javier Pancho Sayas Spring 2014 University of Delaware FJS MATH 612 1 / 38 Plan for this week
More informationExamination paper for TMA4130 Matematikk 4N: SOLUTION
Department of Mathematical Sciences Examination paper for TMA4 Matematikk 4N: SOLUTION Academic contact during examination: Morten Nome Phone: 9 84 97 8 Examination date: December 7 Examination time (from
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationLecture Notes for Math 251: ODE and PDE. Lecture 7: 2.4 Differences Between Linear and Nonlinear Equations
Lecture Notes for Math 51: ODE and PDE. Lecture 7:.4 Differences Between Linear and Nonlinear Equations Shawn D. Ryan Spring 01 1 Existence and Uniqueness Last Time: We developed 1st Order ODE models for
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 informationBarycentric interpolation, Chebfun, rational approximation, Bernstein ellipse,
STABILITY OF BARYCENTRIC INTERPOLATION FORMULAS MARCUS WEBB, LLOYD N. TREFETHEN, AND PEDRO GONNET Abstract. The barycentric interpolation formula defines a stable algorithm for evaluation at points in
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationSymmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota
Symmetry Methods for Differential Equations and Conservation Laws Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Santiago, November, 2010 Symmetry Groups of Differential Equations
More informationExploiting off-diagonal rank structures in the solution of linear matrix equations
Stefano Massei Exploiting off-diagonal rank structures in the solution of linear matrix equations Based on joint works with D. Kressner (EPFL), M. Mazza (IPP of Munich), D. Palitta (IDCTS of Magdeburg)
More informationMATH20411 PDEs and Vector Calculus B
MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationarxiv: v1 [math.na] 6 Dec 2015
arxiv:1512.01803v1 [math.na] 6 Dec 2015 Chopping a Chebyshev series JARED L. AURENTZ and LLOYD N. TREFETHEN December 8, 2015 Abstract Chebfun and related software projects for numerical computing with
More informationGene Golub Day Memorial Workshop HKBU. SVD the Swiss Army Knife of Linear Algebra. Walter Gander. February 7, 2015
1 Gene Golub Day Memorial Workshop HKBU SVD the Swiss Army Knife of Linear Algebra Walter Gander February 7, 15 2 Swiss Army Knife? Dianne P. O'Leary Matrix Factorizations for Information Retrieval (06)
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 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 informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationComputational Multilinear Algebra
Computational Multilinear Algebra Charles F. Van Loan Supported in part by the NSF contract CCR-9901988. Involves Working With Kronecker Products b 11 b 12 b 21 b 22 c 11 c 12 c 13 c 21 c 22 c 23 c 31
More information