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.

A Fast Solver for Generalized Optimal Transport Problems Based on Dynamical System and Algebraic Multigrid

An efficient semismooth Newton-AMG-based inexact primal-dual algorithm for generalized transport problems

A Scalable Deep Learning Approach for Solving High-Dimensional Dynamic Optimal Transport.

A Multiscale Semi-Smooth Newton Method for Optimal Transport

PDOT: a Practical Primal-Dual Algorithm and a GPU-Based Solver for Optimal Transport

Gradient enhanced ADMM Algorithm for dynamic optimal transport on surfaces

Unconditional convergence for discretizations of dynamical optimal transport

Numerical method for feasible and approximately optimal solutions of multi-marginal optimal transport beyond discrete measures

An Accelerated Stochastic Algorithm for Solving the Optimal Transport Problem

A Sparse Smoothing Newton Method for Solving Discrete Optimal Transport Problems

Efficient preconditioners for solving dynamical optimal transport via interior point methods

A fast primal-dual algorithm via dynamical system with variable mass for linearly constrained convex optimization

Duality-based Dynamical Optimal Transport of Discrete Time Systems

Generalized Optimal AMG Convergence Theory for Nonsymmetric and Indefinite Problems

Multilevel Optimal Transport: a Fast Approximation of Wasserstein-1 distances

Structure Preserving Primal Dual Methods for Gradient Flows with Nonlinear Mobility Transport Distances

Asymptotic analysis of domain decomposition for optimal transport

Randomized Methods for Computing Optimal Transport Without Regularization and Their Convergence Analysis

Fast Gradient Computation for Gromov-Wasserstein Distance

Fixed-Time Gradient Flows for Solving Constrained Optimization: A Unified Approach

HOT: An Efficient Halpern Accelerating Algorithm for Optimal Transport Problems