ParaDIAG: Parallel-in-Time Algorithms Based on the Diagonalization Technique

by   Martin J. Gander, et al.

In 2008, Maday and Rønquist introduced an interesting new approach for the direct parallel-in-time (PinT) solution of time-dependent PDEs. The idea is to diagonalize the time stepping matrix, keeping the matrices for the space discretization unchanged, and then to solve all time steps in parallel. Since then, several variants appeared, and we call these closely related algorithms ParaDIAG algorithms. ParaDIAG algorithms in the literature can be classified into two groups: * ParaDIAG-I: direct standalone solvers, * ParaDIAG-II: iterative solvers, We will explain the basic features of each group in this note. To have concrete examples, we will introduce ParaDIAG-I and ParaDIAG-II for the advection-diffusion equation. We will also introduce ParaDIAG-II for the wave equation and an optimal control problem for the wave equation. We show the main known theoretical results in each case, and also provide Matlab codes for testing. The goal of the Matlab codes ( is to help the interested reader understand the key features of the ParaDIAG algorithms, without intention to be highly tuned for efficiency and/or low memory use. We also provide speedup measurements of ParaDIAG algorithms for a 2D linear advection-diffusion equation. These results are obtained on the Tianhe-1 supercomputer in China, which is a multi-array, configurable and cooperative parallel system, and we compare these results to the performance of parareal and MGRiT, two widely used PinT algorithms. In a forthcoming update of this note, we will provide more material on ParaDIAG algorithms, in particular further Matlab codes and parallel computing results,also for more realistic applications.


page 1

page 2

page 3

page 4


Global space-time Trefftz DG schemes for the time-dependent linear nonhomogeneous and anisotropic wave equation

In this paper we are concerned with Trefftz discretizations of the time-...

Fast Iterative Solver for the Optimal Control of Time-Dependent PDEs with Crank-Nicolson Discretization in Time

In this article, we derive a new, fast, and robust preconditioned iterat...

A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

The Helmholtz equation is related to seismic exploration, sonar, antenna...

Time-parallel simulation of the Schrödinger Equation

The numerical simulation of the time-dependent Schrödinger equation for ...

A note on parallel preconditioning for the all-at-once solution of Riesz fractional diffusion equations

The p-step backwards difference formula (BDF) for solving the system of ...

Scalable Algorithms for High Order Approximations on Compact Stencils

The recent development of parallel technologies on modern desktop comput...

Approximation of an optimal control problem for the time-fractional Fokker-Planck equation

In this paper, we study the numerical approximation of a system of PDEs ...

Please sign up or login with your details

Forgot password? Click here to reset