Multi-Implicit Peer Two-Step W-Methods for Parallel Time Integration

Abstract Peer two-step W-methods are designed for integration of stiff initial value problems with parallelism across the method. The essential feature is that in each time step s ‘peer’ approximations are employed having similar properties. In fact, no primary solution variable is distinguished. Parallel implementation of these stages is easy since information from one previous time step is used only and the different linear systems may be solved simultaneously. This paper introduces a subclass having order s−1 where optimal damping for stiff problems is obtained by using different system parameters in different stages. Favourable properties of this subclass are uniform stability for realistic stepsize sequences and a superconvergence property which is proved using a polynomial collocation formulation. Numerical tests on a shared memory computer of a matrix-free implementation with Krylov methods are included.

[1]  Jason Frank,et al.  Parallel iteration of the extended backward differentiation formulas , 1999 .

[2]  L. Shampine,et al.  RKC: an explicit solver for parabolic PDEs , 1998 .

[3]  Rüdiger Weiner,et al.  Parallel Two-Step W-Methods with Peer Variables , 2004, SIAM J. Numer. Anal..

[4]  Helmut Podhaisky,et al.  Two-step W-methods and their application to MOL-systems , 2003 .

[5]  Åke Björck,et al.  Numerical Methods , 2021, Markov Renewal and Piecewise Deterministic Processes.

[6]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[7]  Rüdiger Weiner,et al.  Implicit parallel peer methods for stiff initial value problems , 2005 .

[8]  B. A. Schmitt,et al.  ROWMAP—a ROW-code with Krylov techniques for large stiff ODEs , 1997 .

[9]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[10]  J. Brandts.  [Review of: W. Hundsdorfer, J.G. Verwer (2003) Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations] , 2006 .

[11]  W. Gautschi,et al.  Lower bounds for the condition number of Vandermonde matrices , 1987 .

[12]  Parallel half-block methods for initial value problems , 2000 .

[13]  Helmut Podhaisky,et al.  Parallel Two-Step W-Methods on Singular Perturbation Problems , 2001, PPAM.

[14]  L. Brugnano,et al.  Blended implementation of block implicit methods for ODEs , 2002 .

[15]  John C. Butcher,et al.  General Linear Methods for Stiff Differential Equations , 2001 .

[16]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .