Complex Shifted Laplacian for Large Helmholtz Problems
Large high-frequency free-space Helmholtz equations in, for example, acoustic and RF models can be difficult to solve with the usual iterative solvers such as multigrid and domain decomposition A method to solve such problems is to solve them with a complex wave number as a preconditioner for a Helmholtz equation without damping, thereby modifying the equation for the preconditioner in such a way that the convergence rate is improved. This strategy is called the complex shifted Laplacian (CSL) (even though it is not the Laplacian that is shifted), and it uses the following method:
Consider Helmholtz equation where the linear equation comes from the discretization of the equation
The preconditioner matrix for the CSL preconditioner Mβ is given by the discretization of the equation
where you choose a value of β such that you get evanescent wave solutions. Instead of solving the original system of equation, you now solve
when using left preconditioning (or the corresponding equation for right preconditioning).
You can activate complex shifted Laplacian for the following iterative solvers: Multigrid (GMG, and on multigrid levels for AMG) and Domain Decomposition (Schwarz). For details, see the documentation for those solvers.