Preconditioning and iterative solution of all-at-once systems for evolutionary partial differential equations

McDonald, Eleanor and Pestana, Jennifer and Wathen, Andy (2018) Preconditioning and iterative solution of all-at-once systems for evolutionary partial differential equations. SIAM Journal on Scientific Computing, 40 (2). A1012–A1033. ISSN 1064-8275 (https://doi.org/10.1137/16M1062016)

[thumbnail of McDonald-Pestana-Wathen-2018-SIAM-JOSC-Preconditioning-and-iterative-solution-of-all-at-once-systems]
Preview
 Text. Filename: McDonald_Pestana_Wathen_2018_SIAM_JOSC_Preconditioning_and_iterative_solution_of_all_at_once_systems.pdf
Accepted Author Manuscript
Download (726kB)| Preview

Abstract

Standard Krylov subspace solvers for self-adjoint problems have rigorous convergence bounds based solely on eigenvalues. However, for non-self-adjoint problems, eigenvalues do not determine behavior even for widely used iterative methods. In this paper, we discuss time-dependent PDE problems, which are always non-self-adjoint. We propose a block circulant preconditioner for the all-at-once evolutionary PDE system which has block Toeplitz structure. Through reordering of variables to obtain a symmetric system, we are able to rigorously establish convergence bounds for MINRES which guarantee a number of iterations independent of the number of time-steps for the all-at-once system. If the spatial differential operators are simultaneously diagonalizable, we are able to quickly apply the preconditioner through use of a sine transform, and for those that are not, we are able to use an algebraic multigrid process to provide a good approximation. Results are presented for solution to both the heat and convection diffusion equations.

ORCID iDs

McDonald, Eleanor, Pestana, Jennifer ORCID logoORCID: https://orcid.org/0000-0003-1527-3178 and Wathen, Andy;
CORE (COnnecting REpositories)