Wrong results from Pardiso solver with multi-threading

Hi,

I'm having problems with multi-threading in the Pardiso solver.

For a given sparse matrix (csr format in mat.txt) and right hand side (rhs.txt) I obtain a different and wrong solution with a number of threads equal or greater than 16.

I'm working on windows 7 and I'm compiling with the compiler C++ Visual Studio 2013 and the MKL 11.3 Update 4 (Intel Parallel Studio XE 2016 Update 4). I'm linking with the librairies: mkl_intel_lp64, mkl_intel_thread, mkl_core and libiomp5md.

Could you reproduce the problem? or is there a mistake in my Pardiso usage?

Attached please find my program (./test_mt.exe nb_threads).

Thanks in advance for your help.

Best regards.