A sine transform based preconditioned MINRES method for all-at-once systems from constant and variable-coefficient evolutionary PDEs

Sean Y S Hon*, Po Yin Fung, Jiamei Dong, Stefano Serra-Capizzano

*Corresponding author for this work

Research output: Contribution to journalJournal articlepeer-review

Abstract

In this work, we propose a simple yet generic preconditioned Krylov subspace method for a large class of nonsymmetric block Toeplitz all-at-once systems arising from discretizing evolutionary partial differential equations. Namely, our main result is to propose two novel symmetric positive definite preconditioners, which can be efficiently diagonalized by the discrete sine transform matrix. More specifically, our approach is to first permute the original linear system to obtain a symmetric one and subsequently develop desired preconditioners based on the spectral symbol of the modified matrix. Then, we show that the eigenvalues of the preconditioned matrix sequences are clustered around ±1, which entails rapid convergence when the minimal residual method is devised. Alternatively, when the conjugate gradient method on the normal equations is used, we show that our preconditioner is effective in the sense that the eigenvalues of the preconditioned matrix sequence are clustered around unity. An extension of our proposed preconditioned method is given for high-order backward difference time discretization schemes, which can be applied on a wide range of time-dependent equations. Numerical examples are given, also in the variable-coefficient setting, to demonstrate the effectiveness of our proposed preconditioners, which consistently outperforms an existing block circulant preconditioner discussed in the relevant literature.
Original languageEnglish
JournalNumerical Algorithms
DOIs
Publication statusE-pub ahead of print - 14 Aug 2023

Scopus Subject Areas

  • Applied Mathematics

User-Defined Keywords

  • All-at-once discretization
  • Block Toeplitz systems
  • Krylov subspace methods
  • MINRES
  • Parallel-in-time
  • Sine transform based preconditioners

Fingerprint

Dive into the research topics of 'A sine transform based preconditioned MINRES method for all-at-once systems from constant and variable-coefficient evolutionary PDEs'. Together they form a unique fingerprint.

Cite this