I am using Pardiso, from the MKL that ships with icpc 16.0.0 on Mac OSX.
I have already computed a LU factorization, and I want to use Pardiso as an iterative solver using the previous LU factorization as a preconditionner. For that, I've set iparm = 21 for 2 digits of accuracy. But pardiso_64 returns an error -4 and iparm contains -18.
I don't understand what's going on as the last digit of iparm should be either 1, 2, 3, 4 or 5 according to the MKL documentation : https://software.intel.com/en-us/node/521691
In spite of Pardiso returns error=-4 and iparm=-18 the solution is correct for new matrix with element = 2.1 Could you check this behavior on larger matrix on your side or sent it to us also? It seem that only error output is incorrect in such way of pardiso using.
Best regards, Anton