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

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.

Publication Details

Journal
ArXiv.org
Published
2025-05-08
Primary Topic
Matrix Theory and Algorithms
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

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

ArXiv.org
Matrix Theory and Algorithms
article

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

article en

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.

ArXiv.org
Openalex Percentile: Top 47%
Matrix Theory and Algorithms
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.