Generation of Finite Difference Formulas on Arbitrarily Spaced Grids
|
|
- Clifford Stevens
- 5 years ago
- Views:
Transcription
1 MATHEMATICS OF COMPUTATION VOLUME, NUMBER OCTOBER 9, PAGES 99 Generation of Finite Difference Formulas on Arbitrarily Spaced Grids By Bengt Fornberg Abstract. Simple recursions are derived for calculating the weights in compact finite difference formulas for any order of derivative and to any order of accuracy on onedimensional grids with arbitrary spacing. Tables are included for some special cases (of equispaced grids).. Introduction. Previously published methods to generate finite difference weights (e.g., references [l]-[]) have been of considerable complexity and often been limited to derivatives of low order on equidistantly spaced grids. The most ambitious attempt to tabulate weights for many orders of derivatives and to high orders of accuracy appears to be the work by Keller and Pereyra []. However, their algorithms (limited to equispaced grids) were very involved, and the resulting tables contain both isolated and systematic errors. In the present study we describe two simple recursion relations which give the weights for any order of derivative (including the Oth derivative, corresponding to interpolation), approximated to any order of accuracy on an arbitrary grid in one dimension. Since, in general, only four arithmetic operations are needed to determine each weight, the main anticipated application of the present method is to dynamically changing grids. However, the method is also well suited to generate tables of weights. Such tables (in the special case of equispaced grids, up to the th derivative and up to 9 weights) are included in the cases of one-sided and centered approximations at a grid point and at a 'half-way point' between grid points.. Notation, Algorithm. Given M >, the order of the highest derivative we wish to approximate, and a set of N + grid points (at i-coordinates q,..., a^; N > ), the problem is to find all the weights such that the approximations dmf dxm x=xo n Jèn,uliau), m =,,...,M; n = m,m +,...,N, i/=o become of optimal formal order of accuracy (in general of order n-m +, although it can be higher in special cases). The following algorithm achieves this: Received May, 9; revised December, 9. 9 Mathematics Subject Classification (9 Revision). Primary D. Key words and phrases. Finite difference coefficients, high-order accuracy American Mathematical Society / $. + $. per page
2 BENGT FORNBERG Enter M, N, x, a,ai,a,.,o = cl := for n := to TV do c:=l for v := to n do..,o/v c := an - Oiv c := c c if n < M then < _j v := for m := to min(n, M) do W,u = (K - *o)c-l, - <l!j/c next m next v for m := to min(n, M) do Cn := Càimèn-l,n-l " K-l - *o)c-l, -l) next m cl := c next n Notes.. If the array ^ initially is zero, the statement "if n < M then $i-i,t/ := " can ^e omitted.. In the case of m = (corresponding to interpolation formulas), expressions of the form 'zero*(undefined number)' occur. The result is assumed to be zero.. The order in which the av (all distinct) are given is significant (since the weights corresponding to all leading subsets of the a 's are calculated). There is no restriction on xq coinciding with any a.. Derivation of the Algorithm. For simplicity, assume we seek to approximate the derivatives at the point in =. Let {ao,ot\,...,a ^) be distinct real numbers and denote n (.) wn(x) := J (x-afc). k= The polynomial (-) FnAx) =,, ""(X)-v is the one of minimal degree which takes the value at x = au and at x = ak, < k < n, k t v. For an arbitrary function f(x) and nodes x = a, Lagrange's interpolation polynomial becomes n (.) p(x):=jfnax)f (<*») i/= The desired weights express how the values of [dmp(x)/dxm]x=o vary with changes in f(au). Since only one term in p(x) is influenced by changes in each f(au), we find (-) d" C/ = FnAx) dxr lx=
3 GENERATION OF FINITE DIFFERENCE FORMULAS Therefore, the nth degree polynomial Fn,v(x) can also be expressed as " m (-) Fn, (x) = -^f *m m= From (.) follow (noting that u>(x) = (x an)w _i(i) implies u>'n(x) = (x - an)u'n-lix) + Un-l(x)) (.) FnAx) = ^^Fn-lAx) a - an and to n\ r> I \ wn-l(z) <*>n-(û!n-l),,. (.) F, (x) =-^-=- r^(x-an_i)fn_i, _i(a;) n>l. Wn-l(«n) Wn.i(an) By substituting the expansion (.) into (.) and (.), and by equating powers of x, the desired recursion relations between the weights are obtained: (-) C = -^ KC-,/ - mc-tj and (.9) < = ^(";-;)Kffi-t - a -ic-i, -i). The relation Vf J' (-) «V = (o i/= * ' m >, could be used instead of (.9) to obtain > n. However, this would increase the operation count and might also cause a growth of errors in the case of floatingpoint arithmetic.. Description of the Tables. Special cases which commonly occur are centered and one-sided approximations on equidistant grids. The particular choices of ctu used for Tables - correspond to grid spacings Ax =. For other values of Ax, these coefficients should be divided by (Ax)m (where m, as before, is the order of the derivative).
4 BENGT FORNBERG TABLE Some weights for centered approximations at a grid point (generated by setting M, TV =, xo = and av = {,, -,, -,, -,, -}). Approximations at x = ; i-coordinates at nodes: - - _ _ _ ^ ^ - ^ ^ - _ -9 9 i Si ^-I ^ _
5 GENERATION OF FINITE DIFFERENCE FORMULAS TABLE Some weights for centered approximations at a 'half-way' point (generated by setting M =, N =, x = and av = {/,-/,/,-/,/,-/,/,/}). o dv e a r t v f e O r a ÎC r V?ca f y Approximations at x = ; x-coordinates at nodes: / -/ -/ -/ / / / / _ - - _ _9 9 ZI ^ = -9 9 il ^ - - -s 9
6 BENGT FORNBERG TABLE Some weights for one-sided approximations at a grid point (generated by setting M =, N =, x = and au = {,,,,,,,,}). d v e a r t Approximations at x = ; i-coordinates at nodes: v f e - ^ ^ i\j o - O A Ü = il - 9
7 GENERATION OF FINITE DIFFERENCE FORMULAS TABLE Some weights for one-sided approximations at a 'half-way'' point (generated by setting M=, N=, x= andav = {-l/, /,/, /, /,9/,/,/,/}). d v e a f è O r a ÎÎ r? f y Approximations at x = ; i-coordinates at nodes: -/ / / / / 9/ / / / ^ _ 9 9 _ '.' G ^ ii - -9
8 BENGT FORNBERG Acknowledgment. The comments of the referee have been most helpful in generalizing and simplifying the present algorithm. Exxon Research and Engineering Annandale, New Jersey Company. M. AbramOWITZ k I. A. STEGUN, Handbook of Mathematical Functions, National Bureau of Standards, Appl. Math. Series, no., Washington, D. C., 9.. S. J. CYNAR, "Using Gaussian elimination for computation of the central difference equation coeñcienxs," SIGNUM Newsl. (USA), v., 9, pp B. FORNBERG, "On a Fourier method for the integration of hyperbolic equations," SIAM J. Numer. Anal., v., 9, pp H. B. KELLER & V. PEREYRA, "Symbolic generation of finite difference formulas," Math. Comp., v., 9, pp W. D. LAKIN, "Differentiating matrices for arbitrarily spaced grid points," Internat. J. Numer. Methods Engrg., v., 9, pp. 9-.
Finite Difference Calculus
Chapter 2 Finite Difference Calculus In this chapter we review the calculus of finite differences. The topic is classic and covered in many places. The Taylor series is fundamental to most analysis. A
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationTable for Third-Degree Spline Interpolation Using Equi-Spaced Knots. By W. D. Hoskins
MATHEMATICS OF COMPUTATION, VOLUME 25, NUMBER 116, OCTOBER, 1971 Table for Third-Degree Spline Interpolation Using Equi-Spaced Knots By W. D. Hoskins Abstract. A table is given for the calculation of the
More informationFebruary 23 Math 2335 sec 51 Spring 2016
February 23 Math 2335 sec 51 Spring 2016 Section 4.1: Polynomial Interpolation Interpolation is the process of finding a curve or evaluating a function whose curve passes through a known set of points.
More informationSeptember 18, B Math Test Chapter 1 Name: x can be expressed as: {y y 0, y R}.
September 8, 208 62B Math Test Chapter Name: Part : Objective Questions [ mark each, total 2 marks]. State whether each of the following statements is TRUE or FALSE a) The mapping rule (x, y) (-x, y) represents
More informationNatasha S. Sharma, PhD
Revisiting the function evaluation problem Most functions cannot be evaluated exactly: 2 x, e x, ln x, trigonometric functions since by using a computer we are limited to the use of elementary arithmetic
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationLinear Interpolating Splines
Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 17 Notes Tese notes correspond to Sections 112, 11, and 114 in te text Linear Interpolating Splines We ave seen tat ig-degree polynomial interpolation
More informationThe Total Variation Approach to Approximate Hyperbolic Wave Equations
2 CUMMING The Total Variation Approach to Approximate Hyperbolic Wave Equations Christine Cumming Faculty Sponsor: Dr. Robert H. Hoar, Department of Math ABSTRACT From ultrasound for medical purposes to
More informationSharing Several Secrets based on Lagrange s Interpolation formula and Cipher Feedback Mode
Int. J. Nonlinear Anal. Appl. 5 (2014) No. 2, 60-66 ISSN: 2008-6822 (electronic) http://www.ijnaa.semnan.ac.ir Sharing Several Secrets based on Lagrange s Interpolation formula and Cipher Feedback Mode
More informationNatural Quartic Spline
Natural Quartic Spline Rafael E Banchs INTRODUCTION This report describes the natural quartic spline algorithm developed for the enhanced solution of the Time Harmonic Field Electric Logging problem As
More informationMAT128A: Numerical Analysis Lecture One: Course Logistics and What is Numerical Analysis?
MAT128A: Numerical Analysis Lecture One: Course Logistics and What is Numerical Analysis? September 26, 2018 Lecture 1 September 26, 2018 1 / 19 Course Logistics My contact information: James Bremer Email:
More informationPolynomial Approximation and Interpolation Chapter 4
4.4 LAGRANGE POLYNOMIALS The direct fit polynomial presented in Section 4.3, while quite straightforward in principle, has several disadvantages. It requires a considerable amount of effort to solve the
More informationHyperbola for Curvilinear Interpolation
Applied Mathematical Sciences, Vol. 7, 2013, no. 30, 1477-1481 HIKARI Ltd, www.m-hikari.com Hyperbola for Curvilinear Interpolation G. L. Silver 868 Kristi Lane Los Alamos, NM 87544, USA gsilver@aol.com
More informationGeneral Method for Exponential-Type Equations. for Eight- and Nine-Point Prismatic Arrays
Applied Mathematical Sciences, Vol. 3, 2009, no. 43, 2143-2156 General Method for Exponential-Type Equations for Eight- and Nine-Point Prismatic Arrays G. L. Silver Los Alamos National Laboratory* P.O.
More informationMathematics Scope & Sequence Algebra I
Mathematics Scope & Sequence 2016-17 Algebra I Revised: June 20, 2016 First Grading Period (24 ) Readiness Standard(s) Solving Equations and Inequalities A.5A solve linear equations in one variable, including
More informationIS BINARY ENCODING APPROPRIATE FOR THE PROBLEM-LANGUAGE RELATIONSHIP?
Theoretical Computer Science 19 (1982) 337-341 North-Holland Publishing Company NOTE IS BINARY ENCODING APPROPRIATE FOR THE PROBLEM-LANGUAGE RELATIONSHIP? Nimrod MEGIDDO Statistics Department, Tel Aviv
More informationParallel Implementations of Gaussian Elimination
s of Western Michigan University vasilije.perovic@wmich.edu January 27, 2012 CS 6260: in Parallel Linear systems of equations General form of a linear system of equations is given by a 11 x 1 + + a 1n
More informationAlgorithms and Data Structures
Charles A. Wuethrich Bauhaus-University Weimar - CogVis/MMC June 22, 2017 1/51 Introduction Matrix based Transitive hull All shortest paths Gaussian elimination Random numbers Interpolation and Approximation
More information2 Computation with Floating-Point Numbers
2 Computation with Floating-Point Numbers 2.1 Floating-Point Representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However, real numbers
More informationApproximation of a Fuzzy Function by Using Radial Basis Functions Interpolation
International Journal of Mathematical Modelling & Computations Vol. 07, No. 03, Summer 2017, 299-307 Approximation of a Fuzzy Function by Using Radial Basis Functions Interpolation R. Firouzdor a and M.
More informationPiecewise Polynomial Interpolation, cont d
Jim Lambers MAT 460/560 Fall Semester 2009-0 Lecture 2 Notes Tese notes correspond to Section 4 in te text Piecewise Polynomial Interpolation, cont d Constructing Cubic Splines, cont d Having determined
More informationFriday, 11 January 13. Interpolation
Interpolation Interpolation Interpolation is not a branch of mathematic but a collection of techniques useful for solving computer graphics problems Basically an interpolant is a way of changing one number
More informationInterpolation. TANA09 Lecture 7. Error analysis for linear interpolation. Linear Interpolation. Suppose we have a table x x 1 x 2...
TANA9 Lecture 7 Interpolation Suppose we have a table x x x... x n+ Interpolation Introduction. Polynomials. Error estimates. Runge s phenomena. Application - Equation solving. Spline functions and interpolation.
More informationGeneralised Mean Averaging Interpolation by Discrete Cubic Splines
Publ. RIMS, Kyoto Univ. 30 (1994), 89-95 Generalised Mean Averaging Interpolation by Discrete Cubic Splines By Manjulata SHRIVASTAVA* Abstract The aim of this work is to introduce for a discrete function,
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 19 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More information9-1: Slope NAME: 1. What do you think is meant by the terms rise and run?
9-1: Slope NAME: CUES: PER: DATE: 1. What do you think is meant by the terms rise and run? 2. What is the vertical change between: a. points A and B? b. points A and C? c. points C and D? 3. What is the
More informationspecified or may be difficult to handle, we often have a tabulated data
Interpolation Introduction In many practical situations, for a function which either may not be explicitly specified or may be difficult to handle, we often have a tabulated data where and for In such
More informationLAB 2: Linear Equations and Matrix Algebra. Preliminaries
Math 250C, Section C2 Hard copy submission Matlab # 2 1 Revised 07/13/2016 LAB 2: Linear Equations and Matrix Algebra In this lab you will use Matlab to study the following topics: Solving a system of
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
More informationSect 3.1 Quadratic Functions and Models
Objective 1: Sect.1 Quadratic Functions and Models Polynomial Function In modeling, the most common function used is a polynomial function. A polynomial function has the property that the powers of the
More informationinterpolation, multivariate, Newton polynomial, divided difference, algorithm,
MULTIVARIATE POLYNOMIAL INTERPOLATION IN NEWTON FORMS RICHARD D. NEIDINGER Abstract. Techniques of univariate Newton interpolating polynomials are extended to multivariate data points by different generalizations
More informationWHAT YOU SHOULD LEARN
GRAPHS OF EQUATIONS WHAT YOU SHOULD LEARN Sketch graphs of equations. Find x- and y-intercepts of graphs of equations. Use symmetry to sketch graphs of equations. Find equations of and sketch graphs of
More information1KOd17RMoURxjn2 CSE 20 DISCRETE MATH Fall
CSE 20 https://goo.gl/forms/1o 1KOd17RMoURxjn2 DISCRETE MATH Fall 2017 http://cseweb.ucsd.edu/classes/fa17/cse20-ab/ Today's learning goals Explain the steps in a proof by mathematical and/or structural
More informationNumerical Analysis I - Final Exam Matrikelnummer:
Dr. Behrens Center for Mathematical Sciences Technische Universität München Winter Term 2005/2006 Name: Numerical Analysis I - Final Exam Matrikelnummer: I agree to the publication of the results of this
More informationComputational Mathematics: Models, Methods and Analysis. Zhilin Li
Computational Mathematics: Models, Methods and Analysis Zhilin Li Chapter 1 Introduction Why is this course important (motivations)? What is the role of this class in the problem solving process using
More informationSEQUENCES AND SERIES Sequences CC Standard
N Sequences and Series, Lesson 1, Sequences (r. 2018) SEQUENCES AND SERIES Sequences CC Standard NG Standard F-IF.A.3 Recognize that sequences are functions, sometimes defined recursively, whose domain
More information2 Computation with Floating-Point Numbers
2 Computation with Floating-Point Numbers 2.1 Floating-Point Representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However, real numbers
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationMath 3316, Fall 2016 Due Nov. 3, 2016
Math 3316, Fall 2016 Due Nov. 3, 2016 Project 3 Polynomial Interpolation The first two sections of this project will be checked in lab the week of Oct. 24-26 this completion grade will count for 10% of
More informationEasy way to Find Multivariate Interpolation
International Journal of Emerging Trends in Science Technology IC Value: 7689 (Index Copernicus) Impact Factor: 4219 DOI: https://dxdoiorg/1018535/ijetst/v4i511 Easy way to Find Multivariate Interpolation
More informationIntroduction to Matlab
Technische Universität München WT 21/11 Institut für Informatik Prof Dr H-J Bungartz Dipl-Tech Math S Schraufstetter Benjamin Peherstorfer, MSc October 22nd, 21 Introduction to Matlab Engineering Informatics
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationDOWNLOAD PDF BIG IDEAS MATH VERTICAL SHRINK OF A PARABOLA
Chapter 1 : BioMath: Transformation of Graphs Use the results in part (a) to identify the vertex of the parabola. c. Find a vertical line on your graph paper so that when you fold the paper, the left portion
More informationPrecomposing Equations
Precomposing Equations Let s precompose the function f(x) = x 3 2x + 9 with the function g(x) = 4 x. (Precompose f with g means that we ll look at f g. We would call g f postcomposing f with g.) f g(x)
More informationNumerical Analysis Fall. Numerical Differentiation
Numerical Analysis 5 Fall Numerical Differentiation Differentiation The mathematical definition of a derivative begins with a difference approimation: and as is allowed to approach zero, the difference
More informationA-SSE.1.1, A-SSE.1.2-
Putnam County Schools Curriculum Map Algebra 1 2016-2017 Module: 4 Quadratic and Exponential Functions Instructional Window: January 9-February 17 Assessment Window: February 20 March 3 MAFS Standards
More informationArithmetic in a Finite Field
MATHEMATICS OF COMPUTATION, VOLUME 35, NUMBER 152 OCTOBER 198, PAGES 1353-1359 Arithmetic in a Finite Field By Michael Willett * Abstract. An algorithm for realizing finite field arithmetic is presented.
More informationFor example, the system. 22 may be represented by the augmented matrix
Matrix Solutions to Linear Systems A matrix is a rectangular array of elements. o An array is a systematic arrangement of numbers or symbols in rows and columns. Matrices (the plural of matrix) may be
More information3.1. 3x 4y = 12 3(0) 4y = 12. 3x 4y = 12 3x 4(0) = y = x 0 = 12. 4y = 12 y = 3. 3x = 12 x = 4. The Rectangular Coordinate System
3. The Rectangular Coordinate System Interpret a line graph. Objectives Interpret a line graph. Plot ordered pairs. 3 Find ordered pairs that satisfy a given equation. 4 Graph lines. 5 Find x- and y-intercepts.
More informationMore Ways to Solve & Graph Quadratics The Square Root Property If x 2 = a and a R, then x = ± a
More Ways to Solve & Graph Quadratics The Square Root Property If x 2 = a and a R, then x = ± a Example: Solve using the square root property. a) x 2 144 = 0 b) x 2 + 144 = 0 c) (x + 1) 2 = 12 Completing
More informationPROGRAM EFFICIENCY & COMPLEXITY ANALYSIS
Lecture 03-04 PROGRAM EFFICIENCY & COMPLEXITY ANALYSIS By: Dr. Zahoor Jan 1 ALGORITHM DEFINITION A finite set of statements that guarantees an optimal solution in finite interval of time 2 GOOD ALGORITHMS?
More informationA NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME
J. Appl. Math. & Computing Vol. 18(2005), No. 1-2, pp. 553-562 A NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME YOUCHUN TANG, KWAN PYO KO* AND BYUNG-GOOK LEE Abstract. It is well-known
More informationDaubechies Wavelets and Interpolating Scaling Functions and Application on PDEs
Daubechies Wavelets and Interpolating Scaling Functions and Application on PDEs R. Schneider F. Krüger TUB - Technical University of Berlin November 22, 2007 R. Schneider F. Krüger Daubechies Wavelets
More informationSimplified Calculation of the Exponential Integral
SIMPLIFIED CALCULATION OF THE EXPONENTIAL INTEGRAL 87 Simplified Calculation of the Exponential Integral By James Miller and R. P. Hurst Calculations of the negative exponential integral function [] ()
More information5.1 Introduction to the Graphs of Polynomials
Math 3201 5.1 Introduction to the Graphs of Polynomials In Math 1201/2201, we examined three types of polynomial functions: Constant Function - horizontal line such as y = 2 Linear Function - sloped line,
More informationNumerical Aspects of Special Functions
Numerical Aspects of Special Functions Nico M. Temme In collaboration with Amparo Gil and Javier Segura, Santander, Spain. Nico.Temme@cwi.nl Centrum voor Wiskunde en Informatica (CWI), Amsterdam Numerics
More informationPolynomial Functions Graphing Investigation Unit 3 Part B Day 1. Graph 1: y = (x 1) Graph 2: y = (x 1)(x + 2) Graph 3: y =(x 1)(x + 2)(x 3)
Part I: Polynomial Functions when a = 1 Directions: Polynomial Functions Graphing Investigation Unit 3 Part B Day 1 1. For each set of factors, graph the zeros first, then use your calculator to determine
More informationComputational Physics PHYS 420
Computational Physics PHYS 420 Dr Richard H. Cyburt Assistant Professor of Physics My office: 402c in the Science Building My phone: (304) 384-6006 My email: rcyburt@concord.edu My webpage: www.concord.edu/rcyburt
More informationFloating-point representation
Lecture 3-4: Floating-point representation and arithmetic Floating-point representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However,
More informationAlgorithms of Scientific Computing
Algorithms of Scientific Computing Overview and General Remarks Michael Bader Technical University of Munich Summer 2017 Classification of the Lecture Who is Who? Students of Informatics: Informatics Bachelor
More informationSection 7.2 Characteristics of Quadratic Functions
Section 7. Characteristics of Quadratic Functions A QUADRATIC FUNCTION is a function of the form " # $ N# 1 & ;# & 0 Characteristics Include:! Three distinct terms each with its own coefficient:! An x
More informationKU Leuven Department of Computer Science
Computing with functions on domains with arbitrary shapes Daan Huybrechs Roel Matthysen Report TW 678, June 9, 2017 KU Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium)
More informationSECTION 5.1. Sequences
SECTION 5.1 Sequences Sequences Problem: count number of ancestors one has 2 parents, 4 grandparents, 8 greatgrandparents,, written in a row as 2, 4, 8, 16, 32, 64, 128, To look for pattern of the numbers,
More informationQuasilinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 16 Lecture 3 Quasilinear First-Order PDEs A first order quasilinear PDE is of the form a(x, y, z) + b(x, y, z) x y = c(x, y, z). (1) Such equations
More informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationIndex. C m (Ω), 141 L 2 (Ω) space, 143 p-th order, 17
Bibliography [1] J. Adams, P. Swarztrauber, and R. Sweet. Fishpack: Efficient Fortran subprograms for the solution of separable elliptic partial differential equations. http://www.netlib.org/fishpack/.
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationWhat is a Computer Algorithm?
What is a Computer Algorithm? A computer algorithm is a detailed step-by-step method for solving a problem by using a computer. Problem-Solving (Science and Engineering) Analysis How does it work? Breaking
More informationNumerical methods for volume preserving image registration
Numerical methods for volume preserving image registration SIAM meeting on imaging science 2004 Eldad Haber with J. Modersitzki haber@mathcs.emory.edu Emory University Volume preserving IR p.1/?? Outline
More informationAssignment 1 (concept): Solutions
CS10b Data Structures and Algorithms Due: Thursday, January 0th Assignment 1 (concept): Solutions Note, throughout Exercises 1 to 4, n denotes the input size of a problem. 1. (10%) Rank the following functions
More informationComputation of interpolatory splines via triadic subdivision
Computation of interpolatory splines via triadic subdivision Valery A. Zheludev and Amir Z. Averbuch Abstract We present an algorithm for computation of interpolatory splines of arbitrary order at triadic
More informationRAPID COMPUTATION OF THE DISCRETE FOURIER TRANSFORM*
SIAM J. ScI. COMPUT. Vol. 17, No. 4, pp. 913-919, July 1996 1996 Society for Industrial and Applied Mathematics O08 RAPID COMPUTATION OF THE DISCRETE FOURIER TRANSFORM* CHRIS ANDERSON AND MARIE DILLON
More informationMaths for Signals and Systems Linear Algebra in Engineering. Some problems by Gilbert Strang
Maths for Signals and Systems Linear Algebra in Engineering Some problems by Gilbert Strang Problems. Consider u, v, w to be non-zero vectors in R 7. These vectors span a vector space. What are the possible
More informationDIGITAL TERRAIN MODELLING. Endre Katona University of Szeged Department of Informatics
DIGITAL TERRAIN MODELLING Endre Katona University of Szeged Department of Informatics katona@inf.u-szeged.hu The problem: data sources data structures algorithms DTM = Digital Terrain Model Terrain function:
More informationNew Basis Functions and Their Applications to PDEs
Copyright c 2007 ICCES ICCES, vol.3, no.4, pp.169-175, 2007 New Basis Functions and Their Applications to PDEs Haiyan Tian 1, Sergiy Reustkiy 2 and C.S. Chen 1 Summary We introduce a new type of basis
More informationEvaluating Rational Functions
. Name Date A C T I V I T Y 7 Evaluating Rational Functions A rational function is a function that can be expressed as f(x) = h(x) g(x), where h(x) and g(x) are both polynomials and g(x) is not the zero
More informationHonors Precalculus: Solving equations and inequalities graphically and algebraically. Page 1
Solving equations and inequalities graphically and algebraically 1. Plot points on the Cartesian coordinate plane. P.1 2. Represent data graphically using scatter plots, bar graphs, & line graphs. P.1
More informationContents. Implementing the QR factorization The algebraic eigenvalue problem. Applied 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 informationACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION
Scientific Research of the Institute of Mathematics and Computer Science ACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION Ewa Węgrzyn-Skrzypczak, Tomasz Skrzypczak Institute of Mathematics, Czestochowa
More informationAlgebra 2 Semester 2 Final Exam Study Outline Semester 2 Final Exam Study Tips and Information
Algebra 2 Semester 2 Final Exam Study Outline 2013 Semester 2 Final Exam Study Tips and Information The final exam is CUMULATIVE and will include all concepts taught from Chapter 1 through Chapter 13.
More informationLong time integrations of a convective PDE on the sphere by RBF collocation
Long time integrations of a convective PDE on the sphere by RBF collocation Bengt Fornberg and Natasha Flyer University of Colorado NCAR Department of Applied Mathematics Institute for Mathematics Applied
More informationF.BF.B.3: Graphing Polynomial Functions
F.BF.B.3: Graphing Polynomial Functions 1 Given the graph of the line represented by the equation f(x) = 2x + b, if b is increased by 4 units, the graph of the new line would be shifted 4 units 1) right
More informationarxiv: v1 [math.na] 20 Jun 2014
Iterative methods for the inclusion of the inverse matrix Marko D Petković University of Niš, Faculty of Science and Mathematics Višegradska 33, 18000 Niš, Serbia arxiv:14065343v1 [mathna 20 Jun 2014 Miodrag
More informationlecture 10: B-Splines
9 lecture : -Splines -Splines: a basis for splines Throughout our discussion of standard polynomial interpolation, we viewed P n as a linear space of dimension n +, and then expressed the unique interpolating
More informationON THE NOTION OF «-CARDINALITY
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 69, Number 2, May 1978 ON THE NOTION OF «-CARDINALITY TEODOR C. PRZYMUSIÑSKI Abstract. In this paper we introduce and investigate the notion of n-cardinality,
More informationNumerical Rectification of Curves
Applied Mathematical Sciences, Vol. 8, 2014, no. 17, 823-828 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.39500 Numerical Rectification of Curves B. P. Acharya, M. Acharya and S. B.
More informationMath 213 Exam 2. Each question is followed by a space to write your answer. Please write your answer neatly in the space provided.
Math 213 Exam 2 Name: Section: Do not remove this answer page you will return the whole exam. You will be allowed two hours to complete this test. No books or notes may be used other than a onepage cheat
More informationBlock-based Thiele-like blending rational interpolation
Journal of Computational and Applied Mathematics 195 (2006) 312 325 www.elsevier.com/locate/cam Block-based Thiele-like blending rational interpolation Qian-Jin Zhao a, Jieqing Tan b, a School of Computer
More informationSUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING
Bulletin of Mathematics Vol. 06, No. 0 (20), pp.. SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING Eddy Roflin, Sisca Octarina, Putra B. J Bangun,
More informationVW 1LQH :HHNV 7KH VWXGHQW LV H[SHFWHG WR
PreAP Pre Calculus solve problems from physical situations using trigonometry, including the use of Law of Sines, Law of Cosines, and area formulas and incorporate radian measure where needed.[3e] What
More informationChemnitz Scientific Computing Preprints
Roman Unger Obstacle Description with Radial Basis Functions for Contact Problems in Elasticity CSC/09-01 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN
More informationMAE 384 Numerical Methods for Engineers
MAE 384 Numerical Methods for Engineers Instructor: Huei-Ping Huang office: ERC 359, email: hp.huang@asu.edu (Huei rhymes with way ) Tu/Th 9:00-10:15 PM WGHL 101 Textbook: Numerical Methods for Engineers
More informationIntroduction to Applied and Pre-calculus Mathematics (2008) Correlation Chart - Grade 10 Introduction to Applied and Pre-Calculus Mathematics 2009
Use words and algebraic expressions to describe the data and interrelationships in a table with rows/columns that are not related recursively (not calculated from previous data) (Applied A-1) Use words
More informationUnlabeled equivalence for matroids representable over finite fields
Unlabeled equivalence for matroids representable over finite fields November 16, 2012 S. R. Kingan Department of Mathematics Brooklyn College, City University of New York 2900 Bedford Avenue Brooklyn,
More informationI. Recursive Descriptions A phrase like to get the next term you add 2, which tells how to obtain
Mathematics 45 Describing Patterns in s Mathematics has been characterized as the science of patterns. From an early age students see patterns in mathematics, including counting by twos, threes, etc.,
More informationCAP 6736 Geometric Modeling Team Hamilton B-Spline Basis Functions. Presented by Padmavathi Siddi Vigneshkumar Suresh Jorge Medrano
CAP 6736 Geometric Modeling Team Hamilton B-Spline Basis Functions Presented by Padmavathi Siddi Vigneshkumar Suresh Jorge Medrano 1 Team members contributions Member: Padmavathi Siddi Piecewise Polynomials
More informationImage Smoothing and Segmentation by Graph Regularization
Image Smoothing and Segmentation by Graph Regularization Sébastien Bougleux 1 and Abderrahim Elmoataz 1 GREYC CNRS UMR 6072, Université de Caen Basse-Normandie ENSICAEN 6 BD du Maréchal Juin, 14050 Caen
More informationNumerical Methods for Hyperbolic and Kinetic Equations
Numerical Methods for Hyperbolic and Kinetic Equations Organizer: G. Puppo Phenomena characterized by conservation (or balance laws) of physical quantities are modelled by hyperbolic and kinetic equations.
More information