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

TL;DR

This paper develops a theoretical framework for optimal transfer operators in algebraic two-level methods applicable to nonsymmetric and indefinite problems, extending previous HPD-specific results and providing practical real-valued solutions.

Contribution

It characterizes all inner products for orthogonality of coarse-space correction, develops tight convergence bounds, and constructs optimal real-valued transfer operators for non-HPD matrices.

Findings

Derived tight convergence bounds in specific inner product norms.

Proved the optimality of the constructed transfer operators.

Validated the theory with numerical examples from advection and wave equations.

Abstract

Consider an algebraic two-level method applied to the $n$-dimensional linear system $A \mathbf{x} = \mathbf{b}$ using fine-space preconditioner (i.e., ``relaxation'' or ``smoother'') $M$, with $M \approx A$, restriction and interpolation $R$ and $P$, and algebraic coarse-space operator ${A_c := R^*AP}$. Then, what are the the best possible transfer operators $R$ and $P$ of a given dimension $n_c < 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 $n_c$ 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.