Abstract:Algebraic multigrid (AMG) is known to be an effective solver for many sparse symmetric positive definite (SPD) linear systems. For SPD systems, the convergence theory of AMG is well-understood in terms of the $A$-norm, but in a nonsymmetric setting, such an energy norm is non-existent. For this reason, convergence of AMG for nonsymmetric systems of equations remains an open area of research. A particular aspect missing from theory of nonsymmetric and indefinite AMG is the incorporation of general relaxation schemes. In the SPD setting, the classical form of optimal AMG interpolation provides a useful insight in determining the best possible two-grid convergence rate of a method based on an arbitrary symmetrized relaxation scheme. In this work, we discuss a generalization of the optimal AMG convergence theory targeting nonsymmetric problems, using a certain matrix-induced orthogonality of the left and right eigenvectors of a generalized eigenvalue problem relating the system matrix and relaxation operator. We show that using this generalization of the optimal convergence theory, one can obtain a measure of the spectral radius of the two grid error transfer operator that is mathematically equivalent to the derivation in the SPD setting for optimal interpolation, which instead uses norms. In addition, this generalization of the optimal AMG convergence theory can be further extended for symmetric indefinite problems, such as those arising from saddle point systems so that one can obtain a precise convergence rate of the resulting two-grid method based on optimal interpolation. We provide supporting numerical examples of the convergence theory for nonsymmetric advection-diffusion problems, two-dimensional Dirac equation motivated by $\gamma_5$-symmetry, and the mixed Darcy flow problem corresponding to a saddle point system.

The Convergence of Krylov Subspace Methods for Large Unsymmetric Linear Systems

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

The Modified Solutoin for Large Sparse Symmetric Linear Systems in Electromagnetic Field Analysis

Convergence Properties Of Krylov Subspace Methods For Singular Linear Systems With Arbitrary Index

A Convergence Analysis of the Inexact Rayleigh Quotient Iteration and Simplified Jacobi-Davidson Method for the Large Hermitian Matrix Eigenproblem

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

An Efficient Parallel Krylov-Schur Method for Eigen-Analysis of Large-Scale Power Systems

Approximation Accuracy of the Krylov Subspaces for Linear Discrete Ill-Posed Problems

A semi-conjugate gradient method for solving unsymmetric positive definite linear systems

Fast parallelizable methods for the Hermitian eigenvalue problem

A spectrally preconditioned and initially deflated variant of the restarted block GMRES method for solving multiple right-hand sides linear systems

Accelerating Convergence by Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

Generalized Optimal AMG Convergence Theory for Nonsymmetric and Indefinite Problems

On IGMRES: an Incomplete Generalized Minimal Residual Method for Large Unsymmetric Linear Systems

A global variant of the COCR method for the complex symmetric Sylvester matrix equation AX + XB = C

A Krylov Subspace Method for Large-Scale Second-Order Cone Linear Complementarity Problem

Convergence Analysis of Newton-Like Methods for Inverse Eigenvalue Problems with Multiple Eigenvalues

A block quaternion GMRES method and its convergence analysis

A Cubically Convergent Parallelizable Method for the Hermitian Eigenvalue Problem

On Convergence of the Inexact Rayleigh Quotient Iteration with the Lanczos Method Used for Solving Linear Systems