Round-off Error Analysis of Explicit One-Step Numerical Integration Methods
|
|
- Ashlynn Russell
- 5 years ago
- Views:
Transcription
1 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 et CNRS - Univ. Paris-Saclay 2 U2IS, ÉNSTA ParisTech we thank the IEEE for registration bursary
2 2 / 34 Table of contents 1 Motivations and numerical methods 2 Roundoff errors of RK methods Local roundoff errors Global roundoff errors of classical methods 3 Conclusion and perspectives
3 Ordinary differential equations (ODEs) y (t) = f(y, t). Exact resolution is hard numerical methods. 3 / 34
4 4 / 34 Numerical integration
5 5 / 34 Numerical integration
6 6 / 34 Numerical integration
7 7 / 34 Numerical integration
8 8 / 34 Numerical integration
9 9 / 34 Numerical integration
10 10 / 34 Numerical integration
11 11 / 34 Numerical integration
12 12 / 34 Euler method k 1 = hf(t n, y n ) y n+1 = y n + k 1 + O(h 2 )
13 13 / 34 RK2 method k 1 = h f(t n, y n ) k 2 = h f(t n + h 2, y n + k 1 2 ) y n+1 = y n + k 2 + O(h 3 )
14 14 / 34 Working assumptions FP arithmetic; neither underflow nor overflow, radix 2 double precision, u = 2 53
15 14 / 34 Working assumptions FP arithmetic; neither underflow nor overflow, radix 2 double precision, ODEs; first-order, linear, y R R, u = 2 53 y = λy
16 14 / 34 Working assumptions FP arithmetic; neither underflow nor overflow, radix 2 double precision, ODEs; first-order, linear, y R R, Methods. explicit, one step, constant step. u = 2 53 y = λy
17 15 / 34 RK methods on linear problems: linear stability { y 0 R y n+1 = R(h, λ)y n (R polynomial in hλ)
18 15 / 34 RK methods on linear problems: linear stability { y 0 R y n+1 = R(h, λ)y n (R polynomial in hλ) Stable R(h, λ) < 1:
19 16 / 34 RK methods on linear problems: FP implementation { y 0 R y n+1 = R(h, λ)y n { ỹ0 y 0 ỹ n+1 = R( h, λ, ỹ n )
20 16 / 34 RK methods on linear problems: FP implementation { y 0 R y n+1 = R(h, λ)y n { ỹ0 y 0 ỹ n+1 = R( h, λ, ỹ n ) Euler: R(h, λ) = 1 + hλ; R( h, λ, ỹ n ) = ỹ n h λ ỹ n.
21 RK methods on linear problems: FP implementation { y 0 R y n+1 = R(h, λ)y n { ỹ0 y 0 ỹ n+1 = R( h, λ, ỹ n ) Euler: R(h, λ) = 1 + hλ; R( h, λ, ỹ n ) = ỹ n h λ ỹ n. RK4: R(h, λ) = 1 + hλ (hλ) (hλ) (hλ)4 ; R( h, λ, ỹ n ) = ỹ n h 6 λ ỹ n h 3 λ ỹ n h h 6 λ λ ỹ n h 3 λ ỹ n h h 6 λ λ ỹn h h h 12 λ λ λ ỹ n h 6 λ ỹ n h h 6 λ λ ỹ n h h h 12 λ λ λ ỹ n h h h h 24 λ λ λ λ ỹ n. (> 60 flops!) 16 / 34
22 17 / 34 State-of-the-art Roundoff errors in numerical methods (N = nb of iterations): probabilistic result: error in N [Henrici,1963]; in practice (implicit RK): error in N [Hairer&al, 2008]; interval analysis [Bouissou-Martel, 2006]; numerical integration (fine-grained): Newton-Cotes, Gauss-Legendre,... [Fousse, 2006].
23 State-of-the-art Roundoff errors in numerical methods (N = nb of iterations): probabilistic result: error in N [Henrici,1963]; in practice (implicit RK): error in N [Hairer&al, 2008]; interval analysis [Bouissou-Martel, 2006]; numerical integration (fine-grained): Newton-Cotes, Gauss-Legendre,... [Fousse, 2006]. Our approach: fined-grained analysis; use of mathematical properties of the methods (stability). 17 / 34
24 18 / 34 Table of contents 1 Motivations and numerical methods 2 Roundoff errors of RK methods Local roundoff errors Global roundoff errors of classical methods 3 Conclusion and perspectives
25 19 / 34 Method error vs roundoff error
26 20 / 34 Method error vs roundoff error
27 21 / 34 Method error vs roundoff error
28 Local roundoff error vs global roundoff error Local error: ε 0 = ỹ 0 y 0 n N, ε n = R( h, λ, ỹ n 1 ) R(h, λ)ỹ n 1. Global error: n N, E n = ỹ n y n. 22 / 34
29 23 / 34 From local to global roundoff error Local error: ε 0 = ỹ 0 y 0 n N, ε n = R( h, λ, ỹ n 1 ) R(h, λ)ỹ n 1. Global error: n N, E n = ỹ n y n. Theorem 1: Global absolute error of RK methods Let C R +. Suppose n N, ε n C ỹ n 1. Then, n N, E n (C + R(h, λ) ) n C y 0 (ε 0 + n C + R(h, λ) ).
30 24 / 34 Relative roundoff errors Relative error: n y n C + R(h, λ) C y 0 ỹn ( ) (ε 0 + n y n R(h, λ) C + R(h, λ) ).
31 24 / 34 Relative roundoff errors Relative error: n y n C + R(h, λ) C y 0 ỹn ( ) (ε 0 + n y n R(h, λ) C + R(h, λ) ). If C R(h, λ), then: ỹn y n y n ε 0 + n C y 0 R(h, λ). In practice (Euler, RK2, RK4): C 200u and 200u R(h, λ).
32 25 / 34 Table of contents 1 Motivations and numerical methods 2 Roundoff errors of RK methods Local roundoff errors Global roundoff errors of classical methods 3 Conclusion and perspectives
33 26 / 34 Technical lemma for local roundoff errors Local error of stable Euler s method ( 2 hλ < 0): ε n+1 = ỹ n ( h λ ỹ n ) (1 + hλ)y n.
34 27 / 34 Technical lemma for local roundoff errors Local error of stable Euler s method ( 2 hλ < 0): ε n+1 = ỹ n X1 ( h λ X2 ỹ n y ) ( 1 α1 + hλ α2 ) ỹ n y. Lemma 2 Let y R. Let C 1, C 2 R +. Let α 1, α 2 R. Let X 1, X 2 F s.t. X 1 α 1 y C 1 y and X 2 α 2 C 2. Then: X 1 (X 2 y) (α 1 + α 2 )y y (C 1 + C 2 + u ( α α 2 + C 1 + 2C 2 ) +u 2 (C 2 + α 2 )).
35 27 / 34 Technical lemma for local roundoff errors Local error of stable Euler s method ( 2 hλ < 0): ε n+1 = ỹ n X1 ( h λ X2 ỹ n y ) ( 1 α1 + hλ α2 ) ỹ n y. Lemma 2 Let y R. Let C 1, C 2 R +. Let α 1, α 2 R. Let X 1, X 2 F s.t. X 1 α 1 y C 1 y and X 2 α 2 C 2. Then: X 1 (X 2 y) (α 1 + α 2 )y y (C 1 + C 2 + u ( α α 2 + C 1 + 2C 2 ) +u 2 (C 2 + α 2 )). ỹ n X1 1 α1 ỹ n y 0 C1 ỹ n, h λ hλ 6u y X2 α2 C2 (Gappa [Melquiond]).
36 28 / 34 Local roundoff errors of higher-order methods Stable RK4 method: C ỹ n 0 [ h 6 λ] 1 2u (Gappa) [ỹ n + h 1 6 λỹ n ] = ỹ n h 6 λ ỹ n 4u (Lemma 2) [ h 1 3 λ] 4u (Gappa) [ h h 1 6 λ λ] 12u (Gappa) [ h h h 1 12 λ λ λ] 28u (Gappa) [ h h h h 1 24 λ λ λ λ] 53u (Gappa) (7 Lemma 2) [ỹ n + + h h h h 1 24 λ λ λ λỹ n ] 194u (Lemma 2)
37 29 / 34 Local roundoff errors of classical methods Lemma 3: Local error of Euler method Suppose 2 hλ and 2 60 h 1. Then: n N, ε Euler n u ỹ n. Lemma 4: Local error of RK2 method Suppose 2 hλ and 2 60 h 1. Then: n N, ε RK2 n u ỹ n. Lemma 5: Local error of RK4 method Suppose 3 hλ and 2 60 h 1. Then: n N, ε RK4 n+1 194u ỹ n.
38 30 / 34 Table of contents 1 Motivations and numerical methods 2 Roundoff errors of RK methods Local roundoff errors Global roundoff errors of classical methods 3 Conclusion and perspectives
39 31 / 34 Global roundoff errors of classical methods Bounds on global errors: Euler: C = 11.01u E n (11.01u + R(h, λ) ) n 11.01u y 0 (ε 0 + n 11.01u + R(h, λ) ) ; RK2: C = 28.01u E n (28.01u + R(h, λ) ) n 28.01u y 0 (ε 0 + n 28.01u + R(h, λ) ) ; RK4: C = 194u E n (194u + R(h, λ) ) n 194u y 0 (ε 0 + n 194u + R(h, λ) ). no compensation but reasonable bounds.
40 32 / 34 Table of contents 1 Motivations and numerical methods 2 Roundoff errors of RK methods Local roundoff errors Global roundoff errors of classical methods 3 Conclusion and perspectives
41 33 / 34 General methodology y n+1 = M i=0 α i y n ỹ n+1 = M i=0 α i ỹ n (α i = hλ, hλ 3, hλ 6, h4 λ )
42 33 / 34 General methodology y n+1 = M i=0 α i y n ỹ n+1 = M i=0 α i ỹ n (α i = hλ, hλ 3, hλ 6, h4 λ ) Methodology to bound roundoff errors:
43 33 / 34 General methodology y n+1 = M i=0 α i y n ỹ n+1 = M i=0 α i ỹ n (α i = hλ, hλ 3, hλ 6, h4 λ ) Methodology to bound roundoff errors: 1) bound ε 0 = ỹ 0 y 0 ;
44 33 / 34 General methodology y n+1 = M i=0 α i y n ỹ n+1 = M i=0 α i ỹ n (α i = hλ, hλ 3, hλ 6, h4 λ ) Methodology to bound roundoff errors: 1) bound ε 0 = ỹ 0 y 0 ; 2) bound the error on each term α i (Gappa + stability);
45 33 / 34 General methodology y n+1 = M i=0 α i y n ỹ n+1 = M i=0 α i ỹ n (α i = hλ, hλ 3, hλ 6, h4 λ ) Methodology to bound roundoff errors: 1) bound ε 0 = ỹ 0 y 0 ; 2) bound the error on each term α i (Gappa + stability); 3) bound local errors by M applications of Lemma 2;
46 33 / 34 General methodology y n+1 = M i=0 α i y n ỹ n+1 = M i=0 α i ỹ n (α i = hλ, hλ 3, hλ 6, h4 λ ) Methodology to bound roundoff errors: 1) bound ε 0 = ỹ 0 y 0 ; 2) bound the error on each term α i (Gappa + stability); 3) bound local errors by M applications of Lemma 2; 4) bound the global error by instantiating Theorem 1.
47 34 / 34 Conclusion and perspectives Conclusion: fine and mechanical analysis of roundoff errors; results on useful and classical methods (Euler, RK2, RK4); linear growth of roundoff errors; no compensation (for explicit one-step methods). Perspectives: overflows and underflows; formalization in Coq (proof assistant): based on the Flocq library [Boldo-Melquiond], based on the gappa and interval tactics [Melquiond], more general ODEs: complex, matricial, non-linear; more general methods: multi-step, variable step, implicit.
Round-off Error Analysis of Explicit One-Step Numerical Integration Methods
Round-off Error Analysis of Explicit One-Step Numerical Integration Methods Sylvie Boldo, Florian Faissole, Alexandre Chapoutot To cite this version: Sylvie Boldo, Florian Faissole, Alexandre Chapoutot
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 informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationHow to Compute the Area of a Triangle: a Formal Revisit with a Tighter Error Bound
How to Compute the Area of a Triangle: a Formal Revisit with a Tighter Error Bound Sylvie Boldo To cite this version: Sylvie Boldo. How to Compute the Area of a Triangle: a Formal Revisit with a Tighter
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationInnocuous Double Rounding of Basic Arithmetic Operations
Innocuous Double Rounding of Basic Arithmetic Operations Pierre Roux ISAE, ONERA Double rounding occurs when a floating-point value is first rounded to an intermediate precision before being rounded to
More informationHow to Compute the Area of a Triangle: a Formal Revisit
How to Compute the Area of a Triangle: a Formal Revisit Sylvie Boldo To cite this version: Sylvie Boldo. How to Compute the Area of a Triangle: a Formal Revisit. Alberto Nannarelli and Peter-Michael Seidel
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationA Verified ODE solver and Smale s 14th Problem
A Verified ODE solver and Smale s 14th Problem Fabian Immler Big Proof @ Isaac Newton Institute Jul 6, 2017 = Isabelle λ β α Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for
More informationValidated Explicit and Implicit Runge-Kutta Methods
Validated Explicit and Implicit Runge-Kutta Methods Alexandre Chapoutot joint work with Julien Alexandre dit Sandretto and Olivier Mullier U2IS, ENSTA ParisTech 8th Small Workshop on Interval Methods,
More informationWave Equation Numerical Resolution: a Comprehensive Mechanized Proof of a C Program
Journal of Automated Reasoning manuscript No. (will be inserted by the editor) Wave Equation Numerical Resolution: a Comprehensive Mechanized Proof of a C Program Sylvie Boldo François Clément Jean-Christophe
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE712B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationThe Classical Relative Error Bounds for Computing a2 + b 2 and c/ a 2 + b 2 in Binary Floating-Point Arithmetic are Asymptotically Optimal
-1- The Classical Relative Error Bounds for Computing a2 + b 2 and c/ a 2 + b 2 in Binary Floating-Point Arithmetic are Asymptotically Optimal Claude-Pierre Jeannerod Jean-Michel Muller Antoine Plet Inria,
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI097 TU) Thursday January 014, 18:0-1:0
More informationAutomated design of floating-point logarithm functions on integer processors
23rd IEEE Symposium on Computer Arithmetic Santa Clara, CA, USA, 10-13 July 2016 Automated design of floating-point logarithm functions on integer processors Guillaume Revy (presented by Florent de Dinechin)
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 informationNumerical Analysis II. Problem Sheet 9
H. Ammari W. Wu S. Yu Spring Term 08 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 9 Problem 9. Extrapolation of the Implicit Mid-Point Rule. Taking the implicit mid-point rule as a basis method,
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationOrdinary differential equation II
Ordinary Differential Equations ISC-5315 1 Ordinary differential equation II 1 Some Basic Methods 1.1 Backward Euler method (implicit method) The algorithm looks like this: y n = y n 1 + hf n (1) In contrast
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationAccurate Solution of Triangular Linear System
SCAN 2008 Sep. 29 - Oct. 3, 2008, at El Paso, TX, USA Accurate Solution of Triangular Linear System Philippe Langlois DALI at University of Perpignan France Nicolas Louvet INRIA and LIP at ENS Lyon France
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
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 informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationReducing round-off errors in rigid body dynamics
Reducing round-off errors in rigid body dynamics Gilles Vilmart INRIA Rennes, ENS Cachan Bretagne, avenue Robert Schuman, 357 Bruz, France Université de Genève, Section de mathématiques, -4 rue du Lièvre,
More informationTransformation of a PID Controller for Numerical Accuracy
Replace this file with prentcsmacro.sty for your meeting, or with entcsmacro.sty for your meeting. Both can be found at the ENTCS Macro Home Page. Transformation of a PID Controller for Numerical Accuracy
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 informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
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 informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
More informationTuning PI controllers in non-linear uncertain closed-loop systems with interval analysis
Tuning PI controllers in non-linear uncertain closed-loop systems with interval analysis J. Alexandre dit Sandretto, A. Chapoutot and O. Mullier U2IS, ENSTA ParisTech SYNCOP April 11, 2015 Closed-loop
More informationNumerical Analysis. EE, NCKU Tien-Hao Chang (Darby Chang)
Numerical Analysis EE, NCKU Tien-Hao Chang (Darby Chang) 1 In the previous slide Error (motivation) Floating point number system difference to real number system problem of roundoff Introduced/propagated
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationAccurate Multiple-Precision Gauss-Legendre Quadrature
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 Outline Motivation
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationOn various ways to split a floating-point number
On various ways to split a floating-point number Claude-Pierre Jeannerod Jean-Michel Muller Paul Zimmermann Inria, CNRS, ENS Lyon, Université de Lyon, Université de Lorraine France ARITH-25 June 2018 -2-
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 27 Theme of Last Three
More informationFaithful Horner Algorithm
ICIAM 07, Zurich, Switzerland, 16 20 July 2007 Faithful Horner Algorithm Philippe Langlois, Nicolas Louvet Université de Perpignan, France http://webdali.univ-perp.fr Ph. Langlois (Université de Perpignan)
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationPropagation of Roundoff Errors in Finite Precision Computations: a Semantics Approach 1
Propagation of Roundoff Errors in Finite Precision Computations: a Semantics Approach 1 Matthieu Martel CEA - Recherche Technologique LIST-DTSI-SLA CEA F91191 Gif-Sur-Yvette Cedex, France e-mail : mmartel@cea.fr
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationNewton-Raphson Algorithms for Floating-Point Division Using an FMA
Newton-Raphson Algorithms for Floating-Point Division Using an FMA Nicolas Louvet, Jean-Michel Muller, Adrien Panhaleux Abstract Since the introduction of the Fused Multiply and Add (FMA) in the IEEE-754-2008
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationAccurate 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 informationCONVERGENCE AND STABILITY CONSTANT OF THE THETA-METHOD
Conference Applications of Mathematics 2013 in honor of the 70th birthday of Karel Segeth. Jan Brandts, Sergey Korotov, Michal Křížek, Jakub Šístek, and Tomáš Vejchodský (Eds.), Institute of Mathematics
More informationFixed Points and Contractive Transformations. Ron Goldman Department of Computer Science Rice University
Fixed Points and Contractive Transformations Ron Goldman Department of Computer Science Rice University Applications Computer Graphics Fractals Bezier and B-Spline Curves and Surfaces Root Finding Newton
More informationA Brief Introduction to Numerical Methods for Differential Equations
A Brief Introduction to Numerical Methods for Differential Equations January 10, 2011 This tutorial introduces some basic numerical computation techniques that are useful for the simulation and analysis
More information1 Ordinary Differential Equations
Ordinary Differential Equations.0 Mathematical Background.0. Smoothness Definition. A function f defined on [a, b] is continuous at ξ [a, b] if lim x ξ f(x) = f(ξ). Remark Note that this implies existence
More informationOrdinary Differential Equations: Initial Value problems (IVP)
Chapter Ordinary Differential Equations: Initial Value problems (IVP) Many engineering applications can be modeled as differential equations (DE) In this book, our emphasis is about how to use computer
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationInvestigation on the Most Efficient Ways to Solve the Implicit Equations for Gauss Methods in the Constant Stepsize Setting
Applied Mathematical Sciences, Vol. 12, 2018, no. 2, 93-103 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711340 Investigation on the Most Efficient Ways to Solve the Implicit Equations
More informationThe Definition and Numerical Method of Final Value Problem and Arbitrary Value Problem Shixiong Wang 1*, Jianhua He 1, Chen Wang 2, Xitong Li 1
The Definition and Numerical Method of Final Value Problem and Arbitrary Value Problem Shixiong Wang 1*, Jianhua He 1, Chen Wang 2, Xitong Li 1 1 School of Electronics and Information, Northwestern Polytechnical
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationEnergy-Preserving Runge-Kutta methods
Energy-Preserving Runge-Kutta methods Fasma Diele, Brigida Pace Istituto per le Applicazioni del Calcolo M. Picone, CNR, Via Amendola 122, 70126 Bari, Italy f.diele@ba.iac.cnr.it b.pace@ba.iac.cnr.it SDS2010,
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 32 Theme of Last Few
More informationChapter 10. Initial value Ordinary Differential Equations
Chapter 10 Initial value Ordinary Differential Equations Consider the problem of finding a function y(t) that satisfies the following ordinary differential equation (ODE): dy dt = f(t, y), a t b. The function
More informationSouthern Methodist University.
Title: Continuous extensions Name: Lawrence F. Shampine 1, Laurent O. Jay 2 Affil./Addr. 1: Department of Mathematics Southern Methodist University Dallas, TX 75275 USA Phone: +1 (972) 690-8439 E-mail:
More informationSemantics and Validation Lecture 5. Numerical program analysis.
Semantics and Validation Lecture 5. Numerical program analysis. Laboratoire Modélisation et Analyse de Systèmes en Interaction, CEA-LIST and Ecole Polytechnique Eric Goubault and Sylvie Putot January 14,
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Second Edition Translated by R. Bartels, W. Gautschi, and C. Witzgall With 35 Illustrations Springer Contents Preface to the Second Edition Preface
More informationRunge Kutta Theory and Constraint Programming
Runge Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto julien.alexandre-dit-sandretto@ensta-paristech.fr U2IS, ENSTA ParisTech, Université Paris-Saclay, 828 bd des Maréchaux, 9762
More informationORBIT 14 The propagator. A. Milani et al. Dipmat, Pisa, Italy
ORBIT 14 The propagator A. Milani et al. Dipmat, Pisa, Italy Orbit 14 Propagator manual 1 Contents 1 Propagator - Interpolator 2 1.1 Propagation and interpolation............................ 2 1.1.1 Introduction..................................
More informationDepartment of Applied Mathematics Preliminary Examination in Numerical Analysis August, 2013
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August, 013 August 8, 013 Solutions: 1 Root Finding (a) Let the root be x = α We subtract α from both sides of x n+1 = x
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
More informationLecture Notes Multibody Dynamics B, wb1413
Lecture Notes Multibody Dynamics B, wb43 A. L. Schwab & Guido M.J. Delhaes Laboratory for Engineering Mechanics Mechanical Engineering Delft University of Technolgy The Netherlands April 9, 29 Contents
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationProgram-ing in Coq. Matthieu Sozeau under the direction of Christine Paulin-Mohring
Program-ing in Coq Matthieu Sozeau under the direction of Christine Paulin-Mohring LRI, Univ. Paris-Sud - Démons Team & INRIA Saclay - ProVal Project Foundations of Programming seminar February 15th 2008
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationFloating Point Number Systems. Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le
Floating Point Number Systems Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le 1 Overview Real number system Examples Absolute and relative errors Floating point numbers Roundoff
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationChapter 1 Error Analysis
Chapter 1 Error Analysis Several sources of errors are important for numerical data processing: Experimental uncertainty: Input data from an experiment have a limited precision. Instead of the vector of
More informationNumerical Differential Equations: IVP
Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential
More informationOrdinary Differential Equations
Ordinary Differential Equations Introduction: first order ODE We are given a function f(t,y) which describes a direction field in the (t,y) plane an initial point (t 0,y 0 ) We want to find a function
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
More informationNumerical Program Optimization by Automatic Improvement of the Accuracy of Computations. Nasrine Damouche* and Matthieu Martel
Int. J. of Intelligent Engineering Informatics, Vol. x, No. x, 1 26 1 Numerical Program Optimization by Automatic Improvement of the Accuracy of Computations Nasrine Damouche* and Matthieu Martel Université
More informationAccelerated Shift-and-Add algorithms
c,, 1 12 () Kluwer Academic Publishers, Boston. Manufactured in The Netherlands. Accelerated Shift-and-Add algorithms N. REVOL AND J.-C. YAKOUBSOHN nathalie.revol@univ-lille1.fr, yak@cict.fr Lab. ANO,
More informationConstruction of Real Algebraic Numbers in Coq
Construction of Real Algebraic Numbers in Coq INRIA Saclay Île-de-France LIX École Polytechnique INRIA Microsoft Research Joint Centre cohen@crans.org August 13, 2012 Why algebraic numbers? Field strictly
More informationThe Dynamics of Runge Kutta Methods
The Dynamics of Runge Kutta Methods Julyan H. E. Cartwright & Oreste Piro School of Mathematical Sciences Queen Mary and Westfield College University of London Mile End Road London E1 4NS U.K. Int. J.
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationNumerical Methods for Engineers
Numerical Methods for Engineers SEVENTH EDITION Steven C Chopra Berger Chair in Computing and Engineering Tufts University Raymond P. Canal Professor Emeritus of Civil Engineering of Michiaan University
More informationSolving scalar IVP s : Runge-Kutta Methods
Solving scalar IVP s : Runge-Kutta Methods Josh Engwer Texas Tech University March 7, NOTATION: h step size x n xt) t n+ t + h x n+ xt n+ ) xt + h) dx = ft, x) SCALAR IVP ASSUMED THROUGHOUT: dt xt ) =
More information