Optimal transfer operators in algebraic two-level methods for nonsymmetric and indefinite problems

Abstract

Consider an algebraic two-level method applied to the n-dimensional linear system A x = b using fine-space preconditioner (i.e., ``relaxation'' or ``smoother'') M, with M ≈ A, restriction and interpolation R and P, and algebraic coarse-space operator Ac := R*AP. Then, what are the the best possible transfer operators R and P of a given dimension nc < n? Brannick et al. (2018) showed that when A and M are Hermitian positive definite (HPD), the optimal interpolation is such that its range contains the nc smallest generalized eigenvectors of the matrix pencil (A, M). Recently, in Ali et al. (2025) we generalized this framework to the non-HPD setting, by considering both right (interpolation) and left (restriction) generalized eigenvectors of (A, M) and defining corresponding nonsymmetric transfer operators \R\#,P\#\. Tight convergence bounds for \R\#,P\#\ are derived in spectral radius, as well as a proof of pseudo-optimality. Note, \R\#,P\#\ are typically complex valued, which is not practical for real-valued problems. Here we build on Ali et al. (2025), first characterizing all inner products in which the coarse-space correction defined by \R\#,P\#\ is orthogonal. We then develop tight two-level convergence bounds in these norms, and prove that the underlying transfer operators \R\#,P\#\ are genuinely optimal. As a special case, our theory both recovers and extends the HPD results from Brannick et al. (2018). Finally, we show how to construct optimal, real-valued transfer operators in the case of that A and M are real valued, but are not HPD. Numerical examples arising from discretized advection and wave-equation problems are used to verify and illustrate the theory.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…