I'm a user of MKL 2018 up1 parallel studio.
Usually I handle the large matrix(>50Mx50M) using schur complement.
'iparm = 1' is the my general setting and it shows the schur complement matrix.
However, when I changed the setting 'iparm' from 1 to 2, they shows the error message during factorization (phase=2)
Error code is -4. (not SPD)
Already, I checked iparam and matrix->a[iparam], but there are no zero pivot or wrong conditions for SPD.
I think, if the matrix is factorized with 'iparm = 1', then 'iparm = 2' should be also OK.
I ask your advice.
Thank you very much in advance!!!
Check the Bug list since MKL 2018 u1 at https://software.intel.com/en-us/articles/intel-math-kernel-library-intel-mkl-2018-bug-fixes-list
I haven't seen the exact report about the problem. But what kind of your machine memory and cpu? How about link with ILP64 or when the input size is small , does it work?
or anyway, could you please try MKL 2018 u3 if possible
Thank you for yout quick reply.
As you mentioned, I searched the fix lists, but I can't find any related contents.
Now I'm trying to conduct my program at MKL 2018 u3.
I'll tell you the result ASAP.
Followings are my system info and my program settings
- Xeon E5-2630 v2 and 500GB memory
- Unix redhat, C++
- mtype : 2 (because this matrix is generated from passive circuit design.)
- iparm : =1, =3, =8, =1 (or 2), and [others]=0.
- Schur complement size (MxN) : 3k x 3k
- included library : "cdft_core.a", "intel_ilp64.a", intel_thread.a", "core.a", "blacs_intelmpi_ilp64.a"
Do you have any idea what might be wrong from this condition?