Variable order, directional H2-matrices for Helmholtz problems with complex frequency

Abstract

The sparse approximation of high-frequency Helmholtz-type integral operators has many important physical applications such as problems in wave propagation and wave scattering. The discrete system matrices are huge and densely populated; hence, their sparse approximation is of outstanding importance. In our paper, we will generalize the directional-matrix techniques from the ‘pure’ Helmholtz operator with a purely imaginay frequency to general complex frequencies with non negative real part. In this case, the fundamental solution decreases exponentially for large arguments. We will develop a new admissibility condition that contains the real part of the frequency in an explicit way, and introduces the approximation of the integral kernel function on admissible blocks in terms of frequency-dependent directional expansion functions. We develop an error analysis that is explicit with respect to the expansion order and with respect to both the real and the imaginary part of the frequency. This allows for choosing the variable expansion order in a quasi-optimal way. The complexity analysis shows how higher values of the real part of the frequency reduce the complexity. In certain cases, it even turns out that the discrete matrix can be replaced by its nearfield part. Numerical experiments illustrate the sharpness of the derived estimates and the efficiency of our sparse approximation.