Applying GMRES to the Helmholtz equation with strong trapping: how does the number of iterations depend on the frequency?

by   Pierre Marchand, et al.

We consider GMRES applied to discretisations of the high-frequency Helmholtz equation with strong trapping; recall that in this situation the problem is exponentially ill-conditioned through an increasing sequence of frequencies. Under certain assumptions about the distribution of the eigenvalues, we prove upper bounds on how the number of GMRES iterations grows with the frequency. Our main focus is on boundary-integral-equation formulations of the exterior Dirichlet and Neumann obstacle problems in 2- and 3-d; for these problems, we investigate numerically the sharpness (in terms of dependence on frequency) of both our bounds and various quantities entering our bounds. This paper is therefore the first comprehensive study of the frequency-dependence of the number of GMRES iterations for Helmholtz boundary-integral equations under trapping.



There are no comments yet.


page 1

page 2

page 3

page 4


High-frequency estimates on boundary integral operators for the Helmholtz exterior Neumann problem

We study a commonly-used second-kind boundary-integral equation for solv...

Local absorbing boundary conditions on fixed domains give order-one errors for high-frequency waves

We consider approximating the solution of the Helmholtz exterior Dirichl...

Eigenvalues of the truncated Helmholtz solution operator under strong trapping

For the Helmholtz equation posed in the exterior of a Dirichlet obstacle...

Separability of the kernel function in an integral formulation for anisotropic radiative transfer equation

We study in this work an integral formulation for the radiative transfer...

Spectral Galerkin boundary element methods for high-frequency sound-hard scattering problems

This paper is concerned with the design of two different classes of Gale...

A Fully-Accelerated Surface Integral Equation Method for the Electromagnetic Modeling of Arbitrary Objects

Surface integral equation (SIE) methods are of great interest for the nu...

Solutions of New Potential Integral Equations Using MLFMA Based on the Approximate Stable Diagonalization

We present efficient solutions of recently developed potential integral ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.