Hello all,
I am using the extended eigensolver routines specifically dfeast_sygv function, and I get the following error
==>INFO code =: -4
Intel MKL Extended Eigensolvers Error: Matrix B is not positive definite.
But the matrix B is positive definite!!!
B =
+2.222e-01 +0.000e+00 +1.111e-01 +0.000e+00 +1.111e-01 +0.000e+00 +5.556e-02 +0.000e+00
+0.000e+00 +2.222e-01 +0.000e+00 +1.111e-01 +0.000e+00 +1.111e-01 +0.000e+00 +5.556e-02
+1.111e-01 +0.000e+00 +2.222e-01 +0.000e+00 +5.556e-02 +0.000e+00 +1.111e-01 +0.000e+00
+0.000e+00 +1.111e-01 +0.000e+00 +2.222e-01 +0.000e+00 +5.556e-02 +0.000e+00 +1.111e-01
+1.111e-01 +0.000e+00 +5.556e-02 +0.000e+00 +2.222e-01 +0.000e+00 +1.111e-01 +0.000e+00
+0.000e+00 +1.111e-01 +0.000e+00 +5.556e-02 +0.000e+00 +2.222e-01 +0.000e+00 +1.111e-01
+5.556e-02 +0.000e+00 +1.111e-01 +0.000e+00 +1.111e-01 +0.000e+00 +2.222e-01 +0.000e+00
+0.000e+00 +5.556e-02 +0.000e+00 +1.111e-01 +0.000e+00 +1.111e-01 +0.000e+00 +2.222e-01
using dfeast_syev function give the following eigenvalues.
+5.556e-02
+5.556e-02
+1.667e-01
+1.667e-01
+1.667e-01
+1.667e-01
+5.000e-01
+5.000e-01
all positive.
Here go the matrix A for the general eigenvalues problem
A =
+4.569e+02 +0.000e+00 -3.046e+02 +0.000e+00 +7.616e+01 +0.000e+00 -2.285e+02 +0.000e+00
+0.000e+00 +4.569e+02 +0.000e+00 +7.616e+01 +0.000e+00 -3.046e+02 +0.000e+00 -2.285e+02
-3.046e+02 +0.000e+00 +4.569e+02 +0.000e+00 -2.285e+02 +0.000e+00 +7.616e+01 +0.000e+00
+0.000e+00 +7.616e+01 +0.000e+00 +4.569e+02 +0.000e+00 -2.285e+02 +0.000e+00 -3.046e+02
+7.616e+01 +0.000e+00 -2.285e+02 +0.000e+00 +4.569e+02 +0.000e+00 -3.046e+02 +0.000e+00
+0.000e+00 -3.046e+02 +0.000e+00 -2.285e+02 +0.000e+00 +4.569e+02 +0.000e+00 +7.616e+01
-2.285e+02 +0.000e+00 +7.616e+01 +0.000e+00 -3.046e+02 +0.000e+00 +4.569e+02 +0.000e+00
+0.000e+00 -2.285e+02 +0.000e+00 -3.046e+02 +0.000e+00 +7.616e+01 +0.000e+00 +4.569e+02
can anybody help me here???
I'm getting something wrong???
Im using linux (debian testing).
thanks in advance.
ps. Sorry about the english.