Next Article in Journal
Challenges in Kinetic-Kinematic Driven Musculoskeletal Subject-Specific Infant Modeling
Next Article in Special Issue
On the Convergence of the Damped Additive Schwarz Methods and the Subdomain Coloring
Previous Article in Journal
Reduced Order Modeling Using Advection-Aware Autoencoders
Previous Article in Special Issue
Learning Adaptive Coarse Spaces of BDDC Algorithms for Stochastic Elliptic Problems with Oscillatory and High Contrast Coefficients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Can DtN and GenEO Coarse Spaces Be Sufficiently Robust for Heterogeneous Helmholtz Problems?

by
Niall Bootland
1,* and
Victorita Dolean
1,2
1
Department of Mathematics and Statistics, University of Strathclyde, Glasgow G1 1XH, UK
2
Laboratoire J.A. Dieudonné, Centre National de la Recherche Scientifique (CNRS), University Côte d’Azur, 06000 Nice, France
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2022, 27(3), 35; https://doi.org/10.3390/mca27030035
Submission received: 17 February 2022 / Revised: 19 April 2022 / Accepted: 19 April 2022 / Published: 21 April 2022
(This article belongs to the Special Issue Domain Decomposition Methods)

Abstract

:
Numerical solutions of heterogeneous Helmholtz problems present various computational challenges, with descriptive theory remaining out of reach for many popular approaches. Robustness and scalability are key for practical and reliable solvers in large-scale applications, especially for large wave number problems. In this work, we explore the use of a GenEO-type coarse space to build a two-level additive Schwarz method applicable to highly indefinite Helmholtz problems. Through a range of numerical tests on a 2D model problem, discretised by finite elements on pollution-free meshes, we observe robust convergence, iteration counts that do not increase with the wave number, and good scalability of our approach. We further provide results showing a favourable comparison with the DtN coarse space. Our numerical study shows promise that our solver methodology can be effective for challenging heterogeneous applications.

1. Introduction

Consider solving the heterogeneous Helmholtz problem: for a bounded domain Ω R d , d = 2 , 3 , we wish to find u ( x ) : Ω C such that
Δ u k 2 u = f in Ω ,
C ( u ) = 0 on Ω ,
where C incorporates some appropriate boundary conditions. For the heterogeneous problem, we suppose the wave number k ( x ) > 0 is a function of space, defined by the ratio k = ω / c of the angular frequency ω and the wave speed c ( x ) . In this work, we investigate the use of an overlapping Schwarz preconditioner with a suitably chosen coarse space based on solving local eigenvalue problems on each subdomain.
The ability to compute solutions to the Helmholtz problem (1) is important across many disciplines of science and engineering. As the prototypical model for frequency-domain wave propagation, it features within the fields of optics, acoustics, and seismology, amongst others. Further applications can be found in imaging science, such as through medical imaging techniques and geophysical studies of the Earth’s subsurface within, for instance, the oil industry. Nonetheless, it is challenging to develop efficient computational methods to solve (1), particularly when the wave number k becomes large.
Discretisation of (1) by standard approaches, such as Lagrange finite elements, as we shall use here, results in large linear systems to be solved which are indefinite, non-self-adjoint, and ill-conditioned; see also [1]. These systems present various difficulties to solve, especially in the presence of complex heterogeneities, at high frequencies (large k), or when solutions include many wavelengths in the domain. Classical methods for solving such large systems typically fail for several reasons, as detailed in [2,3], and specialist approaches must be employed for a robust solver. While much progress has been made for symmetric positive definite problems, such techniques cannot be applied out-of-the-box and extensions to tackle indefinite and non-self-adjoint problems may not be clear. This has led to a number of approaches being developed in recent years aiming to bridge this gap. For the Helmholtz problem, this includes parallel direct solvers, such as [4,5], and preconditioned iterative methods that utilise multigrids, such as [6,7], or the so-called “shifted Laplace” approach [8,9]. This latter (complex) shifted Laplace preconditioner has seen much interest into its practical use [10] and further developments through deflation techniques, most recently in [11,12].
Another broad class of solvers are domain decomposition methods, which provide a natural balance between using direct and iterative solvers. Specialist methods are again required for the Helmholtz problem, a popular set of which fall under the heading of “sweeping” methods [13,14]. These are multiplicative domain decomposition methods linked to a variety of other approaches, including optimised Schwarz methods, as detailed in the recent survey [3]. While sweeping is a conceptually serial approach, much work has been done to incorporate parallelism. Of particular note is the “L-Sweeps” method [15], stated to be the first parallel scalable preconditioner for high-frequency Helmholtz problems; a review of developments in this area is provided in the introductions of [15,16]. Other popular approaches are FETI methods; for the Helmholtz problem, these include FETI-H [17] and FETI-DPH [18].
Within the domain decomposition community, there has also been renewed work on additive Schwarz methods, which offer a naturally parallel approach. Following on from the seminal work [19], which utilised Robin (or impedance) transmission conditions to provide a convergent Schwarz method for the Helmholtz problem, a wealth of non-overlapping Schwarz methods have been devised; see the introduction of [20] for a recent overview. In these methods, one has to be careful to either avoid or treat cross points (where three or more subdomains meet), as can be done in the robust treatment of [20]. Many optimised approaches rely on deriving higher-order transmission conditions, such as through second-order impedance operators in [21], absorbing boundary conditions (ABCs) [22], or non-local operators [23]. Ideally, one would use the Dirichlet-to-Neumann (DtN) map (a Poincaré–Steklov operator) to provide transparent transmission conditions, but this is prohibitive in practice, and so these optimised Schwarz methods in essence try to approximate this operator.
Overlapping Schwarz methods for the Helmholtz problem—see for example [24,25,26]—have also received renewed attention in recent years, and it is this type of method we shall consider. A successful approach is to design additive Schwarz methods based on including absorption (a complex shift k 2 k 2 + i ε ), with absorption parameter ε ; see [27,28]. Theoretical work to understand the effectiveness of this approach can be found in [29,30] and for the heterogeneous problem in [31]; see also [32]. To be scalable, such additive Schwarz methods require a second level, known as a coarse space (see [33] for a novel analysis for the absorptive problem). For the two-level methods in [27,29], this is provided through a coarse grid—an approach that is effective also for the time-harmonic Maxwell problem [34].
As well as to provide scalability, coarse spaces have been devised to provide robustness to heterogeneity. This is exemplified by the “Generalised Eigenproblems in the Overlap” (GenEO) approach for symmetric positive definite (SPD) problems [35]. This approach provides a spectral coarse space, where appropriate local eigenvalue problems are solved to provide a two-level method. Another spectral coarse space is the DtN coarse space [36], which has been extended and investigated for the Helmholtz problem in [37,38]. While the standard GenEO theory applies only in the SPD case (see [39]), in this work, we develop and explore a GenEO-type method for the Helmholtz problem and show numerically that, for a 2D model problem of a wave guide, it provides a scalable approach that is robust to heterogeneity and increasing wave number in terms of the iteration count of a preconditioned GMRES method. Companion results for large benchmark problems on coarsely resolved meshes that typically arise in applications, along with comparisons to other methods, are found in [40]. Theoretical results on variants of DtN and GenEO for Helmholtz problems are currently out of reach, but promising numerical results, based on heuristics, were obtained in [37,38], showing the potential of these methods in practice.
The primary aim of this work is to explore the utility of a GenEO-type method for the heterogeneous Helmholtz problem (1); we call this approach H-GenEO. In particular, we highlight the following contributions:
  • We present a range of numerical tests, on pollution-free meshes, comparing our proposed H-GenEO approach with another spectral coarse space applicable to the Helmholtz problem, namely the DtN method.
  • We investigate the use of appropriate thresholding for the required generalised eigenproblems in both the DtN and H-GenEO coarse spaces.
  • We consider robustness to non-uniform decomposition, heterogeneity, and increasing wave number as well as the scalability of the methods. We find that only the H-GenEO approach is scalable and robust to all of these factors for a 2D model problem.
  • We provide both weak and strong scalability tests for H-GenEO applied to high wave number problems.
The remainder of this work is structured as follows. We begin by considering a finite element discretisation of the Helmholtz problem in Section 2.1, before outlining the underlying domain decomposition methodology we use in Section 2.2. The main topic of interest, that of suitable spectral coarse spaces for the heterogeneous Helmholtz problem, is then detailed in Section 2.3. Extensive numerical results on a 2D model problem are provided in Section 3, along with a discussion of our findings. Finally, we draw our conclusions in Section 4.

2. Materials and Methods

2.1. Finite Element Discretisation

The problem we consider in this work is the interior Helmholtz problem (1), for which we must prescribe appropriate boundary conditions. In practical applications, the computational domain Ω is often truncated, and the physically relevant condition, namely the far field Sommerfeld radiation condition, must be approximated on the non-physical boundary of Ω . This allows for appropriate wave behaviour to be modelled in a bounded domain. This simplest approximation that is widely used is that of a Robin (or impedance) condition, and this is what we shall consider; other approaches include ABCs [41] or perfectly matched layers (PML) [42,43]. We also suppose that Dirichlet conditions may be imposed on a boundary Γ D Ω , with the Robin condition on the remaining boundary Γ R = Ω \ Γ D . Thus, in general, we seek the solution of the boundary value problem
Δ u k 2 u = f in Ω ,
u = u Γ D on Γ D ,
u n + i k u = 0 on Γ R ,
where the forcing function f ( x ) incorporates any sources in the domain. Note that if Γ R , the problem is well posed, but if Γ R = , the problem is ill-posed for certain choices of k related to eigenfunctions of the Laplacian.
To discretise (2), we use standard Lagrange finite elements; the details can be found in, for example, [40], and so we provide only an outline here. Defining the relevant trial and test spaces V = u H 1 ( Ω ) : u = u Γ D on Γ D and V 0 = u H 1 ( Ω ) : u = 0 on Γ D , the weak formulation of (2) is to find u V such that
a ( u , v ) = F ( v ) v V 0 ,
where
a ( u , v ) = Ω u · v ¯ k 2 u v ¯ d x + Γ R i k u v ¯ d s and F ( v ) = Ω f v ¯ d x .
Assuming a simplicial mesh T h of Ω with characteristic element diameter h, piecewise polynomial finite element approximation reduces the problem to solving the complex linear system
A u = f ,
with coefficient matrix A C n × n and right-hand side vector f C n stemming from a ( · , · ) and F ( · ) , respectively; see [40]. For accurate discretisation of the Helmholtz problem, the number of degrees of freedom n is required to be large, especially for large k. Indeed, to maintain the same level of accuracy of discrete solutions as k increases, the number of mesh points must increase faster than k, due to the pollution effect [44]. This depends on the polynomial order of the approximation used: for instance, using piecewise linear (P1) finite elements, k 3 h 2 must be bounded, and so h must shrink as O ( k 3 / 2 ) . Higher order finite elements can reduce this restriction on h, but ultimately the interpolation properties of such methods degrade. Here, we utilise standard P1 elements and maintain a discretisation such that k 3 h 2 is fixed to avoid the pollution effect.

2.2. Underlying Domain Decomposition Method

To solve the discrete Helmholtz problem (5), we use GMRES accelerated via a two-level overlapping domain decomposition preconditioner. For the underlying one-level method, we consider the optimised restricted additive Schwarz (ORAS) method [45] (sometimes known as WRAS-H [26]; see discussions in [31,32] and [27,29] for a related IMPHRAS1 method). The key difference in ORAS compared with standard additive Schwarz methods such as RAS is that the local Dirichlet problems on subdomains are replaced by appropriate Robin problems.
To formulate the domain decomposition preconditioner, we first suppose that Ω is decomposed into non-overlapping subdomains Ω s s = 1 N , assumed to be resolved by the mesh T h . To give an overlapping decomposition, a layer of adjoining mesh elements is added to give overlapping subdomains Ω s s = 1 N by way of the extension
Ω s = Int supp ( ϕ j ) Ω s supp ( ϕ j ) ,
where ϕ j j = 1 n are the nodal basis functions of the finite element space, Int ( · ) denotes the interior of a domain, and supp ( · ) is the support of a function. Further layers of elements can be added in a recursive manner in order to obtain subdomains with larger overlap, if desired.
Once we have an overlapping decomposition into subdomains Ω s s = 1 N , we define the required operators for the Schwarz preconditioner. We let R s R n s × n be the discrete form of the restriction operator, restricting functions to the subdomain Ω s , where n s is the number of degrees of freedom in Ω s . The corresponding extension operator, R s T , then acts as an extension by zero outside of Ω s . We also utilise a partition of unity, with the discrete form of a diagonal matrix D s R n s × n s satisfying s = 1 N R s T D s R s = I , which appropriately scales the multiple subdomain contributions in the overlapping regions. Finally, within ORAS, we require the solution of local Robin problems given by
Δ w s k 2 w s = f in Ω s ,
w s n s + i k w s = 0 on Ω s \ Ω ,
C ( w s ) = 0 on Ω s Ω ,
with C representing the underlying problem boundary conditions on Ω , namely (2b) and (2c). Note that the use of the Robin condition (7b) ensures the solvability of these local problems. Defining the equivalent finite element discretisation of (7) to be given by the stiffness matrix A ^ s R n s × n s , the construction of the one-level ORAS preconditioner is given by the sum
M ORAS 1 = s = 1 N R s T D s A ^ s 1 R s .
Note that the local solutions on each subdomain given by A ^ s 1 can be carried out in parallel.
In order to provide robustness and scalability, the ORAS method above must be augmented by the suse of a coarse space to provide a two-level method. A coarse space can be thought of as a collection of linearly independent column vectors Z. The vectors that are incorporated into Z are key to providing scalability, especially for indefinite problems, such as the Helmholtz problems we solve here, where the addition of a coarse space need not improve the performance of the underlying one-level method [46]. There are several ways to incorporate the coarse space; here, we consider an effective approach that is based on deflation. For this, a coarse space operator E = Z A Z is constructed as well as the coarse correction operator Q = Z E 1 Z , which we incorporate to give a two-level ORAS method
M ORAS , 2 1 = M ORAS 1 ( I A Q ) + Q .
We now turn our attention the choice of coarse space.

2.3. Spectral Coarse Spaces

In this work, we consider and explore spectral coarse spaces for the discrete Helmholtz problem (5). These utilise local eigenvalue problems on subdomains in order to build a global coarse space. We review the DtN coarse space [38] before detailing a new coarse space for Helmholtz problems based on GenEO technology [35]. We then show a link between these two approaches.
Remark 1
(Notation). We utilise the following notation for local Dirichlet, Robin, and Neumann matrices: for a variational problem that gives rise to a system matrix B, we denote by B s the corresponding local Dirichlet matrix on Ω s . In the case that Robin conditions are used on internal subdomain interfaces, the local problem matrix is denoted by B ^ s . On the other hand, if Neumann conditions are used on such interfaces, we denote the local matrix by B ˜ s .

2.3.1. The DtN Coarse Space

The Dirichlet-to-Neumann (DtN) coarse space, first studied in [36,47] for elliptic problems, is based on solving local eigenvalue problems on subdomain boundaries related to a DtN map. Harmonic extensions of low-frequency modes on subdomains are then used to provide a coarse space. In order to define this approach for the Helmholtz problem, as explored in [38], we first require the Helmholtz extension operator from the subdomain boundary Ω s .
On each subdomain, let Γ s = Ω s \ Ω , and suppose we have Dirichlet data v Γ s on Γ s . The Helmholtz extension v into Ω s is given by solving
Δ v k 2 v = 0 in Ω s ,
v = v Γ s on Γ s ,
C ( v ) = 0 on Ω s Ω ,
where C ( v ) = 0 represents the original problem boundary conditions, as in (7c). The DtN map takes Dirichlet data v Γ s to the corresponding Neumann data on Γ s , namely
DtN Ω s ( v Γ s ) = v n Γ s
where v is the Helmholtz extension defined by (10). The associated local DtN eigenproblem on subdomain Ω s is given by
DtN Ω s ( u Γ s ) = λ u Γ s ,
for eigenfunctions u Γ s and eigenvalues λ C . In order to build the coarse space, we take the Helmholtz extension of u Γ s in Ω s and extend by zero into the whole domain Ω using the partition of unity; see [38].
To formulate the discrete version of the eigenproblems to be solved, we require the coefficient matrices A ˜ s , corresponding to local Neumann problems on Ω s with boundary conditions C = 0 on Ω s Ω , similar to that supplying the local Robin problems in (7). Furthermore, we must distinguish between degrees of freedom on the boundary and the interior of the subdomain Ω s , and so we let Γ s and I s be the set of indices on the boundary and interior, respectively. Recalling that ϕ j are our nodal basis functions, we also define
M Γ s = Γ s ϕ j ϕ i i , j Γ s
to be the mass matrix on the subdomain interface. Using standard block notation to denote submatrices of A s and A ˜ s , the discrete DtN eigenproblem can be written as
A ˜ Γ s , Γ s A Γ s , I s A I s , I s 1 A I s , Γ s u Γ s = λ M Γ s u Γ s .
We then make use of the Helmholtz extension of u Γ s to degrees of freedom in I s given by u I s = A I s , I s 1 A I s , Γ s u Γ s . Letting u s denote the complete local vector representing the Helmholtz extension, the corresponding global vector that enters the coarse space Z is R s T D s u s . Further motivation and details on the DtN eigenproblems can be found in [38].
What remains is to determine which eigenvectors of (14) should be incorporated into the coarse space. A variety of selection criteria were investigated in [38], which made it clear that the best choice was to select eigenvectors corresponding to eigenvalues with the smallest real part. That is, a threshold on the abscissa η = Re ( λ ) should be used, namely
η < η max ,
where η max depends on k s = max x Ω s k ( x ) . The choice η max = k s is advocated in [38]; however, we recently showed that taking a slightly larger threshold η max = k s 4 / 3 can be beneficial in certain cases in order to gain robustness to the wave number [37]. Unfortunately, this only occurs for the homogeneous problem with sufficiently uniform subdomains. To construct a more robust coarse space, we build upon the GenEO approach.

2.3.2. The GenEO Coarse Space

The Generalised Eigenproblems in the Overlap (GenEO) coarse space was derived in [35] to provide a rigorously robust approach for symmetric positive definite problems even in the presence of heterogeneities. In recent years, this approach has been extended and used within various settings and applications—for example, [39,45,48,49]; see also the discussion on developments for other spectral coarse spaces in [50].
Within the original derivation [35], the generalised eigenproblems are defined in a variational framework on the assumption that a ( · , · ) is a symmetric and coercive bilinear form. On a subdomain Ω s , an overlapping zone Ω s is defined as the parts of Ω s that overlap with another subdomain, and the local eigenproblem
a Ω s ( u , v ) = λ a Ω s Ξ s ( u ) , Ξ s ( v ) v V ( Ω s ) ,
is solved for small λ , where Ξ s represents the action of the partition of unity operator on Ω s . To be clear, a D ( · , · ) represents the underlying variational problem on the domain D with problem boundary conditions on Ω and natural (“do nothing”) conditions on D \ Ω . The eigenproblem (16) provides an appropriate link in order to bound the condition number of the preconditioned operator independently of the heterogeneity and number of subdomains. This bound depends on the smallest eigenvalue λ whose corresponding eigenfunction is not incorporated into the GenEO coarse space. Hence, to achieve a desired rate of convergence, all eigenfunctions corresponding to eigenvalues smaller than a threshold (say, λ < λ max ) must be computed. Alternatively, one may opt in practice to compute a fixed number of eigenfunctions per subdomain and use these in order to accelerate convergence.
The restriction within the right-hand side of (16) to the overlapping zone is not an essential requirement, and alternative formulations can be used. In order to remove the need to track overlap regions, one possibility is to replace Ω s with the whole of Ω s , as in [45,48]. When formulating the discrete eigenproblem, this requires only the local Neumann matrix A ˜ s to be constructed. In this case, we must solve
A ˜ s u = λ D s A s D s u ,
where A s = R s A R s T is the local Dirichlet matrix—a sub-matrix of A. As with the DtN method, the vectors which then go into the coarse space Z are R s T D s u . It is this form of the GenEO eigenproblem in (17) that we shall build upon to develop an approach tailored to the heterogeneous Helmholtz problem.

2.3.3. H-GenEO: A GenEO-Type Coarse Space for Helmholtz Problems

In consideration of GenEO approaches for the Helmholtz problem, a key hurdle is the loss of operators that are definite or self-adjoint. While some progress has been made at the theoretical level for closely related problems [39,51], currently available analysis has yet to overcome all the challenges present for the Helmholtz problem, and rigorous justification for a choice of spectral coarse space remains out of reach. The approach taken in [39,51] is to formulate the GenEO eigenproblem for a nearby symmetric positive definite problem, here corresponding to a Laplace problem and hence given the name Δ -GenEO. This problem can be given by setting k = 0 in (2) and (4). Letting L s be the local Dirichlet matrix for this problem in Ω s and L ˜ s be the equivalent Neumann matrix, the Δ -GenEO eigenproblem is given by
L ˜ s u = λ D s L s D s u .
This is a positive (semi-)definite eigenproblem with real non-negative eigenvalues λ and, as such, eigenvectors can be chosen in the standard way using a threshold λ < λ max . Unfortunately, for the Helmholtz problem, the Δ -GenEO approach can perform rather poorly when k becomes large (see Table 1 in Section 3), as might be anticipated from the fact that solutions to the Laplace problem differ considerably to those of the Helmholtz problem in this range.
In order to provide an appropriate spectral coarse space for the Helmholtz problem, it stands to reason that a Helmholtz operator must be included. If we try to apply the GenEO eigenproblem (17) as is, with matrices stemming from the Helmholtz bilinear form in (4), we must first note that the problem is non-self-adjoint and, as such, eigenvalues λ are no longer real in general. As a threshold criterion, as with the DtN approach, we can consider the abscissa η = Re ( λ ) instead and seek eigenvectors corresponding to η < η max . Unfortunately, this formulation fails to be robust, and numerically, we have found that the eigensolver we use can often break down. Given that, without any relevant theory, it is no longer clear that (17) provides an appropriate eigenproblem; thus, we develop a different approach that yields a robust method in our numerical experiments.
While the above approaches consider only Helmholtz, or only Laplace, operators in their formulation, we instead link the underlying Helmholtz problem to the positive definite Laplace problem. Since this GenEO-type method targets the Helmholtz problem, we call it “H-GenEO”. The local eigenproblem utilised is given by
A ˜ s u = λ D s L s D s u .
Since eigenvalues are complex (though, in our experience, they tend to cluster close to the real line) we threshold based on the abscissa η < η max . We see the effectiveness of the H-GenEO coarse space for a 2D model problem in Section 3. When used with a threshold (typically, we use η max = 1 2 ), we observe that the number of GMRES iterations for the H-GenEO approach is independent of the wave number k, which is not true of the DtN method in general. Before continuing to our numerical results, we first show a link between the DtN and GenEO eigenproblems and how the H-GenEO method is related.

2.3.4. A Link between DtN and GenEO

While the DtN and GenEO eigenproblems look rather different, here we show a link between the two approaches. In particular, consider the GenEO method when we remove the partition of unity matrices D s from (17), so that we are left with the Neumann matrix on the left and the Dirichlet matrix on the right, namely A ˜ s u = λ A s u . Further, let us use subscripts O s and I s to denote overlap and interior degrees of freedom in Ω s , respectively. Then, moving all terms to the left-hand side, the GenEO eigenproblem can be written as
A ˜ O s , O s λ A O s , O s ( 1 λ ) A O s , I s ( 1 λ ) A I s , O s ( 1 λ ) A I s , I s u O s u I s = 0 0 ,
where we have used the fact that the Neumann matrix is equal to Dirichlet matrix except in the block associated with the overlap degrees of freedom. Forming the Schur complement of (20) with respect to the A I s , I s block and dividing by 1 λ , we obtain
A ˜ O s , O s A O s , I s A I s , I s 1 A I s , O s u O s = λ 1 λ A O s , O s A ˜ O s , O s u O s .
Now, if the overlap degrees of freedom O s are precisely those used as boundary degrees of freedom Γ s in the DtN method, we see that (21) resembles the DtN eigenproblem (14). The primary difference stems from the right-hand side, where we now have the difference between the Dirichlet and Neumann matrices on the boundary degrees of freedom as opposed to a mass matrix, though we note that these can coincide for certain choices of discretisation; for example, a simple finite difference scheme on a Cartesian grid where forward or backward differences are used to approximate the Neumann condition. The other difference, aside from the partition of unity in the true GenEO approach, stems from a transformation of the eigenvalues, namely the comparative eigenvalues for the DtN method are μ = λ 1 λ .
Now consider the H-GenEO eigenproblem (19) without the partition of unity matrices. Splitting for the moment A = L k 2 M , we have
L ˜ O s , O s λ L O s , O s k 2 M O s , O s ( 1 λ ) L O s , I s k 2 M O s , I s ( 1 λ ) L I s , O s k 2 M I s , O s ( 1 λ ) L I s , I s k 2 M I s , I s u O s u I s = 0 0 .
Now, dividing by 1 λ and defining κ 2 = k 2 1 λ and the corresponding Helmholtz matrices with wave number κ as B = L κ 2 M , we have
B ˜ O s , O s + λ 1 λ L ˜ O s , O s L O s , O s B O s , I s B I s , O s B I s , I s u O s u I s = 0 0 .
Using the fact that L ˜ O s , O s L O s , O s = B ˜ O s , O s B O s , O s , the resulting Schur complement system is then given as
B ˜ O s , O s B O s , I s B I s , I s 1 B I s , O s u O s = λ 1 λ B O s , O s B ˜ O s , O s u O s .
Thus, we see that H-GenEO corresponds to solving GenEO problems based on a wave number κ that varies with the eigenvalue. Since we are predominantly looking for small eigenvalues λ , this wave number κ = k ( 1 λ ) 1 / 2 will then be close to k.
To further exhibit the link between these spectral coarse spaces for the Helmholtz problem, we consider example eigenfunctions for both the DtN (14) and H-GenEO (19) eigenproblems on the central subdomain of a 5 × 5 decomposition into square subdomains for a homogeneous model problem with k = 46.5 (see Section 3 for details). Figure 1 displays a selection of eigenfunctions, and in the top two rows, we plot DtN (top row) and H-GenEO (middle row) eigenfunctions, which show the same features, visually being very similar. Note that, since the central subdomain does not touch the Robin boundary Γ R , both eigenproblems are real and symmetric, albeit indefinite, and so all eigenvalues λ are real. With DtN, as λ increases, variation in the eigenfunctions tends to be restricted to the boundary, as can be seen in Figure 1c,e. Such behaviour is observed for many H-GenEO eigenfunction as well; however, we also obtain distinct eigenfunctions which are not found amongst DtN eigenfunctions: examples are given in the bottom row of Figure 1, and we note that they tend to exhibit large variation in the interior of the subdomain. Nonetheless, there is a clear link between many of the DtN and H-GenEO eigenfunctions.

3. Results and Discussion

In this section, we present and discuss the numerical results computed using FreeFEM [52], in particular through the functionality of ffddm, which handles the underlying domain decomposition data structures. As a model problem, we consider the case of a wave guide in 2D, defined on the unit square Ω = ( 0 , 1 ) 2 . We impose homogeneous Dirichlet conditions on two opposite sides, namely (2b) with u Γ D = 0 on Γ D = { 0 , 1 } × [ 0 , 1 ] , and Robin conditions on the two remaining sides, that is (2c) on Γ R = [ 0 , 1 ] × { 0 , 1 } . A point source is located in the centre of the domain at ( 1 2 , 1 2 ) and provides the forcing function f. A schematic of this model problem is found in Figure 2.
To discretise the problem, we triangulate Ω using a Cartesian grid with spacing h and alternating diagonals to form a simplicial mesh (see Figure 2 where h = 1 4 ). The discrete problem (5) is then built using a P1 finite element approximation on this mesh. In order to avoid the pollution effect, we choose k and h simultaneously so that k 3 h 2 = 2 π 10 is fixed. The large sparse linear system (5) is solved using right-preconditioned GMRES, with the preconditioner given by the two-level ORAS method (9) and the choice of coarse space as stated. We terminate the GMRES iteration once a relative residual tolerance of 10 6 is reached. Unless otherwise stated, within the domain decomposition preconditioner, we use minimal overlap: that is, one layer of adjoining mesh elements are added to the non-overlapping subdomains via the extension in (6). To solve the local eigenvalue problems in the two-level methods, we make use of ARPACK [53], while both the subdomain solves and coarse space operator solves are given by MUMPS [54].
We first compare coarse spaces in the simplest case of a homogeneous problem, investigating the choice of eigenvalue threshold used. We then show results on how the methods perform in a variety of settings; for instance, with non-uniform subdomains, heterogeneity, and additional overlap. Finally, both weak and strong scalability tests are performed for H-GenEO applied to high wave number problems with timings reported.

3.1. A Comparison of Methods for the Homogeneous Problem with Uniform Partitioning

In Table 1, we give some benchmark results for the simplest problem of a homogeneous wave guide using a uniform decomposition into 25 square subdomains. We see that the one-level ORAS method (8) performs relatively poorly as the wave number k increases. The standard DtN coarse space (14) (with η max = k ) is able to reduce iteration counts with a relatively small coarse space size; however, there is still a clear increase in iterations as k increases. The Δ -GenEO method (18), with λ max = 1 2 , performs poorly here, often performing worse than the one-level method despite a larger coarse space than the DtN approach; this may be because the impedance conditions from the wave guide problem are not included in the definition of the Δ -GenEO coarse space and so the eigenfunctions are not appropriate here. Finally, the standard H-GenEO method (19) (with η max = 1 2 ) performs well and significantly reduces the iteration counts, by a factor of 10 for the largest wave number, and provides robustness to increasing wave number k (in fact, iteration counts tend to decrease with k). We note that the size of the H-GenEO coarse space is larger than the DtN coarse space, and we now explore this further.
In Table 2, we provide results for both the DtN and H-GenEO methods with differing eigenvalue thresholds η max . For DtN, we use the standard threshold η max = k , the suggested threshold from [37] η max = k 4 / 3 , and the larger threshold of η max = k 3 / 2 . For H-GenEO, we use the standard threshold η max = 1 2 as well as the weaker thresholds of 1 4 and 1 8 , the latter giving coarse space sizes more comparable to the standard DtN approach. To differentiate between these methods, we use the notation DtN( η max ) and H-GenEO( η max ) where η max is as specified. We notice in Table 2 that increasing the DtN threshold to η max = k 4 / 3 significantly improves the iteration counts for this problem so that they are almost independent of k, albeit very slightly growing. Increasing the threshold further only marginally improves the iteration counts, which still grow slightly with k, but at the expense of a coarse space almost twice the size. On the other hand, if we relax the H-GenEO threshold, we start to see higher iteration counts and lose some robustness, but generally iteration counts do not increase with k as they do for the standard DtN method. We note that, for this homogeneous problem, roughly comparable coarse space sizes give approximately similar iteration counts, and so it is primarily the thresholds used in DtN and H-GenEO that dictate the different growth behaviour we observe.
To explore the growth in the size of the coarse space further, in Figure 3a, we plot coarse space size against the wave number k. From this, we can see that growth for DtN( k 4 / 3 ) is approximately proportional to k 4 / 3 , while for H-GenEO( 1 2 ), it is around k 5 / 3 for our model problem. When the thresholds are relaxed, the precise relationship becomes less clear, but we note that the coarse space sizes grows more slowly with a weaker threshold, especially for the DtN approach. The faster growth seen for DtN( k 4 / 3 ) and H-GenEO( 1 2 ) may help to accommodate the stronger robustness to k observed in the iteration counts of Table 2. These two approaches appear to provide the best trade-off for obtaining a well-behaved method, and so we focus primarily on these approaches, but first we consider the question of scalability.

3.2. Scalability of DtN and H-GenEO for the Homogeneous Problem with Uniform Partitioning

We now investigate the scalability of the DtN and H-GenEO methods. This will depend on the threshold used, and so we compare results for DtN(k) and H-GenEO( 1 8 ) as well as DtN( k 4 / 3 ) and H-GenEO( 1 2 ), with each pair of approaches giving broadly similar iteration counts. Results for the homogeneous problem with k = 73.8 and h 1 = 800 are given in Table 3 for an increasing number of uniform square subdomains N. We see that the DtN(k) approach does not exhibit scalability here, with iteration counts that noticeably increase with N. Similarly, with the weaker threshold, H-GenEO( 1 8 ) also fails to be scalable. On the other hand, both DtN( k 4 / 3 ) and H-GenEO( 1 2 ) are scalable here, with low iteration counts that vary little with N.
Comparing the size of the coarse spaces employed, we see that the DtN( k 4 / 3 ) coarse space grows faster with N and becomes larger than the H-GenEO( 1 2 ) coarse space, which may account for its particularly strong robustness to N here. From Figure 3b, we see that for both DtN approaches, the coarse space size grows approximately proportional to N 0.6 , while for the H-GenEO method, it is around N 0.4 . This suggests that H-GenEO may be advantageous when N becomes large due to the smaller coarse space required. For each method, the average number of eigenvectors taken per subdomain decreases as N increases; this means that, as well as smaller subdomain solves, we also benefit from requiring fewer eigenvectors to be computed per subdomain, even if the global coarse space size increases. We note that the size of the DtN(k) or H-GenEO( 1 8 ) coarse space sometimes shrinks as we increase N, and in these cases, the iteration counts tend to be particularly poor, suggesting that the thresholds of η max = k for DtN and η max = 1 8 are not doing a suitable job in capturing the eigenfunctions required for scalability. As such, we now narrow our focus to the DtN( k 4 / 3 ) and H-GenEO( 1 2 ) approaches.

3.3. Robustness of DtN and H-GenEO for the Homogeneous Problem with METIS Decomposition

We now consider utilising non-uniform subdomains, as provided through the software METIS [55]. In Table 4, we compare DtN( k 4 / 3 ) and H-GenEO( 1 2 ) in this situation, again for the homogeneous problem. We observe that both methods retain their robustness, meaning that iteration counts only depend mildly on the wave number k and the number of subdomains N. For H-GenEO( 1 2 ), the scalability becomes more favourable for larger k, and we can see that for the smallest wave numbers with large N, the average number of eigenvectors per subdomain becomes very small; in fact, on many subdomains, only a single eigenvector is taken, and so the achieved tolerance on the eigenvalue may be somewhat weaker than 1 2 , which may explain the slightly poorer performance. One way this could be overcome is by always taking a minimum number of eigenvectors per subdomain; for instance, using at least 5 eigenvectors per subdomain when k = 18.5 gives iteration counts bounded by 19. Since we are primarily interested in approaches that remain effective for increasingly large wave numbers, where this issue diminishes, we do not worry further about this and assert that for problems of interest, H-GenEO( 1 2 ) provides good scalability. On the other hand, for DtN( k 4 / 3 ), the mild increase in iteration count is seen as k increases here, with scalability observed for all wave numbers.
We note that the coarse space sizes for each method tend to be slightly larger with the more general decompositions used by METIS, but otherwise the same trends are seen. As such, we conclude that non-uniform decompositions can be well handled by the spectral coarse spaces employed here.

3.4. The Effect of Heterogeneity

We now turn our attention to the key property of robustness to heterogeneities. For this, we consider layered media within the wave guide. Three configurations, each having 10 layers, are used and are detailed in Figure 4. The heterogeneity is introduced in the wave speed c ( x ) , and in each case, c takes values from 1 to ρ , where ρ is a contrast parameter determining the strength of the heterogeneity. The wave number is then given by k = ω / c where ω is the angular frequency; we vary both ω and ρ in our tests. Note that for the DtN method, the eigenvalue threshold now depends on k s = max x Ω s k ( x ) , which may be different for different subdomains Ω s . To avoid notational clutter, we omit the subscript when referring to the method, namely retaining the name DtN( k 4 / 3 ).
Results for DtN( k 4 / 3 ) and H-GenEO( 1 2 ) for the increasing layers problem (Figure 4a) are provided in Table 5. Unfortunately, we see that the DtN( k 4 / 3 ) loses robustness with the heterogeneity present in this problem. In particular, we lose any robustness to the wave number k, and for the largest wave number used, we also see that changes in the contrast, given by ρ , can begin to have a sizeable impact on the iteration counts despite otherwise being relatively stable to changes in ρ . To a lesser extent, we also lose scalability with DtN( k 4 / 3 ) as the iteration counts now slowly increase with N. On the other hand, H-GenEO( 1 2 ) has strong robustness throughout, both with respect to the wave number k, the contrast in the heterogeneity ρ , and scalability as N increases. As such, we see a clear preference for H-GenEO( 1 2 ) as a stable and reliable method for heterogeneous problems.
We also consider the case of a diagonal layers problem (Figure 4c) in Table 6. Here, the issues with DtN( k 4 / 3 ) are reduced, but there is still some increase in iteration counts, especially for higher wave numbers k and a larger number of subdomains N. We note, in results not shown here, that DtN( k 3 / 2 ) also suffers from the same lack of robustness. For H-GenEO( 1 2 ), however, we still have good robustness to the parameters of the problem.
We further consider the heterogeneous problem when making use of METIS for more general non-uniform subdomain decompositions in order to ensure that H-GenEO is able to handle both difficulties together. For this, we consider the alternating layers problem (Figure 4b) and provide results in Table 7. We find that H-GenEO( 1 2 ) performs very well and continues to provide a rather robust method, even in the presence of heterogeneity on non-uniform subdomains. This further evidences the strength of the H-GenEO( 1 2 ) approach, and we briefly study this more closely, dropping reference to the eigenvalue tolerance and simply denoting the method as H-GenEO. First, however, we provide a brief examination into the use of higher order discretisation.

3.5. Higher Order Finite Elements

We briefly investigate the use of higher order finite elements. In particular, we consider the use of P2 elements as opposed to P1 elements. To give a direct comparison, we utilise the same meshes, and in Table 8 we give results for the heterogeneous diagonal layers problem with ρ = 10 ; equivalent results for P1 elements are given in Table 6. We observe that both iteration counts and coarse space sizes remain rather similar to the case of P1 elements. This suggests that, by itself, the order of the underlying finite elements used does not strongly affect performance. However, the typical meshes employed for higher order elements may be coarser, and further studies would be required to observe how this affects the utility of the spectral coarse spaces presented here. Further investigation into which method and order of finite element approximation provides the most efficient choice is beyond the scope of the present study.

3.6. The Effect of Boundary Conditions within the H-GenEO Eigenproblem

We now consider the choice of boundary conditions within H-GenEO in light of the fact that, for wave propagation problems, impedance conditions can often prove more practical within overlapping Schwarz methods. To this end, we consider the H-GenEO eigenproblem where the Neumann boundary condition is replaced by an impedance condition instead (i.e., the Robin condition in (7b)). Results for this impedance–H-GenEO method are given in Table 9 for the homogeneous problem with uniform square subdomains. We see that, while the use of this eigenproblem retains the good behaviour of H-GenEO as k increases, it lacks scalability as we increase the number of subdomains N. We note that the size of the coarse space is very similar to that when the standard Neumann condition is used (a direct comparison can be made for k = 73.8 with results in Table 3), and so this is not simply an artefact of a smaller coarse space. This shows that the Neumann condition within the eigenproblem is an important aspect of H-GenEO.

3.7. The Effect of More Overlap When Using H-GenEO

So far, all our results use minimal overlap. Here, we consider the case of increasing the overlap between subdomains—this is done by adding on layers of adjoining elements to each subdomain in a symmetric way so that minimal overlap is given by an overlap parameter of 2; that is, the overlapping region has a width of 2 elements. In Table 10, we report results for increasing overlap when using H-GenEO for the homogeneous problem with k = 46.5 and h 1 = 400 and where a uniform decomposition is used. We see that adding on a small amount of overlap can slightly decrease the iteration counts, but increasing the overlap width further can give much poorer results, especially when using a large number of subdomains N. One possible explanation could be an increase in the “colouring constant” (see, e.g., [45] (Definition 5.5)). We note that the size of the coarse space decreases somewhat as the overlap is increased; however, the extra computational effort required to deal with the larger subdomains will hamper any gains from this, along with the increased iteration counts. From these results, we determine that the H-GenEO coarse space is best suited to the case of minimal overlap, as we have used elsewhere throughout this work.

3.8. Weak Scalability and Timing Results for H-GenEO

In this section, we consider the weak scalability of the H-GenEO method. To approach this, we consider a growing wave guide domain where fixed size subdomains are added, each with the same number of degrees of freedom (dofs). This is done by repeatedly adding a unit square, which is split into 25 non-overlapping square subdomains, to the right of the existing domain L times to give Ω = ( 0 , L ) × ( 0 , 1 ) . Along the long edges of the global domain, we prescribe homogeneous Dirichlet boundary conditions, while Robin conditions are used at each end of the wave guide; a schematic for this weak scaling test is given in Figure 5. Heterogeneity is given by the alternating layers problem (see Figure 4b) with the 10 layers extending across the length of the wave guide.
To deal with the large problem sizes (reaching up to 10 , 253 , 601 dofs) and provide appropriate timing results, we assign one core per subdomain and solve using the ARCHIE-WeSt supercomputing facility on up to 400 cores (the machine uses Intel Xeon Gold 6138 processors at 2.0 GHz with 4.8 GB RAM per core).
In Table 11, we give results for a length of domain L = 2 up to L = 16 ; that is, from N = 50 to N = 400 subdomains, using k = 73.8 , h 1 = 800 and ρ = 100 . We observe that the iteration counts remain almost constant, increasing only very mildly with N. On the other hand, the coarse space size grows linearly with N, as expected given that subdomains have a fixed size, and we see that the time taken to solve the eigenproblems stays constant. However, due the increasing size of the coarse space, the run times slowly increase, meaning that the efficiency degrades as we solve increasingly large problems. Further, the setup of the decomposition and partitioning is performed sequentially, and so this setup time increases with N. Nonetheless, when removing this initial setup, the efficiency still slowly reduces as we increase the problem size, as can be seen in the final row of Table 11. To overcome this, the factorisation of the coarse space operator and associated solves must be made scalable. While a multi-level approach is therefore an attractive option, it is not yet clear how to formulate such a strategy for Helmholtz problems. Additionally, the coarse problem is only assigned to one process here, and so spreading it over more processes as N increases may stop the coarse solve from becoming a bottleneck.

3.9. High Wave Number Strong Scalability and Timing Results for H-GenEO

To conclude our numerical results, we consider the use of H-GenEO within a high wave number example and explore timings and the overall strong scalability of the approach. To this end, we consider the original homogeneous wave guide, as outlined in Figure 2, with a wave number of k = 186.0 and mesh width h 1 = 3200 , giving a total problem size of 10 , 246 , 401 dofs. To solve, we use METIS to give a non-uniform decomposition into subdomains and again utilise the ARCHIE-WeSt supercomputing facility to assign one core per subdomain.
Results are tabulated in Table 12, in which we detail the run time (in seconds) and the percentage of time spent by the eigensolver and coarse factorisation. The run times are further used to determine the parallel efficiency based on the smallest run on N = 80 cores. Timing results are also displayed graphically in Figure 6, where we show that the bulk of the computation is spent in the eigensolves and factorisation of the coarse problem; as such, while the METIS decomposition is sequential, we do not separate out the setup time here, given that it is comparatively rather small. From Table 12, we see that the iteration counts show good scalability, increasing only very mildly as we increase the number of subdomains fivefold. This is also seen in the run times, which decrease as we use more cores, and hence more subdomains, to solve the problem. In particular, we see that the parallel efficiency remains over 100%, showing the strong scalability of the approach, in part due to the fact that as N increases, we have to solve smaller eigenproblems, and so the percentage of time spent by the eigensolver drops significantly as we increase N. Overall, these provide promising results that the H-GenEO method can be effective for the solution of high wave number problems in 2D.

4. Conclusions

In this work, we have developed and explored a GenEO-type coarse space for additive Schwarz methods that is appropriate for the heterogeneous Helmholtz problem. We have conducted extensive numerical tests to show how this approach behaves on a 2D model test problem of a wave guide discretised using finite elements on a pollution-free mesh, comparing our method with the DtN coarse space. We find that only our H-GenEO approach is robust to heterogeneity and increasing wave number k and further provides good scaling such that iteration counts of right-preconditioned GMRES are only mildly dependent on the number of subdomains. This dependence is strongest for lower wave numbers on many subdomains and reduces as k grows. Furthermore, convergence does not deteriorate with increasing wave number, albeit at the cost of a coarse space that grows as k increases. This is achieved consistently for non-uniform partitioning into subdomains and in the presence of strong heterogeneity. These results show promise that H-GenEO can be used as an effective coarse space for challenging heterogeneous problems.
Finally, we discuss how these findings differ from that in the companion paper [40], where none of the approaches (including the spectral coarse spaces explored in detail here) are seen to be clearly favourable over a wide range of problem settings. A key difference is that well-resolved pollution-free meshes are used in the present work, while the more realistic benchmark problems examined in [40] use under-resolved meshes with a fixed number of points per wavelength, as is typical in engineering practice. As well as the more complex test cases considered, another contributing factor to the difference in studies is that in [40], a fixed number of eigenvectors was taken per subdomain (limiting the size of the coarse space), whereas here, eigenvalue thresholding was employed to provide a more robust approach, at the cost of a potentially larger coarse space. As such, the present study can be thought of as investigating the extent of what can be achieved in terms of a robust method in the ideal case, while [40] presents a viewpoint within the regime of more challenging practical problems. Identifying coarse spaces that can bridge this gap and provide a consistently robust approach even for the most demanding real-world applications remains an open area of research.

Author Contributions

Conceptualization, N.B. and V.D.; methodology, N.B.; software, N.B.; validation, N.B.; formal analysis, N.B. and V.D.; investigation, N.B.; resources, V.D.; writing—original draft preparation, N.B.; writing—review and editing, N.B. and V.D.; visualization, N.B.; supervision, V.D.; project administration, V.D.; funding acquisition, V.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by EPSRC grant number EP/S004017/1.

Acknowledgments

Numerical results were obtained using the ARCHIE-WeSt High Performance Computer (www.archie-west.ac.uk accessed on 16 February 2022) based at the University of Strathclyde.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moiola, A.; Spence, E.A. Is the Helmholtz equation really sign-indefinite? SIAM Rev. 2014, 56, 274–312. [Google Scholar] [CrossRef]
  2. Ernst, O.G.; Gander, M.J. Why it is difficult to solve Helmholtz problems with classical iterative methods. In Numerical Analysis of Multiscale Problems; Graham, I.G., Hou, T.Y., Lakkis, O., Scheichl, R., Eds.; Springer: Berlin, Germany, 2012; pp. 325–363. [Google Scholar]
  3. Gander, M.J.; Zhang, H. A class of iterative solvers for the Helmholtz equation: Factorizations, sweeping preconditioners, source transfer, single layer potentials, polarized traces, and optimized Schwarz methods. SIAM Rev. 2019, 60, 3–76. [Google Scholar] [CrossRef] [Green Version]
  4. Gillman, A.; Barnett, A.H.; Martinsson, P.-G. A spectrally accurate direct solution technique for frequency-domain scattering problems with variable media. BIT 2015, 55, 141–170. [Google Scholar] [CrossRef] [Green Version]
  5. Wang, S.; de Hoop, M.V.; Xia, J. On 3D modeling of seismic wave propagation via a structured parallel multifrontal direct Helmholtz solver. Geophys. Prospect. 2011, 59, 857–873. [Google Scholar] [CrossRef]
  6. Calandra, H.; Gratton, S.; Pinel, X.; Vasseur, X. An improved two-grid preconditioner for the solution of three-dimensional Helmholtz problems in heterogeneous media. Numer. Linear Algebra Appl. 2013, 20, 663–688. [Google Scholar] [CrossRef]
  7. Hu, Q.; Zhang, H. Substructuring preconditioners for the systems arising from plane wave discretization of Helmholtz equations. SIAM J. Sci. Comput. 2016, 38, A2232–A2261. [Google Scholar] [CrossRef]
  8. Erlangga, Y.; Vuik, C.; Oosterlee, C.W. On a class of preconditioners for solving the Helmholtz equation. Appl. Numer. Math. 2004, 50, 409–425. [Google Scholar] [CrossRef] [Green Version]
  9. Erlangga, Y.; Oosterlee, C.W.; Vuik, C. A novel multigrid based preconditioner for heterogeneous Helmholtz problems. SIAM J. Sci. Comput. 2006, 27, 1471–1492. [Google Scholar] [CrossRef]
  10. Cocquet, P.-H.; Gander, M.J. How large a shift is needed in the shifted Helmholtz preconditioner for its effective inversion by multigrid? SIAM J. Sci. Comput. 2017, 39, A438–A478. [Google Scholar] [CrossRef]
  11. Dwarka, V.; Vuik, C. Scalable convergence using two-level deflation preconditioning for the Helmholtz equation. SIAM J. Sci. Comput. 2020, 42, A901–A928. [Google Scholar] [CrossRef]
  12. Lahaye, D.; Vuik, C. How to choose the shift in the shifted Laplace preconditioner for the Helmholtz equation combined with deflation. In Modern Solvers for Helmholtz Problems; Lahaye, D., Tang, J., Vuik, K., Eds.; Birkhäuser: Cham, Switzerland, 2017; pp. 85–112. [Google Scholar]
  13. Engquist, B.; Ying, L. Sweeping preconditioner for the Helmholtz equation: Hierarchical matrix representation. Comm. Pure Appl. Math. 2011, 64, 697–735. [Google Scholar] [CrossRef] [Green Version]
  14. Engquist, B.; Ying, L. Sweeping preconditioner for the Helmholtz equation: Moving perfectly matched layers. Multiscale. Model. Simul. 2011, 9, 686–710. [Google Scholar] [CrossRef] [Green Version]
  15. Taus, M.; Zepeda-Núñez, L.; Hewett, R.J.; Demanet, L. L-Sweeps: A scalable, parallel preconditioner for the high-frequency Helmholtz equation. J. Comput. Phys. 2020, 420, 109706. [Google Scholar] [CrossRef]
  16. Dai, R.; Modave, A.; Remacle, J.-F.; Geuzaine, C. Multidirectional sweeping preconditioners with non-overlapping checkerboard domain decomposition for Helmholtz problems. J. Comput. Phys. 2022, 453, 110887. [Google Scholar] [CrossRef]
  17. Farhat, C.; Macedo, A.; Lesoinne, M. A two-level domain decomposition method for the iterative solution of high frequency exterior Helmholtz problems. Numer. Math. 2000, 85, 283–308. [Google Scholar] [CrossRef]
  18. Farhat, C.; Avery, P.; Tezaur, R.; Li, J. FETI-DPH: A dual-primal domain decomposition method for acoustic scattering. J. Comput. Acoust. 2005, 13, 499–524. [Google Scholar] [CrossRef]
  19. Després, B. Domain decomposition method for the Helmholtz problem. C R Math. Acad. Sci. Paris I Math. 1990, 311, 313–316. [Google Scholar]
  20. Claeys, X.; Parolin, E. Robust treatment of cross points in Optimized Schwarz Methods. arXiv 2020, arXiv:2003.06657. [Google Scholar]
  21. Gander, M.J.; Magoules, F.; Nataf, F. Optimized Schwarz methods without overlap for the Helmholtz equation. SIAM J. Sci. Comput. 2002, 24, 38–60. [Google Scholar] [CrossRef] [Green Version]
  22. Boubendir, Y.; Antoine, X.; Geuzaine, C. A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation. J. Comput. Phys. 2012, 231, 262–280. [Google Scholar] [CrossRef] [Green Version]
  23. Collino, F.; Ghanemi, S.; Joly, P. Domain decomposition method for harmonic wave propagation: A general presentation. Comput. Methods Appl. Mech. Engrg. 2000, 184, 171–211. [Google Scholar] [CrossRef] [Green Version]
  24. Cai, X.-C.; Casarin, M.A.; Elliott, F.W., Jr.; Widlund, O.B. Overlapping Schwarz algorithms for solving Helmholtz’s equation. Contemp. Math. 1998, 218, 391–399. [Google Scholar]
  25. Gander, M.J.; Zhang, H. Optimized Schwarz methods with overlap for the Helmholtz equation. SIAM J. Sci. Comput. 2016, 38, A3195–A3219. [Google Scholar] [CrossRef] [Green Version]
  26. Kimn, J.-H.; Sarkis, M. Restricted overlapping balancing domain decomposition methods and restricted coarse problems for the Helmholtz problem. Comput. Methods Appl. Mech. Engrg. 2007, 196, 1507–1514. [Google Scholar]
  27. Graham, I.G.; Spence, E.A.; Vainikko, E. Recent results on domain decomposition preconditioning for the high-frequency Helmholtz equation using absorption. In Modern Solvers for Helmholtz Problems; Lahaye, D., Tang, J., Vuik, K., Eds.; Birkhäuser: Cham, Switzerland, 2017; pp. 3–26. [Google Scholar]
  28. Kimn, J.-H.; Sarkis, M. Shifted Laplacian RAS solvers for the Helmholtz equation. In Domain Decomposition Methods in Science and Engineering XX; Bank, R., Holst, M., Widlund, O., Xu, J., Eds.; Springer: Berlin, Germany, 2013; pp. 151–158. [Google Scholar]
  29. Graham, I.G.; Spence, E.A.; Vainikko, E. Domain decomposition preconditioning for high-frequency Helmholtz problems with absorption. Math. Comp. 2017, 86, 2089–2127. [Google Scholar] [CrossRef] [Green Version]
  30. Graham, I.G.; Spence, E.A.; Zou, J. Domain Decomposition with local impedance conditions for the Helmholtz equation with absorption. SIAM J. Numer. Anal. 2020, 58, 2515–2543. [Google Scholar] [CrossRef]
  31. Gong, S.; Graham, I.G.; Spence, E.A. Domain decomposition preconditioners for high-order discretizations of the heterogeneous Helmholtz equation. IMA J. Numer. Anal. 2021, 41, 2139–2185. [Google Scholar] [CrossRef]
  32. Gong, S.; Gander, M.J.; Graham, I.G.; Lafontaine, D.; Spence, E.A. Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation. arXiv 2021, arXiv:2106.05218. [Google Scholar]
  33. Bootland, N.; Dolean, V.; Kyriakis, A.; Pestana, J. Analysis of parallel Schwarz algorithms for time-harmonic problems using block Toeplitz matrices. Electron. Trans. Numer. Anal. 2022, 55, 112–141. [Google Scholar] [CrossRef]
  34. Bonazzoli, M.; Dolean, V.; Graham, I.G.; Spence, E.A.; Tournier, P.-H. Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption. Math. Comp. 2019, 88, 2559–2604. [Google Scholar] [CrossRef] [Green Version]
  35. Spillane, N.; Dolean, V.; Hauret, P.; Nataf, F.; Pechstein, C.; Scheichl, R. Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps. Numer. Math. 2014, 126, 741–770. [Google Scholar] [CrossRef] [Green Version]
  36. Nataf, F.; Xiang, H.; Dolean, V.; Spillane, N. A coarse space construction based on local Dirichlet-to-Neumann maps. SIAM J. Sci. Comput. 2011, 33, 1623–1642. [Google Scholar] [CrossRef] [Green Version]
  37. Bootland, N.; Dolean, V. On the Dirichlet-to-Neumann coarse space for solving the Helmholtz problem using domain decomposition. In Numerical Mathematics and Advanced Applications ENUMATH 2019; Vermolen, F.J., Vuik, C., Eds.; Springer: Cham, Switzerland, 2021; pp. 175–184. [Google Scholar]
  38. Conen, L.; Dolean, V.; Krause, R.; Nataf, F. A coarse space for heterogeneous Helmholtz problems based on the Dirichlet-to-Neumann operator. J. Comput. Appl. Math. 2014, 271, 83–99. [Google Scholar] [CrossRef]
  39. Bootland, N.; Dolean, V.; Graham, I.G.; Ma, C.; Scheichl, R. Overlapping Schwarz methods with GenEO coarse spaces for indefinite and non-self-adjoint problems. arXiv 2021, arXiv:2110.13537. [Google Scholar]
  40. Bootland, N.; Dolean, V.; Jolivet, P.; Tournier, P.-H. A comparison of coarse spaces for Helmholtz problems in the high frequency regime. Comput. Math. Appl. 2021, 98, 239–253. [Google Scholar] [CrossRef]
  41. Zarmi, A.; Turkel, E. A general approach for high order absorbing boundary conditions for the Helmholtz equation. J. Comput. Phys. 2013, 242, 387–404. [Google Scholar] [CrossRef]
  42. Beriot, H.; Modave, A. An automatic perfectly matched layer for acoustic finite element simulations in convex domains of general shape. Int. J. Numer. Methods Engergy 2021, 122, 1239–1261. [Google Scholar] [CrossRef]
  43. Harari, I.; Slavutin, M.; Turkel, E. Analytical and numerical studies of a finite element PML for the Helmholtz equation. J. Comput. Acoust. 2000, 8, 121–137. [Google Scholar] [CrossRef]
  44. Babuska, I.M.; Sauter, S.A. Is the pollution effect of the FEM avoidable for the Helmholtz equation considering high wave numbers? SIAM J. Numer. Anal. 1997, 34, 2392–2423. [Google Scholar] [CrossRef] [Green Version]
  45. Dolean, V.; Jolivet, P.; Nataf, F. An Introduction to Domain Decomposition Methods: Algorithms, Theory, and Parallel Implementation; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2015. [Google Scholar]
  46. Fish, J.; Qu, Y. Global-basis two-level method for indefinite systems. Part 1: Convergence studies. Int. J. Numer. Methods Engergy 2000, 49, 439–460. [Google Scholar] [CrossRef]
  47. Nataf, F.; Xiang, H.; Dolean, V. A two level domain decomposition preconditioner based on local Dirichlet-to-Neumann maps. C R Math. Acad. Sci. Paris I Math. 2010, 348, 1163–1167. [Google Scholar] [CrossRef]
  48. Haferssas, R.; Jolivet, P.; Nataf, F. An additive Schwarz method type theory for Lions’s algorithm and a symmetrized optimized restricted additive Schwarz method. SIAM J. Sci. Comput. 2017, 39, A1345–A1365. [Google Scholar] [CrossRef]
  49. Nataf, F.; Tournier, P.-H. A GenEO domain decomposition method for saddle point problems. arXiv 2019, arXiv:1911.01858. [Google Scholar]
  50. Spillane, N. An abstract theory of domain decomposition methods with coarse spaces of the GenEO family. arXiv 2021, arXiv:2104.00280. [Google Scholar]
  51. Bootland, N.; Dolean, V.; Graham, I.G.; Ma, C.; Scheichl, R. GenEO coarse spaces for heterogeneous indefinite elliptic problems. In Domain Decomposition Methods in Science and Engineering XXVI; Brenner, S., Chung, E.T.S., Klawonn, A., Kwok, F., Xu, J., Zou, J., Eds.; Springer: Cham, Switzerland, 2017; accepted. [Google Scholar]
  52. Hecht, F. New development in FreeFem++. J. Numer. Math. 2012, 20, 251–266. [Google Scholar] [CrossRef]
  53. Lehoucq, R.B.; Sorensen, D.C.; Yang, C. ARPACK Users’ Guide: Solution of Large-Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 1998. [Google Scholar]
  54. Amestoy, P.R.; Duff, I.S.; L’Excellent, J.-Y.; Koster, J. A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM J. Matrix Anal. Appl. 2001, 23, 15–41. [Google Scholar] [CrossRef] [Green Version]
  55. Karypis, G.; Kumar, V. A fast and high quality multilevel scheme for partitioning irregular graphs. SIAM J. Sci. Comput. 1998, 20, 359–392. [Google Scholar] [CrossRef]
Figure 1. Local eigenfunctions for k = 46.5 . Top row: Examples using DtN (14). Middle row: Equivalent examples using H-GenEO (19). Bottom row: Examples using H-GenEO which are not found amongst the DtN eigenfunctions.
Figure 1. Local eigenfunctions for k = 46.5 . Top row: Examples using DtN (14). Middle row: Equivalent examples using H-GenEO (19). Bottom row: Examples using H-GenEO which are not found amongst the DtN eigenfunctions.
Mca 27 00035 g001
Figure 2. Schematic of the 2D wave guide model problem with example triangular mesh.
Figure 2. Schematic of the 2D wave guide model problem with example triangular mesh.
Mca 27 00035 g002
Figure 3. The size of coarse space utilised for the homogeneous problem when using ORAS with the DtN and H-GenEO coarse spaces. A uniform decomposition into N × N square subdomains is used. (a) Varying the wave number k for N = 25 , (b) Varying the number of subdomains N for k = 73.8 .
Figure 3. The size of coarse space utilised for the homogeneous problem when using ORAS with the DtN and H-GenEO coarse spaces. A uniform decomposition into N × N square subdomains is used. (a) Varying the wave number k for N = 25 , (b) Varying the number of subdomains N for k = 73.8 .
Mca 27 00035 g003
Figure 4. Piecewise constant layer profiles for the wave speed c ( x ) . For the darkest shade c ( x ) = 1 , while for the lightest shade c ( x ) = ρ , with ρ being the contrast factor.
Figure 4. Piecewise constant layer profiles for the wave speed c ( x ) . For the darkest shade c ( x ) = 1 , while for the lightest shade c ( x ) = ρ , with ρ being the contrast factor.
Mca 27 00035 g004
Figure 5. Schematic of the growing 2D wave guide model problem used in a weak scaling test on N = 25 L fixed size subdomains, with the underlying non-overlapping subdomains shown in grey.
Figure 5. Schematic of the growing 2D wave guide model problem used in a weak scaling test on N = 25 L fixed size subdomains, with the underlying non-overlapping subdomains shown in grey.
Mca 27 00035 g005
Figure 6. Timings for the homogeneous problem when using ORAS with H-GenEO ( 1 2 ) and a varying number of subdomains for k = 186.0 and h 1 = 3200 , giving a total of 10 , 246 , 401 dofs. A non-uniform decomposition into N subdomains is used, given by METIS.
Figure 6. Timings for the homogeneous problem when using ORAS with H-GenEO ( 1 2 ) and a varying number of subdomains for k = 186.0 and h 1 = 3200 , giving a total of 10 , 246 , 401 dofs. A non-uniform decomposition into N subdomains is used, given by METIS.
Mca 27 00035 g006
Table 1. Preconditioned GMRES iteration counts and size of coarse space (in parentheses) for the homogeneous problem when using ORAS and various coarse spaces. A uniform decomposition into 5 × 5 square subdomains is used, giving 25 subdomains in total.
Table 1. Preconditioned GMRES iteration counts and size of coarse space (in parentheses) for the homogeneous problem when using ORAS and various coarse spaces. A uniform decomposition into 5 × 5 square subdomains is used, giving 25 subdomains in total.
kh−1One-LevelDtNΔ-GenEOH-GenEO
18.51007319 (147)53 (135)21 (164)
29.32009726 (218)100 (271)18 (370)
46.540012535 (303)148 (560) 17 (779)
73.880015642 (502)220 (1120)15 (1712)
Table 2. Preconditioned GMRES iteration counts and size of coarse space (in parentheses) for the homogeneous problem when using ORAS and the DtN and H-GenEO coarse spaces with varying eigenvalue thresholds. A uniform decomposition into 5 × 5 square subdomains is used, giving 25 subdomains in total.
Table 2. Preconditioned GMRES iteration counts and size of coarse space (in parentheses) for the homogeneous problem when using ORAS and the DtN and H-GenEO coarse spaces with varying eigenvalue thresholds. A uniform decomposition into 5 × 5 square subdomains is used, giving 25 subdomains in total.
kh−1DtN (k)DtN (k4/3)DtN (k3/2)H-GenEO ( 1 8 )H-GenEO ( 1 4 )H-GenEO ( 1 2 )
18.510019 (147)13 (260)11 (403)46 (80)31 (105)21 (164)
29.320026 (218)14 (483)13 (759)53 (139)33 (189)18 (370)
46.540035 (303)14 (868)12 (1479)56 (245)35 (378)17 (779)
73.880042 (502)16 (1588)15 (2925)40 (546)25 (800)15 (1712)
Table 3. Preconditioned GMRES iteration counts (above), size of coarse space (middle), and average number of eigenvectors taken per subdomain (below) for the homogeneous problem when using ORAS with the DtN and H-GenEO coarse spaces and a varying number of subdomains N for k = 73.8 and h 1 = 800 . A uniform decomposition into N × N square subdomains is used.
Table 3. Preconditioned GMRES iteration counts (above), size of coarse space (middle), and average number of eigenvectors taken per subdomain (below) for the homogeneous problem when using ORAS with the DtN and H-GenEO coarse spaces and a varying number of subdomains N for k = 73.8 and h 1 = 800 . A uniform decomposition into N × N square subdomains is used.
N49162536496481100121144169196
DtN (k)283240425176499490363796154
DtN (k4/3)15161916161615161515161717
H-GenEO ( 1 8 )2631364071706512781116247194138
H-GenEO ( 1 2 )13151515161616181618181819
DtN (k)12425136250260573684310009461329155415291327
DtN (k4/3)39279011751588199423662753317636113976436949555188
H-GenEO ( 1 8 )2003054085465366007887339369277809741264
H-GenEO ( 1 2 )852111614281712190322612444262931203204348238823816
DtN (k)31.027.922.620.116.815.013.212.39.511.010.89.06.8
DtN (k4/3)98.087.873.463.555.448.343.039.236.132.930.329.326.1
H-GenEO ( 1 8 )50.033.925.521.814.912.212.39.09.47.75.45.86.4
H-GenEO ( 1 2 )213.0124.089.368.552.946.138.232.531.226.524.223.019.5
Table 4. Preconditioned GMRES iteration counts (above), size of coarse space (middle), and average number of eigenvectors taken per subdomain (below) for the homogeneous problem when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A non-uniform decomposition into N subdomains is used, given by METIS.
Table 4. Preconditioned GMRES iteration counts (above), size of coarse space (middle), and average number of eigenvectors taken per subdomain (below) for the homogeneous problem when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A non-uniform decomposition into N subdomains is used, given by METIS.
Number of Subdomains N
DtN (k4/3) H-GenEO ( 1 2 )
kh−1204080120160200204080120160200
18.5100101010101010151719222727
29.3200121511121212151719202223
46.5400121315131313151616182020
73.8800151514161416151617171719
117.21600141516171516141515161616
18.510028142265284310051157201285383471524589
29.3200477758113014101693192240057478395810971245
46.54009591466213226773151355386911931670200822532507
73.8800169525633751467254866199186324563433414747495338
117.216003049469568318486989611,092 423856807575904910,27311,305
18.510014.110.68.27.06.35.810.17.14.83.93.32.9
29.320023.918.914.111.810.69.620.014.39.88.06.96.2
46.540048.036.626.622.319.717.843.529.820.916.714.112.5
73.880084.864.146.938.934.331.093.261.442.934.629.726.7
117.21600152.4117.485.470.761.955.5211.9142.094.775.4164.256.5
Table 5. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the heterogeneous increasing layers problem when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Table 5. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the heterogeneous increasing layers problem when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Number of Subdomains N
DtN (k4/3) H-GenEO ( 1 2 )
ωh−1ρ163664100144196163664100144196
29.320010293741525558151619181819
1000444450585252151517181717
46.540010323841666573151616191818
1000636974847371141516181717
73.880010354342405869151716171817
10008993107111114109141515161716
29.320010116173234363399467224354452662679754
100084111136285329371222350446642679741
46.540010208317405600704812458706990123415231678
1000144176202421496554450693990121615121666
73.88001037955769311421217140493014252074258430603553
100025429432674878483891414092058257230593534
Table 6. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the heterogeneous diagonal layers problem when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Table 6. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the heterogeneous diagonal layers problem when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Number of Subdomains N
DtN (k4/3) H-GenEO ( 1 2 )
ωh−1ρ163664100144196163664100144196
29.320010131413142125161820182325
1000131414142225161820182325
46.540010151414162531161717262122
1000151415162534161718272222
73.880010141816152026161717171920
1000151816153239161717171920
29.320010336593866109013761390260376499689737828
1000336594866109013751390259375499687737826
46.540010621107515401910237026225437891095138415991825
1000621107515391907236826145417901093138115961824
73.880010116419472692359241454608114516362243282332333681
1000116319462693359241314569114116332239282232323671
Table 7. Preconditioned GMRES iteration counts for the heterogeneous alternating layers problem with ρ = 10 / 100 / 1000 when using ORAS with H-GenEO( 1 2 ) and a varying number of subdomains. A non-uniform decomposition into N subdomains is used, given by METIS.
Table 7. Preconditioned GMRES iteration counts for the heterogeneous alternating layers problem with ρ = 10 / 100 / 1000 when using ORAS with H-GenEO( 1 2 ) and a varying number of subdomains. A non-uniform decomposition into N subdomains is used, given by METIS.
Number of Subdomains N with Sub-Columns for ρ = 10/100/1000
ωh−1204080120160200
18.5100171717191919212121272727313131333333
29.3200161616171717191919202020212121232323
46.5400171818181818222323252626272828282929
73.8800161616171717181818181919192020232323
117.21600151515151616161616161616161616161616
Table 8. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for P2 finite element discretisation of the heterogeneous diagonal layers problem with ρ = 10 when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Table 8. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for P2 finite element discretisation of the heterogeneous diagonal layers problem with ρ = 10 when using ORAS with DtN( k 4 / 3 ) or H-GenEO( 1 2 ) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Number of Subdomains N
DtN (k4/3) H-GenEO ( 1 2 )
ωh−1ρ163664100144196163664100144196
18.510010131110101018151619182324
29.320010141212131925151718182325
46.540010151212152330151617202122
73.880010171614131825151616171820
18.510010151326510706898937125186231346398519
29.320010300608892108614441516260377507686733824
46.540010589110815722069249126385407941100140315941820
73.88001091919162862361444094748114416452239280533253695
Table 9. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the homogeneous problem when using ORAS with impedance-H-GenEO (the eigenproblem (19) is altered to have impedance as opposed to Neumann boundary conditions on the left-hand side) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Table 9. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the homogeneous problem when using ORAS with impedance-H-GenEO (the eigenproblem (19) is altered to have impedance as opposed to Neumann boundary conditions on the left-hand side) and a varying number of subdomains. A uniform decomposition into N × N square subdomains is used.
Number of Subdomains N
kh−149162536496481100121144169196
18.510017192327303642454358616167
29.320017192225343341383549606265
46.540015181922262526394347525154
73.880015191920252726353339434451
18.510068102140158204217236287334341404477504
29.3200148215296370392521504576720768725793908
46.54003604926287549179881236117614681550174018071930
73.8800848110614201696187722182432257429603180344338343732
Table 10. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the homogeneous problem when using ORAS with H-GenEO, varying the amount of overlap (in terms of element width, with 2 representing minimal overlap) and number of subdomains for k = 46.5 and h 1 = 400 . A uniform decomposition into N × N square subdomains is used.
Table 10. Preconditioned GMRES iteration counts (above) and size of coarse space (below) for the homogeneous problem when using ORAS with H-GenEO, varying the amount of overlap (in terms of element width, with 2 representing minimal overlap) and number of subdomains for k = 46.5 and h 1 = 400 . A uniform decomposition into N × N square subdomains is used.
Number of Subdomains N
Overlap49162536496481100121144169196
214151517161616202619192221
410111112121312171416172120
88101010131312202322263127
16162126263777617586109178157164
236849264477993810301248119514761558175818452016
43524726006998719471088112412961449168916971690
8336436538650799863102498111321395151114251512
1631641750061073371792094211081239108612121280
Table 11. Weak scaling results and timings for the alternating layers problem when using ORAS with H-GenEO ( 1 2 ) and a varying number of subdomains for k = 73.8 , h 1 = 800 and ρ = 100 . A uniform decomposition into N = 25 L subdomains is used, as depicted in Figure 5. Note that setup refers to the initial decomposition and partitioning, which is performed sequentially, while the local problems and eigensolves are carried out in parallel.
Table 11. Weak scaling results and timings for the alternating layers problem when using ORAS with H-GenEO ( 1 2 ) and a varying number of subdomains for k = 73.8 , h 1 = 800 and ρ = 100 . A uniform decomposition into N = 25 L subdomains is used, as depicted in Figure 5. Note that setup refers to the initial decomposition and partitioning, which is performed sequentially, while the local problems and eigensolves are carried out in parallel.
N50100150200250300350400
Iteration count1718181919202121
Coarse space size30106150929012,43015,57018,71021,85024,990
Total run time (s)45.848.653.058.763.570.079.788.1
Weak scaling efficiency94.2%86.4%78.0%72.1%65.4%57.5%52.0%
Eigensolve time (s)37.137.937.938.337.837.937.937.7
Setup time (s)5.57.712.916.519.923.827.430.8
Efficiency without setup-98.5%100.5%95.5%92.4%87.2%77.1%70.3%
Table 12. Strong scaling results and timings for the homogeneous problem when using ORAS with H-GenEO ( 1 2 ) and a varying number of subdomains for k = 186.0 and h 1 = 3200 , giving a total of 10 , 246 , 401 dofs. A non-uniform decomposition into N subdomains is used, given by METIS. The average local eigenproblem size is given approximately as the number of dofs divided by N.
Table 12. Strong scaling results and timings for the homogeneous problem when using ORAS with H-GenEO ( 1 2 ) and a varying number of subdomains for k = 186.0 and h 1 = 3200 , giving a total of 10 , 246 , 401 dofs. A non-uniform decomposition into N subdomains is used, given by METIS. The average local eigenproblem size is given approximately as the number of dofs divided by N.
N80120160200240280320360400
Iteration count141615161717181919
Coarse space size16,01419,01821,34823,74725,56027,27028,79330,35731,773
Total run time (s)1214.4614.6404.4279.3217.3195.0159.4154.0147.6
Parallel efficiency132%150%174%186%178%190%175%165%
Eigenproblem size 128,08085,38764,04051,23242,69336,59432,02028,46225,616
Average no. of eigenvectors200.2158.5133.4118.7106.597.490.084.379.4
Eigensolve time68.6%70.1%66.8%62.7%53.7%52.3%41.6%39.3%37.3%
Coarse factorisation time29.6%27.8%30.8%34.4%42.8%44.3%54.0%56.3%58.3%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bootland, N.; Dolean, V. Can DtN and GenEO Coarse Spaces Be Sufficiently Robust for Heterogeneous Helmholtz Problems? Math. Comput. Appl. 2022, 27, 35. https://doi.org/10.3390/mca27030035

AMA Style

Bootland N, Dolean V. Can DtN and GenEO Coarse Spaces Be Sufficiently Robust for Heterogeneous Helmholtz Problems? Mathematical and Computational Applications. 2022; 27(3):35. https://doi.org/10.3390/mca27030035

Chicago/Turabian Style

Bootland, Niall, and Victorita Dolean. 2022. "Can DtN and GenEO Coarse Spaces Be Sufficiently Robust for Heterogeneous Helmholtz Problems?" Mathematical and Computational Applications 27, no. 3: 35. https://doi.org/10.3390/mca27030035

APA Style

Bootland, N., & Dolean, V. (2022). Can DtN and GenEO Coarse Spaces Be Sufficiently Robust for Heterogeneous Helmholtz Problems? Mathematical and Computational Applications, 27(3), 35. https://doi.org/10.3390/mca27030035

Article Metrics

Back to TopTop