Analyzing the wave number dependency of the convergence rate of a multigrid preconditioned Krylov method for the Helmholtz equation with an absorbing layer

This paper analyzes the Krylov convergence rate of a Helmholtz problem preconditioned with Multigrid. The multigrid method is applied to the Helmholtz problem formulated on a complex contour and uses GMRES as a smoother substitute at each level. A one-dimensional model is analyzed both in a continuous and discrete way. It is shown that the Krylov convergence rate of the continuous problem is independent of the wave number. The discrete problem, however, can deviate significantly from this bound due to a pitchfork in the spectrum. It is further shown in numerical experiments that the convergence rate of the Krylov method approaches the continuous bound as the grid distance $h$ gets small

Local Fourier Analysis of the Complex Shifted Laplacian preconditioner for Helmholtz problems

In this paper we solve the Helmholtz equation with multigrid preconditioned Krylov subspace methods. The class of Shifted Laplacian preconditioners are known to significantly speed-up Krylov convergence. However, these preconditioners have a parameter beta, a measure of the complex shift. Due to contradictory requirements for the multigrid and Krylov convergence, the choice of this shift parameter can be a bottleneck in applying the method. In this paper, we propose a wavenumber-dependent minimal complex shift parameter which is predicted by a rigorous k-grid Local Fourier Analysis (LFA) of the multigrid scheme. We claim that, given any (regionally constant) wavenumber, this minimal complex shift parameter provides the reader with a parameter choice that leads to efficient Krylov convergence. Numerical experiments in one and two spatial dimensions validate the theoretical results. It appears that the proposed complex shift is both the minimal requirement for a multigrid V-cycle to converge, as well as being near-optimal in terms of Krylov iteration count.Comment: 20 page

Fourier Method for Approximating Eigenvalues of Indefinite Stekloff Operator

We introduce an efficient method for computing the Stekloff eigenvalues associated with the Helmholtz equation. In general, this eigenvalue problem requires solving the Helmholtz equation with Dirichlet and/or Neumann boundary condition repeatedly. We propose solving the related constant coefficient Helmholtz equation with Fast Fourier Transform (FFT) based on carefully designed extensions and restrictions of the equation. The proposed Fourier method, combined with proper eigensolver, results in an efficient and clear approach for computing the Stekloff eigenvalues.Comment: 12 pages, 4 figure

Matrix Reordering Methods for Table and Network Visualization

International audienceThis survey provides a description of algorithms to reorder visual matrices of tabular data and adjacency matrix of networks. The goal of this survey is to provide a comprehensive list of reordering algorithms published in different fields such as statistics, bioinformatics, or graph theory. While several of these algorithms are described in publications and others are available in software libraries and programs, there is little awareness of what is done across all fields. Our survey aims at describing these reordering algorithms in a unified manner to enable a wide audience to understand their differences and subtleties. We organize this corpus in a consistent manner, independently of the application or research field. We also provide practical guidance on how to select appropriate algorithms depending on the structure and size of the matrix to reorder, and point to implementations when available