topic MKL PARDISO produces FPE in Intel® oneAPI Math Kernel Library https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096477#M23593 <P>With the attached general matrix and rhs, PARDISO solution contains NaN.</P> <P>Major version: 11<BR /> Minor version: 1<BR /> Update version: 1<BR /> Product status:&nbsp; Product<BR /> Build: n20131010<BR /> Processor optimization: Intel(R) Advanced Vector Extensions (Intel(R) AVX) Enabled Processor</P> <P>linux 64 bit</P> <P>PARDISO_matrix.txt format:<BR /> n<BR /> index base (0 means 0 based)<BR /> nz<BR /> i j a (nonzeros)</P> <P>PARDISO_rhs.txt and PARDISO_sol.txt format:<BR /> n<BR /> dense/sparse (1 for dense and 0 for sparse)<BR /> index base(0 means 0 based)<BR /> nz (for dense format, nz=0)<BR /> i,b (nonzeros)</P> <P>&nbsp;</P> Thu, 11 Aug 2016 19:52:27 GMT xian-zhong_guous_cd- 2016-08-11T19:52:27Z MKL PARDISO produces FPE https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096477#M23593 <P>With the attached general matrix and rhs, PARDISO solution contains NaN.</P> <P>Major version: 11<BR /> Minor version: 1<BR /> Update version: 1<BR /> Product status:&nbsp; Product<BR /> Build: n20131010<BR /> Processor optimization: Intel(R) Advanced Vector Extensions (Intel(R) AVX) Enabled Processor</P> <P>linux 64 bit</P> <P>PARDISO_matrix.txt format:<BR /> n<BR /> index base (0 means 0 based)<BR /> nz<BR /> i j a (nonzeros)</P> <P>PARDISO_rhs.txt and PARDISO_sol.txt format:<BR /> n<BR /> dense/sparse (1 for dense and 0 for sparse)<BR /> index base(0 means 0 based)<BR /> nz (for dense format, nz=0)<BR /> i,b (nonzeros)</P> <P>&nbsp;</P> Thu, 11 Aug 2016 19:52:27 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096477#M23593 xian-zhong_guous_cd- 2016-08-11T19:52:27Z Have you evaluated MKL 11.3? https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096478#M23594 <P>Have you evaluated MKL 11.3? Such sort of problems has been fixed into this version.</P> Fri, 12 Aug 2016 04:21:50 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096478#M23594 Gennady_F_Intel 2016-08-12T04:21:50Z Same error using mkl 11.3.3. https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096479#M23595 <P>Same error using mkl 11.3.3.</P> Fri, 12 Aug 2016 17:42:30 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096479#M23595 xian-zhong_guous_cd- 2016-08-12T17:42:30Z Hi, https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096480#M23596 <P>Hi,</P> <P>Can you sent me your set of iparm and other pardiso parameters? With default parameters I have no Nans in result:</P> <P>i &nbsp; x<I> &nbsp; b<I></I></I></P> <P>0 -4.959544e+12 1.000000e+00<BR /> 1 6.805074e+14 -4.585740e-07<BR /> 2 -4.959372e+12 2.000000e+00<BR /> 3 6.805046e+14 -2.292070e-07<BR /> 4 -4.959751e+12 3.000000e+00<BR /> 5 6.805064e+14 -4.585740e-07<BR /> 6 2.673112e+13 4.000000e+00<BR /> 7 2.673106e+13 -2.292070e-07<BR /> 8 -4.959035e+12 5.000000e+00<BR /> 9 6.805016e+14 -4.585740e-07<BR /> 10 2.673123e+13 6.000000e+00</P> <P>Thanks,</P> <P>Alex&nbsp;</P> Sun, 14 Aug 2016 02:47:58 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096480#M23596 Alexander_K_Intel2 2016-08-14T02:47:58Z MKL PARDISO is called https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096481#M23597 <P>MKL PARDISO is called multiple times. Enclosed are the iparm for each call and solution is bad at last call. I notice some of iparam values becomes big numbers after the symbolic which are suspicious.</P> <P><BR /> iparm for phase = 11<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 22<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 11<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 22<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 11<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 22<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 319964<BR /> 220367<BR /> 205251<BR /> 51849976<BR /> 12857<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 319964<BR /> 144327<BR /> 1568<BR /> 103911000<BR /> iparm for phase = 11<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 22<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 11<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> -1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> iparm for phase = 22<BR /> 1<BR /> 2<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 13<BR /> 1<BR /> 0<BR /> 1<BR /> 0<BR /> 44089<BR /> 32155<BR /> 30783<BR /> 7666270<BR /> 1146<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> -1<BR /> 1<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 0<BR /> 1<BR /> 44089<BR /> 27467<BR /> 1028<BR /> 103911000</P> <P>A floating point exception has occurred: floating point exception [Invalid operation].</P> Mon, 15 Aug 2016 17:29:10 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096481#M23597 xian-zhong_guous_cd- 2016-08-15T17:29:10Z Hi, https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096482#M23598 <P>Hi,</P> <P>Are you call pardiso in loop with same matrix? I just try to implement reproducer of your issue. And if yes do you call pardiso with ido=-1 to free the memory?</P> <P>Thanks,</P> <P>Alex</P> Tue, 16 Aug 2016 04:23:42 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096482#M23598 Alexander_K_Intel2 2016-08-16T04:23:42Z PARDISO is called multiple https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096483#M23599 <P>PARDISO is called multiple time, each time with different matrix. Enclosed are my output message. cleanUpPARDISO calls PARDISO with phase = -1.</P> <P><BR /> Arithmetic type=float; matrix type=general; N=2727; NZ(A)=2727.<BR /> MKLPARDISO::init starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::init ends Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::symbolicFact starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::symbolicFact ends Tue Aug 16 09:53:29 2016<BR /> Arithmetic type=float; matrix type=general; N=2727; NZ(A)=2727.<BR /> Memory Estimate(MB)<BR /> --------------------------------------------------------------------------------<BR /> mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(A)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(fact)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> In-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> Out-of-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> --------------------------------------------------------------------------------<BR /> MKLPARDISO::numericalFact starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::numericalFact ends Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::solve starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::solve ends Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::cleanUpPARDISO starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::cleanUpPARDISO ends Tue Aug 16 09:53:29 2016<BR /> Arithmetic type=float; matrix type=general; N=81; NZ(A)=81.<BR /> MKLPARDISO::init starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::init ends Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::symbolicFact starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::symbolicFact ends Tue Aug 16 09:53:29 2016<BR /> Arithmetic type=float; matrix type=general; N=81; NZ(A)=81.<BR /> Memory Estimate(MB)<BR /> --------------------------------------------------------------------------------<BR /> mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(A)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(fact)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> In-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> Out-of-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> --------------------------------------------------------------------------------<BR /> MKLPARDISO::numericalFact starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::numericalFact ends Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::solve starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::solve ends Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::cleanUpPARDISO starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::cleanUpPARDISO ends Tue Aug 16 09:53:29 2016<BR /> Arithmetic type=float; matrix type=general; N=231360; NZ(A)=5144569.<BR /> MKLPARDISO::init starts Tue Aug 16 09:53:29 2016<BR /> MKLPARDISO::init ends Tue Aug 16 09:53:30 2016<BR /> MKLPARDISO::symbolicFact starts Tue Aug 16 09:53:30 2016<BR /> MKLPARDISO::symbolicFact ends Tue Aug 16 09:53:34 2016<BR /> Arithmetic type=float; matrix type=general; N=231360; NZ(A)=5144569.<BR /> Memory Estimate(MB)<BR /> --------------------------------------------------------------------------------<BR /> mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(A)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(fact)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> In-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 39&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 415&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 454&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> Out-of-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 39&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 312&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 351&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> --------------------------------------------------------------------------------<BR /> MKLPARDISO::numericalFact starts Tue Aug 16 09:53:34 2016<BR /> MKLPARDISO::numericalFact ends Tue Aug 16 09:53:36 2016<BR /> MKLPARDISO::solve starts Tue Aug 16 09:53:36 2016<BR /> MKLPARDISO::solve ends Tue Aug 16 09:53:36 2016<BR /> MKLPARDISO::cleanUpPARDISO starts Tue Aug 16 09:53:36 2016<BR /> MKLPARDISO::cleanUpPARDISO ends Tue Aug 16 09:53:36 2016<BR /> Arithmetic type=float; matrix type=symmetric positive definite; N=81; NZ(A)=81.<BR /> MKLPARDISO::init starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::init ends Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::symbolicFact starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::symbolicFact ends Tue Aug 16 09:53:38 2016<BR /> Arithmetic type=float; matrix type=symmetric positive definite; N=81; NZ(A)=81.<BR /> Memory Estimate(MB)<BR /> --------------------------------------------------------------------------------<BR /> mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(A)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(fact)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> In-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> Out-of-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> --------------------------------------------------------------------------------<BR /> MKLPARDISO::numericalFact starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::numericalFact ends Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::solve starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::solve ends Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::cleanUpPARDISO starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::cleanUpPARDISO ends Tue Aug 16 09:53:38 2016<BR /> Arithmetic type=float; matrix type=general; N=58320; NZ(A)=520078.<BR /> MKLPARDISO::init starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::init ends Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::symbolicFact starts Tue Aug 16 09:53:38 2016<BR /> MKLPARDISO::symbolicFact ends Tue Aug 16 09:53:39 2016<BR /> Arithmetic type=float; matrix type=general; N=58320; NZ(A)=520078.<BR /> Memory Estimate(MB)<BR /> --------------------------------------------------------------------------------<BR /> mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(A)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem(fact)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mem&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> In-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 61&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 65&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> Out-of-core&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 43&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 47&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<BR /> --------------------------------------------------------------------------------<BR /> MKLPARDISO::numericalFact starts Tue Aug 16 09:53:39 2016<BR /> MKLPARDISO::numericalFact ends Tue Aug 16 09:53:39 2016<BR /> MKLPARDISO::solve starts Tue Aug 16 09:53:39 2016<BR /> MKLPARDISO::solve ends Tue Aug 16 09:53:39 2016<BR /> MKLPARDISO::cleanUpPARDISO starts Tue Aug 16 09:53:39 2016<BR /> MKLPARDISO::cleanUpPARDISO ends Tue Aug 16 09:53:39 2016<BR /> A floating point exception has occurred: floating point exception [Invalid operation].</P> Tue, 16 Aug 2016 16:57:53 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096483#M23599 xian-zhong_guous_cd- 2016-08-16T16:57:53Z Can I ask you to provide full https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096484#M23600 <P>Can I ask you to provide full reproducer that call set of matrix to play with? I cannot reproduce issue with one matrix</P> <P>Thanks,</P> <P>Alex</P> Tue, 16 Aug 2016 17:53:31 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096484#M23600 Alexander_K_Intel2 2016-08-16T17:53:31Z I've attached the all https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096485#M23601 <P>I've attached the all matrices, rhs and sol for 5 calls.</P> Tue, 16 Aug 2016 21:48:26 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096485#M23601 xian-zhong_guous_cd- 2016-08-16T21:48:26Z Any update? Thanks. https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096486#M23602 <P>Any update? Thanks.</P> Fri, 19 Aug 2016 17:36:55 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096486#M23602 xian-zhong_guous_cd- 2016-08-19T17:36:55Z Hi, https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096487#M23603 <P>Hi,</P> <P>Could you also attach code that was aborted as you mentioned before? I implemented my version of program and it passed correctly (please see attached output file)</P> <P>Thanks,</P> <P>Alex&nbsp;</P> Mon, 22 Aug 2016 09:22:58 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096487#M23603 Alexander_K_Intel2 2016-08-22T09:22:58Z I solve in single precision https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096488#M23604 <P>I solve in single precision yet you are solving in double precision.</P> Tue, 23 Aug 2016 00:01:36 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096488#M23604 xian-zhong_guous_cd- 2016-08-23T00:01:36Z HI, https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096489#M23605 <P>HI,<BR /> I attach file with float precision run. Looks like the devil is in the details - can we connect directly to look at the code?<BR /> Thanks,</P> <P>Alex</P> Tue, 23 Aug 2016 04:41:19 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-produces-FPE/m-p/1096489#M23605 Alexander_K_Intel2 2016-08-23T04:41:19Z