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

08/03/2021
by   Jun Liu, et al.
0

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.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/23/2022

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

Explicit time-marching schemes are popular for solving time-dependent pa...
research
07/30/2019

Exponential integrators for large-scale stiff matrix Riccati differential equations

Matrix Riccati differential equations arise in many different areas and ...
research
08/02/2021

C1 Triangular Isogeometric Analysis of the von Karman Equations

In this paper, we report the use of rational Triangular Bezier Splines (...
research
02/04/2020

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

In [McDonald, Pestana and Wathen, SIAM J. Sci. Comput., 40 (2018), pp. A...
research
10/21/2022

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

In many applications, linear systems arise where the coefficient matrix ...
research
02/23/2021

Explicit high-order generalized-α methods for isogeometric analysis of structural dynamics

We propose a new family of high-order explicit generalized-α methods for...
research
05/10/2017

A performance spectrum for parallel computational frameworks that solve PDEs

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

Please sign up or login with your details

Forgot password? Click here to reset