topic Strange (and seemingly wrong) result from Pardiso in IntelĀ® oneAPI Math Kernel Library
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Strange-and-seemingly-wrong-result-from-Pardiso/m-p/1009002#M19139
<P>I obtained strange (and seemingly wrong) result from Pardiso. I used Pardiso to solve a very ill-conditioned complex linear system arising from a circuital analysis of a three-dimensional conducting structure. I have no error from the analysis and numerical factorization phases.</P>
<P>Data and results are reported in the attached document. For instance, the Pardiso parameters and log are reported in the table, for three cases (three frequencies of analysis). When not present, parameters have zero value.</P>
<P>The right hand vector is shown in Figure 1: there are only two not-zero entries (corresponding to the two point excited on the structure). Figure 2 reports the magnitude of the solution vector. The case relative to the green curve (1500 Hz) exhibits unexpected values and leads to wrong results. There is no physical reason to have values with this magnitude. The other two cases (two frequencies very near to the wrong one) give the right result.</P>
<P>Do you suggest some changes in the parameters (the matrix is very ill-conditioned), or there is a way to recognize the wrong case?</P>
<P>Thank you,</P>
<P>Alessandro M.</P>Fri, 03 Apr 2015 15:03:48 GMTAlessandro_M_2015-04-03T15:03:48ZStrange (and seemingly wrong) result from Pardiso
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Strange-and-seemingly-wrong-result-from-Pardiso/m-p/1009002#M19139
<P>I obtained strange (and seemingly wrong) result from Pardiso. I used Pardiso to solve a very ill-conditioned complex linear system arising from a circuital analysis of a three-dimensional conducting structure. I have no error from the analysis and numerical factorization phases.</P>
<P>Data and results are reported in the attached document. For instance, the Pardiso parameters and log are reported in the table, for three cases (three frequencies of analysis). When not present, parameters have zero value.</P>
<P>The right hand vector is shown in Figure 1: there are only two not-zero entries (corresponding to the two point excited on the structure). Figure 2 reports the magnitude of the solution vector. The case relative to the green curve (1500 Hz) exhibits unexpected values and leads to wrong results. There is no physical reason to have values with this magnitude. The other two cases (two frequencies very near to the wrong one) give the right result.</P>
<P>Do you suggest some changes in the parameters (the matrix is very ill-conditioned), or there is a way to recognize the wrong case?</P>
<P>Thank you,</P>
<P>Alessandro M.</P>Fri, 03 Apr 2015 15:03:48 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Strange-and-seemingly-wrong-result-from-Pardiso/m-p/1009002#M19139Alessandro_M_2015-04-03T15:03:48ZHello,Ā Alessandro
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Strange-and-seemingly-wrong-result-from-Pardiso/m-p/1009003#M19140
<P>Hello, <SPAN style="font-size: 12px; line-height: 18px;">Alessandro</SPAN></P>
<P>I see that you use parallel reordering (iparm(2)=3), so the results may be different for diffrent calls. Please use iparm(2)=2, it has to produce stable results.</P>
<P>For <SPAN style="font-size: 12px; line-height: 18px;">very ill-conditioned matrix direct solver can produce poor solution. In such case the default value of iparm(10)=13 (pivoting pertrubation) may be </SPAN>inappropriate. Could you play with iparm(10) = from 1 to 15. It is possible that one of the value resolve the problem.</P>
<P>Best regards, Anton</P>Wed, 08 Apr 2015 08:32:39 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Strange-and-seemingly-wrong-result-from-Pardiso/m-p/1009003#M19140Anton_A_Intel2015-04-08T08:32:39Z