DeepAI AI Chat
Log In Sign Up

Parallel-in-time preconditioners for the Sinc-Nyström method

by   Jun Liu, et al.

The Sinc-Nyström method in time is a high-order spectral method for solving evolutionary differential equations and it has wide applications in scientific computation. But in this method we have to solve all the time steps implicitly at one-shot, which may results in a large-scale nonsymmetric dense system that is expensive to solve. In this paper, we propose and analyze a parallel-in-time (PinT) preconditioner for solving such Sinc-Nyström systems, where both the parabolic and hyperbolic PDEs are investigated. Attributed to the special Toeplitz-like structure of the Sinc-Nyström systems, the proposed PinT preconditioner is indeed a low-rank perturbation of the system matrix and we show that the spectrum of the preconditioned system is highly clustered around one, especially when the time step size is refined. Such a clustered spectrum distribution matches very well with the numerically observed mesh-independent GMRES convergence rates in various examples. Several linear and nonlinear ODE and PDE examples are presented to illustrate the convergence performance of our proposed PinT preconditioners, where the achieved exponential order of accuracy are especially attractive to those applications in need of high accuracy.


page 1

page 2

page 3

page 4


A boundary-penalized isogeometric analysis for second-order hyperbolic equations

Explicit time-marching schemes are popular for solving time-dependent pa...

Exponential integrators for large-scale stiff matrix Riccati differential equations

Matrix Riccati differential equations arise in many different areas and ...

C1 Triangular Isogeometric Analysis of the von Karman Equations

In this paper, we report the use of rational Triangular Bezier Splines (...

An all-at-once preconditioner for evolutionary partial differential equations

In [McDonald, Pestana and Wathen, SIAM J. Sci. Comput., 40 (2018), pp. A...

Analysis of GMRES for Low-Rank and Small-Norm Perturbations of the Identity Matrix

In many applications, linear systems arise where the coefficient matrix ...

ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems

ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for s...

A performance spectrum for parallel computational frameworks that solve PDEs

Important computational physics problems are often large-scale in nature...