How to install pardiso on matlab with `pthread` on windows?

11 次查看（过去 30 天）

wei zhang 2020-6-30

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/557149-how-to-install-pardiso-on-matlab-with-pthread-on-windows

评论： Rylan 2021-8-11

I am trying to use PARDISO in Matlab on windows. I notice that some compile options may not availible for windows system. This is my building script for pardisofree.mexw64 as an example.

mex  -largeArrayDims -lpardiso -lmwlapack -lmwblas ...
    -LD:\IntelSWTools\parallel_studio_xe_2019.0.045\compilers_and_libraries_2019\windows\compiler\lib\intel64_win -lifcoremt ...
    -LD:\zw\matlab\PARDISO\pthreads-w32-2-9-1-release\Pre-built.2\lib\x64 -lpthread ...
    -lm -output pardisofree ...
    common.cpp matlabmatrix.cpp sparsematrix.cpp pardisoinfo.cpp pardisofree.cpp

I have installed the intel parallel XE to get the windows version of gfortran, as lifcoremt.lib.

I also download pthread for windows version from sourceforge https://sourceforge.net/projects/pthreads4w/

So the command above runs ok for me. But I still have any way to activate the multithreads. I got no examples for that. How to build a matlab file to check it? Do I need to compile OPENMP into it? Any suggestion would be appreciated. I am really confused in this process.

14 个评论
显示 12更早的评论隐藏 12更早的评论

wei zhang 2020-7-8

编辑：wei zhang 2020-7-8

@Jeremy Hu I also tryied to add "openmp" into compile. This my entire build code for all the files. You may need this.

filename = {'pardisoinit.cpp','pardisoreorder.cpp','pardisofactor.cpp','pardisosolve.cpp','pardisofree.cpp'};
intel_lib_path = '-LD:\IntelSWTools\parallel_studio_xe_2019.0.045\compilers_and_libraries_2019\windows\compiler\lib\intel64_win';
pthread_lib_path = '-LD:\zw\matlab\PARDISO\pthreads-w32-2-9-1-release\Pre-built.2\lib\x64';
for i=1:length(filename)
mex('COMPFLAGS="$COMPFLAGS /openmp"',...
    '-largeArrayDims','-lpardiso','-lmwlapack',...
    intel_lib_path,'-lifcoremt','-liomp5md',...
    pthread_lib_path,'-lpthread',...
    '-lm', '-output',filename{i},...
    'common.cpp', 'matlabmatrix.cpp', 'sparsematrix.cpp', 'pardisoinfo.cpp', filename{i})
end

Unfortunately, Matlab still crash when I set the iparm[2] = 2 in the pardisoinfo.cpp, no matter before or after pardisoinit_ line. If I set the iparm[2] = 1, I could run the example m file welll. I guess the problem may from OMP_NUM_THREADS. But I don't know how to set it up.

Jeremy Hu 2021-8-11

Hi Rylan, I wasn't able to get PARDISO to work using MEX, but I did get it working using:

addpath('C:\Program Files (x86)\IntelSWTools\compilers_and_libraries_2020.1.216\windows\mkl\include'); % header files
addpath('C:\Program Files (x86)\IntelSWTools\parallel_studio_xe_2020.1.086\compilers_and_libraries_2020\windows\redist\intel64_win\mkl'); % dll
dllname = 'mkl_rt.dll';
global libraryname
libraryname = 'pardiso_mkl';
headername = 'mkl_pardiso.h';
headername2 = 'mkl_dss.h';
if ~libisloaded(libraryname)
    [notfound, warnings] = loadlibrary(dllname, headername, 'addheader', headername2, 'alias', libraryname);
end