Improvements to the theoretical estimates of the Schwarz preconditioner with $Δ$-GenEO coarse space for the indefinite Helmholtz problem

Victorita Dolean,Mark Fry,Matthias Langer
2024-06-18
Abstract:The purpose of this work is to improve the estimates for the $\Delta$-GenEO method from the paper "Overlapping Schwarz methods with GenEO coarse spaces for indefinite and nonself-adjoint problems" by N. Bootland, V. Dolean, I. G Graham, C. Ma, R. Scheichl (<a class="link-external link-https" href="https://doi.org/10.1093/imanum/drac036" rel="external noopener nofollow">this https URL</a>) when applied to the indefinite Helmholtz equation. We derive k-dependent estimates of quantities of interest ensuring the robustness of the method.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the theoretical estimates of the ∆-GenEO coarse - space Schwarz preconditioner applied to the indefinite Helmholtz equation. Specifically, the author aims to derive estimates that depend on the wavenumber \(k\) to ensure the robustness of the method. ### Description of the Main Problem in the Paper When solving the indefinite Helmholtz equation, in the case of dealing with highly variable coefficients and increasing wavenumbers, the Schwarz preconditioner usually requires an additional global coarse - grid solution to enhance its scalability and robustness. For this purpose, the author uses the local generalized eigenvalue problem and obtains a symmetric but indefinite linear system through finite - element discretization. Due to the indefiniteness of the system, the iterative solver GMRES is adopted, and its convergence analysis depends on Elman's theory, which requires an upper bound on the preconditioned matrix norm and a lower bound on the distance from its numerical range to the origin. ### Main Improvements The main theoretical results of the paper provide strict and explicitly \(k\)-dependent upper bounds on the coarse - grid diameter \(H\) and the eigenvalue tolerance \(\tau\) to ensure robust and mesh - independent convergence of GMRES when applying the preconditioned problem. Specifically, the author improves the estimates in previous literature: - The original estimates are: \[ H\lesssim k^{- 2}\quad\text{and}\quad(1 + C_{\text{stab}})^2k^8\lesssim\tau \] - The improved estimates are: \[ H\lesssim k^{-1}\quad\text{and}\quad(1 + C_{\text{stab}})^2k^4\lesssim\tau \] Although these improvements still show a strong dependence on the wavenumber \(k\), they are more optimistic than the previous estimates. ### Conclusion Through these improvements, the paper enhances the theoretical estimates of the preconditioner, thereby ensuring better robustness and mesh - independence of the GMRES method when dealing with indefinite Helmholtz problems. This provides theoretical support for more efficiently solving such problems in practical applications.