Accurate Multiple-Precision Gauss-Legendre Quadrature
|
|
- Ralph Flowers
- 5 years ago
- Views:
Transcription
1 Accurate Multiple-Precision Gauss-Legendre Quadrature Laurent Fousse 1,2 1 Université Henri-Poincaré Nancy 1 2 INRIA Rocquencourt 18th IEEE International Symposium on Computer Arithmetic
2 Outline Motivation Numerical Quadrature Certified arithmetics Our Results Gauss-Legendre Quadrature Complete example Conclusion
3 Quadrature is ubiquitous In physics (crystallography [Jézéquel Chesneaux 2004]): g(a) = + 0 [(exp(x) + exp( x)) a exp(ax) exp( ax)] dx where 0 < a < 2 and g(a) is the crystal energy. probabilities [Trefethen, Kern 2002] 2 = 1 π π π dφ 3 4 cos φ + cos 2 φ + 16ɛ 2 (parameterized random walk in Z 2 ).
4 Quadrature is hard At least you can fail spectacularly. One example with Pari/GP 2.3.0: 10 0 x 2 sin(x 3 )dx?default(realprecision,115);intnum(x=0,10,x^2*sin(x^3)) %2 = [...]
5 Quadrature is hard At least you can fail spectacularly. One example with Pari/GP 2.3.0: 10 0 x 2 sin(x 3 )dx?default(realprecision,115);intnum(x=0,10,x^2*sin(x^3)) %2 = [...]?default(realprecision,116);intnum(x=0,10,x^2*sin(x^3)) %2 = [...] Which result do you trust?
6 Quadrature is hard At least you can fail spectacularly. One example with Pari/GP 2.3.0: 10 0 x 2 sin(x 3 )dx?default(realprecision,115);intnum(x=0,10,x^2*sin(x^3)) %2 = [...]?default(realprecision,116);intnum(x=0,10,x^2*sin(x^3)) %2 = [...] Which result do you trust? 10 0 x 2 sin(x 3 )dx = 1 (1 cos(1000)) 3 0,
7 Quadrature is hard A famous computer algebra software: \^/ Maple 10 (IBM INTEL LINUX)._ \ / _. Copyright (c) Maplesoft, a division of Waterloo Maple Inc \ MAPLE / All rights reserved. Maple is a trademark of < > Waterloo Maple Inc. Type? for help. > evalf(int(exp(-x^2)*ln(x), x=17..42));
8 Quadrature is hard A famous computer algebra software: \^/ Maple 10 (IBM INTEL LINUX)._ \ / _. Copyright (c) Maplesoft, a division of Waterloo Maple Inc \ MAPLE / All rights reserved. Maple is a trademark of < > Waterloo Maple Inc. Type? for help. > evalf(int(exp(-x^2)*ln(x), x=17..42)); > Digits:=20: evalf(int(exp(-x^2)*ln(x), x=17..42));
9 Quadrature is hard A famous computer algebra software: \^/ Maple 10 (IBM INTEL LINUX)._ \ / _. Copyright (c) Maplesoft, a division of Waterloo Maple Inc \ MAPLE / All rights reserved. Maple is a trademark of < > Waterloo Maple Inc. Type? for help. > evalf(int(exp(-x^2)*ln(x), x=17..42)); > Digits:=20: evalf(int(exp(-x^2)*ln(x), x=17..42)); > Digits:=50: evalf(int(exp(-x^2)*ln(x), x=17..42));
10 Quadrature is hard A famous computer algebra software: \^/ Maple 10 (IBM INTEL LINUX)._ \ / _. Copyright (c) Maplesoft, a division of Waterloo Maple Inc \ MAPLE / All rights reserved. Maple is a trademark of < > Waterloo Maple Inc. Type? for help. > evalf(int(exp(-x^2)*ln(x), x=17..42)); > Digits:=20: evalf(int(exp(-x^2)*ln(x), x=17..42)); > Digits:=50: evalf(int(exp(-x^2)*ln(x), x=17..42)); > Digits:=100: evalf(int(exp(-x^2)*ln(x), x=17..42)); [...] 10
11 Quadrature is hard exp( x 2 ) log xdx 0, Precision Maple answer Correct digits 10 0, , , , , , , , 7
12 Miscomputations are bad Ariane 5
13 Miscomputations are bad Ariane 5 Patriot anti-missile
14 Miscomputations are bad Ariane 5 Patriot anti-missile Vancouver Stock Exchange
15 Miscomputations are bad Every time you miscompute a kitten gets killed. Please, think of the kitten! = mimic the IEEE 754 standard as closely as possible for numerical integration.
16 Gauss-Legendre Quadrature in short rule of interpolatory type: b a f (x)dx b n 1 a x xj n 1 f (x i ) = w i f (x i ) xi x j i=0 the x i s are the roots of the n-th Legendre polynomial P n the weights w i are easily computed w i = mathematical error bounded by 2 (1 x 2 i )P 2 n (x i ). (b a) 2n+1 (n!) 4 (2n + 1)[(2n)!] 3 f (2n). i=0
17 Gauss-Legendre Quadrature in CRQ Correctly Rounded Quadrature Compute P n via (n + 1)P n+1 (X) = (2n + 1)XP n (X) np n 1 (X) Root isolation (Uspensky) and refinement (interval Newton) Forward error analysis to derive bounds and their proof
18 Hypothesis y f(x) b a f(x)dx a b x f is a C black-box function in arbitrary precision. bounds M k on f (k) are known.
19 Quadrature Algorithm INPUTS: â, b a, (ŵ i ), f, ( v i ), n where v i = ( 1+x i 2 ). OUTPUT: Î b a f (x)dx 1: for i 0 to n 1 do 2: t (( b a) v i ) 3: x i (t + â) 4: fi (f ( x i )) 5: ŷ i ( f i ŵ i ) 6: end for 7: Ŝ sum(ŷ i, i = 0... n 1) 8: D ( b a)/2 9: return ( DŜ) = Î
20 Theorem Theorem Let δ byi = 11 4 ulp(ŷ i ) + 6M 1 ŵ i ulp( x i ), p 2 the working precision. With n integration points the error Î I is bounded by: B total = 21 4 ulp(î) + 5n 4 D max(δ byi ) + (b a)2n+1 (n!) 4 (2n + 1)[(2n)!] 3 M 2n.
21 A glimpse of the proof Let v i = 1+x i 2 Then v i v i 1 2 ulp( v i ) x i = ( ( v i ( b a)) + â). ( v i b a) v i (b a) 5 2 ulp( ( v i b a))
22 A glimpse of the proof Let v i = 1+x i 2 Then v i v i 1 2 ulp( v i ) x i = ( ( v i ( b a)) + â). x i x i 1 2 ulp( x i ) ulp( ( v i b a)) ulp(â) 17 4 ulp( x i ).
23 Let v i = 1+x i 2 Then A glimpse of the proof v i v i 1 2 ulp( v i ) x i = ( ( v i ( b a)) + â). x i x i 1 2 ulp( x i ) ulp( ( v i b a)) ulp(â) 17 4 ulp( x i ). (f ( ˆx i )) f (x i ) (f ( ˆx i )) f ( ˆx i ) + f ( ˆx i ) f (x i )
24 A glimpse of the proof Let v i = 1+x i 2 x i x i 1 2 ulp( x i ) ulp( ( v i b a)) ulp(â) 17 4 ulp( x i ). (f ( ˆx i )) f (x i ) (f ( ˆx i )) f ( ˆx i ) + f ( ˆx i ) f (x i ) ulp( (f ( ˆx i ))) + M 1 ˆx i x i
25 Let v i = 1+x i 2 A glimpse of the proof x i x i 1 2 ulp( x i ) ulp( ( v i b a)) ulp(â) 17 4 ulp( x i ). (f ( ˆx i )) f (x i ) (f ( ˆx i )) f ( ˆx i ) + f ( ˆx i ) f (x i ) ulp( (f ( ˆx i ))) + M 1 ˆx i x i ulp( (f ( ˆx i ))) M 1ulp( ˆx i )
26 Our favorite example I = g(x) = e x 2 i=0 e x 2 log xdx h(x) = log x n ( ) n d f (n) i dn i (x) = g(x) h(x) i dx i dx n i After some computations: f (n) n n!e 172 ( (n + 1)42 n log 42 + (n 1)42 n 2).
27 B math vs. B rounding (p = 1000) math rounding log(error) n
28 Measured error vs. computed error bound (p = 1000) -200 measured total log(error) n
29 Quadrature coefficients computation time 3500 Computation time Time (ms) n
30 CRQ vs. Rest of the World I = 1 0 max(sin(x), cos(x))dx = 2 cos 1 Precision CRQ Time 26ms 110ms 1s 11s 108s 3138s Error 1 2.6e-2 9.2e Maple Time 6s 38s 371s 2618s Error 2.1e-1 1.7e-1 1.7e-1 2.9e-1 PARI/GP Time 20ms 80ms 480ms 3s 22s 301s Error 3.1e25 2.0e54 2.7e e e e1496
31 Summary Multiple-precision quadrature can be fast...
32 Summary Multiple-precision quadrature can be fast... and certified at the same time.
33 Summary Multiple-precision quadrature can be fast... and certified at the same time. Code available at
34 Summary Multiple-precision quadrature can be fast... and certified at the same time. Code available at TODOs: Speed-up Gauss-Legendre coefficients computation (parallel Newton iteration?), Look at other methods in the Gauss family, Automatic computation of bounds for specific functions, Optimal composition strategy.
Accurate Multiple-Precision Gauss-Legendre Quadrature
Accurate Multiple-Precision Gauss-Legendre Quadrature Laurent Fousse Université Henri-Poincaré Nancy 1 laurent@komite.net Abstract Numerical integration is an operation that is frequently available in
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationNotes on floating point number, numerical computations and pitfalls
Notes on floating point number, numerical computations and pitfalls November 6, 212 1 Floating point numbers An n-digit floating point number in base β has the form x = ±(.d 1 d 2 d n ) β β e where.d 1
More informationElements of Floating-point Arithmetic
Elements of Floating-point Arithmetic Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Floating-point Numbers Representations IEEE Floating-point Standards Underflow
More informationElements of Floating-point Arithmetic
Elements of Floating-point Arithmetic Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Floating-point Numbers Representations IEEE Floating-point Standards Underflow
More informationComputation of the error functions erf and erfc in arbitrary precision with correct rounding
Computation of the error functions erf and erfc in arbitrary precision with correct rounding Sylvain Chevillard Arenaire, LIP, ENS-Lyon, France Sylvain.Chevillard@ens-lyon.fr Nathalie Revol INRIA, Arenaire,
More informationODE Background: Differential (1A) Young Won Lim 12/29/15
ODE Background: Differential (1A Copyright (c 2011-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationDeterministic Global Optimization for Dynamic Systems Using Interval Analysis
Deterministic Global Optimization for Dynamic Systems Using Interval Analysis Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame,
More informationACM 106a: Lecture 1 Agenda
1 ACM 106a: Lecture 1 Agenda Introduction to numerical linear algebra Common problems First examples Inexact computation What is this course about? 2 Typical numerical linear algebra problems Systems of
More informationCan You Count on Your Computer?
Can You Count on Your Computer? Professor Nick Higham School of Mathematics University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ p. 1/33 p. 2/33 Counting to Six I asked my computer
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
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 informationMath Numerical Analysis
Math 541 - Numerical Analysis Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State University
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationINTRODUCTION TO COMPUTATIONAL MATHEMATICS
INTRODUCTION TO COMPUTATIONAL MATHEMATICS Course Notes for CM 271 / AMATH 341 / CS 371 Fall 2007 Instructor: Prof. Justin Wan School of Computer Science University of Waterloo Course notes by Prof. Hans
More informationIntroduction to Scientific Computing Languages
1 / 21 Introduction to Scientific Computing Languages Prof. Paolo Bientinesi pauldj@aices.rwth-aachen.de Numerical Representation 2 / 21 Numbers 123 = (first 40 digits) 29 4.241379310344827586206896551724137931034...
More informationRound-off Error Analysis of Explicit One-Step Numerical Integration Methods
Round-off Error Analysis of Explicit One-Step Numerical Integration Methods 24th IEEE Symposium on Computer Arithmetic Sylvie Boldo 1 Florian Faissole 1 Alexandre Chapoutot 2 1 Inria - LRI, Univ. Paris-Sud
More informationIntegration. Antiderivatives and Indefinite Integration 3/9/2015. Copyright Cengage Learning. All rights reserved.
Integration Copyright Cengage Learning. All rights reserved. Antiderivatives and Indefinite Integration Copyright Cengage Learning. All rights reserved. 1 Objectives Write the general solution of a differential
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationHigher order derivatives of the inverse function
Higher order derivatives of the inverse function Andrej Liptaj Abstract A general recursive and limit formula for higher order derivatives of the inverse function is presented. The formula is next used
More information1 Finding and fixing floating point problems
Notes for 2016-09-09 1 Finding and fixing floating point problems Floating point arithmetic is not the same as real arithmetic. Even simple properties like associativity or distributivity of addition and
More informationIntroduction to Scientific Computing Languages
1 / 19 Introduction to Scientific Computing Languages Prof. Paolo Bientinesi pauldj@aices.rwth-aachen.de Numerical Representation 2 / 19 Numbers 123 = (first 40 digits) 29 4.241379310344827586206896551724137931034...
More informationRigorous Polynomial Approximations and Applications
Rigorous Polynomial Approximations and Applications Mioara Joldeș under the supervision of: Nicolas Brisebarre and Jean-Michel Muller École Normale Supérieure de Lyon, Arénaire Team, Laboratoire de l Informatique
More informationDRAFT - Math 102 Lecture Note - Dr. Said Algarni
Math02 - Term72 - Guides and Exercises - DRAFT 7 Techniques of Integration A summery for the most important integrals that we have learned so far: 7. Integration by Parts The Product Rule states that if
More informationPower Series Solutions to the Legendre Equation
Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre s equation. When α Z +, the equation has polynomial
More informationNUMERICAL ANALYSIS PROBLEMS
NUMERICAL ANALYSIS PROBLEMS JAMES KEESLING The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.. Solving Equations Problem.
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationThe problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.
NUMERICAL ANALYSIS PRACTICE PROBLEMS JAMES KEESLING The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.. Solving Equations
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationNumerical Integration with the double exponential method
Numerical Integration with the double exponential method Pascal Molin Université Paris 7 December 7 2016 Numerical integration Number-theory context, want fast quadrature high precision (100 1000 digits)
More informationFundamentals of Numerical Mathematics. Charles University Prague Faculty of Mathematics and Physics Czech Republic
Fundamentals of Numerical Mathematics Vít Dolejší Charles University Prague Faculty of Mathematics and Physics Czech Republic dolejsi@karlin.mff.cuni.cz December 18, 2017 Preface These lecture notes more
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More information2-4 Zeros of Polynomial Functions
Write a polynomial function of least degree with real coefficients in standard form that has the given zeros. 33. 2, 4, 3, 5 Using the Linear Factorization Theorem and the zeros 2, 4, 3, and 5, write f
More information13. Nonlinear least squares
L. Vandenberghe ECE133A (Fall 2018) 13. Nonlinear least squares definition and examples derivatives and optimality condition Gauss Newton method Levenberg Marquardt method 13.1 Nonlinear least squares
More informationMATH1190 CALCULUS 1 - NOTES AND AFTERNOTES
MATH90 CALCULUS - NOTES AND AFTERNOTES DR. JOSIP DERADO. Historical background Newton approach - from physics to calculus. Instantaneous velocity. Leibniz approach - from geometry to calculus Calculus
More informationERROR BOUNDS ON COMPLEX FLOATING-POINT MULTIPLICATION
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 005-5718(XX)0000-0 ERROR BOUNDS ON COMPLEX FLOATING-POINT MULTIPLICATION RICHARD BRENT, COLIN PERCIVAL, AND PAUL ZIMMERMANN In memory of
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationIntroduction to ODE's (0A) Young Won Lim 3/9/15
Introduction to ODE's (0A) Copyright (c) 2011-2014 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationExercises * on Linear Algebra
Exercises * on Linear Algebra Laurenz Wiskott Institut für Neuroinformatik Ruhr-Universität Bochum, Germany, EU 4 February 7 Contents Vector spaces 4. Definition...............................................
More informationNumerical integration in arbitrary-precision ball arithmetic
/ 24 Numerical integration in arbitrary-precision ball arithmetic Fredrik Johansson (LFANT, Bordeaux) Journées FastRelax, Inria Sophia Antipolis 7 June 208 Numerical integration in Arb 2 / 24 New code
More informationBinary floating point
Binary floating point Notes for 2017-02-03 Why do we study conditioning of problems? One reason is that we may have input data contaminated by noise, resulting in a bad solution even if the intermediate
More informationClassroom Tips and Techniques: Electric Field from Distributed Charge
Classroom Tips and Techniques: Electric Field from Distributed Charge Introduction Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft This past summer I was asked if Maple could
More informationOn Newton-type methods with cubic convergence
Journal of Computational and Applied Mathematics 176 (2005) 425 432 www.elsevier.com/locate/cam On Newton-type methods with cubic convergence H.H.H. Homeier a,b, a Science + Computing Ag, IT Services Muenchen,
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 informationRomberg Integration and Gaussian Quadrature
Romberg Integration and Gaussian Quadrature P. Sam Johnson October 17, 014 P. Sam Johnson (NITK) Romberg Integration and Gaussian Quadrature October 17, 014 1 / 19 Overview We discuss two methods for integration.
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechanics for Scientists and Engineers David Miller Vector spaces, operators and matrices Vector spaces, operators and matrices Vector space Vector space We need a space in which our vectors exist
More informationStudy 4.10 #465, 471, , 487, , , 515, 517, 521, 523
Goals: 1. Understand that antiderivatives are the functions from which the present derivative was found. 2. The process of finding an antiderivative or indefinite integral requires the reverse process
More informationComputing the Monodromy Group of a Plane Algebraic Curve Using a New Numerical-modular Newton-Puiseux Algorithm
Computing the Monodromy Group of a Plane Algebraic Curve Using a New Numerical-modular Newton-Puiseux Algorithm Poteaux Adrien XLIM-DMI UMR CNRS 6172 Université de Limoges, France SNC'07 University of
More informationCompute the behavior of reality even if it is impossible to observe the processes (for example a black hole in astrophysics).
1 Introduction Read sections 1.1, 1.2.1 1.2.4, 1.2.6, 1.3.8, 1.3.9, 1.4. Review questions 1.1 1.6, 1.12 1.21, 1.37. The subject of Scientific Computing is to simulate the reality. Simulation is the representation
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationAutomating the Verification of Floating-point Algorithms
Introduction Interval+Error Advanced Gappa Conclusion Automating the Verification of Floating-point Algorithms Guillaume Melquiond Inria Saclay Île-de-France LRI, Université Paris Sud, CNRS 2014-07-18
More information1 Backward and Forward Error
Math 515 Fall, 2008 Brief Notes on Conditioning, Stability and Finite Precision Arithmetic Most books on numerical analysis, numerical linear algebra, and matrix computations have a lot of material covering
More informationƒ f(x)dx ~ X) ^i,nf(%i,n) -1 *=1 are the zeros of P«(#) and where the num
ZEROS OF THE HERMITE POLYNOMIALS AND WEIGHTS FOR GAUSS' MECHANICAL QUADRATURE FORMULA ROBERT E. GREENWOOD AND J. J. MILLER In the numerical integration of a function ƒ (x) it is very desirable to choose
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationComparison of Clenshaw-Curtis and Gauss Quadrature
WDS'11 Proceedings of Contributed Papers, Part I, 67 71, 211. ISB 978-8-7378-184-2 MATFYZPRESS Comparison of Clenshaw-Curtis and Gauss Quadrature M. ovelinková Charles University, Faculty of Mathematics
More informationAutomatic determination of numerical properties of software and systems
Automatic determination of numerical properties of software and systems Eric Goubault and Sylvie Putot Modelling and Analysis of Interacting Systems, CEA LIST MASCOT-NUM 2012 Meeting, March 21-23, 2012
More informationSection 0.2 & 0.3 Worksheet. Types of Functions
MATH 1142 NAME Section 0.2 & 0.3 Worksheet Types of Functions Now that we have discussed what functions are and some of their characteristics, we will explore different types of functions. Section 0.2
More informationCore Mathematics 3 Algebra
http://kumarmathsweeblycom/ Core Mathematics 3 Algebra Edited by K V Kumaran Core Maths 3 Algebra Page Algebra fractions C3 The specifications suggest that you should be able to do the following: Simplify
More informationChange of Variables: Indefinite Integrals
Change of Variables: Indefinite Integrals Mathematics 11: Lecture 39 Dan Sloughter Furman University November 29, 2007 Dan Sloughter (Furman University) Change of Variables: Indefinite Integrals November
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationSection 4.4. Using the Fundamental Theorem. Difference Equations to Differential Equations
Difference Equations to Differential Equations Section 4.4 Using the Fundamental Theorem As we saw in Section 4.3, using the Fundamental Theorem of Integral Calculus reduces the problem of evaluating a
More informationNumerical integration - I. M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V
Numerical integration - I M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V deterministic methods in 1D equispaced points (trapezoidal, Simpson...), others...
More informationStep 1. Step 2. Step 4. The corrected trial solution y with evaluated coefficients d 1, d 2,..., d k becomes the particular solution y p.
The Basic Trial Solution Method Outlined here is the method for a second order differential equation ay + by + cy = f(x). The method applies unchanged for nth order equations. Step 1. Step 2. Repeatedly
More informationLecture 28 The Main Sources of Error
Lecture 28 The Main Sources of Error Truncation Error Truncation error is defined as the error caused directly by an approximation method For instance, all numerical integration methods are approximations
More informationMATH 2400: Calculus III, Fall 2013 FINAL EXAM
MATH 2400: Calculus III, Fall 2013 FINAL EXAM December 16, 2013 YOUR NAME: Circle Your Section 001 E. Angel...................... (9am) 002 E. Angel..................... (10am) 003 A. Nita.......................
More informationNeural Networks in Structured Prediction. November 17, 2015
Neural Networks in Structured Prediction November 17, 2015 HWs and Paper Last homework is going to be posted soon Neural net NER tagging model This is a new structured model Paper - Thursday after Thanksgiving
More informationMAT 460: Numerical Analysis I. James V. Lambers
MAT 460: Numerical Analysis I James V. Lambers January 31, 2013 2 Contents 1 Mathematical Preliminaries and Error Analysis 7 1.1 Introduction............................ 7 1.1.1 Error Analysis......................
More informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
More informationNonlinear Equations and Continuous Optimization
Nonlinear Equations and Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationCalculus I Announcements
Slide 1 Calculus I Announcements Read sections 4.2,4.3,4.4,4.1 and 5.3 Do the homework from sections 4.2,4.3,4.4,4.1 and 5.3 Exam 3 is Thursday, November 12th See inside for a possible exam question. Slide
More informationMath 480 The Vector Space of Differentiable Functions
Math 480 The Vector Space of Differentiable Functions The vector space of differentiable functions. Let C (R) denote the set of all infinitely differentiable functions f : R R. Then C (R) is a vector space,
More informationMath 155 Prerequisite Review Handout
Math 155 Prerequisite Review Handout August 23, 2010 Contents 1 Basic Mathematical Operations 2 1.1 Examples...................................... 2 1.2 Exercises.......................................
More informationRound-off error propagation and non-determinism in parallel applications
Round-off error propagation and non-determinism in parallel applications Vincent Baudoui (Argonne/Total SA) vincent.baudoui@gmail.com Franck Cappello (Argonne/INRIA/UIUC-NCSA) Georges Oppenheim (Paris-Sud
More informationMathematics Notes for Class 12 chapter 7. Integrals
1 P a g e Mathematics Notes for Class 12 chapter 7. Integrals Let f(x) be a function. Then, the collection of all its primitives is called the indefinite integral of f(x) and is denoted by f(x)dx. Integration
More informationOutline. 1 Numerical Integration. 2 Numerical Differentiation. 3 Richardson Extrapolation
Outline Numerical Integration Numerical Differentiation Numerical Integration Numerical Differentiation 3 Michael T. Heath Scientific Computing / 6 Main Ideas Quadrature based on polynomial interpolation:
More informationfunction independent dependent domain range graph of the function The Vertical Line Test
Functions A quantity y is a function of another quantity x if there is some rule (an algebraic equation, a graph, a table, or as an English description) by which a unique value is assigned to y by a corresponding
More informationExperience in Factoring Large Integers Using Quadratic Sieve
Experience in Factoring Large Integers Using Quadratic Sieve D. J. Guan Department of Computer Science, National Sun Yat-Sen University, Kaohsiung, Taiwan 80424 guan@cse.nsysu.edu.tw April 19, 2005 Abstract
More informationTest 2 Review Math 1111 College Algebra
Test 2 Review Math 1111 College Algebra 1. Begin by graphing the standard quadratic function f(x) = x 2. Then use transformations of this graph to graph the given function. g(x) = x 2 + 2 *a. b. c. d.
More informationCh 7 Summary - POLYNOMIAL FUNCTIONS
Ch 7 Summary - POLYNOMIAL FUNCTIONS 1. An open-top box is to be made by cutting congruent squares of side length x from the corners of a 8.5- by 11-inch sheet of cardboard and bending up the sides. a)
More informationIsolating critical cases for reciprocals using integer factorization. John Harrison Intel Corporation ARITH-16 Santiago de Compostela 17th June 2003
0 Isolating critical cases for reciprocals using integer factorization John Harrison Intel Corporation ARITH-16 Santiago de Compostela 17th June 2003 1 result before the final rounding Background Suppose
More information1 Strict local optimality in unconstrained optimization
ORF 53 Lecture 14 Spring 016, Princeton University Instructor: A.A. Ahmadi Scribe: G. Hall Thursday, April 14, 016 When in doubt on the accuracy of these notes, please cross check with the instructor s
More informationPrinciples of Scientific Computing Local Analysis
Principles of Scientific Computing Local Analysis David Bindel and Jonathan Goodman last revised January 2009, printed February 25, 2009 1 Among the most common computational tasks are differentiation,
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationSolving Triangular Systems More Accurately and Efficiently
17th IMACS World Congress, July 11-15, 2005, Paris, France Scientific Computation, Applied Mathematics and Simulation Solving Triangular Systems More Accurately and Efficiently Philippe LANGLOIS and Nicolas
More informationI II III IV V VI VII VIII IX Total
DEPARTMENT OF MATHEMATICS AND STATISTICS QUEEN S UNIVERSITY AT KINGSTON MATH 121 - DEC 2014 CDS/Section 700 Students ONLY INSTRUCTIONS: Answer all questions, writing clearly in the space provided. If you
More informationModule III: Partial differential equations and optimization
Module III: Partial differential equations and optimization Martin Berggren Department of Information Technology Uppsala University Optimization for differential equations Content Martin Berggren (UU)
More informationMathematical Models with Maple
Algebraic Biology 005 151 Mathematical Models with Maple Tetsu YAMAGUCHI Applied System nd Division, Cybernet Systems Co., Ltd., Otsuka -9-3, Bunkyo-ku, Tokyo 11-001, Japan tetsuy@cybernet.co.jp Abstract
More informationNonlinear Stationary Subdivision
Nonlinear Stationary Subdivision Michael S. Floater SINTEF P. O. Box 4 Blindern, 034 Oslo, Norway E-mail: michael.floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research Center
More informationAccelerating Convergence
Accelerating Convergence MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation We have seen that most fixed-point methods for root finding converge only linearly
More informationStep 1. Step 2. Step 4. The corrected trial solution y with evaluated coefficients d 1, d 2,..., d k becomes the particular solution y p.
Definition Atoms A and B are related if and only if their successive derivatives share a common atom. Then x 3 is related to x and x 101, while x is unrelated to e x, xe x and x sin x. Atoms x sin x and
More informationExercises given in lecture on the day in parantheses.
A.Miller M22 Fall 23 Exercises given in lecture on the day in parantheses. The ɛ δ game. lim x a f(x) = L iff Hero has a winning strategy in the following game: Devil plays: ɛ > Hero plays: δ > Devil plays:
More information