Semi-Automatic Floating-Point Implementation of Special Functions
|
|
- Edith Franklin
- 5 years ago
- Views:
Transcription
1 Semi-Automatic Floating-Point Implementation of Special Functions Christoph Lauter 1 Marc Mezzarobba 1,2 Pequan group 1 Université Paris 6 2 CNRS ARITH 22, Lyon,
2 }main() { int temp; float celsius; char repeat; char flag; do { flag='n"; do { if(flag=='n') printf("input a valid temperature :"); else printf("input a valid temperature,stupid:"); scanf("%d",&temp); flag='y'; } while (temp<0 temp >100); celsius=(5.0/9.0)*(temp-32); printf("%d degrees F is %6.2f degrees celsius\n",temp,cel printf("do you have another temperature?"); repeat=getchar(); putchar('\n'); Code Generation for Mathematical Functions explicit function expression equation black box f = log(x) f = e sin(x) Φ(f ) = 0 x f (x) double fun(double x) {... }
3 }main() { int temp; float celsius; char repeat; char flag; do { flag='n"; do { if(flag=='n') printf("input a valid temperature :"); else printf("input a valid temperature,stupid:"); scanf("%d",&temp); flag='y'; } while (temp<0 temp >100); celsius=(5.0/9.0)*(temp-32); printf("%d degrees F is %6.2f degrees celsius\n",temp,cel printf("do you have another temperature?"); repeat=getchar(); putchar('\n'); Code Generation for Mathematical Functions explicit function expression equation black box f = log(x) f = e sin(x) Φ(f ) = 0 x f (x) double fun(double x) {... }
4 Our Focus: Special Functions Ai(x) Ci(x) erf(x) J 0 (x) D-finite Functions = sol ns of linear ODEs with polynomial coeff ts p r (x) f (r) (x) + + p 1 (x) f (x) + p 0 (x) f (x) = 0 p 0, p 1,..., p r R[x] Example: f (x) = arctan(x) (x 2 + 1) f (x) + 2x f (x) = 0, [ ] f (0) f = (0) [ ] 0 1
5 }main() { int temp; float celsius; char repeat; char flag; do { flag='n"; do { if(flag=='n') printf("input a valid temperature :"); else printf("input a valid temperature,stupid:"); scanf("%d",&temp); flag='y'; } while (temp<0 temp >100); celsius=(5.0/9.0)*(temp-32); printf("%d degrees F is %6.2f degrees celsius\n",temp,cel printf("do you have another temperature?"); repeat=getchar(); putchar('\n'); Problem Statement { pr f (r) + + p 0 f = 0, f (0),..., f (r 1) (0) D R, ε > 0 double fun(double x); x D double, fun(x) f (x) f (x) ε Semi-automatic : may require implementation hints Future goal: full automation Rigorous in principle, some shortcuts in current prototype
6 differential equation Overview sets of polynomial approximations NumGfun (Maple) rigorous multiple precision solver for ODE Sollya Metalibm-lutetia (cf. previous talk by O. Kupriianova)
7 differential equation The : Taylor Series p r (x) f (r) (x) + p 0 (x) f (x) = 0 Ansatz: f (x) = f (x) = x f (x) = n=0 c n x n n=0 (n + 1)c n+1 x n c n 1 x n n=1 c n = b 0 (n)c n 1 + b 1 (n)c n 2 + f (x) = c 0 + c 1 x + + c d 1 x d 1 +
8 differential equation The : Error Bounds f (x) = x ρ d 1 n=0 c n x n + n=d c n x n } {{ }? Majorizing series : using the ODE, find a simple ˆf (x) = ĉ n x n such that c n ĉ n for all n. f (x) d 1 c n x n ĉ n ρ n! ε. n=0 n=d
9 The : Analytic Continuation differential equation What if D [ ρ, ρ]? i i
10 The : Analytic Continuation differential equation What if D [ ρ, ρ]? i [ ] f (x1 ) f (x 1 ) i x 1 = [ ] [ ] f (0) f (0)
11 The : Analytic Continuation differential equation What if D [ ρ, ρ]? i [ ] f (x1 ) f (x 1 ) i x 1 = [ ] [ ] f (0) f (0)
12 The : Analytic Continuation differential equation What if D [ ρ, ρ]? i [ ] f (x1 ) f (x 1 ) [ ] f (x2 ) f (x 2 ) i x 1 x 2 = [ = [ ] [ ] f (0) f (0) ] [ ] f (x 1 ) f (x 1 )
13 The : Analytic Continuation differential equation What if D [ ρ, ρ]? i [ ] f (x1 ) f (x 1 ) [ ] f (x2 ) f (x 2 )... i x 1 x 2 = [ = [ ] [ ] f (0) f (0) ] [ ] f (x 1 ) f (x 1 )
14 differential equation The : Economization Have: f (x) (c c d 1 x d 1 ) ε 2 x C, x 1 c c 7 x 7 + c 8 x 8 + c 9 x 9 c (1) c (1) 7 x7 + c (1) 8 x8 c (2) c (2) 7 x7 c T 8(x) c T 9(x). T n 1 over [ 1, 1] f (x) ( c (k) c (k) d 1 k xk ) ε x [ 1, 1]
15 The : Tight Approximation differential equation : D = D i x D i p i (x) f (x) ε deg p i = O(500) : D = Di x Di pi (x) f (x) 1 ε deg pi = O(10)
16 The : Tight Approximation differential equation : leveraging Metalibm Reuse domain splitting algorithm sketched in last talk D = D i Reuse minimax implementation with relative error bounds x Di pi (x) f (x) 1 ε Addressing issues Pure black-box interface too slow Zeros of f in domain D?
17 differential equation The : FP Polynomials Implementation need: p i F[x] Leverage existing FP minimax heuristics
18 differential equation The : FP Polynomials Implementation need: p i F[x] Leverage existing FP minimax heuristics unless f has a zero in the domain
19 differential equation The : Zeros of f Suppose f (c) = 0 for c D We want pi (x) f (x) 1 ε We need pi (c) = 0 No if pi F[x] but c R\F
20 differential equation The : Zeros of f Suppose f (c) = 0 for c D Actually we want x F. pi (x) f (x) 1 ε p i (x) never needs to be 0 exactly Okay but we need to compute pi Express c as a symbol, evaluated with Newton-Raphson on f Give f (x c)/x k to minimax algorithm, yielding q Deduce pi from x k q(x)
21 differential equation The : Evaluation scheme Reuse Metalibm-Lutetia to generate Horner scheme We could also use Metalibm-Lugdunum here
22 Examples: erfc Consider erfc(x) = 1 2 x e t2 dt π over D = [ 2; 2] at ε = x Describe erfc with f (x) + 2xf (x) = 0, f (0) = 1, f (0) = 2 π Generated has domain split into 16 subdomains Code runs in 110 to 350 cycles, libm exponential in 80 cycles Code is bitwise the same if we take MPFR s erfc instead of ODE 1e-19 8e-20 6e-20 4e-20 2e e-20-4e-20-6e-20-8e-20-1e e
23 Examples: J I I I x 30 I Consider Bessel function J0 over D = [0.5; 42] at ε = 2 45 I J0 is given as 40 xf 00 (x) + f 0 (x) + xf (x) = 0, Code generation takes about 30 minutes Generator has to cope with 13 zeros in the domain Code runs in 60 to 500 cycles, libm exponential in 80 cycles f (x) 1, x 0 1.5e-14 1e-14 5e e-15-1e e-14-2e
24 Conclusion and Outlook Don t. Have your special functions generated! Process starts from the very basic definition Only small domains handled fully automatically Full range implementions require manual intervention -backend interface not satisfactory Code generation performance unpredictable
25 Thanks! Thank you for your attention! Questions?
Computation 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 informationDesigning a Correct Numerical Algorithm
Intro Implem Errors Sollya Gappa Norm Conc Christoph Lauter Guillaume Melquiond March 27, 2013 Intro Implem Errors Sollya Gappa Norm Conc Outline 1 Introduction 2 Implementation theory 3 Error analysis
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 informationEfficient and accurate computation of upper bounds of approximation errors
Laboratoire de l Informatique du Parallélisme École Normale Supérieure de Lyon Unité Mixte de Recherche CNRS-INRIA-ENS LYON-UCBL n o 5668 Efficient and accurate computation of upper bounds of approximation
More informationA Branch-and-Bound Algorithm for Unconstrained Global Optimization
SCAN 2010, Lyon, September 27 30, 2010 1 / 18 A Branch-and-Bound Algorithm for Unconstrained Global Optimization Laurent Granvilliers and Alexandre Goldsztejn Université de Nantes LINA CNRS Interval-based
More informationTable-Based Polynomials for Fast Hardware Function Evaluation
ASAP 05 Table-Based Polynomials for Fast Hardware Function Evaluation Jérémie Detrey Florent de Dinechin Projet Arénaire LIP UMR CNRS ENS Lyon UCB Lyon INRIA 5668 http://www.ens-lyon.fr/lip/arenaire/ CENTRE
More informationSolutions to Exercises, Section 2.5
Instructor s Solutions Manual, Section 2.5 Exercise 1 Solutions to Exercises, Section 2.5 For Exercises 1 4, write the domain of the given function r as a union of intervals. 1. r(x) 5x3 12x 2 + 13 x 2
More informationComputing Machine-Efficient Polynomial Approximations
Computing Machine-Efficient Polynomial Approximations N. Brisebarre, S. Chevillard, G. Hanrot, J.-M. Muller, D. Stehlé, A. Tisserand and S. Torres Arénaire, LIP, É.N.S. Lyon Journées du GDR et du réseau
More informationAlgorithmic Tools for the Asymptotics of Linear Recurrences
Algorithmic Tools for the Asymptotics of Linear Recurrences Bruno Salvy Inria & ENS de Lyon Computer Algebra in Combinatorics, Schrödinger Institute, Vienna, Nov. 2017 Motivation p 0 (n)a n+k + + p k (n)a
More informationSecond Order Function Approximation Using a Single Multiplication on FPGAs
FPL 04 Second Order Function Approximation Using a Single Multiplication on FPGAs Jérémie Detrey Florent de Dinechin Projet Arénaire LIP UMR CNRS ENS Lyon UCB Lyon INRIA 5668 http://www.ens-lyon.fr/lip/arenaire/
More informationMA 8019: Numerical Analysis I Solution of Nonlinear Equations
MA 8019: Numerical Analysis I Solution of Nonlinear Equations Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationNumerical Methods. V. Leclère May 15, x R n
Numerical Methods V. Leclère May 15, 2018 1 Some optimization algorithms Consider the unconstrained optimization problem min f(x). (1) x R n A descent direction algorithm is an algorithm that construct
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationFinding Hyperexponential Solutions of Linear ODEs by Numerical Evaluation
Finding Hyperexponential Solutions of Linear ODEs by Numerical Evaluation Fredrik Johansson* 1 Manuel Kauers* 1 Marc Mezzarobba* 1,2 *RISC-Linz ISSAC 2013, Boston, MA 1 Supported by the Austrian Science
More information1 Lecture 8: Interpolating polynomials.
1 Lecture 8: Interpolating polynomials. 1.1 Horner s method Before turning to the main idea of this part of the course, we consider how to evaluate a polynomial. Recall that a polynomial is an expression
More informationNonlinear Control as Program Synthesis (A Starter)
Nonlinear Control as Program Synthesis (A Starter) Sicun Gao MIT December 15, 2014 Preliminaries Definition (L RF ) L RF is the first-order language over the reals that allows arbitrary numerically computable
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
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 informationarxiv: v1 [cs.sc] 7 Jul 2016
Rigorous Multiple-Precision Evaluation of D-Finite Functions in SageMath Marc Mezzarobba CNRS, LIP6, Université Pierre et Marie Curie, Paris, France marc@mezzarobba.net, http://marc.mezzarobba.net/ arxiv:1607.01967v1
More informationSection 5.8. Taylor Series
Difference Equations to Differential Equations Section 5.8 Taylor Series In this section we will put together much of the work of Sections 5.-5.7 in the context of a discussion of Taylor series. We begin
More informationMAT 311 Midterm #1 Show your work! 1. The existence and uniqueness theorem says that, given a point (x 0, y 0 ) the ODE. y = (1 x 2 y 2 ) 1/3
MAT 3 Midterm # Show your work!. The existence and uniqueness theorem says that, given a point (x 0, y 0 ) the ODE y = ( x 2 y 2 ) /3 has a unique (local) solution with initial condition y(x 0 ) = y 0
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 9, 2018 Outline 1 Extremal Values 2
More informationWebAssign Lesson 6-3 Taylor Series (Homework)
WebAssign Lesson 6-3 Taylor Series (Homework) Current Score : / 56 Due : Tuesday, August 5 204 0:59 AM MDT Jaimos Skriletz Math 75, section 3, Summer 2 204 Instructor: Jaimos Skriletz. /4 points Consider
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 2017 Outline Extremal Values The
More informationTOPIC 3. Taylor polynomials. Mathematica code. Here is some basic mathematica code for plotting functions.
TOPIC 3 Taylor polynomials Main ideas. Linear approximating functions: Review Approximating polynomials Key formulas: P n (x) =a 0 + a (x x )+ + a n (x x ) n P n (x + x) =a 0 + a ( x)+ + a n ( x) n where
More informationQuestions Pool. Amnon Ta-Shma and Dean Doron. January 2, Make sure you know how to solve. Do not submit.
Questions Pool Amnon Ta-Shma and Dean Doron January 2, 2017 General guidelines The questions fall into several categories: (Know). (Mandatory). (Bonus). Make sure you know how to solve. Do not submit.
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationTwo hours UNIVERSITY OF MANCHESTER. 21 January
Two hours MATH20111 UNIVERSITY OF MANCHESTER REAL ANALYSIS 21 January 2015 09.45 11.45 Answer ALL SIX questions in Section A (50 marks in total). Answer TWO of the THREE questions in Section B (30 marks
More informationAlgebra II: Strand 5. Power, Polynomial, and Rational Functions; Topic 3. Rational Functions; Task 5.3.1
TASK 5.3.: RATIONAL FUNCTIONS AND THEIR RECIPROCALS Solutions Rational functions appear frequently in business, science, engineering, and medical applications. This activity explores some aspects of rational
More informationSolution of Algebric & Transcendental Equations
Page15 Solution of Algebric & Transcendental Equations Contents: o Introduction o Evaluation of Polynomials by Horner s Method o Methods of solving non linear equations o Bracketing Methods o Bisection
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationSolving Third Order Linear Differential Equations in Terms of Second Order Equations
Solving Third Order Linear Differential Equations in Terms of Second Order Equations Mark van Hoeij (Florida State University) ISSAC 2007 Talk presented by: George Labahn (University of Waterloo) Notations.
More information18.01 EXERCISES. Unit 3. Integration. 3A. Differentials, indefinite integration. 3A-1 Compute the differentials df(x) of the following functions.
8. EXERCISES Unit 3. Integration 3A. Differentials, indefinite integration 3A- Compute the differentials df(x) of the following functions. a) d(x 7 + sin ) b) d x c) d(x 8x + 6) d) d(e 3x sin x) e) Express
More informationReview: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.
Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:
More informationPower Series Solutions We use power series to solve second order differential equations
Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations
More informationCHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0:
TAYLOR SERIES Construct a polynomial with the following behavior at x = 0: CHALLENGE! P( x) = a + ax+ ax + ax + ax 2 3 4 0 1 2 3 4 P(0) = 1 P (0) = 2 P (0) = 3 P (0) = 4 P (4) (0) = 5 Sounds hard right?
More informationMa 530 Power Series II
Ma 530 Power Series II Please note that there is material on power series at Visual Calculus. Some of this material was used as part of the presentation of the topics that follow. Operations on Power Series
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 3 Lecture 3 3.1 General remarks March 4, 2018 This
More informationINTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES PHILIP FOTH 1. Cauchy s Formula and Cauchy s Theorem 1. Suppose that γ is a piecewise smooth positively ( counterclockwise ) oriented simple closed
More informationConvergence of sequences and series
Convergence of sequences and series A sequence f is a map from N the positive integers to a set. We often write the map outputs as f n rather than f(n). Often we just list the outputs in order and leave
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More information1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy
.. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution
More informationPolynomial Review Problems
Polynomial Review Problems 1. Find polynomial function formulas that could fit each of these graphs. Remember that you will need to determine the value of the leading coefficient. The point (0,-3) is on
More informationMath WW09 Solutions November 24, 2008
Math 352- WW09 Solutions November 24, 2008 Assigned problems: 8.7 0, 6, ww 4; 8.8 32, ww 5, ww 6 Always read through the solution sets even if your answer was correct. Note that like many of the integrals
More informationThe final is cumulative, but with more emphasis on chapters 3 and 4. There will be two parts.
Math 141 Review for Final The final is cumulative, but with more emphasis on chapters 3 and 4. There will be two parts. Part 1 (no calculator) graphing (polynomial, rational, linear, exponential, and logarithmic
More informationImplementation of an αbb-type underestimator in the SGO-algorithm
Implementation of an αbb-type underestimator in the SGO-algorithm Process Design & Systems Engineering November 3, 2010 Refining without branching Could the αbb underestimator be used without an explicit
More informationMath 227 Sample Final Examination 1. Name (print) Name (sign) Bing ID number
Math 227 Sample Final Examination 1 Name (print) Name (sign) Bing ID number (Your instructor may check your ID during or after the test) No books, notes, or electronic devices (calculators, cell phones,
More information3 Interpolation and Polynomial Approximation
CHAPTER 3 Interpolation and Polynomial Approximation Introduction A census of the population of the United States is taken every 10 years. The following table lists the population, in thousands of people,
More informationAlgorithmic Tools for the Asymptotics of Diagonals
Algorithmic Tools for the Asymptotics of Diagonals Bruno Salvy Inria & ENS de Lyon Lattice walks at the Interface of Algebra, Analysis and Combinatorics September 19, 2017 Asymptotics & Univariate Generating
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
More informationWeek 1: need to know. November 14, / 20
Week 1: need to know How to find domains and ranges, operations on functions (addition, subtraction, multiplication, division, composition), behaviors of functions (even/odd/ increasing/decreasing), library
More informationMath 2: Algebra 2, Geometry and Statistics Ms. Sheppard-Brick Chapter 4 Test Review
Chapter 4 Test Review Students will be able to (SWBAT): Write an explicit and a recursive function rule for a linear table of values. Write an explicit function rule for a quadratic table of values. Determine
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationHomework 1 Solutions
Homework 1 Solutions 1. Consider the function f(x) = e x. (a) Derive the nth Taylor polynomial p n (x) as well as the remainder term R n (x) for the function f(x), expanded about the point x = 0. Let s
More informationReducibility in polynomialtime computable analysis. Akitoshi Kawamura (U Tokyo) Wadern, Germany, September 24, 2015
Reducibility in polynomialtime computable analysis Akitoshi Kawamura (U Tokyo) Wadern, Germany, September 24, 2015 Computable Analysis Applying computability / complexity theory to problems involving real
More information1. Pace yourself. Make sure you write something on every problem to get partial credit. 2. If you need more space, use the back of the previous page.
***THIS TIME I DECIDED TO WRITE A LOT OF EXTRA PROBLEMS TO GIVE MORE PRACTICE. The actual midterm will have about 6 problems. If you want to practice something with approximately the same length as the
More informationMath 4310 Solutions to homework 7 Due 10/27/16
Math 4310 Solutions to homework 7 Due 10/27/16 1. Find the gcd of x 3 + x 2 + x + 1 and x 5 + 2x 3 + x 2 + x + 1 in Rx. Use the Euclidean algorithm: x 5 + 2x 3 + x 2 + x + 1 = (x 3 + x 2 + x + 1)(x 2 x
More informationMATH 103 Pre-Calculus Mathematics Test #3 Fall 2008 Dr. McCloskey Sample Solutions
MATH 103 Pre-Calculus Mathematics Test #3 Fall 008 Dr. McCloskey Sample Solutions 1. Let P (x) = 3x 4 + x 3 x + and D(x) = x + x 1. Find polynomials Q(x) and R(x) such that P (x) = Q(x) D(x) + R(x). (That
More informationODE Final exam - Solutions
ODE Final exam - Solutions May 3, 018 1 Computational questions (30 For all the following ODE s with given initial condition, find the expression of the solution as a function of the time variable t You
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 informationNewton-Raphson. Relies on the Taylor expansion f(x + δ) = f(x) + δ f (x) + ½ δ 2 f (x) +..
2008 Lecture 7 starts here Newton-Raphson When the derivative of f(x) is known, and when f(x) is well behaved, the celebrated (and ancient) Newton- Raphson method gives the fastest convergence of all (
More informationContinued Fractions for Special Functions: Handbook and Software
Continued Fractions for Special Functions: Handbook and Software Team: A. Cuyt, V.B. Petersen, J. Van Deun, H. Waadeland, F. Backeljauw, S. Becuwe, M. Colman, W.B. Jones Universiteit Antwerpen Sør-Trøndelag
More information, applyingl Hospital s Rule again x 0 2 cos(x) xsinx
Lecture 3 We give a couple examples of using L Hospital s Rule: Example 3.. [ (a) Compute x 0 sin(x) x. To put this into a form for L Hospital s Rule we first put it over a common denominator [ x 0 sin(x)
More information11.10a Taylor and Maclaurin Series
11.10a 1 11.10a Taylor and Maclaurin Series Let y = f(x) be a differentiable function at x = a. In first semester calculus we saw that (1) f(x) f(a)+f (a)(x a), for all x near a The right-hand side of
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 informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationLinear Differential Equations as a Data-Structure
Linear Differential Equations as a Data-Structure Bruno Salvy Inria & ENS de Lyon FoCM, July 14, 2017 Computer Algebra Effective mathematics: what can we compute exactly? And complexity: how fast? (also,
More information1. First-order ODE s
18.03 EXERCISES 1. First-order ODE s 1A. Introduction; Separation of Variables 1A-1. Verif that each of the following ODE s has the indicated solutions (c i,a are constants): a) 2 + = 0, = c 1 e x +c 2
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
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 informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationTheory and Practice of Symbolic Integration in Maple
Theory and Practice of Symbolic Integration in Maple John May Senior Developer Mathematical Software September 2010 Introduction When the talking version of a popular doll said Math is Hard!, she was almost
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationMath 12: Discrete Dynamical Systems Homework
Math 12: Discrete Dynamical Systems Homework Department of Mathematics, Harvey Mudd College Several of these problems will require computational software to help build our insight about discrete dynamical
More informationGlobal optimization, polynomial optimization, polynomial system solving, real
PROBABILISTIC ALGORITHM FOR POLYNOMIAL OPTIMIZATION OVER A REAL ALGEBRAIC SET AURÉLIEN GREUET AND MOHAB SAFEY EL DIN Abstract. Let f, f 1,..., f s be n-variate polynomials with rational coefficients of
More informationNumerical Analysis Homework, Problems
Numerical Analysis Homework, Problems. Modify the following program to return the sum, the product and the difference of the two arguments #include #include mainint argc, char **argv)
More informationOrder of convergence
Order of convergence Linear and Quadratic Order of convergence Computing square root with Newton s Method Given a > 0, p def = a is positive root of equation Newton s Method p k+1 = p k p2 k a 2p k = 1
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 informationProving Unsatisfiability in Non-linear Arithmetic by Duality
Proving Unsatisfiability in Non-linear Arithmetic by Duality [work in progress] Daniel Larraz, Albert Oliveras, Enric Rodríguez-Carbonell and Albert Rubio Universitat Politècnica de Catalunya, Barcelona,
More informationTaylor Series. richard/math230 These notes are taken from Calculus Vol I, by Tom M. Apostol,
Taylor Series Professor Richard Blecksmith richard@math.niu.edu Dept. of Mathematical Sciences Northern Illinois University http://math.niu.edu/ richard/math230 These notes are taken from Calculus Vol
More informationNumerical Algorithms for ODEs/DAEs (Transient Analysis)
Numerical Algorithms for ODEs/DAEs (Transient Analysis) Slide 1 Solving Differential Equation Systems d q ( x(t)) + f (x(t)) + b(t) = 0 dt DAEs: many types of solutions useful DC steady state: state no
More informationIowa State University. Instructor: Alex Roitershtein Summer Homework #5. Solutions
Math 50 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 205 Homework #5 Solutions. Let α and c be real numbers, c > 0, and f is defined
More informationD-MATH Algebra I HS18 Prof. Rahul Pandharipande. Solution 6. Unique Factorization Domains
D-MATH Algebra I HS18 Prof. Rahul Pandharipande Solution 6 Unique Factorization Domains 1. Let R be a UFD. Let that a, b R be coprime elements (that is, gcd(a, b) R ) and c R. Suppose that a c and b c.
More informationComputing Real Roots of Real Polynomials
Computing Real Roots of Real Polynomials and now For Real! Alexander Kobel Max-Planck-Institute for Informatics, Saarbrücken, Germany Fabrice Rouillier INRIA & Université Pierre et Marie Curie, Paris,
More informationMath Practice Exam 3 - solutions
Math 181 - Practice Exam 3 - solutions Problem 1 Consider the function h(x) = (9x 2 33x 25)e 3x+1. a) Find h (x). b) Find all values of x where h (x) is zero ( critical values ). c) Using the sign pattern
More informationLecture 7: Minimization or maximization of functions (Recipes Chapter 10)
Lecture 7: Minimization or maximization of functions (Recipes Chapter 10) Actively studied subject for several reasons: Commonly encountered problem: e.g. Hamilton s and Lagrange s principles, economics
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationSection 10.7 Taylor series
Section 10.7 Taylor series 1. Common Maclaurin series 2. s and approximations with Taylor polynomials 3. Multiplication and division of power series Math 126 Enhanced 10.7 Taylor Series The University
More informationCHEE 222: PROCESS DYNAMICS AND NUMERICAL METHODS
CHEE 222: PROCESS DYNAMICS AND NUMERICAL METHODS Winter 2017 Implementation of Numerical Methods via MATLAB Instructor: Xiang Li 1 Outline 1. Introduction - Command, script and function - MATLAB function
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 informationMathematics 102 Fall 1999 The formal rules of calculus The three basic rules The sum rule. The product rule. The composition rule.
Mathematics 02 Fall 999 The formal rules of calculus So far we have calculated the derivative of each function we have looked at all over again from scratch, applying what is essentially the definition
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 informationMath 106: Review for Exam II - SOLUTIONS
Math 6: Review for Exam II - SOLUTIONS INTEGRATION TIPS Substitution: usually let u a function that s inside another function, especially if du (possibly off by a multiplying constant) is also present
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 informationFourier and Partial Differential Equations
Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the
More informationSection Example Determine the Maclaurin series of f (x) = e x and its the interval of convergence.
Example Determine the Maclaurin series of f (x) = e x and its the interval of convergence. Example Determine the Maclaurin series of f (x) = e x and its the interval of convergence. f n (0)x n Recall from
More informationTaylor Series and Numerical Approximations
Taylor Series and Numerical Approximations Hilary Weller h.weller@reading.ac.uk August 7, 05 An introduction to the concept of a Taylor series and how these are used in numerical analysis to find numerical
More informationTS Method Summary. T k (x,y j 1 ) f(x j 1,y j 1 )+ 2 f (x j 1,y j 1 ) + k 1
TS Method Summary Let T k (x,y j 1 ) denote the first k +1 terms of the Taylor series expanded about the discrete approximation, (x j 1,y j 1 ), and ẑ k,j (x) be the polynomial approximation (to y(x))
More information