Fast Computation of Fourier Integral Operators
Emmanuel Candes, Laurent Demanet, Lexing Ying
Abstract
We introduce a general purpose algorithm for rapidly computing certain types of oscillatory integrals which frequently arise in problems connected to wave propagation and general hyperbolic equations. The problem is to evaluate numerically a so-called Fourier integral operator (FIO) of the form ∫ e2πi Φ(x,ξ) a(x,ξ) f(ξ) dξ at points given on a Cartesian grid. Here, ξ is a frequency variable, f(ξ) is the Fourier transform of the input f, a(x,ξ) is an amplitude and Φ(x,ξ) is a phase function, which is typically as large as |ξ|; hence the integral is highly oscillatory at high frequencies. Because an FIO is a dense matrix, a naive matrix vector product with an input given on a Cartesian grid of size N by N would require O(N4) operations. This paper develops a new numerical algorithm which requires O(N2.5 N) operations, and as low as O(N) in storage space. It operates by localizing the integral over polar wedges with small angular aperture in the frequency plane. On each wedge, the algorithm factorizes the kernel e2 πi Φ(x,ξ) a(x,ξ) into two components: 1) a diffeomorphism which is handled by means of a nonuniform FFT and 2) a residual factor which is handled by numerical separation of the spatial and frequency variables. The key to the complexity and accuracy estimates is that the separation rank of the residual kernel is provably independent of the problem size. Several numerical examples demonstrate the efficiency and accuracy of the proposed methodology. We also discuss the potential of our ideas for various applications such as reflection seismology.
Create a lesson
Related papers
Solving the Incompressible Navier-Stokes Equations on Oriented Curved Surfaces Discretized by Point Clouds
Alejandra Foggia, Ivo F. Sbalzarini
A coercive space-time variational approach to fractional diffusion problems
Herbert Egger, Marvin Fritz, Barbara Wohlmuth
Model reduction of port-Hamiltonian systems via neural networks
Silke Glas, Alexander Heinlein, Harald Monsuur et al.
Integral Invariants of Vasudeva Murthy's Relaxation Systems: Analysis and Numerical Validation
Sudipta Sahu, Rathan Samala
Efficient primal--dual splitting methods for a Poisson-constrained JKO scheme for Poisson-Nernst-Planck models
Wei Wu, Jin Zeng, Zhen Zhang et al.
Magnetic Field Conforming Multiscale Formulations for Locally-Confined Nonlinear Eddy Current Problems Using the FE-HMM Method
Innocent Niyonzima, Gérard Meunier, Antoine Marteau et al.