Overlapping Schwarz methods with GenEO coarse spaces for indefinite and nonself-adjoint problems

Bootland, Niall and Dolean, Victorita and Graham, Ivan G and Ma, Chupeng and Scheichl, Robert (2023) Overlapping Schwarz methods with GenEO coarse spaces for indefinite and nonself-adjoint problems. IMA Journal of Numerical Analysis, 43 (4). pp. 1899-1936. ISSN 0272-4979 (https://doi.org/10.1093/imanum/drac036)

[thumbnail of Bootland-etal-JNA-2022-Overlapping-Schwarz-methods-with-GenEO-coarse-spaces-for-indefinite]
Text. Filename: Bootland_etal_JNA_2022_Overlapping_Schwarz_methods_with_GenEO_coarse_spaces_for_indefinite.pdf
Accepted Author Manuscript
License: Strathprints license 1.0

Download (1MB)| Preview


Generalized eigenvalue problems on the overlap(GenEO) is a method for computing an operator-dependent spectral coarse space to be combined with local solves on subdomains to form a robust parallel domain decomposition preconditioner for elliptic PDEs. It has previously been proved, in the self-adjoint and positive-definite case, that this method, when used as a preconditioner for conjugate gradients, yields iteration numbers that are completely independent of the heterogeneity of the coefficient field of the partial differential operator. We extend this theory to the case of convection–diffusion–reaction problems, which may be nonself-adjoint and indefinite, and whose discretizations are solved with preconditioned GMRES. The GenEO coarse space is defined here using a generalized eigenvalue problem based on a self-adjoint and positive-definite subproblem. We prove estimates on GMRES iteration counts that are independent of the variation of the coefficient of the diffusion term in the operator and depend only very mildly on variations of the other coefficients. These are proved under the assumption that the subdomain diameter is sufficiently small and the eigenvalue tolerance for building the coarse space is sufficiently large. While the iteration number estimates do grow as the nonself-adjointness and indefiniteness of the operator increases, practical tests indicate the deterioration is much milder. Thus, we obtain an iterative solver that is efficient in parallel and very effective for a wide range of convection–diffusion–reaction problems.