Search

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.

In this paper, a few dual least-squares finite element methods and their application to scalar linear hyperbolic problems are studied. The purpose is to obtain $L^2$-norm approximations on finite element spaces of the exact solutions to hyperbolic partial differential equations of interest. This is approached by approximating the generally infeasible quadratic minimization, that defines the $L^2$-orthogonal projection of the exact solution, by feasible least-squares principles using the ideas of the original $\mathcal{L}\mathcal{L}^*$ method proposed in the context of elliptic equations. All methods in this paper are founded upon and extend the $\mathcal{L}\mathcal{L}^*$ approach which is rather general and applicable beyond the setting of elliptic problems. Error bounds are shown that point to the factors affecting the convergence and provide conditions that guarantee optimal rates. Furthermore, the preconditioning of the resulting linear systems is discussed. Numerical results are provided to illustrate the behavior of the methods on common finite element spaces.

Algebraic multigrid (AMG) solvers and preconditioners are some of the fastest numerical methods to solve linear systems, particularly in a parallel environment, scaling to hundreds of thousands of cores. Most AMG methods and theory assume a symmetric positive definite operator. This paper presents a new variation on classical AMG for nonsymmetric matrices (denoted lAIR), based on a local approximation to the ideal restriction operator, coupled with F-relaxation. A new block decomposition of the AMG error-propagation operator is used for a spectral analysis of convergence, and the efficacy of the algorithm is demonstrated on systems arising from the discrete form of the advection-diffusion-reaction equation. lAIR is shown to be a robust solver for various discretizations of the advection-diffusion-reaction equation, including time-dependent and steady-state, from purely advective to purely diffusive. Convergence is robust for discretizations on unstructured meshes and using higher-order finite elements, and is particularly effective on upwind discontinuous Galerkin discretizations. Although the implementation used here is not parallel, each part of the algorithm is highly parallelizable, avoiding common multigrid adjustments for strong advection such as line-relaxation and K- or W-cycles that can be effective in serial, but suffer from high communication costs in parallel, limiting their scalability., Comment: SIAM Journal on Scientific Computing, 40(6), A4105-A4130