Abstract:We propose a matrix-free parallel two-level-deflation preconditioner combined with the Complex Shifted Laplacian preconditioner(CSLP) for the two-dimensional Helmholtz problems. The Helmholtz equation is widely studied in seismic exploration, antennas, and medical imaging. It is one of the hardest problems to solve both in terms of accuracy and convergence, due to scalability issues of the numerical solvers. Motivated by the observation that for large wavenumbers, the eigenvalues of the CSLP-preconditioned system shift towards zero, deflation with multigrid vectors, and further high-order vectors were incorporated to obtain wave-number-independent convergence. For large-scale applications, high-performance parallel scalable methods are also indispensable. In our method, we consider the preconditioned Krylov subspace methods for solving the linear system obtained from finite-difference discretization. The CSLP preconditioner is approximated by one parallel geometric multigrid V-cycle. For the two-level deflation, the matrix-free Galerkin coarsening as well as high-order re-discretization approaches on the coarse grid are studied. The results of matrix-vector multiplications in Krylov subspace methods and the interpolation/restriction operators are implemented based on the finite-difference grids without constructing any coefficient matrix. These adjustments lead to direct improvements in terms of memory consumption. Numerical experiments of model problems show that wavenumber independence has been obtained for medium wavenumbers. The matrix-free parallel framework shows satisfactory weak and strong parallel scalability.

A robust and adaptive GenEO-type domain decomposition preconditioner for $\mathbf{H}(\mathbf{curl})$ problems in three-dimensional general topologies

Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption

Substructuring preconditioners with novel interface solvers for general elliptic-type equations in three dimensions

Scalable solvers for complex electromagnetics problems

Nonoverlapping Domain Decomposition Methods with A Simple Coarse Space for Elliptic Problems

Schwarz preconditioner with $H_k$-GenEO coarse space for the indefinite Helmholtz problem

Robust domain decomposition methods for high-contrast multiscale problems on irregular domains with virtual element discretizations

A substructuring preconditioner with vertex-related interface solvers for elliptic-type equations in three dimensions

Space-time balancing domain decomposition

A scalable multilevel domain decomposition preconditioner with a subspace-based coarsening algorithm for the neutron transport calculations

Coarse spaces for non-symmetric two-level preconditioners based on local generalized eigenproblems

Inexact subdomain solves using deflated GMRES for Helmholtz problems

Auxiliary Splines Space Preconditioning for B-Splines Finite Elements: The case of $\bm{H}(\bm{curl},Ω)$ and $\bm{H}(div,Ω)$ elliptic problems

Two-level hybrid Schwarz Preconditioners for The Helmholtz Equation with high wave number

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

Hybrid discontinuous Galerkin discretisation and domain decomposition preconditioners for the Stokes problem

A Novel Preconditioner Based On Csl Operator For Solving The Helmholtz Equations

BDDC and FETI-DP Preconditioners with Adaptive Coarse Spaces for Three-Dimensional Elliptic Problems with Oscillatory and High Contrast Coefficients.

Efficient Parallel Preconditioners for High-Order Finite Element Discretizations of H(grad) and H(curl) Problems

A PRECONDITIONER BASED ON NONOVERLAPPING DOMAIN DECOMPOSITION FOR SOLVING H(curl) ELLIPTIC PROBLEM

Neumann-Neumann type domain decomposition of elliptic problems on metric graphs