An adaptive delaminating Levin method in two dimensions

Abstract: We present an adaptive delaminating Levin method for evaluating bivariate oscillatory integrals over rectangular domains. Whereas previous analyses of Levin methods impose non-resonance conditions that exclude stationary and resonance points, we rigorously establish the existence of a slowly-varying, approximate solution to the Levin PDE across all frequency regimes, even when the non-resonance condition is violated. This allows us to derive error estimates for the numerical solution of the Levin PDE via the Chebyshev spectral collocation method, and for the evaluation of the corresponding oscillatory integrals, showing that high accuracy can be achieved regardless of whether or not stationary and resonance points are present. We then present a Levin method incorporating adaptive subdivision in both two and one dimensions, as well as delaminating Chebyshev spectral collocation, which is effective in both the presence and absence of stationary and resonance points. We demonstrate the effectiveness of our algorithm with a number of numerical experiments.