Parallelizing spectral deferred corrections across the method

In this paper we present two strategies to enable “parallelization across the method” for spectral deferred corrections (SDC). Using standard low-order time-stepping methods in an iterative fashion, SDC can be seen as preconditioned Picard iteration for the collocation problem. Typically, a serial Gauß–Seidel-like preconditioner is used, computing updates for each collocation node one by one. The goal of this paper is to show how this process can be parallelized, so that all collocation nodes are updated simultaneously. The first strategy aims at finding parallel preconditioners for the Picard iteration and we test three choices using four different test problems. For the second strategy we diagonalize the quadrature matrix of the collocation problem directly. In order to integrate non-linear problems we employ simplified and inexact Newton methods. Here, we estimate the speed of convergence depending on the time-step size and verify our results using a non-linear diffusion problem.

[1]  Rolf Krause,et al.  A multi-level spectral deferred correction method , 2013, BIT Numerical Mathematics.

[2]  Rolf Krause,et al.  Inexact spectral deferred corrections , 2016 .

[3]  C. Kelley Iterative Methods for Linear and Nonlinear Equations , 1987 .

[4]  Kenneth R. Jackson,et al.  The potential for parallelism in Runge-Kutta methods. Part 1: RK formulas in standard form , 1995 .

[5]  Laurent O. Jay,et al.  Inexact Simplified Newton Iterations for Implicit Runge-Kutta Methods , 2000, SIAM J. Numer. Anal..

[6]  Kevin Burrage,et al.  Parallel methods for initial value problems , 1993 .

[7]  Martin J. Gander,et al.  A Direct Solver for Time Parallelization , 2016 .

[8]  Jingfang Huang,et al.  Accelerating the convergence of spectral deferred correction methods , 2006, J. Comput. Phys..

[9]  Elizabeth L. Bouzarth,et al.  A multirate time integrator for regularized Stokeslets , 2010, J. Comput. Phys..

[10]  Ben P. Sommeijer,et al.  Analysis of parallel diagonally implicit iteration of Runge-Kutta methods , 1993 .

[11]  P. Houwen,et al.  Parallel iteration of high-order Runge-Kutta methods with stepsize control , 1990 .

[12]  Michael L. Minion,et al.  TOWARD AN EFFICIENT PARALLEL IN TIME METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS , 2012 .

[13]  Michael L. Minion,et al.  Semi-implicit projection methods for incompressible flow based on spectral deferred corrections , 2004 .

[14]  Tao Tang,et al.  High-Order Convergence of Spectral Deferred Correction Methods on General Quadrature Nodes , 2013, J. Sci. Comput..

[15]  John C. Butcher,et al.  The Use of Butcher Series in the Analysis of Newton-Like Iterations in Runge-Kutta Formulas , 2007 .

[16]  L. Greengard,et al.  Spectral Deferred Correction Methods for Ordinary Differential Equations , 2000 .

[17]  Colin B. Macdonald,et al.  Parallel High-Order Integrators , 2010, SIAM J. Sci. Comput..

[18]  Robert Speck,et al.  Spectral deferred corrections with fast-wave slow-wave splitting , 2016, SIAM J. Sci. Comput..

[19]  Michael L. Minion,et al.  Conservative multi-implicit spectral deferred correction methods for reacting gas dynamics , 2004 .

[20]  Robert Speck,et al.  A high-order Boris integrator , 2014, J. Comput. Phys..

[21]  M. Minion Semi-implicit spectral deferred correction methods for ordinary differential equations , 2003 .

[22]  Martin Weiser,et al.  Faster SDC convergence on non-equidistant grids by DIRK sweeps , 2015 .

[23]  Kevin Burrage Parallel methods for ODEs , 1997, Adv. Comput. Math..

[24]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[25]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[26]  Zhaosheng Feng,et al.  Traveling wave behavior for a generalized fisher equation , 2008 .

[27]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[28]  Yinhua Xia,et al.  Efficient time discretization for local discontinuous Galerkin methods , 2007 .