topic NRHS has to be one. The in Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038546#M20580
<P>NRHS has to be one. The calls are part of an Arnoldi iteration to find eigenvalues, so the rhs vectors are generated on the fly one at a time.</P>
<P>I tried installing the MKL for Composer XE 2016, but it would not work with XE 2013. I should point out that I cannot upgrade my version of Composer XE because I am not the administrator of the license.</P>Tue, 22 Sep 2015 02:54:35 GMTBrian_Murphy2015-09-22T02:54:35Zneed help understanding error=-4 from pardiso
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038542#M20576
<P>I am calling pardiso with mtype=11 (real nonsymmetric).</P>
<P>I first call it with phase=11, iparm(1)=0, and the result is error=0. So all is good.</P>
<P>I then call it with phase=22, and error=0. So all is still good.</P>
<P>I then call it multiple times with phase=33 with different rhs's. The first two calls return error=0, but the third call returns error=-4. What does this mean? Does -4 mean I have a bug somewhere, or can -4 be generated by something it doesn't like about the rhs? The rhs for the third call looks normal to me (i.e. no outlier elements).</P>
<P> </P>Mon, 21 Sep 2015 22:05:58 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038542#M20576Brian_Murphy2015-09-21T22:05:58ZI ran an unsymmetric problem
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038543#M20577
<P>I ran an unsymmetric problem with n=576, n_nz=2088 (matrix data file in Matrix Market format) with three different right hand sides, and there were no cases in which the error argument came back with a non-zero value. The code to read the data files and call Pardiso was in C. Conclusion: if there is a problem, it is not generic, but may surface only with specific inputs to Pardiso.</P>
<P>Can you make up an example (code + data) with the failure that you described? Did you ask for iterative refinement? If not, I don't see why you should get that specific error (-4).</P>Mon, 21 Sep 2015 23:42:00 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038543#M20577mecej42015-09-21T23:42:00ZThanks for the reply, mecej4.
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038544#M20578
<P>Thanks for the reply, mecej4. The main thing I am trying to find out are the implications of the -4 error code. I read the documentation about it, and looked at the iparam values returned by the call, but I cannot figure out what pardiso is trying to tell me with this error code. I don't want to burden anyone by sending the matrix and the three rhs's. If someone (pardiso's developers?) can tell me what sort of rhs can cause -4, that will give me a clue as to what to look for. Or if someone can say that the rhs cannot be the reason for -4, that would also help me.</P>
<P>For the time being, I am having my code just ignore the -4 error and continue running. Results agree with what I get using umfpack. Umfpack doesn't generate any error codes for this case.</P>
<P>The code+data that you mentioned did not come from me. I have run lots of other test cases and they run fine and do not generate any nonzero error codes.</P>Tue, 22 Sep 2015 01:04:45 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038544#M20578Brian_Murphy2015-09-22T01:04:45ZHi Brian,
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038545#M20579
<P>Hi Brian, </P>
<P>how about you define the nrhs from the begin of <SPAN style="font-size: 12px; line-height: 18px;"> phase=11 call. I recalled there is issue about nrhs change in phase =33 step in early MKL version. </SPAN></P>
<P><SPAN style="font-size: 12px; line-height: 18px;">There are two workaround, 1. define the nrhs=max of possible value at the begin of program. or use latest version. for example, the MKL 2016 in </SPAN><A href="http://software.intel.com/sites/campaigns/nest/">http://software.intel.com/sites/campaigns/nest/</A></P>
<P>Best Regards,</P>
<P>Ying</P>
<P> </P>Tue, 22 Sep 2015 02:25:36 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038545#M20579Ying_H_Intel2015-09-22T02:25:36ZNRHS has to be one. The
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038546#M20580
<P>NRHS has to be one. The calls are part of an Arnoldi iteration to find eigenvalues, so the rhs vectors are generated on the fly one at a time.</P>
<P>I tried installing the MKL for Composer XE 2016, but it would not work with XE 2013. I should point out that I cannot upgrade my version of Composer XE because I am not the administrator of the license.</P>Tue, 22 Sep 2015 02:54:35 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038546#M20580Brian_Murphy2015-09-22T02:54:35ZHi Brian,
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038547#M20581
<P>Hi Brian,<BR />
<BR />
Thanks for your question.<BR />
Can you please provide additional information, so that we can investigate this problem and find the root cause:<BR />
Set parameter msglvl to 1 and show me the output info and print number of pivots after factorization step (phase = 22) which value is iparm(30) ( iparm[29] in C ).</P>
<P><BR />
Best regards,<BR />
Maria.</P>Tue, 22 Sep 2015 12:32:25 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/need-help-understanding-error-4-from-pardiso/m-p/1038547#M20581MariaZh2015-09-22T12:32:25Z