On Compatible Transfer Operators in Nonsymmetric Algebraic Multigrid

arXiv2307.05900MaRDI QIDQ6511028FDOQ6511028

Authors: Ben S. Southworth, Thomas A. Manteuffel 

Abstract: The standard goal for an effective algebraic multigrid (AMG) algorithm is to develop relaxation and coarse-grid correction schemes that attenuate complementary error modes. In the nonsymmetric setting, coarse-grid correction $Pi$ will almost certainly be nonorthogonal (and divergent) in any known inner product, meaning $|Pi|>1$. This introduces a new consideration, that one wants coarse-grid correction to be as close to orthogonal as possible, in an appropriate norm. In addition, due to non-orthogonality, $Pi$ may actually amplify certain error modes that are in the range of interpolation. Relaxation must then not only be complementary to interpolation, but also rapidly eliminate any error amplified by the non-orthogonal correction, or the algorithm may diverge. This note develops analytic formulae on how to construct ``compatible transfer operators in nonsymmetric AMG such that $|Pi|=1$ in any standard matrix-induced norm. Discussion is provided on different options for norm in the nonsymmetric setting, the relation between ``ideal transfer operators in different norms, and insight into the convergence of nonsymmetric reduction-based AMG.

This page was built for publication: On Compatible Transfer Operators in Nonsymmetric Algebraic Multigrid

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6511028)