pith. sign in

arxiv: 2605.31552 · v1 · pith:PTSMKZBBnew · submitted 2026-05-29 · 🧮 math.NA · cs.NA

Spectral coarse spaces based on indefinite operators: the H_k-GenEO method

Th\'eophile Chaumont-Frelet , Victorita Dolean , Mark Fry , Ivan G. Graham , Matthias Langer This is my paper

Pith reviewed 2026-06-28 21:09 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords GenEOcoarse spacesindefinite operatorsdomain decompositionpreconditioningGMRESspectral methods
0
0 comments X

The pith

The H_k-GenEO method builds coarse spaces from local indefinite eigenvalue problems to precondition highly indefinite PDEs robustly.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a spectral coarse space method called H_k-GenEO for preconditioning iterative solvers applied to discrete PDEs that are highly indefinite, characterized by a large parameter k. It constructs the coarse space by solving local eigenvalue problems that are copies of the global indefinite operator, without restricting the diameter of the subdomains. This approach is shown to be more robust than previous GenEO methods that use positive semi-definite local problems as k grows. Sufficient conditions are given for the robustness of the preconditioned GMRES method in terms of the tolerance used in the local eigenproblems and the size of the subdomains. Experiments indicate the method works well even under weaker conditions and with variable coefficients.

Core claim

The H_k-GenEO method forms local eigenvalue problems from indefinite copies of the global problem on subdomains of arbitrary diameter. Sufficient conditions on the eigenproblem tolerance and subdomain size ensure robustness of the preconditioned GMRES iterative method for the global solve.

What carries the argument

The H_k-GenEO coarse space, constructed by combining selected modes from local indefinite eigenproblems based on the global operator.

Load-bearing premise

Local eigenvalue problems formed from indefinite copies of the global operator remain well-posed and computationally tractable on subdomains of arbitrary diameter.

What would settle it

Observation that the number of GMRES iterations grows without bound as k increases, even when local eigenproblem tolerance and subdomain sizes satisfy the sufficient conditions.

read the original abstract

GenEO (`Generalised Eigenvalue problems on the Overlap') is a method for constructing coarse spaces used in the preconditioning of iterative solvers for discrete PDEs. This method combines a (small) number of modes of local PDE eigenproblems to obtain a global coarse space. A coarse solve is then combined with local solves of the global PDE to obtain the preconditioner. A substantial theory for GenEO has been developed for the case when the local elgenproblems are positive semi-definite. This has been applied mostly to positive definite global PDEs, but also recently extended to the case of convection--diffusion--reaction problems, which may be neither self-adjoint, nor positive definite. However, when the global problem is highly indefinite, coarse spaces built from positive semi-definite local eigenproblems fail to be robust in practice. In this paper we consider highly indefinite global PDE problems, characterised by a large parameter $k$ (allowing also highly variable coefficients), and we develop a new spectral coarse space built from solving eigenvalue problems based on \textit{local copies of the global problem}. We put no constraint on the diameters of the local domains, thus allowing the local eigenvalue problems to be indefinite. The new method (which we call $H_k$-GenEO) is seen to be much more robust as $k$ increases than methods based on positive semi-definite eigenproblems. We provide sufficient conditions for robustness of the preconditioned GMRES iterative method, in terms of the tolerance of the local eigenproblems and the size of the subdomains for the local PDE solves. In practice the method is observed to be robust with respect to $k$ under even weaker conditions on the local eigenproblem tolerance. The experiments also suggest the method can be resilient to high variation in PDE coefficients.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper introduces the H_k-GenEO method for constructing spectral coarse spaces in domain-decomposition preconditioners for highly indefinite discrete PDEs characterized by a large parameter k (and possibly highly variable coefficients). Unlike prior GenEO variants that rely on positive semi-definite local eigenproblems, H_k-GenEO solves local eigenproblems formed from indefinite copies of the global operator on subdomains of arbitrary diameter. The central claims are that the resulting coarse space yields a preconditioner for which GMRES is robust as k grows, that sufficient conditions guaranteeing this robustness can be stated in terms of the local eigenproblem tolerance τ and subdomain diameter H, and that experiments confirm practical robustness even under weaker conditions on τ together with resilience to coefficient variation.

Significance. If the sufficient conditions are shown to be independent of k for fixed H and the experimental evidence is reproducible, the work would meaningfully extend the GenEO framework to a class of problems (high-frequency indefinite operators) where existing positive semi-definite coarse-space constructions are known to lose robustness. The absence of any diameter restriction on the subdomains is a potentially useful feature for practical implementations.

major comments (2)
  1. [analysis section deriving the sufficient conditions] The abstract states that sufficient conditions for GMRES robustness are provided in terms of the local eigenproblem tolerance and subdomain size. The load-bearing question is whether the hidden constants in those conditions remain independent of k when H is held fixed; if the analysis involves factors that grow with k (e.g., via the number of negative eigenvalues or kH-type terms), the claimed robustness for arbitrary diameters would not follow. The manuscript must make this independence explicit in the derivation of the sufficient conditions.
  2. [numerical experiments] The experimental section reports that the method is observed to be robust under weaker conditions than the sufficient ones. To support the central claim, the experiments must demonstrate that, for fixed τ and H independent of k, the number of retained eigenmodes does not explode and GMRES iteration counts remain bounded as k increases; otherwise the practical robustness may still rely on k-dependent tuning.
minor comments (1)
  1. Notation for the local indefinite operator and the tolerance parameter τ should be introduced with a clear reference to the global problem statement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and for the constructive comments. We address each major comment below.

read point-by-point responses

  1. Referee: [analysis section deriving the sufficient conditions] The abstract states that sufficient conditions for GMRES robustness are provided in terms of the local eigenproblem tolerance and subdomain size. The load-bearing question is whether the hidden constants in those conditions remain independent of k when H is held fixed; if the analysis involves factors that grow with k (e.g., via the number of negative eigenvalues or kH-type terms), the claimed robustness for arbitrary diameters would not follow. The manuscript must make this independence explicit in the derivation of the sufficient conditions.

    Authors: The analysis derives the sufficient conditions using the spectral information from the indefinite local eigenproblems (with tolerance τ) and the subdomain diameter H. The resulting bounds on the GMRES convergence factor depend only on τ and H; no factors involving the number of negative eigenvalues or explicit kH growth appear because the local operators are direct copies of the global indefinite operator. To make this independence fully explicit, we will add a short remark immediately after the main theorem stating that all constants are independent of k for fixed H and τ. revision: partial

  2. Referee: [numerical experiments] The experimental section reports that the method is observed to be robust under weaker conditions than the sufficient ones. To support the central claim, the experiments must demonstrate that, for fixed τ and H independent of k, the number of retained eigenmodes does not explode and GMRES iteration counts remain bounded as k increases; otherwise the practical robustness may still rely on k-dependent tuning.

    Authors: The numerical experiments already include series of tests performed with fixed τ and fixed H while k is increased over several orders of magnitude. In these runs the number of retained eigenmodes grows only mildly and the GMRES iteration counts remain essentially bounded, consistent with the observed practical robustness. We will revise the experimental section to label these fixed-(τ,H) sequences more prominently and to add a short paragraph summarizing the observed independence from k. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation is a direct construction with independent sufficient conditions

full rationale

The paper introduces the H_k-GenEO method as an explicit construction using local copies of the indefinite global operator on subdomains of arbitrary diameter. It states sufficient conditions on eigenproblem tolerance and subdomain size for GMRES robustness, and reports that experiments succeed under weaker conditions. No quoted step equates a prediction to a fitted input by construction, renames a known result, or reduces the central claim to a self-citation chain. Self-citations to prior GenEO theory are present but not load-bearing for the new indefinite-operator extension or the k-robustness claim. The analysis remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on abstract only; no explicit free parameters, axioms, or invented entities are identifiable beyond standard assumptions of linear algebra and PDE theory.

pith-pipeline@v0.9.1-grok · 5878 in / 1080 out tokens · 23705 ms · 2026-06-28T21:09:13.588920+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 9 canonical work pages

  1. [1]

    AGULLO, L

    E. AGULLO, L. GIRAUD ANDL. POIREL,Robust preconditioners via generalized eigenproblems for hybrid sparse linear solvers, SIAM J. Matrix Anal. Appl. 40, 417–439, 2019

    2019

  2. [2]

    BABUSKA ANDR

    I. BABUSKA ANDR. LIPTON,Optimal local approximation spaces for generalized finite element methods with application to multiscale problems, Multiscale Model. Simul. 9, 373–406, 2011

    2011

  3. [3]

    BASTIAN, R

    P. BASTIAN, R. SCHEICHL, L. SEELINGER ANDA. STREHLOW,Multilevel spectral domain decomposition, SIAM J. Sci. Comput. 45, S1–S26, 2023

    2023

  4. [4]

    BOOTLAND ANDV

    N. BOOTLAND ANDV. DOLEAN,On the Dirichlet-to-Neumann coarse space for solving the Helmholtz problem using domain decomposition. InNumerical Mathematics and Advanced Applications ENUMATH 2019; Vermolen, F. J., Vuik, C., Eds.; Springer: Cham, Switzerland, 2021; pp. 175–184

    2019

  5. [5]

    BOOTLAND, V

    N. BOOTLAND, V. DOLEAN, I.G. GRAHAM, C. MA ANDR. SCHEICHL,GenEO coarse spaces for heterogeneous indefinite elliptic problems. InDomain Decomposition Methods in Science and Engineering XXVI; S. Brenner et al., Eds.; Springer: Cham, Switzerland, 2022; pp. 117–125

    2022

  6. [6]

    BOOTLAND, V

    N. BOOTLAND, V. DOLEAN, I.G. GRAHAM, C. MA ANDR. SCHEICHL,Overlapping Schwarz methods with GenEO coarse spaces for indefinite and non-self-adjoint problems, IMA J. Numer. Anal. 43, 1899–1936, 2023

    1936

  7. [7]

    BOOTLAND, V

    N. BOOTLAND, V. DOLEAN, P. JOLIVET ANDP.-H. TOURNIER,A comparison of coarse spaces for Helmholtz problems in the high frequency regime, Comput. Math. Appl. 98, 239–253, 2021

    2021

  8. [8]

    CIARLET,The Finite Element Method for Elliptic Problems, SIAM, 2002

    P.G. CIARLET,The Finite Element Method for Elliptic Problems, SIAM, 2002

    2002

  9. [9]

    CONEN, V

    L. CONEN, V. DOLEAN, R. KRAUSE ANDF. NATAF,A coarse space for heterogeneous Helmholtz problems based on the Dirichlet-to-Neumann operator, J. Comput. Appl. Math. 271, 83–99, 2014

    2014

  10. [10]

    Statistical inverse problems: discretization, model reduction and inverse crimes

    V. DOLEAN, M. FRY ANDM. LANGER,Can Symmetric Positive Definite (SPD) coarse spaces perform well for indefinite Helmholtz problems?, J. of Comput. Appl. Math. 484, 117403, https://doi.org/10.1016/j.cam. 2026.117403, 2026

    work page doi:10.1016/j.cam 2026

  11. [11]

    DOLEAN, M

    V. DOLEAN, M. FRY, M. LANGER, E. PAROLIN ANDP.-H. TOURNIER,Achieving wavenumber robustness in domain decomposition for heterogeneous Helmholtz equation: an overview of spectral coarse spaces, https: //arxiv.org/abs/2509.02131, 2025

    arXiv 2025

  12. [12]

    EFENDIEV, J

    Y. EFENDIEV, J. GALVIS ANDT.-Y. HOU,Generalized multiscale finite element methods (GMsFEM), J. Comput. Phys. 251, 116–135, 2013

    2013

  13. [13]

    EISENSTAT, H.C

    S.C. EISENSTAT, H.C. ELMAN ANDM.H. SCHULTZ,Variational iterative methods for nonsymmetric systems of linear equations, SIAM J. Numer. Anal. 20, 345–357, 1983

    1983

  14. [14]

    S. FU, S. GONG, G. LI ANDY. WANG,On edge multiscale space based hybrid Schwarz preconditioner for Helmholtz problems with large wavenumbers, https://arxiv.org/abs/2408.08198, 2024

    arXiv 2024

  15. [15]

    GALKOWSKI ANDE.A

    J. GALKOWSKI ANDE.A. SPENCE,Convergence theory for two-level hybrid Schwarz preconditioners for high-frequency Helmholtz problems, SIAM J. Numer. Anal. 64, 29–54, https://doi.org/10.1137/25M1726972, 2026

    work page doi:10.1137/25m1726972 2026

  16. [16]

    GALVIS, E.T

    J. GALVIS, E.T. CHUNG, Y. EFENDIEV ANDW.T. LEUNG,On Overlapping Domain Decomposition Methods for High-Contrast Multiscale Problems. In: P. Bjørstad et al. (eds), Domain Decomposition Methods in Science and Engineering XXIV ,Lect. Notes Comput. Sci. Eng., vol 125, pp. 45–57. Springer, Cham, 2018. COARSE SPACES BASED ON INDEFINITE OPERATORS41

    2018

  17. [17]

    GALVIS ANDY

    J. GALVIS ANDY. EFENDIEV,Domain decomposition preconditioners for multiscale flows in high-contrast media, Multiscale Model. Simul. 8, 1461–1483, 2010

    2010

  18. [18]

    GALVIS ANDY

    J. GALVIS ANDY. EFENDIEV,Domain decomposition preconditioners for multiscale flows in high contrast media: reduced dimension coarse spaces, Multiscale Model. Simul. 8, 1621–1644, 2010

    2010

  19. [19]

    GRAHAM ANDE.A

    I.G. GRAHAM ANDE.A. SPENCE,Theory of two-level Schwarz preconditioners with piecewise-polynomial coarse spaces for the high-frequency Helmholtz equation, https://arxiv.org/abs/2501.15976, 2025

    Pith/arXiv arXiv 2025

  20. [20]

    GRAHAM, E.A

    I.G. GRAHAM, E.A. SPENCE ANDE. VAINIKKO,Domain decomposition preconditioning for high-frequency Helmholtz problems with absorption, Math. Comput. 86, 2089-2127, 2017

    2089

  21. [21]

    GRAHAM, E.A

    I.G. GRAHAM, E.A. SPENCE ANDJ. ZOU,Domain Decomposition with local impedance conditions for the Helmholtz equation with absorption, SIAM J. Numer. Anal. 58, 2515–2543, 2020

    2020

  22. [22]

    HEINLEIN, A

    A. HEINLEIN, A. KLAWONN, J. KNEPPER ANDO. RHEINBACH,Adaptive GDSW coarse spaces for overlapping Schwarz methods in three dimensions, SIAM J. Sci. Comput. 41, A3045–A3072, 2019

    2019

  23. [23]

    Q. HU ANDZ. LI,A novel coarse space applying to the weighted Schwarz method for Helmholtz equations, https://arxiv.org/abs/2402.06905, 2024

    arXiv 2024

  24. [24]

    SARKIS,Restricted overlapping balancing domain decomposition methods and restricted coarse problems for the Helmholtz problem, Comput

    J-H KIMN ANDM. SARKIS,Restricted overlapping balancing domain decomposition methods and restricted coarse problems for the Helmholtz problem, Comput. Methods Appl. Mech. Engrg. 196, 1507–1514, 2007

    2007

  25. [25]

    LEONI, A First Course in Sobolev Spaces, American Mathematical Soc., 2017

    G. LEONI, A First Course in Sobolev Spaces, American Mathematical Soc., 2017

    2017

  26. [26]

    P. LU, X. XU, B. ZHENG ANDJ. ZOU,Two-level hybrid Schwarz Preconditioners for the Helmholtz Equation with high wave number, SIAM J. Numer. Anal. 63, 2187–2220, https://doi.org/10.1137/24M168533X, 2025

    work page doi:10.1137/24m168533x 2025

  27. [27]

    C. MA, C. ALBER ANDR. SCHEICHL,Wavenumber explicit convergence of a multiscale generalized finite element method for heterogeneous Helmholtz problems, SIAM J. Numer. Anal. 61, 1546–1584, https://doi. org/10.1137/21M1466748, 2023

    work page doi:10.1137/21m1466748 2023

  28. [28]

    C. MA, C. ALBER, R. SCHEICHL ANDY. ZHANG,Two-level restricted additive Schwarz preconditioner based on multiscale spectral generalized FEM for heterogeneous Helmholtz problems, J. Sci. Comput. 105, 99, https://doi.org/10.1007/s10915-025-03138-y, 2025

    work page doi:10.1007/s10915-025-03138-y 2025

  29. [29]

    MELENK ANDS

    J.M. MELENK ANDS. SAUTER,Convergence analysis for finite element discretizations of the Helmholtz equation with Dirichlet-to-Neumann boundary conditions, Math. Comp. 79, 1871–1914, https://doi.org/10. 1090/S0025-5718-10-02362-8, 2010

    1914

  30. [30]

    PETERSEIM.Eliminating the pollution effect in Helmholtz problems by local subscale correction, Math

    D. PETERSEIM.Eliminating the pollution effect in Helmholtz problems by local subscale correction, Math. Comp. 86, 1005–1036, https://doi.org/10.1090/mcom/3156, 2017

    work page doi:10.1090/mcom/3156 2017

  31. [31]

    SCHATZ,An observation concerning Ritz–Galerkin methods with indefinite bilinear forms, Math

    A.H. SCHATZ,An observation concerning Ritz–Galerkin methods with indefinite bilinear forms, Math. Comput. 28, 959–962, https://doi.org/10.1090/S0025-5718-96-00649-7, 1974

    work page doi:10.1090/s0025-5718-96-00649-7 1974

  32. [32]

    SCHATZ ANDJ.P

    A.H. SCHATZ ANDJ.P. WANG,Some new error estimates for Ritz–Galerkin methods with minimal regularity assumptions, Math. Comput., 213(65), 19–27, 1996

    1996

  33. [33]

    SPILLANE,Toward a new fully algebraic preconditioner for symmetric positive definite problems, In Domain Decomposition Methods in Science and Engineering XXVI, S

    N. SPILLANE,Toward a new fully algebraic preconditioner for symmetric positive definite problems, In Domain Decomposition Methods in Science and Engineering XXVI, S. Brenner et al., Eds.; Springer: Cham, Switzerland, 2022; 745–752, https://doi.org/10.1007/978-3-030-95025-5 81

    work page doi:10.1007/978-3-030-95025-5 2022

  34. [34]

    SPILLANE, F

    N. SPILLANE, F. NATAF, V. DOLEAN, P. HAURET, C. PECHSTEIN ANDR. SCHEICHL,Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps, Numer. Math. 4, 741–740, https://doi.org/10.1007/s00211-013-0576-y, 2014

    work page doi:10.1007/s00211-013-0576-y 2014

  35. [35]

    TOSELLI ANDO

    A. TOSELLI ANDO. WIDLUND,Domain Decomposition Methods: Algorithms and Theory, Springer Series in Computational Mathematics, Springer-Verlag, 2005

    2005