Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Gaure__Simen
Beginner
276 Views

MKL's FFTW wrappers block FFTW on linux (at least)

Hi,

This problem popped up downstream in R (https://bugs.r-project.org/bugzilla/show_bug.cgi?id=17443) when it is linked with MKL on Linux (I don't know about Windows/OSX). R is a an extensible program, one may write dlls/shared libraries in C which is linked into R at runtime. Now, here is a problem. MKL contains wrappers for fftw, some of them do not work, for instance multidimensional r2r transforms, or strided r2r transforms (which can easily be utilized to make multidimensional ones).  I make a shared library for R which uses these parts of fftw3:

# gcc -shared -o foo.so foo.c -lfftw3

Everything seems fine. I start R, it is linked with MKL to have a nice optimized, parallel BLAS, say -lmkl_rt, and loads symbols from there. When I from R load (dlopen) my foo.so, the fftw-references in there are resolved to MKL which is already loaded in R, and they do not work. There is no easy way for my foo.so to get hold of the real fftw3-routines. A workaround is a little bit cumbersome, I can either link R with parts of MKL static (like -Bstatic,-lmkl_gf_lp64,-Bdynamic -lmkl_gnu_threads ...), or link my foo.so with a static libfftw3.a, or preload R with LD_PRELOAD=libfftw3.so.

My question is, would it be possible to ship MKL's fftw-wrappers separate from the working parts of MKL so that it would be possible to avoid the poor fftw masquerade?

0 Kudos
1 Reply
Ying_H_Intel
Employee
276 Views

Hi Simen,

​Thank you a lot for report the issue here and R,  also share your experience for workaround the problem.

For your questions, 
would it be possible to ship MKL's fftw-wrappers separate from the working parts of MKL so that it would be possible to avoid the poor fftw masquerade?
​Yes,  I can understand two FFTW could be problem.  and  it is possible to separate the MKL Blas , MKL FFT etc .
The tool is to use the MKL custom DLL,  which is supposed for that purpose.
https://software.intel.com/en-us/mkl-linux-developer-guide-building-custom-shared-objects
​for example, under Linux:
>cp -r /opt/intel/compilers_and_libraries_2018/linux/mkl/tools/builder .
​>cd builder

>make libintel64 parallel=gnu export=blas_example_list name=libRblas
Which will build one Rblas library based on selected BLAS API. Thus no MKL FFTW inference is mixed.   

​I will consider to add it in our publish articles and notify R if possible.

https://software.intel.com/en-us/articles/using-intel-mkl-with-r​

Best Regards,
​Ying

Reply