3.3NADec 6, 2012
A Gaussian quadrature rule for oscillatory integrals on a bounded intervalAndreas Asheim, Alfredo Deaño, Daan Huybrechs et al.
We investigate a Gaussian quadrature rule and the corresponding orthogonal polynomials for the oscillatory weight function $e^{iωx}$ on the interval $[-1,1]$. We show that such a rule attains high asymptotic order, in the sense that the quadrature error quickly decreases as a function of the frequency $ω$. However, accuracy is maintained for all values of $ω$ and in particular the rule elegantly reduces to the classical Gauss-Legendre rule as $ω\to 0$. The construction of such rules is briefly discussed, and though not all orthogonal polynomials exist, it is demonstrated numerically that rules with an even number of points are always well defined. We show that these rules are optimal both in terms of asymptotic order as well as in terms of polynomial order.
1.2NAJun 21, 2015
Asymptotic solvers for ordinary differential equations with multiple frequenciesMarissa Condon, Alfredo Deano, Jing Gao et al.
We construct asymptotic expansions for ordinary differential equations with highly oscillatory forcing terms, focussing on the case of multiple, non-commensurate frequencies. We derive an asymptotic expansion in inverse powers of the oscillatory parameter and use its truncation as an exceedingly effective means to discretize the differential equation in question. Numerical examples illustrate the effectiveness of the method.
1.2CAJan 13, 2010
Asymptotic zero distribution of complex orthogonal polynomials associated with Gaussian quadratureA. Deano, D. Huybrechs, A. B. J. Kuijlaars
In this paper we study the asymptotic behavior of a family of polynomials which are orthogonal with respect to an exponential weight on certain contours of the complex plane. The zeros of these polynomials are the nodes for complex Gaussian quadrature of an oscillatory integral on the real axis with a high order stationary point, and their limit distribution is also analyzed. We show that the zeros accumulate along a contour in the complex plane that has the S-property in an external field. In addition, the strong asymptotics of the orthogonal polynomials is obtained by applying the nonlinear Deift--Zhou steepest descent method to the corresponding Riemann--Hilbert problem.