DeepAI AI Chat
Log In Sign Up

A well-conditioned direct PinT algorithm for first- and second-order evolutionary equations

by   Jun Liu, et al.

In this paper, we propose a direct parallel-in-time (PinT) algorithm for time-dependent problems with first- or second-order derivative. We use a second-order boundary value method as the time integrator that leads to a tridiagonal time discretization matrix. Instead of solving the corresponding all-at-once system iteratively, we diagonalize the time discretization matrix, which yields a direct parallel implementation across all time levels. A crucial issue on this methodology is how the condition number of the eigenvector matrix V grows as n is increased, where n is the number of time levels. A large condition number leads to large roundoff error in the diagonalization procedure, which could seriously pollute the numerical accuracy. Based on a novel connection between the characteristic equation and the Chebyshev polynomials, we present explicit formulas for computing V and V^-1, by which we prove that Cond_2(V)=𝒪(n^2). This implies that the diagonalization process is well-conditioned and the roundoff error only increases moderately as n grows and thus, compared to other direct PinT algorithms, a much larger n can be used to yield satisfactory parallelism. Numerical results on parallel machine are given to support our findings, where over 60 times speedup is achieved with 256 cores.


page 1

page 2

page 3

page 4


Substitution Method for Fractional Differential Equations

Numerical solving differential equations with fractional derivatives req...

Subdiffusion with Time-Dependent Coefficients: Improved Regularity and Second-Order Time Stepping

This article concerns second-order time discretization of subdiffusion e...

A pseudo-spectral Strang splitting method for linear dispersive problems with transparent boundary conditions

The present work proposes a second-order time splitting scheme for a lin...

A Direct Parallel-in-Time Quasi-Boundary Value Method for Inverse Space-Dependent Source Problems

Inverse source problems arise often in real-world applications, such as ...

A Second-Order TGV Discretization with Some Invariance Properties

In this work, we propose a new discretization for second-order total gen...

Full operator preconditioning and the accuracy of solving linear systems

Unless special conditions apply, the attempt to solve ill-conditioned sy...

Efficient computation of the density matrix with error control on distributed computer systems

The recursive polynomial expansion for construction of a density matrix ...