https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardisoinit/m-p/1098986#M23743 <P>One of the purposes of the existence of pardisoinit() is stated in the MKL manual as &nbsp;"...set default iparm values in accordance with the matrix type". You do not need to call pardisoinit() at all; if you do, do not set any iparm[] values before calling pardisoinit. As you found out, many of those values will be set in pardisoinit().</P> Tue, 08 Dec 2015 00:32:30 GMT mecej4 2015-12-08T00:32:30Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardisoinit/m-p/1098985#M23742 <P>I guess this is another bug: the function pardisoinit resets iparm[34] to zero.</P> <P>I must reset iparm[34] = 1 between calling pardisoinit() and pardiso().</P> Mon, 07 Dec 2015 18:03:53 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardisoinit/m-p/1098985#M23742 Roman_F_ 2015-12-07T18:03:53Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardisoinit/m-p/1098986#M23743 <P>One of the purposes of the existence of pardisoinit() is stated in the MKL manual as &nbsp;"...set default iparm values in accordance with the matrix type". You do not need to call pardisoinit() at all; if you do, do not set any iparm[] values before calling pardisoinit. As you found out, many of those values will be set in pardisoinit().</P> Tue, 08 Dec 2015 00:32:30 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardisoinit/m-p/1098986#M23743 mecej4 2015-12-08T00:32:30Z