Separation-of-variables as a preconditioner for an iterative Helmholtz solver

A preconditioned iterative method based on separation-of-variables for solving the Helmholtz equation in an inhomogeneous medium is tested. The preconditioner is constructed by approximating the wavenumber by a sum of two terms, one depending only on one spatial coordinate, say x, and the other depending on the remaining coordinates. The Helmholtz equation can be solved efficiently if the wavenumber has such a separable form. First, an eigenvalue-eigenvector decomposition is applied in the x-direction. With these eigenvectors, a change of variables is performed in order to obtain a set of independent systems with one dimension less than the original one. For smooth models and low frequencies, the convergence rate with this preconditioner is satisfactory. Unfortunately, it rapidly deteriorates when the roughness of the model or the frequency increases. Examples from seismic modeling are given to illustrate this behaviour. Moreover, numerical evidence is presented that suggests that the decomposition of the wavenumber in the sum of two terms cannot be improved with this approach.