topic How to use DFT to compute a cosine Fourier Tranform in Intel® oneAPI Math Kernel Library

How to use DFT to compute a cosine Fourier Tranform

Chang_Lei — Fri, 01 Oct 2010 04:11:37 GMT

As the title.

Best regards!

How to use DFT to compute a cosine Fourier Tranform

Dmitry_B_Intel — Fri, 01 Oct 2010 04:30:41 GMT

Hi Chang,

Intel MKL 10.2 and later versions support major part of FFTW3 interface. In particular 1D trigonometric transforms are supported. Example:

#include "fftw3.h" /* compile with -I$mkl/include/fftw */
fftw_plan plan = fftw_plan_r2r_1d(N,x,y,FFTW_REDFT01,FFTW_MEASURE);
...

For details see fftw documentation.

Thanks
Dima

How to use DFT to compute a cosine Fourier Tranform

Chang_Lei — Sun, 03 Oct 2010 06:09:38 GMT

Thanks.

Current version of MKL I'm using is CMKL 9.0, does that mean I must update my MKL to 10.2 or later versions?

And another question is: what's themain benefits of using MKL FFTwwrapperscomparing with using FFTwdirectly? Can I run it with muti-threads(OpenMP threads)?

With Best Regards!

Chang lei

How to use DFT to compute a cosine Fourier Tranform

Dmitry_B_Intel — Mon, 04 Oct 2010 07:30:48 GMT

Hi Chang,

Intel MKL 9.0 doesn't provide wrappers for trigonometric transforms. If you want to use MKL FFTW3 interfaces for this, you should upgrade to MKL 10.2 or later versions. Alternatively, with MKL 9.0 you can use trigonometric transforms with a special interface (look into MKL 9.0 documentation, section "Partial Differential Equations Support" -> "Trigonometric Transforms Routines"). Examples of using these can be found in directory examples/pdett of your MKL installation.

To the question aboutbenefits,Intel MKL provides performance optimized functions, and this is the main benefit.However, if youare concerned withportability,then a more popular interface may be the choice, such as FFTW.

OpenMP threads should be harnessedautomatically, provided that you linked your application with proper (threaded) library,perform recommended call sequence to employ threads [e.g. FFTW interface requires a call tofftw_plan_with_nthreads() be made in advance], and/or set environment [e.g. OMP_NUM_THREADS=4]. Details should be found inthe product documentation.

Thanks
Dima

How to use DFT to compute a cosine Fourier Tranform

Chang_Lei — Tue, 05 Oct 2010 07:07:26 GMT

It's very helpful and clear for me.

Thanks for yourpatiently reply. That's so nice of you.