Constructing spatial discretizations for sparse multivariate trigonometric polynomials that allow for a fast discrete Fourier transform

Abstract

The paper discusses the construction of high dimensional spatial discretizations for arbitrary multivariate trigonometric polynomials, where the frequency support of the trigonometric polynomial is known. We suggest a construction based on the union of several rank-1 lattices as sampling scheme. We call such schemes multiple rank-1 lattices. This approach automatically makes available a fast discrete Fourier transform (FFT) on the data. The key objective of the construction of spatial discretizations is the unique reconstruction of the trigonometric polynomial using the sampling values at the sampling nodes. We develop different construction methods for multiple rank-1 lattices that allow for this unique reconstruction. The symbol M denotes the total number of sampling nodes within the multiple rank-1 lattice. In addition, we assume that the multivariate trigonometric polynomial is a linear combination of T trigonometric monomials. The ratio of the number M of sampling points that are sufficient for the unique reconstruction to the number T of distinct monomials is called oversampling factor in this context. The presented construction methods for multiple rank-1 lattices allow for estimates of this number M . Roughly speaking, the oversampling factor M / T is independent of the spatial dimension and, with high probability, only logarithmic in T , which is much better than the oversampling factor that is expected for a sampling method that uses one single rank-1 lattice. The newly developed approaches for the construction of spatial discretizations are probabilistic methods. The arithmetic complexity of these algorithms depend only linearly on the spatial dimension and, with high probability, only linearly on T up to some logarithmic factors. Furthermore, we analyze the computational complexities of the resulting FFT algorithms, that exploits the structure of the suggested multiple rank-1 lattice spatial discretizations, in detail and obtain upper bounds in O ( M log ⁡ M ) , where the constants depend only linearly on the spatial dimension. With high probability, we construct spatial discretizations where M / T ≤ C log ⁡ T holds, which implies that the complexity of the corresponding FFT converts to O ( T log 2 ⁡ T ) .