Abstract:Solving time-harmonic wave propagation problems by iterative methods is a difficult task, and over the last two decades an important research effort has gone into developing preconditioners for the simplest representative of such wave propagation problems, the Helmholtz equation. A specific class of these new preconditioners is considered here. They were developed by researchers with various backgrounds using formulations and notations that are very different, and all are among the most promising preconditioners for the Helmholtz equation. \indent The goal of the present article is to show that this class of preconditioners is based on a common mathematical principle, and that they can all be formulated in the context of domain decomposition methods known as optimized Schwarz methods. This common formulation allows us to explain in detail how and why all these methods work. The domain decomposition formulation also allows us to avoid technicalities in the implementation description we give of these recent methods. \indent The equivalence of these methods with optimized Schwarz methods translates at the discrete level into equivalence with approximate block LU decomposition preconditioners, and in each case we give the algebraic version, including a detailed description of the approximations used. While we choose to use the Helmholtz equation for which these methods were developed, our notation is completely general and the algorithms we give are written for an arbitrary second-order elliptic operator. The algebraic versions are even more general, assuming only a connectivity pattern in the discretization matrix. \indent All the new methods studied here are based on sequential decomposition of the problem in space into a sequence of subproblems, and they have in their optimal form the property to lead to nilpotent iterations, like an exact block LU factorization. Using our domain decomposition formulation, we finally present an algorithm for two-dimensional decompositions, i.e., decompositions that contain cross points, which is still nilpotent in its optimal form. Its approximation is currently an active area of research, and it would have been difficult to discover such an algorithm without the domain decomposition framework.

A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

A matrix-free parallel two-level deflation preconditioner for the two-dimensional Helmholtz problems

A Parallel Sweeping Preconditioner for Heterogeneous 3D Helmholtz Equations.

Adaptive complex frequency with V-cycle GMRES for preconditioning 3D Helmholtz equation

Parallel Chebyshev spectral element method for Helmholtz equation in 3D domain

An anisotropic multilevel preconditioner for solving the Helmholtz equation with unequal directional sampling intervals

A Neural Multigrid Solver for Helmholtz Equations with High Wavenumber and Heterogeneous Media

Exponentially convergent multiscale methods for high frequency heterogeneous Helmholtz equations

An accurate and efficient multiscale finite-difference frequency-domain method for the scalar Helmholtz equation

A novel direct Helmholtz solver in inhomogeneous media based on the operator Fourier transform functional calculus

Stand-alone Multigrid for Helmholtz Revisited: Towards Convergence Using Standard Components

Parallel Preconditioned GMRES Solvers for 3-D Helmholtz Equations in Regional Non-hydrostatic Atmosphere Model

A parallel shared-memory implementation of a high-order accurate solution technique for variable coefficient Helmholtz problems

Matrix-Free Parallel Scalable Multilevel Deflation Preconditioning for Heterogeneous Time-Harmonic Wave Problems

Parallel block multigrid preconditioner for 3D navier-stokes equations on unstructured grids

A Class of Iterative Solvers for the Helmholtz Equation: Factorizations, Sweeping Preconditioners, Source Transfer, Single Layer Potentials, Polarized Traces, and Optimized Schwarz Methods

A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations

Multigrid-Augmented Deep Learning Preconditioners for the Helmholtz Equation using Compact Implicit Layers