Comparison of Two-Level Preconditioners Derived from Deflation, Domain Decomposition and Multigrid Methods

For various applications, it is well-known that a multi-level, in particular two-level, preconditioned CG (PCG) method is an efficient method for solving large and sparse linear systems with a coefficient matrix that is symmetric positive definite. The corresponding two-level preconditioner combines traditional and projection-type...

Theoretical and numerical comparison of various projection methods derived from deflation, domain decomposition and multigrid methods

Convergence bounds for preconditioned GMRES using element-by-element estimates of the field of values

By combining element-by-element estimates for the field of values of a preconditioned matrix with GMRES-convergence estimates it is possible to derive an easily computable upper bound on the GMRES-residual norm. This method can be applied to general finite element systems, but the preconditioner has to be Hermitian and positive definite. The...

A preconditioner for the Helmholtz equation with perfectly matched layer

This paper discusses an iterative method for solving the Helmholtz equation with the perfectly matched layer (PML). The method consists of an outer and inner iteration process. The inner iteration is used to approximately solve a preconditioner, which in this case is based on a modified PML equation. The outer iteration is a Krylov subspace...

Convergence bounds for preconditioned GMRES using element-by-element estimates of the field of values

By combining element-by-element estimates for the field of values of a preconditioned matrix with GMRES-convergence estimates it is possible to derive an easily computable upper bound on the GMRES-residual norm. This method can be applied to general finite element systems, but the preconditioner has to be Hermitian and positive definite. The...

Preconditioning a Finite Element Solver of the Exterior Helmholtz Equation

We consider acoustic scattering about a general body. This is described by the Helmholtz equation exterior to the body. In order to truncate the infinite domain we use the BGT absorbing boundary condition. The resultant problem in a finite domain is solved by a finite element procedure. This yields a large sparse system of linear equations which...

A new iterative solver for the time-harmonic wave equation

The time-harmonic wave equation, also known as the Helmholtz equation, is obtained if the constant-density acoustic wave equation is transformed from the time domain to the frequency domain. Its discretization results in a large, sparse, linear system of equations. In two dimensions, this system can be solved efficiently by a direct method. In...

Spectral analysis of the discrete Helmholtz operator preconditioned with a shifted Laplacian

Shifted Laplace preconditioners have attracted considerable attention as a technique to speed up convergence of iterative solution methods for the Helmholtz equation. In this paper we present a comprehensive spectral analysis of the Helmholtz operator preconditioned with a shifted Laplacian. Our analysis is valid under general conditions. The...

A robust and efficient iterative method for the numerical solution of the Helmholtz equation

A novel multigrid based preconditioner for heterogeneous Helmholtz problems

An iterative solution method, in the form of a preconditioner for a Krylov subspace method, is presented for the Helmholtz equation. The preconditioner is based on a Helmholtz type differential operator with a complex term. A multigrid iteration is used for approximately inverting the preconditioner. The choice of multigrid components for the...

Shifted Laplace preconditioners for the Helmholtz equations

On a class of preconditioners for solving the Helmholtz equation

Some numerical aspects for solving sparse large linear systems derived from the Helmholtz equation

In this report, several numerical aspects and diculties for solving a linear system derived from the time-harmonic wave equations are overviewed. The presentation begins with the derivation of the governing equation for waves propagating in general inhomogeneous media. Due to the need of numerical solutions, various discretizations based on...