Developer Reference for Intel® oneAPI Math Kernel Library for C
Fourier Transform Functions
The general form of the discrete Fourier transform is
\[z_{k_1,k_2,\dots,k_d} = \sigma \times \sum_{j_d = 0}^{n_d - 1} \dots \sum_{j_2 = 0}^{n_2 - 1} \sum_{j_1 = 0}^{n_1 - 1} w_{j_1,j_2,\dots,j_d} exp \left( \delta i 2 \pi \sum_{l = 1}^{d} \frac{j_l k_l}{n_l} \right)\]
where
\(k_{l} = 0, \dots, n_l - 1\) for \(l = 1, \dots, d\),
\(\sigma\) is a scale factor,
\(\delta = -1\) for the forward transform, and \(\delta = 1\) for the inverse (backward) transform.
In the forward transform, the input (periodic) sequence \(w_{j_{1}, j_{2}, \dots, j_{d}}\) belongs to the set of complex-valued sequences and real-valued sequences. Respective domains for the backward transform are represented by complex-valued sequences and complex-valued conjugate-even sequences.
The Intel® oneAPI Math Kernel Library (oneMKL) provides an interface for computing a discrete Fourier transform through the fast Fourier transform algorithm. Prefixes Dfti in function names and DFTI in the names of configuration parameters stand for Discrete Fourier Transform Interface.
The manual describes the following implementations of the fast Fourier transform functions available in Intel® oneMKL :
Fast Fourier transform (FFT) functions for single-processor or shared-memory systems (see FFT Functions)
Cluster FFT functions for distributed-memory architectures (available only for Intel® 64 architectures)
Both FFT and Cluster FFT functions compute an FFT in five steps:
Allocate a fresh descriptor for the problem with a call to the DftiCreateDescriptor or DftiCreateDescriptorDM function. The descriptor captures the configuration of the transform, such as the dimensionality (or rank), sizes, number of transforms, memory layout of the input/output data (defined by strides), and scaling factors. Many of the configuration settings are assigned default values in this call which you might need to modify in your application.
Optionally adjust the descriptor configuration with a call to the DftiSetValue or DftiSetValueDM function as needed. Typically, you must carefully define the data storage layout for an FFT or the data distribution among processes for a Cluster FFT. The configuration settings of the descriptor, such as the default values, can be obtained with the DftiGetValue or DftiGetValueDM function.
Commit the descriptor with a call to the DftiCommitDescriptor or DftiCommitDescriptorDM function, that is, make the descriptor ready for the transform computation. Once the descriptor is committed, the parameters of the transform, such as the type and number of transforms, strides and distances, the type and storage layout of the data, and so on, are “frozen” in the descriptor.
Compute the transform with a call to the DftiComputeForward / DftiComputeBackward or DftiComputeForwardDM / DftiComputeBackwardDM functions as many times as needed. Because the descriptor is defined and committed separately, all that the compute functions do is take the input and output data and compute the transform as defined. To modify any configuration parameters for another call to a compute function, use DftiSetValue followed by DftiCommitDescriptor (DftiSetValueDM followed by DftiCommitDescriptorDM) or create and commit another descriptor.
Deallocate the descriptor with a call to the DftiFreeDescriptor or DftiFreeDescriptorDM function. This returns the memory internally consumed by the descriptor to the operating system.
All the above functions return an integer status value, which is zero upon successful completion of the operation. You can interpret a non-zero status with the help of the DftiErrorClass or DftiErrorMessage function.
The FFT functions support lengths with arbitrary factors. You can improve performance of the Intel® oneMKL FFT if the length of your data vector permits factorization into powers of optimized radices. See the Intel® oneAPI Math Kernel Library (oneMKL) Developer Guide Linux or Windows for specific radices supported efficiently.
Product and Performance Information |
|---|
Performance varies by use, configuration and other factors. Learn more at www.Intel.com/PerformanceIndex . Notice revision #20201201 |