Saved in:
Bibliographic Details
Main Authors: Southworth, Ben S., Manteuffel, Thomas A.
Format: Preprint
Published: 2023
Subjects:
Online Access:https://arxiv.org/abs/2307.05900
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909072553410560
author Southworth, Ben S.
Manteuffel, Thomas A.
author_facet Southworth, Ben S.
Manteuffel, Thomas A.
contents 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 $Π$ will almost certainly be nonorthogonal (and divergent) in any known standard product, meaning $\|Π\| > 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, $Π$ 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 paper develops analytic formulae on how to construct ``compatible'' transfer operators in nonsymmetric AMG such that $\|Π\| = 1$ in some standard matrix-induced norm. Discussion is provided on different options for the norm in the nonsymmetric setting, the relation between ``ideal'' transfer operators in different norms, and insight into the convergence of nonsymmetric reduction-based AMG.
format Preprint
id arxiv_https___arxiv_org_abs_2307_05900
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle On Compatible Transfer Operators in Nonsymmetric Algebraic Multigrid
Southworth, Ben S.
Manteuffel, Thomas A.
Numerical Analysis
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 $Π$ will almost certainly be nonorthogonal (and divergent) in any known standard product, meaning $\|Π\| > 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, $Π$ 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 paper develops analytic formulae on how to construct ``compatible'' transfer operators in nonsymmetric AMG such that $\|Π\| = 1$ in some standard matrix-induced norm. Discussion is provided on different options for the norm in the nonsymmetric setting, the relation between ``ideal'' transfer operators in different norms, and insight into the convergence of nonsymmetric reduction-based AMG.
title On Compatible Transfer Operators in Nonsymmetric Algebraic Multigrid
topic Numerical Analysis
url https://arxiv.org/abs/2307.05900