luisgarciar / faber-helmholtz

Implementation of a polynomial preconditioner for the Helmholtz equation based on Faber series

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Polynomial preconditioners for the Helmholtz equation based on Faber polynomials

This repository contains an implementation of the polynomial preconditioner for the solution of Helmholtz linear systems resulting from Helmholtz PDE problems described in the paper [GRSN21. It also contains standalone implementations of the multigrid method and the shifted Laplace preconditioner for Helmholtz problems. In addition to solving the problems presented in the paper and the dissertation, the code in the repository can be adapted to solve other Helmholtz problems in 1D and 2D.

Installation

The software has been developed in MATLAB 2017b and tested with versions up to MATLAB 2020a. To clone this repository, navigate using the terminal to your desired location and type git clone https://github.com/luisgarciar/faber-helmholtz.git

Next, add the path to faber-helmholtz to your MATLAB path, which can be done using the graphical interface in MATLAB (File -> Set Path -> Add with Subfolders).

Usage

The folder faber contains the implementations of the
Faber preconditioner and the folder helmholtz the code necessary for constructing the matrices and a multigrid solver.

The folder new_experiments contains a variety of numerical experiments for Helmholtz problems in 1D and 2D using the GMRES method and the Faber preconditioner, including those in the paper [GRSN21]. These can be adapted to solve more general problems

References

[GRSN21] L. Garc’a Ramos, O. Sète, R. Nabben, Preconditioning the Helmholtz equation with the shifted Laplacian and Faber polynomials, Electronic Transactions on Numerical Analysis 54, pp. 534-557, 2021.

About

Implementation of a polynomial preconditioner for the Helmholtz equation based on Faber series


Languages

Language:MATLAB 62.1%Language:HTML 28.6%Language:Jupyter Notebook 8.4%Language:TeX 0.7%Language:Python 0.1%Language:Roff 0.0%Language:Rich Text Format 0.0%Language:M 0.0%