Introduction to Scientific Computing Languages
|
|
- Hugh Ellis
- 6 years ago
- Views:
Transcription
1 1 / 19 Introduction to Scientific Computing Languages Prof. Paolo Bientinesi pauldj@aices.rwth-aachen.de
2 Numerical Representation 2 / 19 Numbers 123 = (first 40 digits) π = In general: Infinite number of digits
3 Numerical Representation 2 / 19 Numbers 123 = (first 40 digits) π = In general: Infinite number of digits Computers Finite memory
4 Computers: Inexact Numbers 3 / 19 Infinite numbers vs. finite memory Approximated numbers
5 Computers: Inexact Numbers 3 / 19 Infinite numbers vs. finite memory Approximated numbers How many digits? A pre-determined amount π=
6 Computers: Inexact Numbers 3 / 19 Infinite numbers vs. finite memory Approximated numbers How many digits? A pre-determined amount π= Using 4 digits: π = Modern computers: normally 8 or 16 digits, single / double precision.
7 Computers: Inexact Numbers 3 / 19 Infinite numbers vs. finite memory Approximated numbers How many digits? A pre-determined amount π= Using 4 digits: π = Modern computers: normally 8 or 16 digits, single / double precision. Alternatively? Extended precision Variable precision Symbolic representation
8 4 / 19 Computers: Approximated Computations 4-digit representation Inexact Arithmetic = =
9 4 / 19 Computers: Approximated Computations 4-digit representation Inexact Arithmetic = =
10 4 / 19 Computers: Approximated Computations 4-digit representation Inexact Arithmetic = = Truncated Rounded 124.0
11 4 / 19 Computers: Approximated Computations 4-digit representation Inexact Arithmetic = = Truncated Rounded Associativity? Exact arithmetic: ( ) = ( )
12 4 / 19 Computers: Approximated Computations 4-digit representation Inexact Arithmetic = = Truncated Rounded Associativity? Exact arithmetic: ( ) = ( ) Inexact arithmetic: ( ) = 124.4
13 4 / 19 Computers: Approximated Computations 4-digit representation Inexact Arithmetic = = Truncated Rounded Associativity? No! Exact arithmetic: ( ) = ( ) Inexact arithmetic: ( ) = ( ) = 124.5
14 Error Analysis 5 / 19 f : A B, y = f(x) Es.: f(x) = x 2 + sin(2 x) x = π, f(x) =? 123
15 Error Analysis 5 / 19 f : A B, y = f(x) Es.: f(x) = x 2 + sin(2 x) x = π, f(x) =? 123 Exact arithmetic: ( π ) 2 ( + sin 2 π ) =...
16 Error Analysis 5 / 19 f : A B, y = f(x) Es.: f(x) = x 2 + sin(2 x) x = π, f(x) =? 123 Exact arithmetic: ( π ) 2 ( + sin 2 π ) =... Inexact arithmetic: x ˆx, f ˆf ˆf(ˆx) instead of f(x)
17 Errors 6 / 19 Representation Errors Roundoff Errors Algorithmic Errors
18 Known Disasters 7 / 19
19 Known Disasters 7 / 19 Patriot Missile, 1991 Scud launched from Iraq againt US military base in South Arabia. US Patriot s missile missed the incoming Scud. 28 casualties. Cancellation
20 Known Disasters 7 / 19 Patriot Missile, 1991 Scud launched from Iraq againt US military base in South Arabia. US Patriot s missile missed the incoming Scud. 28 casualties. Cancellation Spaceship Ariane 5 launched in 1996, destroyed 37 seconds after liftoff. Overflow
21 Floating Point Numbers 8 / 19 y = ±d 0.d 1 d 2... d t 1 β e
22 Floating Point Numbers 8 / 19 β: base (radix) y = ±d 0.d 1 d 2... d t 1 β e
23 Floating Point Numbers 8 / 19 β: base (radix) y = ±d 0.d 1 d 2... d t 1 β e t: precision = # slots for the mantissa (d 0.d 1... d t 1 )
24 Floating Point Numbers 8 / 19 y = ±d 0.d 1 d 2... d t 1 β e β: base (radix) t: precision = # slots for the mantissa (d 0.d 1... d t 1 ) 0 d i β 1: digits
25 Floating Point Numbers 8 / 19 y = ±d 0.d 1 d 2... d t 1 β e β: base (radix) t: precision = # slots for the mantissa (d 0.d 1... d t 1 ) 0 d i β 1: digits e min e e max : exponent range
26 Floating Point Numbers 8 / 19 y = ±d 0.d 1 d 2... d t 1 β e β: base (radix) t: precision = # slots for the mantissa (d 0.d 1... d t 1 ) 0 d i β 1: digits e min e e max : exponent range Normalization: d 0 = 1; d 0 used for the sign ±
27 Floating Point Numbers 8 / 19 y = ±d 0.d 1 d 2... d t 1 β e β: base (radix) t: precision = # slots for the mantissa (d 0.d 1... d t 1 ) 0 d i β 1: digits e min e e max : exponent range Normalization: d 0 = 1; d 0 used for the sign ± Arithmetic β t e min e max Single precision Double precision
28 Floating Point Numbers 8 / 19 y = ±d 0.d 1 d 2... d t 1 β e β: base (radix) t: precision = # slots for the mantissa (d 0.d 1... d t 1 ) 0 d i β 1: digits e min e e max : exponent range Normalization: d 0 = 1; d 0 used for the sign ± Arithmetic β t e min e max Single precision Double precision HP calculator IBM Setun 3 Quadruple prec
29 Floating Point Numbers (2) 9 / 19 What are the nonnegative points in β = 2, t = 3, e min = 1, e max = 3?
30 Floating Point Numbers (2) 9 / 19 What are the nonnegative points in β = 2, t = 3, e min = 1, e max = 3? Floating point numbers are non-equidistant!
31 Floating Point Numbers (2) 9 / 19 What are the nonnegative points in β = 2, t = 3, e min = 1, e max = 3? Floating point numbers are non-equidistant! How to represent 0?
32 Floating Point Numbers (2) 9 / 19 What are the nonnegative points in β = 2, t = 3, e min = 1, e max = 3? Floating point numbers are non-equidistant! How to represent 0? underflow? overflow? NaN?
33 IEEE single precision FYI 10 / 19
34 IEEE double precision FYI 11 / 19
35 Machine Precision 12 / 19 Machine Precision u Smallest positive number such that [1 + u] 1
36 Machine Precision 12 / 19 Machine Precision u Smallest positive number such that [1 + u] 1 Largest positive number such that [1 + u] = 1
37 Machine Precision 12 / 19 Machine Precision u Smallest positive number such that [1 + u] 1 Largest positive number such that [1 + u] = β1 t
38 Machine Precision 12 / 19 Machine Precision u Smallest positive number such that [1 + u] 1 Largest positive number such that [1 + u] = β1 t Distance between 1 and the next floating point number Machine epsilon, ɛ M, u
39 Representation Error 13 / 19 f min = smallest positive floating point number f max = largest positive floating point number x = [x] = floating point representation of x Theorem: Let x R and x [f min, f max ] Then x = x(1 + δ 1 ) where δ 1 u
40 Representation Error 13 / 19 f min = smallest positive floating point number f max = largest positive floating point number x = [x] = floating point representation of x Theorem: Let x R and x [f min, f max ] Then x = x(1 + δ 1 ) where δ 1 u Also, x = x/(1 + δ 2 ) where δ 2 u Note: δ 1 and δ 2 are functions of x
41 Smallest distance 0 14 / 19
42 Smallest distance 0 14 / 19 x and y are floating point numbers, x y; how small can z be? z := x y
43 Smallest distance 0 14 / 19 x and y are floating point numbers, x y; how small can z be? z := x y vs. z := x y max( x, y )
44 Smallest distance 0 14 / 19 x and y are floating point numbers, x y; how small can z be? z := x y vs. z := x y max( x, y ) z := x y 1) = 2 1 t 2) emin = 2 emin t+1 3) emin+t 1 = 2 emin
45 Smallest distance 0 14 / 19 x and y are floating point numbers, x y; how small can z be? z := x y vs. z := x y max( x, y ) z := x y 1) = 2 1 t 2) emin = 2 emin t+1 3) emin+t 1 = 2 emin
46 Smallest distance 0 14 / 19 x and y are floating point numbers, x y; how small can z be? z := x y vs. z := x y max( x, y ) z := x y 1) = 2 1 t 2) emin = 2 emin t+1 3) emin+t 1 = 2 emin z := x y max( x, y ) 1) = 2 1 t 2) emin = 2 emin t+1 3) emin+t 1 = 2 emin
47 Smallest distance 0 14 / 19 x and y are floating point numbers, x y; how small can z be? z := x y vs. z := x y max( x, y ) z := x y 1) = 2 1 t 2) emin = 2 emin t+1 3) emin+t 1 = 2 emin z := x y max( x, y ) 1) = 2 1 t 2) emin = 2 emin t+1 3) emin+t 1 = 2 emin
48 Roundoff Error 15 / 19 Notation: [exp] denotes the evaluation of exp in floating point arithmetic. Assuming a left-to-right evaluation, it holds [ [[x] ] [ ] ] [x + y + z/w] = + [y] + [z]/[w]
49 Roundoff Error 15 / 19 Notation: [exp] denotes the evaluation of exp in floating point arithmetic. Assuming a left-to-right evaluation, it holds [ [[x] ] [ ] ] [x + y + z/w] = + [y] + [z]/[w] Floating Point Arithmetic Theorem: (Standard and Alternative Computational Models) Let x and y be floating point numbers Then [x op y] = (x op y)(1 + ɛ 1 ), where ɛ 1 u, and op {+,,, /}
50 Roundoff Error 15 / 19 Notation: [exp] denotes the evaluation of exp in floating point arithmetic. Assuming a left-to-right evaluation, it holds [ [[x] ] [ ] ] [x + y + z/w] = + [y] + [z]/[w] Floating Point Arithmetic Theorem: (Standard and Alternative Computational Models) Let x and y be floating point numbers Then [x op y] = (x op y)(1 + ɛ 1 ), where ɛ 1 u, and op {+,,, /} Also, [x op y] = x op y (1 + ɛ 2 ), where ɛ 2 u, and op {+,,, /} Note: ɛ 1 and ɛ 2 are functions of x, y and op
51 Example: Dot Product 16 / 19 x, y R n ; κ := x T y ( ((χ0 κ := ψ 0 + χ 1 ψ 1 ) + ) ) + χ n 2 ψ n 2 + χ n 1 ψ n 1 ˇκ = = ( ((χ0 ψ 0 (1 + ɛ (0) ) + χ 1 ψ 1 (1 + ɛ (1) ) ) (1 + ɛ (1) + ) + ) +χ n 1 ψ n 1 (1 + ɛ (n 1) ) (1 + ɛ (n 1) + ) n 1 n 1 χ i ψ i (1 + ɛ (i) ) (1 + ɛ (j) i=0 j=i + ) ) (1 + ɛ (n 2) + ) where ɛ (0) + = 0 and ɛ (0), ɛ (j), ɛ (j) + u for j = 1,..., n 1
52 Case Studies Bisection Recursive implementation Base case? 17 / 19
53 Case Studies Bisection Recursive implementation Base case? Π 2 / i=1 1 i 2 vs. 1 i= i 2 17 / 19
54 Case Studies Bisection Recursive implementation Base case? Π 2 / i=1 1 i 2 vs. 1 i= i 2 sqsqrt vs. sqrtsq 17 / 19
55 Case Studies Bisection Recursive implementation Base case? Π 2 / i=1 1 i 2 vs. 1 i= i 2 sqsqrt vs. sqrtsq Ax = b A = [ ɛ ] 17 / 19
56 Backward Stability 18 / 19 Let f : D R be a map from the domain D to the range R. Let ˆf : D R represent the execution in floating point arithmetic of a given algorithm A that computes f. A is said to be backward stable if for all x D there exists a perturbed input x D, close to x, such that ˆf(x) = f( x).
57 Backward Stability 18 / 19 Let f : D R be a map from the domain D to the range R. Let ˆf : D R represent the execution in floating point arithmetic of a given algorithm A that computes f. A is said to be backward stable if for all x D there exists a perturbed input x D, close to x, such that ˆf(x) = f( x). I.e., the result computed in floating point arithmetic ( ˆf(x)) equals the result obtained when the mathematically exact function (f) is applied to slightly perturbed data ( x). The difference between x and x, is the perturbation to the original input x.
58 References 19 / 19 IEEE and IEEE : Standard for Floating-Point Arithmetic Book: Accuracy and Stability of Numerical Algorithms, by Nick Higham Article: What every computer scientist should know about floating-point arithmetic, by David Goldberg Book: Numerical Computing with IEEE Floating Point Arithmetic, by Michael Overton
Introduction 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 informationLecture 7. Floating point arithmetic and stability
Lecture 7 Floating point arithmetic and stability 2.5 Machine representation of numbers Scientific notation: 23 }{{} }{{} } 3.14159265 {{} }{{} 10 sign mantissa base exponent (significand) s m β e A floating
More informationECS 231 Computer Arithmetic 1 / 27
ECS 231 Computer Arithmetic 1 / 27 Outline 1 Floating-point numbers and representations 2 Floating-point arithmetic 3 Floating-point error analysis 4 Further reading 2 / 27 Outline 1 Floating-point numbers
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 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 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 informationNotes for Chapter 1 of. Scientific Computing with Case Studies
Notes for Chapter 1 of Scientific Computing with Case Studies Dianne P. O Leary SIAM Press, 2008 Mathematical modeling Computer arithmetic Errors 1999-2008 Dianne P. O'Leary 1 Arithmetic and Error What
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 informationArithmetic and Error. How does error arise? How does error arise? Notes for Part 1 of CMSC 460
Notes for Part 1 of CMSC 460 Dianne P. O Leary Preliminaries: Mathematical modeling Computer arithmetic Errors 1999-2006 Dianne P. O'Leary 1 Arithmetic and Error What we need to know about error: -- how
More informationFloating-point Computation
Chapter 2 Floating-point Computation 21 Positional Number System An integer N in a number system of base (or radix) β may be written as N = a n β n + a n 1 β n 1 + + a 1 β + a 0 = P n (β) where a i are
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 5. Ax = b.
CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 5 GENE H GOLUB Suppose we want to solve We actually have an approximation ξ such that 1 Perturbation Theory Ax = b x = ξ + e The question is, how
More informationMathematical preliminaries and error analysis
Mathematical preliminaries and error analysis Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan September 12, 2015 Outline 1 Round-off errors and computer arithmetic
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More information1 Floating point arithmetic
Introduction to Floating Point Arithmetic Floating point arithmetic Floating point representation (scientific notation) of numbers, for example, takes the following form.346 0 sign fraction base exponent
More informationLecture Notes 7, Math/Comp 128, Math 250
Lecture Notes 7, Math/Comp 128, Math 250 Misha Kilmer Tufts University October 23, 2005 Floating Point Arithmetic We talked last time about how the computer represents floating point numbers. In a floating
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 informationNumerical Analysis. Yutian LI. 2018/19 Term 1 CUHKSZ. Yutian LI (CUHKSZ) Numerical Analysis 2018/19 1 / 41
Numerical Analysis Yutian LI CUHKSZ 2018/19 Term 1 Yutian LI (CUHKSZ) Numerical Analysis 2018/19 1 / 41 Reference Books BF R. L. Burden and J. D. Faires, Numerical Analysis, 9th edition, Thomsom Brooks/Cole,
More informationChapter 4 Number Representations
Chapter 4 Number Representations SKEE2263 Digital Systems Mun im/ismahani/izam {munim@utm.my,e-izam@utm.my,ismahani@fke.utm.my} February 9, 2016 Table of Contents 1 Fundamentals 2 Signed Numbers 3 Fixed-Point
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationApplied Mathematics 205. Unit 0: Overview of Scientific Computing. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit 0: Overview of Scientific Computing Lecturer: Dr. David Knezevic Scientific Computing Computation is now recognized as the third pillar of science (along with theory and experiment)
More informationChapter 1 Mathematical Preliminaries and Error Analysis
Numerical Analysis (Math 3313) 2019-2018 Chapter 1 Mathematical Preliminaries and Error Analysis Intended learning outcomes: Upon successful completion of this chapter, a student will be able to (1) list
More information1 ERROR ANALYSIS IN COMPUTATION
1 ERROR ANALYSIS IN COMPUTATION 1.2 Round-Off Errors & Computer Arithmetic (a) Computer Representation of Numbers Two types: integer mode (not used in MATLAB) floating-point mode x R ˆx F(β, t, l, u),
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 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 information4.2 Floating-Point Numbers
101 Approximation 4.2 Floating-Point Numbers 4.2 Floating-Point Numbers The number 3.1416 in scientific notation is 0.31416 10 1 or (as computer output) -0.31416E01..31416 10 1 exponent sign mantissa base
More informationFLOATING POINT ARITHMETHIC - ERROR ANALYSIS
FLOATING POINT ARITHMETHIC - ERROR ANALYSIS Brief review of floating point arithmetic Model of floating point arithmetic Notation, backward and forward errors 3-1 Roundoff errors and floating-point arithmetic
More informationFLOATING POINT ARITHMETHIC - ERROR ANALYSIS
FLOATING POINT ARITHMETHIC - ERROR ANALYSIS Brief review of floating point arithmetic Model of floating point arithmetic Notation, backward and forward errors Roundoff errors and floating-point arithmetic
More informationIntroduction of Computer-Aided Nano Engineering
Introduction of Computer-Aided Nano Engineering Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Computational
More informationNumerical Methods - Preliminaries
Numerical Methods - Preliminaries Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Preliminaries 2013 1 / 58 Table of Contents 1 Introduction to Numerical Methods Numerical
More informationContents Experimental Perturbations Introduction to Interval Arithmetic Review Questions Problems...
Contents 2 How to Obtain and Estimate Accuracy 1 2.1 Basic Concepts in Error Estimation................ 1 2.1.1 Sources of Error.................... 1 2.1.2 Absolute and Relative Errors............. 4
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 informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationNumerical Algorithms. IE 496 Lecture 20
Numerical Algorithms IE 496 Lecture 20 Reading for This Lecture Primary Miller and Boxer, Pages 124-128 Forsythe and Mohler, Sections 1 and 2 Numerical Algorithms Numerical Analysis So far, we have looked
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Computer Representation of Numbers Counting numbers (unsigned integers) are the numbers 0,
More informationLecture 1. Introduction to Numerical Methods. T. Gambill. Department of Computer Science University of Illinois at Urbana-Champaign.
Lecture 1 Introduction to Numerical Methods T. Gambill Department of Computer Science University of Illinois at Urbana-Champaign June 10, 2013 T. Gambill (UIUC) CS 357 June 10, 2013 1 / 52 Course Info
More informationErrors. Intensive Computation. Annalisa Massini 2017/2018
Errors Intensive Computation Annalisa Massini 2017/2018 Intensive Computation - 2017/2018 2 References Scientific Computing: An Introductory Survey - Chapter 1 M.T. Heath http://heath.cs.illinois.edu/scicomp/notes/index.html
More informationNumerical Analysis and Computing
Numerical Analysis and Computing Lecture Notes #02 Calculus Review; Computer Artihmetic and Finite Precision; and Convergence; Joe Mahaffy, mahaffy@math.sdsu.edu Department of Mathematics Dynamical Systems
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationHomework 2 Foundations of Computational Math 1 Fall 2018
Homework 2 Foundations of Computational Math 1 Fall 2018 Note that Problems 2 and 8 have coding in them. Problem 2 is a simple task while Problem 8 is very involved (and has in fact been given as a programming
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 1 Chapter 4 Roundoff and Truncation Errors PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationComputer Arithmetic. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Computer Arithmetic
Computer Arithmetic MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Machine Numbers When performing arithmetic on a computer (laptop, desktop, mainframe, cell phone,
More informationCS412: Introduction to Numerical Methods
CS412: Introduction to Numerical Methods MIDTERM #1 2:30PM - 3:45PM, Tuesday, 03/10/2015 Instructions: This exam is a closed book and closed notes exam, i.e., you are not allowed to consult any textbook,
More informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationIntroduction CSE 541
Introduction CSE 541 1 Numerical methods Solving scientific/engineering problems using computers. Root finding, Chapter 3 Polynomial Interpolation, Chapter 4 Differentiation, Chapter 4 Integration, Chapters
More informationApplied Numerical Analysis (AE2220-I) R. Klees and R.P. Dwight
Applied Numerical Analysis (AE0-I) R. Klees and R.P. Dwight February 018 Contents 1 Preliminaries: Motivation, Computer arithmetic, Taylor series 1 1.1 Numerical Analysis Motivation..........................
More informationAn Introduction to Numerical Analysis
An Introduction to Numerical Analysis Department of Mathematical Sciences, NTNU 21st august 2012 Practical issues webpage: http://wiki.math.ntnu.no/tma4215/2012h/start Lecturer:, elenac@math.ntnu.no Lecures:
More informationChapter 1 Computer Arithmetic
Numerical Analysis (Math 9372) 2017-2016 Chapter 1 Computer Arithmetic 1.1 Introduction Numerical analysis is a way to solve mathematical problems by special procedures which use arithmetic operations
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 informationRound-off Errors and Computer Arithmetic - (1.2)
Round-off Errors and Comuter Arithmetic - (.). Round-off Errors: Round-off errors is roduced when a calculator or comuter is used to erform real number calculations. That is because the arithmetic erformed
More informationIntroduction, basic but important concepts
Introduction, basic but important concepts Felix Kubler 1 1 DBF, University of Zurich and Swiss Finance Institute October 7, 2017 Felix Kubler Comp.Econ. Gerzensee, Ch1 October 7, 2017 1 / 31 Economics
More informationInteractions between parallelism and numerical stability, accuracy
Interactions between parallelism and numerical stability, accuracy Prof. Richard Vuduc Georgia Institute of Technology CSE/CS 8803 PNA: Parallel Numerical Algorithms [L.13] Tuesday, February 19, 2008 1
More informationChapter 1. Numerical Errors. Module No. 1. Errors in Numerical Computations
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-73209 email: anita.buie@gmail.com . Chapter Numerical Errors Module
More information1.1 COMPUTER REPRESENTATION OF NUM- BERS, REPRESENTATION ERRORS
Chapter 1 NUMBER REPRESENTATION, ERROR ANALYSIS 1.1 COMPUTER REPRESENTATION OF NUM- BERS, REPRESENTATION ERRORS Floating-point representation x t,r of a number x: x t,r = m t P cr, where: P - base (the
More informationAccurate polynomial evaluation in floating point arithmetic
in floating point arithmetic Université de Perpignan Via Domitia Laboratoire LP2A Équipe de recherche en Informatique DALI MIMS Seminar, February, 10 th 2006 General motivation Provide numerical algorithms
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 informationChapter 1: Introduction and mathematical preliminaries
Chapter 1: Introduction and mathematical preliminaries Evy Kersalé September 26, 2011 Motivation Most of the mathematical problems you have encountered so far can be solved analytically. However, in real-life,
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 informationMath 128A: Homework 2 Solutions
Math 128A: Homework 2 Solutions Due: June 28 1. In problems where high precision is not needed, the IEEE standard provides a specification for single precision numbers, which occupy 32 bits of storage.
More informationQUADRATIC PROGRAMMING?
QUADRATIC PROGRAMMING? WILLIAM Y. SIT Department of Mathematics, The City College of The City University of New York, New York, NY 10031, USA E-mail: wyscc@cunyvm.cuny.edu This is a talk on how to program
More informationIntroduction to Scientific Computing
(Lecture 2: Machine precision and condition number) B. Rosić, T.Moshagen Institute of Scientific Computing General information :) 13 homeworks (HW) Work in groups of 2 or 3 people Each HW brings maximally
More informationChapter 1 Mathematical Preliminaries and Error Analysis
Chapter 1 Mathematical Preliminaries and Error Analysis Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128A Numerical Analysis Limits and Continuity
More informationLecture notes to the course. Numerical Methods I. Clemens Kirisits
Lecture notes to the course Numerical Methods I Clemens Kirisits November 8, 08 ii Preface These lecture notes are intended as a written companion to the course Numerical Methods I taught from 06 to 08
More informationHomework 2. Matthew Jin. April 10, 2014
Homework Matthew Jin April 10, 014 1a) The relative error is given by ŷ y y, where ŷ represents the observed output value, and y represents the theoretical output value. In this case, the observed output
More informationALU (3) - Division Algorithms
HUMBOLDT-UNIVERSITÄT ZU BERLIN INSTITUT FÜR INFORMATIK Lecture 12 ALU (3) - Division Algorithms Sommersemester 2002 Leitung: Prof. Dr. Miroslaw Malek www.informatik.hu-berlin.de/rok/ca CA - XII - ALU(3)
More informationLecture 2 MODES OF NUMERICAL COMPUTATION
1. Diversity of Numbers Lecture 2 Page 1 It is better to solve the right problem the wrong way than to solve the wrong problem the right way. The purpose of computing is insight, not numbers. Richard Wesley
More informationAn Introduction to Numerical Analysis. James Brannick. The Pennsylvania State University
An Introduction to Numerical Analysis James Brannick The Pennsylvania State University Contents Chapter 1. Introduction 5 Chapter 2. Computer arithmetic and Error Analysis 7 Chapter 3. Approximation and
More informationHow do computers represent numbers?
How do computers represent numbers? Tips & Tricks Week 1 Topics in Scientific Computing QMUL Semester A 2017/18 1/10 What does digital mean? The term DIGITAL refers to any device that operates on discrete
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 informationAn Introduction to the Quality of Computed Solutions 1
An Introduction to the Quality of Computed Solutions 1 Sven Hammarling The Numerical Algorithms Group Ltd Wilkinson House Jordan Hill Road Oxford, OX2 8DR, UK sven@nag.co.uk October 12, 2005 1 Based on
More informationExercises MAT-INF1100. Øyvind Ryan
Exercises MAT-INF1100 Øyvind Ryan February 19, 2013 1. Formulate an algorithm for adding two three-digit numbers. You may assume that it is known how to sum one-digit numbers. Answer: We represent the
More informationCHAPTER 11. A Revision. 1. The Computers and Numbers therein
CHAPTER A Revision. The Computers and Numbers therein Traditional computer science begins with a finite alphabet. By stringing elements of the alphabet one after another, one obtains strings. A set of
More information8/13/16. Data analysis and modeling: the tools of the trade. Ø Set of numbers. Ø Binary representation of numbers. Ø Floating points.
Data analysis and modeling: the tools of the trade Patrice Koehl Department of Biological Sciences National University of Singapore http://www.cs.ucdavis.edu/~koehl/teaching/bl5229 koehl@cs.ucdavis.edu
More informationESO 208A: Computational Methods in Engineering. Saumyen Guha
ESO 208A: Computational Methods in Engineering Introduction, Error Analysis Saumyen Guha Department of Civil Engineering IIT Kanpur What is Computational Methods or Numerical Methods in Engineering? Formulation
More informationBASIC COMPUTER ARITHMETIC
BASIC COMPUTER ARITHMETIC TSOGTGEREL GATUMUR Abstract. First, we consider how integers and fractional numbers are represented and manipulated internally on a computer. Then we develop a basic theoretical
More informationCOMPUTER SCIENCE & ENGINEERING DEPARTMENT Harmattan Semester, 2013/2014 Session. CSC 307: Numerical Computations I [PRACTICAL LAB ONE]
ỌBÁFÉ. MI AWÓLÓ. WÒ. UNIVERSITY COMPUTER SCIENCE & ENGINEERING DEPARTMENT Harmattan Semester, 2013/2014 Session. CSC 307: Numerical Computations I [PRACTICAL LAB ONE] ỌDÉ. JỌBÍ Ọdé.túnjí Àjàdí THIS DOCUMENT
More informationOne-Sided Difference Formula for the First Derivative
POLYTECHNIC UNIVERSITY Department of Computer and Information Science One-Sided Difference Formula for the First Derivative K. Ming Leung Abstract: Derive a one-sided formula for the first derive of a
More informationLineare Algebra. Endliche Arithmetik. Walter Gander. ETH Zürich
Lineare Algebra Endliche Arithmetik Walter Gander ETH Zürich Contents Chapter 1. Finite Arithmetic................... 1 1.1 Introductory Example.................... 1 1. Real Numbers and Machine Numbers............
More informationOutline. Math Numerical Analysis. Errors. Lecture Notes Linear Algebra: Part B. Joseph M. Mahaffy,
Math 54 - Numerical Analysis Lecture Notes Linear Algebra: Part B Outline Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences
More informationTu: 9/3/13 Math 471, Fall 2013, Section 001 Lecture 1
Tu: 9/3/13 Math 71, Fall 2013, Section 001 Lecture 1 1 Course intro Notes : Take attendance. Instructor introduction. Handout : Course description. Note the exam days (and don t be absent). Bookmark the
More informationProgram 1 Foundations of Computational Math 1 Fall 2018
Program 1 Foundations of Computational Math 1 Fall 2018 Due date: 11:59PM on Friday, 28 September 2018 Written Exercises Problem 1 Consider the summation σ = n ξ i using the following binary fan-in tree
More informationCombining Static Analysis and Testing for Overflow and Roundoff Error Detection DO THI BICH NGOC
Combining Static Analysis and Testing for Overflow and Roundoff Error Detection by DO THI BICH NGOC submitted to Japan Advanced Institute of Science and Technology in partial fulfillment of the requirements
More informationGetting tight error bounds in floating-point arithmetic: illustration with complex functions, and the real x n function
Getting tight error bounds in floating-point arithmetic: illustration with complex functions, and the real x n function Jean-Michel Muller collaborators: S. Graillat, C.-P. Jeannerod, N. Louvet V. Lefèvre
More informationIntroduction to Finite Di erence Methods
Introduction to Finite Di erence Methods ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: Introduction to Finite Di erence Approximation
More informationNumerical Mathematical Analysis
Numerical Mathematical Analysis Numerical Mathematical Analysis Catalin Trenchea Department of Mathematics University of Pittsburgh September 20, 2010 Numerical Mathematical Analysis Math 1070 Numerical
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 information1 What is numerical analysis and scientific computing?
Mathematical preliminaries 1 What is numerical analysis and scientific computing? Numerical analysis is the study of algorithms that use numerical approximation (as opposed to general symbolic manipulations)
More informationUNIT V FINITE WORD LENGTH EFFECTS IN DIGITAL FILTERS PART A 1. Define 1 s complement form? In 1,s complement form the positive number is represented as in the sign magnitude form. To obtain the negative
More informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationBinary Floating-Point Numbers
Binary Floating-Point Numbers S exponent E significand M F=(-1) s M β E Significand M pure fraction [0, 1-ulp] or [1, 2) for β=2 Normalized form significand has no leading zeros maximum # of significant
More informationOptimizing MPC for robust and scalable integer and floating-point arithmetic
Optimizing MPC for robust and scalable integer and floating-point arithmetic Liisi Kerik * Peeter Laud * Jaak Randmets * * Cybernetica AS University of Tartu, Institute of Computer Science January 30,
More informationINTRODUCTION TO SCIENTIFIC COMPUTING - DRAFT -](DRAFT)
INTRODUCTION TO SCIENTIFIC COMPUTING - DRAFT -](DRAFT) Abstract These are notes used in the upper level undergraduate courses on numerical analysis. Since the students taking this course sequence come
More informationDSP Design Lecture 2. Fredrik Edman.
DSP Design Lecture Number representation, scaling, quantization and round-off Noise Fredrik Edman fredrik.edman@eit.lth.se Representation of Numbers Numbers is a way to use symbols to describe and model
More informationUniform Random Binary Floating Point Number Generation
Uniform Random Binary Floating Point Number Generation Prof. Dr. Thomas Morgenstern, Phone: ++49.3943-659-337, Fax: ++49.3943-659-399, tmorgenstern@hs-harz.de, Hochschule Harz, Friedrichstr. 57-59, 38855
More informationROUNDOFF ERRORS; BACKWARD STABILITY
SECTION.5 ROUNDOFF ERRORS; BACKWARD STABILITY ROUNDOFF ERROR -- error due to the finite representation (usually in floatingpoint form) of real (and complex) numers in digital computers. FLOATING-POINT
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 informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Notions of error Notes for 2016-09-02 The art of numerics is finding an approximation with a fast algorithm, a form that is easy to analyze, and an error bound Given a task, we want to engineer an approximation
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 2
2.29 Numerical Fluid Mechanics Fall 2011 Lecture 2 REVIEW Lecture 1 1. Syllabus, Goals and Objectives 2. Introduction to CFD 3. From mathematical models to numerical simulations (1D Sphere in 1D flow)
More informationON REPRESENTATIONS OF COEFFICIENTS IN IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 10(1), 57-71 (2004) ON REPRESENTATIONS OF COEFFICIENTS IN IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE ANDRZEJ MARCINIAK 1,2 AND BARBARA SZYSZKA 3 1 Poznań
More information2. Motivation and Introduction: Numerical Algorithms in CSE Basics and Applications
2. Motivation and Introduction: Numerical Algorithms in CSE Basics and Applications Numerical Programming I (for CSE), Hans-Joachim Bungartz page 1 of 46 What is Numerics? Numerical Mathematics: Part of
More information