Skip to main content

Fast‚ numerically stable computation of oscillatory integrals with stationary points

Sheehan Olver

Abstract

We present a numerically stable way to compute oscillatory integrals of the form ∫1-1 ƒ(x)eg(x) dx. For each additional frequency, only a small, well-conditioned linear system with a Hessenberg matrix must be solved, and the amount of work needed decreases as the frequency increases. Moreover, we can modify the method for computing oscillatory integrals with stationary points. This is the first stable algorithm for oscillatory integrals with stationary points which does not lose accuracy as the frequency increases and does not require deformation into the complex plane.

Institution
Oxford University Computing Laboratory
Month
May
Number
NA−09/02
Pages
24
Series
NA Group technical reports
Year
2009
[an error occurred while processing this directive]