Parallel general linear methods for stiff ordinary differential and differential algebraic equations
|
|
- Silvester Spencer
- 5 years ago
- Views:
Transcription
1 ~-~H-~-~ ~ APPLIED ~ NUMERICAL MATHEMATICS ELSEVIER Applied Numerical Mathematics 17 (1995) Parallel general linear methods for stiff ordinary differential and differential algebraic equations J.C. Butcher ~"*, P. Chartier b,l Department of Mathematics The Unicersity of Auckland, Auckland, New Zealand b I.R.I.S.A., Rennes, France Abstract For numerical methods of Runge-Kutta type, the value of thc stability function at infinity, R(~), is a key factor to the solution of stiff ODEs and DAEs. In particular, methods with R(~) = 0 are especially attractive. The aim of this paper is to make some preliminary comparisons of two recently introduced approaches to the construction of L-stable parallel general linear methods, that differ essentially by the behaviour they require at infinity. 1. Introduction Since parallelism across the method makes sense only in a multi-stage environment, it is natural to consider Runge-Kutta methods. However, Iserles and Ncrsett [4] have shown that the order of a Runge-Kutta method cannot exceed s + 1 where s is the number of sequential stages. In an attempt to circumvent this barrier, van der Houwen and Sommeijer [6] have proposed a parallel iterative scheme for the nonlinear system to be solved at each step. Another possibility is to consider the wider class of "general linear methods" (GLMs). Two recent approaches to parallel GLMs are the "Generalized Backward Differentiation Formulas" [3] and DIMSIM methods of type 4 [1]. If we denote by yen-~l the r-block vector of output values at step number n - 1 (these are also the input values to step number n) and by Y the s-block vector of stage values, then one step of these method can be written as Yi[n-l] = h ~aijf(yj[n-1])jr_ ~ u i l y j,[n 1] i = l,...,s, (1) j - I ) 1 y}n]: h ~ bijf(yj[n I])_~_ ~ "'ijyj"[n- l], i = I,..., r, (2) j - 1 j = l where the initial value problem under consideration is y'(x)=f(y(x)). b) * Corresponding author. Telephone: Fax: butcher(c mat.auckland.ac.nz. l Telephone: (33) Fax: (33) chartier(a'irisa.fr /95/$ Elsevier Science B.V. All righls reserved SSDI (95 )
2 214 J.C. Butcher, P. Chartier /Applied Numerical Mathematics 17 (1995) The r output values have here a natural interpretation in the sense that y!'-11, i = 1,..., r, is assumed to be an approximation at point x, _ 1 of order p to the quantity p ;i(x) = E k = 0 In Section 2, we will present methods that are an obvious extension to type-4 methods (see [1]) and that furthermore satisfy B = Vil. Section 3 will consider "pure" type-4 methods, characterized by A = AI. In particular, we shall construct L-stable examples. Since for both methods the evaluation of the stages are independent, they are natural candidates for parallel implementation. In Section 4, preliminary numerical results will be given. Finally, further possible developments will be sketched. (3) 2. Pure A-V methods We consider here methods for which A is a nonconstant diagonal matrix, A = diag(a 1..., As). In this case, the amplification matrix, obtained by a standard linear stability analysis, is of the form M ( z ) = V + z B ( I - z A ) - I u. While the condition p ( M ( ~ ) ) = 0 is sufficient to achieve a good damping of the errors, it may seem even more desirable to require M(~) = 0. This ensures that, in the situation of differential algebraic equations, the numerical solution lies on the constraint manifold. If we arbitrarily choose U = I, this is equivalent to B = VA whenever A is non-singular, and the order conditions take a particularly simple form: Theorem 1. Let us consider a generalized type-4 method of the form [Avil / ] (4) and let p = r - 1 = s - 1. Denote by l:(x) the jth Lagrange polynomial based on the abscissae c:, j = 1...., s; the method is of order p and stage-order p iff V = L - AL', where the (i, j) elements of L and 12 are respectively given by lj(1 + c i ) and Ii(1 + c i ). (6) (5) Proof. Since the order and the stage-order do not differ, the conditions for order p = r - 1 can be written as V e cz = ( I - z A ) e <]l+c,z --t- O ( z r ), (7) where c = (c c,) T, 1[ = ( ) T, e <z = (exp(clz),..., exp(csz)) T and e (~ +<)z = (exp((1 + el)z)..., exp((1 + cs)z))) v. This implies that VP(c) = P ( 1 + c ) - A P ' ( ~ + c ) (8) for all polynomials of degree less than r - 1. Taking P ( x ) = lj(x) for j = 1..., s then gives the coefficients of V. []
3 J.C. Butcher, P. Chartier /Applied Numerical Mathematics 17 (1995) It can be noted that the usual B D F methods (Backward Differentiation Formulas) form a special subset of the set of methods described by the previous theorem. They can be obtained by taking c i = - ( s - 1) + i, i = 1,..., s, and A = d i a g ( 0,..., 0, As), where A, is chosen so as to achieve stage-order p = r and order p. However, it is a well-known fact that BDFs are not A-stable for orders greater than 2. Other choices of the ci's seem to offer only tiny improvements on the usual BDFs (see [2]). As an illustration of this fact, we give below a m e t h o d based on the abscissae c = [ - 3, ~, 1] v, with stage-order 3 and order 3, r = s = 3: l! ( {) (9) In a parallel environment and for a large problem, this m e t h o d has a computational cost similar to a B D F m e t h o d but better properties (its error constant is C = 0.22 and it is A(86.7 )-stable 1 compared to respectively - ~ and 86 for the B D F of order 3). However, it can definitely not be considered as a major competitor of BDF. This seems to suggest that A-stable m ethods of order p --- r = s and stage-order p might not exist and that good methods within that class are difficult to construct. On the other hand, it has been shown [3], that other choices of the A/'s may lead to E-stable methods up to very high orders (at least 12). The optimal choice of the Ai's and c/'s in the general situation still remains an open question. As an example, we take c = [0, 1 34 ' 1ft. A numerical search for A-stable methods with small coefficients has led to the m e t h o d l~/ _ ) (18(/ ( l (10) Numerical evaluation of the coefficients in the V matrix (B = VA) gives V = (11)
4 216 J.C. Butcher, P. Chartier /Applied Numerical Mathematics 17 (1995) If the vector of input values y["-q lies on the exact solution 37(x,_ 1) at the point xn_l, then the error after one step of method (10) is of the form with Now let h4(a y(4)(xn_l) ) + O(hS), (12) A = [ I T ~ ~ ~ J " [ T U ~ [ ~ ~ , ~ J be the normalized left eigenvector of V associated with the eigenvalue 1. The local truncation error 6 of method (10) is given by "~h4~(4)(,," (~ = h 4 ( L : T A ) y ( 4 ' ( X n - 1 ) + O(kS) = ( 1 ~ 2 1 ' " -Y ~,~n-l) + O(h5) (13) A natural way to get an estimation of this local error is to consider a linear combination of the scaled derivatives of the stage values 4 (14) E Oihf(Y,[~-lJ). i = l It is then easily seen that 4 E Oihf(Yii In-l]) = ( 1 ]'" ~ 2 ] h 4 w ( 4 ) ( y* y, ~ n - 1) + O(hS), ( 1 5 ) i = 1 if and only if E,4:10i = O, E4: loici = O, E4_10,c2 = 0 ' (16) Et loic3 = " " Solving the linear system (16) leads to I i],[i:l = ( ~ 0 (17) Thus, an asymptotically correct error estimator is 15216;90706 {hf( Y1 ["- 1]) _ 2hf(YLzn-1]) + 2hf(y~n-,]) _ hf(y4 ~"- 1])}. (18) It can be computed without any additional evaluation of f.
5 J.C. Butcher, P. Chattier/Applied Numerical Mathematics 17 (1995) Pure type-4 methods By a "pure type-4 method", we will mean a method in which r = s and the matrix A is of the form AI. For reasons of simplicity and for good stability at 0, we will assume that the matrix V is of the form V = ]u x, where vt]i = 1. This assumption will guarantee that o-(v)= {1, 0}, as for many methods such as those of Adams type, but will also go a little further. In fact the Jordan canonical form for V is diag(1, 0,..., 0). It is convenient to arbitrarily choose U as the unit matrix. This assumption means that the r quantities computed in a step and on which the stages of the following step depend are precisely equal to these later stages minus h times the corresponding stage derivatives. A standard linear stability analysis of a DIMSIM leads to the stability matrix M = V + z B ( I - z A ) ~U. With the special choice we have made for A and U, this simplifies to M ( z ) = V B, 1 - A z Z with characteristic polynomial of the form d e t ( I w - M ( z ) ) = P ( w, z) = (1 -Az)--~fi(w, z), where the polynomial /3(w, z) is given by / 5 ( w, z ) = ~ C k ( z ) ( 1 - A z ) ~ ~w ~-k, k = () and the polynomials C k have degrees k for k = 1, 2...., s, with Co(z) = 1. On the assumption that V has rank 1, many coefficients in these polynomials vanish. In fact, C~(z) = O(z~-l), for k = 1, 2,..., s. The remaining coefficients in C 1, C2....,C s are functions of the elements in the vector v and can be adjusted to fine tune the coefficients in the stability polynomial. One natural choice amongst the many available is to require that C k has degree only k - 1, for k = 1, 2,..., s. This will imply that the spectral radius of M is zero when z = ~. Thus C k = - a k z k ', where ak, k = 1, 2,..., s, as well as A are completely determined by the order conditions. Having determined the stability polynomial in this way, there is still some freedom in the 1 choice of the vector of abscissae c. If for r = s = 3 we choose c = [0, 3, 1] T, the method takes the form [) 0 1 (I 0 A o (~ 1 o ~ ( JA + 276A 2 ) [) A ~( A 2 ) J,(83 208A + I(12A 2) ~,(79 174A +72A 2) ~( A - 72A 2) 6~( A +72A 2) " ] ~ ( _ A + I I O 4 A 2 ) 1( 199~497A 240A 2) /4( A_816A2) 1(79 174A~72A2) 1( A 72A 2) 61( A A 2 ) ] [ ~(ii0 280A+138A -~) ~a( A 241)A 2 ) J6( A-+2()4A 2 ) 1(79 174A+72A2) 1( A 72A 2) 6~( A A 2 ) J where A is a zero of the polynomial A - 18A 2 + 3A 3. An analysis of the stability in this (19)
6 2 ] ~ J.C. Butcher, P. Chartier/Applied Numerical Mathematics 17 (1995) shows that this method is A-stable (and hence L-stable) if A is chosen as the zero approximately equal to The local truncation error 6 of this scheme can be estimated as described in the previous section (see Eqs. (12) and (13)) = 2~4( A - 168A2)h4y )( x,,_ l) + O(hS). (20) Getting an error estimator is however slightly more complicated, since only three stages are at our disposal. Hence, we seek a linear combination of the scaled derivatives hf(y}~-ll), i = 1, 2, 3, of the stage values and of the incoming values y!~-l], i = 1, 2, 3, such that 3 3 y ~ - - [n O, Yi 1] [n - + ~_,Oihf(Y~ q ) = ~ ( A A 2 ) h a y ( 4 ) ( X n _ l ) + O ( h 5 ). i=1 i=1 (21) Assuming that the local truncation error vector has been chosen so as to lie along the preconsistency vector, a Taylor expansion of the left-hand side gives the following system of equations El} l~)i = O, E;'= O, Jr- E~ lljktei= J i,j = O, 2~_ loici -~- 2 i = I(~)zE) = lbi,j(cj -- 1) = O, z~= ~o~: + z ~ _, ~,, z ~ =, b,,~(~ - 1) 2 = o, (22) zl,ei~;~ + z~=,,~,zi~= ~b/j(~- 1) ~ = 3!. 2 4 ( A - 168A2). In order to get a reliable estimation for small values of h, we furthermore impose that 3 0 = Y~j= 1 j 0. A few calculations eventually lead to the following values of the ~bi's and 0 i s, A A 2 3(1-2 A ) _, _168 &2 = 3(1-2A) ' 4'3-1 ¼(-1-4 A ) ] 2A ¼(1-4 A ) (23) so that an asymptotically correct error estimator reads as A A 2 ) 3 ( 1-2A) {-Yln 1]+2y~,, q_y~n l]+¼(_l_4a)hf(y~n 1]) + 2Ahf(Y[2 n 1]) + ¼(1-4A)hf(Y[3n-a])}. (24)
7 J.C. Butcher, P. Chartier /Applied Numerical Mathematics 1 7 (1995) Unfortunately, this method is unattractive from another point of view. Numerical evaluation of the coefficients in the B and V matrices gives their approximate values as ] B = /, (25) ] ] V = /. (26) ] The large values of many of the coefficients in these matrices suggest that this method is unsatisfactory because of the large roundoff errors that might result from its use. It is possible that this might be remedied by other choices of the c vector. We now consider an alternative method based on the choice of the c vector c = [ - 1, 0, 1] T. With the same value of A, the new m e t h o d has the following tableau 0 0,~ A " ~ ( A + 156A 2) {( A 6(Ia 2) 1~(65-166A + 84A2) ] 16(57-149A + 78A 2) ~( A 60A 2) ~ ( A A 2) [ t ( a a 2) ½( a - 60a z) ~ ( ,~ a z) Numerical values for this case are as follows {( A A 2) ½( A 18A 2) I~( 16 39A + 18A2) ~3( A 18h 2) ~ ( A + 18A 2 ) ½ ( A A 2 ) 0 ~6(34 57A + 18A 2) 1 t~,(34 57A + 18A 2) ~(34-57A + 18A 2) (27) B = V = ] , ] To estimate errors for this method, one can use the following formula, obtained by solving the system of equations (22) where the right-hand side of the last equation has been replaced by the appropriate value, (28) (29) A - 168A 2 ) { _ y ~ 11+2Y~2 n 11_y~n-11+3( _ 1-2 A ) f i f ( y l P ' ) + 2ahf(Y2 t"-'l) + ½(1-2a)hf(Y~"-l])}. (30) 4. N u m e r i c a l t e s t s Methods (10) and (27) were tested on the following problems:
8 220 J.C. Butcher, P. Chartier /Applied Numerical Mathematics 17 (1995) Prothero-Robinson Problem (Stiff ODE). y ' ( x ) = ( y ( x ) - s i n ( x ) ) / E + c o s ( x ), Y0-- 0, x ~ [0, 11], (31) with exact solution y(x) = sin(x). We chose here e = Modified Kaps Problem (Index-2 DAE). y ' i ( x ) = - ( 2 + e - 1 ) Y l ( X ) + e - l y ~ ( x ), y l ( O ) = l y'z(x) = - e 1-z'~'x), Y2(0) = 1, 0 =Y2(X) 2 - y, ( x ) e ~(1-c s~y~-y2(~)2)), z(o) = 1 with exact solution y, ( x ) = e -2~ y 2 ( x ) = e -~, z ( x ) = ~ + x In our experiments, e was set to 10-2 and a to Nonlinear Index-3 Problem (see [6]). y',(x) = 2yl(x)yz(x)z,(X)Zz(X), y,(o) = 1, y'2(x) = - Y l ( X ) Y 2 ( X ) Z 2 ( X ), Y2(O) = 1, Z'I(X)= ( Y l ( X ) Y 2 ( X ) +Z1(X)Z2(X))U(X), zi(o)= 1, Z'z(X ) = -Yl(X)y2(x)z~(x)u2(x), z2(0 ) = 1, O = y, ( x ) y 2 ( x ) - 1, u(o) = 1, x ~ [0, 1], (32) (33) x ~ [0, 11, (34) with exact solution y, ( x ) = e 2x, Y 2 ( x ) = e ~, z ~ ( x ) = e 2~, z 2 ( x ) = e -x, u ( x ) = e x (35) For comparison purposes, we have also included in our tests the R u n g e - K u t t a method of order 3 based on the Radau IIA abscissae. The systems have been integrated with fixed h. Numerical results for these problems are shown graphically in Fig. 1. The observed orders of convergence for the R u n g e - K u t t a method (2 for the y- and z-components and 1 for the u-component) are in agreement with those predicted by [5, T h e o r e m 6.1, Chapter IV]. 5. Conclusion Both methods solve the three test problems without difficulty and apparently without the order reduction found for the Radau IIA method. Our hope is to understand how the quality of the method depends on the spacing of the ci's and the spectral properties of M(0) = V and M(oo) = V - B A - 1 U. Our preliminary results, in which (10) performs especially well, suggests that M(oo) -- 0 is a desirable property. Method (10) also seems to benefit by having s > p, which results into an accuracy comparable to the Radau IIA method. This has the additional
9 2. C. Butcher, P. Chartier /Applied Numerical Mathematics 17 (1995) Prothero-Robineon Problem -~/ i '! i '!,. y ~ 2[ Modified Kaps Problem : y-cemlaenerit i o - > =,. F ~. m, ~ a..... ~ :..... ~ -4 =, ~ a d a ~ l ~ _>,, +-" ~ = ;..'-* [ ~"'- 5. r I... ~ "..... i '"-,,.,t,-,'., /2 -i -&.0',-&.6 o 02,2.,.3, -3' -2.,.,'5 Ioglo(h) leglo(h) Modified Kaps Problem : z-component Nonlinear index 3 problem : y-component 2 i -1!, ~ > F~sl m o t t ~ d -3 o => F+~t method ~ - "=> RaeautlA3 o :, r " ; ~ : : / :\ / : o / w 4" ~ " 2 ", ~ :..,K i ~ -3.5, -3 i -25 i -2 i -15 4, i i -2 I -1.5 I leg1d{h) ~ l o l h ) Nominear index 3 problem : z-componem Nonlinear index 3 problem : u-component -1 2 [ I -3" O => FirSl method ~ : => Second fnelilod ~. ~ / /, / -3 >= N ~ ~"..... ~' i / / 4,, i i I -3 ~2.5-2 iogl0(h), Fig. 1., J i i L Iogl0(h)
10 222 Z C. Butcher, P. Chartier /Applied Numerical Mathematics 17 (1995) advantage to allow an easy and reliable error estimate. On the other hand, method (27) shows a more regular behaviour and this is likely to result from having V of rank one. Further work on these questions is proceeding. References [1] J.C. Butcher, Diagonally-implicit multi-stage integration methods, Appl. Numer. Math. 11 (1993) [2] J.C. Butcher and P. Chartier, The construction of DIMSIMs for stiff ODEs and DAEs, Report Series, University of Auckland, Auckland, New Zealand (1994). [3] P. Chartier, L-stable parallel One-block methods for ordinary differential equations, SlAM J. Numer. Anal 31 (1994) [4] A. Iserles and S.P. Ncrsett, On the theory of parallel Runge-Kutta methods, IMA J. Numer. Anal. 10 (1990) [5] L.O. Jay, M6thodes du type Runge-Kutta pour des 6quations diff6rentielles alg6briques d'index 3 avec des applications aux syst~mes hamiltoniens, Thesis, University of Geneva, Geneva, Switzerland (1994). [6] P.J. van der Houwen and B.P. Sommeijer, Parallel iteration of high-order Runge-Kutta methods with stepsize control, J. Comput. Appl. Math. 29 (1990)
Scientific 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 informationImproved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3
Improved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3 J.H. Verner August 3, 2004 Abstract. In [5], Jackiewicz and Verner derived formulas for, and tested the implementation of two-step
More informationNumerical solution of stiff ODEs using second derivative general linear methods
Numerical solution of stiff ODEs using second derivative general linear methods A. Abdi and G. Hojjati University of Tabriz, Tabriz, Iran SciCADE 2011, Toronto 0 5 10 15 20 25 30 35 40 1 Contents 1 Introduction
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 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 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 informationOn the Diagonal Approximation of Full Matrices
On the Diagonal Approximation of Full Matrices Walter M. Lioen CWI P.O. Box 94079, 090 GB Amsterdam, The Netherlands ABSTRACT In this paper the construction of diagonal matrices, in some
More informationA New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations
A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations Rei-Wei Song and Ming-Gong Lee* d09440@chu.edu.tw, mglee@chu.edu.tw * Department of Applied Mathematics/
More informationarxiv: v1 [math.na] 15 Nov 2008
arxiv:0811.2481v1 [math.na] 15 Nov 2008 A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions Abstract D. F. Papadopoulos a, Z. A. Anastassi
More informationParallel Methods for ODEs
Parallel Methods for ODEs Levels of parallelism There are a number of levels of parallelism that are possible within a program to numerically solve ODEs. An obvious place to start is with manual code restructuring
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
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 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 informationEfficient Simulation of Physical System Models Using Inlined Implicit Runge Kutta Algorithms
Efficient Simulation of Physical System Models Using Inlined Implicit Runge Kutta Algorithms by Vicha Treeaporn A Thesis Submitted to the Faculty of the Department of Electrical and Computer Engineering
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationExplicit General Linear Methods with quadratic stability
Explicit with quadratic stability Angelamaria Cardone a, Zdzislaw Jackiewicz b ancardone@unisa.it, jackiewi@math.la.asu.edu a Dipartimento di Matematica e Informatica, Università degli Studi di Salerno,
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 information3 (Maths) Linear Algebra
3 (Maths) Linear Algebra References: Simon and Blume, chapters 6 to 11, 16 and 23; Pemberton and Rau, chapters 11 to 13 and 25; Sundaram, sections 1.3 and 1.5. The methods and concepts of linear algebra
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationJim Lambers MAT 772 Fall Semester Lecture 21 Notes
Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving
More informationSOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS
SOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS ERNST HAIRER AND PIERRE LEONE Abstract. We prove that to every rational function R(z) satisfying R( z)r(z) = 1, there exists a symplectic Runge-Kutta method
More informationA NOTE ON EXPLICIT THREE-DERIVATIVE RUNGE-KUTTA METHODS (ThDRK)
BULLETIN OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN 303-4874 (p), ISSN (o) 303-4955 www.imvibl.org / JOURNALS / BULLETIN Vol. 5(015), 65-7 Former BULLETIN OF THE SOCIETY OF MATHEMATICIANS
More informationNumerical solution of stiff systems of differential equations arising from chemical reactions
Iranian Journal of Numerical Analysis and Optimization Vol 4, No. 1, (214), pp 25-39 Numerical solution of stiff systems of differential equations arising from chemical reactions G. Hojjati, A. Abdi, F.
More informationCS 450 Numerical Analysis. Chapter 5: Nonlinear 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 informationFifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations
Australian Journal of Basic and Applied Sciences, 6(3): 9-5, 22 ISSN 99-88 Fifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations Faranak Rabiei, Fudziah Ismail Department
More informationJordan Canonical Form
Jordan Canonical Form Massoud Malek Jordan normal form or Jordan canonical form (named in honor of Camille Jordan) shows that by changing the basis, a given square matrix M can be transformed into a certain
More informationHere is an example of a block diagonal matrix with Jordan Blocks on the diagonal: J
Class Notes 4: THE SPECTRAL RADIUS, NORM CONVERGENCE AND SOR. Math 639d Due Date: Feb. 7 (updated: February 5, 2018) In the first part of this week s reading, we will prove Theorem 2 of the previous class.
More informationStarting Methods for Two-Step Runge Kutta Methods of Stage-Order 3 and Order 6
Cambridge International Science Publishing Cambridge CB1 6AZ Great Britain Journal of Computational Methods in Sciences and Engineering vol. 2, no. 3, 2, pp. 1 3 ISSN 1472 7978 Starting Methods for Two-Step
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationPRECONDITIONING AND PARALLEL IMPLEMENTATION OF IMPLICIT RUNGE-KUTTA METHODS.
PRECONDITIONING AND PARALLEL IMPLEMENTATION OF IMPLICIT RUNGE-KUTTA METHODS. LAURENT O. JAY Abstract. A major problem in obtaining an efficient implementation of fully implicit Runge- Kutta IRK methods
More informationDownloaded 12/15/15 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 40, No. 4, pp. 56 537 c 2002 Society for Industrial and Applied Mathematics PREDICTOR-CORRECTOR METHODS OF RUNGE KUTTA TYPE FOR STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE
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 informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
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 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 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 informationVirtual University of Pakistan
Virtual University of Pakistan File Version v.0.0 Prepared For: Final Term Note: Use Table Of Content to view the Topics, In PDF(Portable Document Format) format, you can check Bookmarks menu Disclaimer:
More informationAlmost Runge-Kutta Methods Of Orders Up To Five
Almost Runge-Kutta Methods Of Orders Up To Five ABRAHAM OCHOCHE Department of Information Technology, Federal University of Technology, Minna, Niger State PMB 65, Nigeria abochoche@gmailcom PETER NDAJAH
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 informationEfficient implementation of Radau collocation methods
Efficient implementation of Radau collocation methods Luigi Brugnano Felice Iavernaro Cecilia Magherini arxiv:1302.1037v1 [math.na] 5 Feb 2013 December 29, 2012 Abstract In this paper we define an efficient
More information3 Matrix Algebra. 3.1 Operations on matrices
3 Matrix Algebra A matrix is a rectangular array of numbers; it is of size m n if it has m rows and n columns. A 1 n matrix is a row vector; an m 1 matrix is a column vector. For example: 1 5 3 5 3 5 8
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 informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More 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 informationParallel Numerical Methods for Ordinary Differential Equations: a Survey
Parallel Numerical Methods for Ordinary Differential Equations: a Survey Svyatoslav I. Solodushkin 1,2 and Irina F. Iumanova 1 1 Ural Federal University, Yekaterinburg, Russia 2 Krasovskii Institute of
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 informationREPORTRAPPORT. Centrum voor Wiskunde en Informatica. Analysis of Approximate Factorization in Iteration Methods
Centrum voor Wiskunde en Informatica REPORTRAPPORT Analysis of Approximate Factorization in Iteration Methods C. Eichler-Liebenow, P.J. van der Houwen, B.P. Sommeijer Modelling, Analysis and Simulation
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
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 informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
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 informationA Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital Problems
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal A Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital
More informationA New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution of Oscillatory Problems
Applied Mathematical Sciences, Vol. 1, 16, no. 44, 157-178 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/ams.16.64146 A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution
More informationSPECIALIZED RUNGE-KUTTA METHODS FOR INDEX 2 DIFFERENTIAL-ALGEBRAIC EQUATIONS
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 641 654 S 0025-5718(05)01809-0 Article electronically published on December 19, 2005 SPECIALIZED RUNGE-KUTTA METHODS FOR INDEX 2 DIFFERENTIAL-ALGEBRAIC
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-71309 email: anita.buie@gmail.com 1 . Chapter 8 Numerical Solution of
More informationOn the reliability and stability of direct explicit Runge-Kutta integrators
Global Journal of Pure and Applied Mathematics. ISSN 973-78 Volume, Number 4 (), pp. 3959-3975 Research India Publications http://www.ripublication.com/gjpam.htm On the reliability and stability of direct
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More informationON THE MATRIX EQUATION XA AX = X P
ON THE MATRIX EQUATION XA AX = X P DIETRICH BURDE Abstract We study the matrix equation XA AX = X p in M n (K) for 1 < p < n It is shown that every matrix solution X is nilpotent and that the generalized
More informationA NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN FILM FLOW PROBLEM
Indian J. Pure Appl. Math., 491): 151-167, March 218 c Indian National Science Academy DOI: 1.17/s13226-18-259-6 A NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN
More informationBlended implementation of block implicit methods for ODEs
Applied Numerical Mathematics 42 (2002) 29 45 www.elsevier.com/locate/apnum Blended implementation of block implicit methods for ODEs Luigi Brugnano, Cecilia Magherini Dipartimento di Matematica U. Dini,
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 informationAN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)
AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
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 informationMartin Luther Universität Halle Wittenberg Institut für Mathematik
Martin Luther Universität Halle Wittenberg Institut für Mathematik Implicit peer methods for large stiff ODE systems St. Beck, R. Weiner, H. Podhaisky and B. A. Schmitt Report No. 07 (2010) Editors: Professors
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
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 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 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 informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationA Note on Trapezoidal Methods for the Solution of Initial Value Problems. By A. R. Gourlay*
MATHEMATICS OF COMPUTATION, VOLUME 24, NUMBER 111, JULY, 1970 A Note on Trapezoidal Methods for the Solution of Initial Value Problems By A. R. Gourlay* Abstract. The trapezoidal rule for the numerical
More informationIntegrating-Factor-Based 2-Additive Runge Kutta methods for Advection-Reaction-Diffusion Equations
Integrating-Factor-Based 2-Additive Runge Kutta methods for Advection-Reaction-Diffusion Equations A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements
More informationSums of diagonalizable matrices
Linear Algebra and its Applications 315 (2000) 1 23 www.elsevier.com/locate/laa Sums of diagonalizable matrices J.D. Botha Department of Mathematics University of South Africa P.O. Box 392 Pretoria 0003
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 informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationA SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 3, Fall 2009 A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS SERGEY KHASHIN ABSTRACT. A new approach based on the use of new
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 informationAdditional exercises with Numerieke Analyse
Additional exercises with Numerieke Analyse March 10, 017 1. (a) Given different points x 0, x 1, x [a, b] and scalars y 0, y 1, y, z 1, show that there exists at most one polynomial p P 3 with p(x i )
More informationSymplectic Runge-Kutta methods satisfying effective order conditions
/30 Symplectic Runge-Kutta methods satisfying effective order conditions Gulshad Imran John Butcher (supervisor) University of Auckland 5 July, 20 International Conference on Scientific Computation and
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The
More informationIntroduction to Applied Linear Algebra with MATLAB
Sigam Series in Applied Mathematics Volume 7 Rizwan Butt Introduction to Applied Linear Algebra with MATLAB Heldermann Verlag Contents Number Systems and Errors 1 1.1 Introduction 1 1.2 Number Representation
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
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 informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
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 informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationTwo Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. (8) 7-77 ISSN - Two Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation T.V. Triantafyllidis,
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
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 informationControl Systems. Linear Algebra topics. L. Lanari
Control Systems Linear Algebra topics L Lanari outline basic facts about matrices eigenvalues - eigenvectors - characteristic polynomial - algebraic multiplicity eigenvalues invariance under similarity
More information