Results 11  20
of
161
Accelerating the nonuniform Fast Fourier Transform
 SIAM REVIEW
, 2004
"... The nonequispaced Fourier transform arises in a variety of application areas, from medical imaging to radio astronomy to the numerical solution of partial differential equations. In a typical problem, one is given an irregular sampling of N data in the frequency domain and one is interested in recon ..."
Abstract

Cited by 39 (2 self)
 Add to MetaCart
(Show Context)
The nonequispaced Fourier transform arises in a variety of application areas, from medical imaging to radio astronomy to the numerical solution of partial differential equations. In a typical problem, one is given an irregular sampling of N data in the frequency domain and one is interested in reconstructing the corresponding function in the physical domain. When the sampling is uniform, the fast Fourier transform (FFT) allows this calculation to be computed in O(N log N) operations rather than O(N 2) operations. Unfortunately, when the sampling is nonuniform, the FFT does not apply. Over the last few years, a number of algorithms have been developed to overcome this limitation and are often referred to as nonuniform FFTs (NUFFTs). These rely on a mixture of interpolation and the judicious use of the FFT on an oversampled grid [A. Dutt and V. Rokhlin, SIAM J. Sci. Comput., 14 (1993), pp. 1368–1383]. In this paper, we observe that one of the standard interpolation or “gridding ” schemes, based on Gaussians, can be accelerated by a significant factor without precomputation and storage of the interpolation weights. This is of particular value in two and threedimensional settings, saving either 10dN in storage in d dimensions or a factor of about 5–10 in CPUtime (independent of dimension).
Random sampling of multivariate trigonometric polynomials
 SIAM J. Math. Anal
, 2004
"... We investigate when a trigonometric polynomial p of degree M in d variables is uniquely determined by its sampled values p(xj) on a random set of points xj in the unit cube (the “sampling problem for trigonometric polynomials”) and estimate the probability distribution of the condition number for th ..."
Abstract

Cited by 31 (3 self)
 Add to MetaCart
We investigate when a trigonometric polynomial p of degree M in d variables is uniquely determined by its sampled values p(xj) on a random set of points xj in the unit cube (the “sampling problem for trigonometric polynomials”) and estimate the probability distribution of the condition number for the associated Vandermondetype and Toeplitzlike matrices. The results provide a solid theoretical foundation for some efficient numerical algorithms that are already in use.
Some applications of generalized FFTs
 In Proceedings of DIMACS Workshop in Groups and Computation
, 1997
"... . Generalized FFTs are efficient algorithms for computing a Fourier transform of a function defined on finite group, or a bandlimited function defined on a compact group. The development of such algorithms has been accompanied and motivated by a growing number of both potential and realized applicat ..."
Abstract

Cited by 31 (5 self)
 Add to MetaCart
(Show Context)
. Generalized FFTs are efficient algorithms for computing a Fourier transform of a function defined on finite group, or a bandlimited function defined on a compact group. The development of such algorithms has been accompanied and motivated by a growing number of both potential and realized applications. This paper will attempt to survey some of these applications. Appendices include some more detailed examples. 1. A brief history The now "classical" Fast Fourier Transform (FFT) has a long and interesting history. Originally discovered by Gauss, and later made famous after being rediscovered by Cooley and Tukey [21], it may be viewed as an algorithm which efficiently computes the discrete Fourier transform or DFT. In between Gauss and CooleyTukey others developed special cases of the algorithm, usually motivated by the need to make efficient data analysis of one sort or another. To cite but a few examples, Gauss was interested in efficiently interpolating the orbits of asteroids [43...
Efficient Algorithms for DiffusionGenerated Motion by Mean Curvature
 J. Comput. Phys
, 1996
"... We accept this thesis as conforming to the required standard ..."
Abstract

Cited by 23 (5 self)
 Add to MetaCart
(Show Context)
We accept this thesis as conforming to the required standard
A framework for discrete integral transformations II – the 2D 31 Radon transform
"... This paper is dedicated to the memory of Professor Moshe Israeli 19402007, who passed away on February 18. Computing the Fourier transform of a function in polar coordinates is an important building block in many scientific disciplines and numerical schemes. In this paper we present the pseudopola ..."
Abstract

Cited by 21 (10 self)
 Add to MetaCart
(Show Context)
This paper is dedicated to the memory of Professor Moshe Israeli 19402007, who passed away on February 18. Computing the Fourier transform of a function in polar coordinates is an important building block in many scientific disciplines and numerical schemes. In this paper we present the pseudopolar Fourier transform that samples the Fourier transform on the pseudopolar grid, also known as the concentric squares grid. The pseudopolar grid consists of equally spaced samples along rays, where different rays are equally spaced and not equally angled. The pseudopolar Fourier transform Fourier transform is shown to be fast (the same complexity as the FFT), stable, invertible, requires only
Fast and accurate Polar Fourier transform
 Appl. Comput. Harmon. Anal.
, 2006
"... In a wide range of applied problems of 2D and 3D imaging a continuous formulation of the problem places great emphasis on obtaining and manipulating the Fourier transform in Polar coordinates. However, the translation of continuum ideas into practical work with data sampled on a Cartesian grid is pr ..."
Abstract

Cited by 19 (1 self)
 Add to MetaCart
(Show Context)
In a wide range of applied problems of 2D and 3D imaging a continuous formulation of the problem places great emphasis on obtaining and manipulating the Fourier transform in Polar coordinates. However, the translation of continuum ideas into practical work with data sampled on a Cartesian grid is problematic. In this article we develop a fast high accuracy Polar FFT. For a given twodimensional signal of size N × N, the proposed algorithm’s complexity is O(N^2 log N), just like in a Cartesian 2DFFT. A special feature of our approach is that it involves only 1D equispaced FFT’s and 1D interpolations. A central tool in our method is the pseudoPolar FFT, an FFT where the evaluation frequencies lie in an oversampled set of nonangularly equispaced points. We describe the concept of pseudoPolar domain, including fast forward and inverse transforms. For those interested primarily in Polar FFT’s, the pseudoPolar FFT plays the role of a halfway point—a nearlyPolar system from which conversion to Polar coordinates uses processes relying purely on 1D FFT’s and interpolation operations. We describe the conversion process, and give an error analysis of it. We compare accuracy results obtained by a Cartesianbased unequallysampled FFT method to ours, both algorithms using a smallsupport interpolation and no precompensating, and show marked advantage to the use of the pseudoPolar initial grid.
FAST COMPUTATION OF FOURIER INTEGRAL OPERATORS
, 2007
"... We introduce a general purpose algorithm for rapidly computing certain types of oscillatory integrals which frequently arise in problems connected to wave propagation, general hyperbolic equations, and curvilinear tomography. The problem is to numerically evaluate a socalled Fourier integral operat ..."
Abstract

Cited by 18 (6 self)
 Add to MetaCart
We introduce a general purpose algorithm for rapidly computing certain types of oscillatory integrals which frequently arise in problems connected to wave propagation, general hyperbolic equations, and curvilinear tomography. The problem is to numerically evaluate a socalled 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, ξ) isan amplitude, and Φ(x, ξ) is a phase function, which is typically as large as ξ; hence the integral is highly oscillatory. Because a 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(N 4) operations. This paper develops a new numerical algorithm which requires O(N 2.5 log N) operations and as low as O ( √ N) in storage space (the constants in front of these estimates are small). 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 the fact that the separation rank of the residual kernel is provably independent of the problem size. Several numerical examples demonstrate the numerical accuracy and low computational complexity of the proposed methodology. We also discuss the potential of our ideas for various applications such as reflection seismology.
A Fast Spherical Filter with Uniform Resolution
, 1997
"... This paper introduces a fast algorithm for obtaining a uniform resolution representation of a function known at a latitudelongitude grid on the surface of a sphere, equivalent to a triangular, isotropic truncation of the spherical harmonic coe#cients for the function. The proposed spectral truncati ..."
Abstract

Cited by 18 (0 self)
 Add to MetaCart
This paper introduces a fast algorithm for obtaining a uniform resolution representation of a function known at a latitudelongitude grid on the surface of a sphere, equivalent to a triangular, isotropic truncation of the spherical harmonic coe#cients for the function. The proposed spectral truncation method, which is based on the fast multipole method and the fast Fourier transform, projects the function to a space with uniform resolution while avoiding surface harmonic transformations. The method requires O#N 2 log N# operations for O#N 2 # grid points, as opposed to O#N 3 # operations for the standard spectral transform method, providing a reducedcomplexity spectral method obviating the pole problem in the integration of timedependent partial di#erential equations on the sphere. The #lter's performance is demonstrated with numerical examples. Key Words. fast multipole method, spectral transform method, spectral truncation method, spherical harmonics, pole problem AMS#MOS# s...
On Approximation of Functions by Exponential Sums
, 2005
"... We introduce a new approach, and associated algorithms, for the efficient approximation of functions and sequences by short linear combinations of exponential functions with complexvalued exponents and coefficients. These approximations are obtained for a finite but arbitrary accuracy and typically ..."
Abstract

Cited by 17 (1 self)
 Add to MetaCart
We introduce a new approach, and associated algorithms, for the efficient approximation of functions and sequences by short linear combinations of exponential functions with complexvalued exponents and coefficients. These approximations are obtained for a finite but arbitrary accuracy and typically have significantly fewer terms than Fourier representations. We present several examples of these approximations and discuss applications to fast algorithms. In particular, we show how to obtain a short separated representation (sum of products of onedimensional functions) of certain multidimensional Green’s functions.
Combinatorial sublineartime fourier algorithms,” Submitted. Available at http://www.ima.umn.edu/∼iwen/index.html
, 2008
"... We study the problem of estimating the best k term Fourier representation for a given frequencysparse signal (i.e., vector) A of length N ≫ k. More explicitly, we investigate how to deterministically identify k of the largest magnitude frequencies of Â, and estimate their coefficients, in polynomia ..."
Abstract

Cited by 17 (5 self)
 Add to MetaCart
(Show Context)
We study the problem of estimating the best k term Fourier representation for a given frequencysparse signal (i.e., vector) A of length N ≫ k. More explicitly, we investigate how to deterministically identify k of the largest magnitude frequencies of Â, and estimate their coefficients, in polynomial(k, log N) time. Randomized sublinear time algorithms which have a small (controllable) probability of failure for each processed signal exist for solving this problem [24, 25]. In this paper we develop the first known deterministic sublinear time sparse Fourier Transform algorithm which is guaranteed to produce accurate results. As an added bonus, a simple relaxation of our deterministic Fourier result leads to a new Monte Carlo Fourier algorithm with similar runtime/sampling bounds to the current best randomized Fourier method [25]. Finally, the Fourier algorithm we develop here implies a simpler optimized version of the deterministic compressed sensing method previously developed in [30]. 1