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

Tools

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]
Preview
 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

Abstract

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.

ORCID iDs

Bootland, Niall ORCID logoORCID: https://orcid.org/0000-0002-3207-5395, Dolean, Victorita ORCID logoORCID: https://orcid.org/0000-0002-5885-1903, Graham, Ivan G, Ma, Chupeng and Scheichl, Robert;
  • Item type:Article
    ID code:81938
    Dates:
    Date
    Event
    3 August 2023
    Published
    2 August 2022
    Published Online
    18 June 2022
    Accepted
    24 March 2022
    Submitted
    Subjects:Science > Mathematics
    Department:Faculty of Science > Mathematics and Statistics
    Depositing user: Pure Administrator
    Date deposited:19 Aug 2022 10:52
    Last modified:12 Dec 2024 13:40
    URI:https://strathprints.strath.ac.uk/id/eprint/81938
CORE (COnnecting REpositories)