Accuracy of Eigenvalue decomposition

Accuracy of Eigenvalue decomposition

I have a 3x3 matrix for which i want to calculate eigenvalues/-vectors. Using the function "ippmEigenValuesSym_m_64f" i get an ippStsSingularErr error meaning the matrix is singular, which is not the case.

The function "ippmLUDecomp_m_64f" however performs without problems on the same matrix.

The matrix elements are in the range of 1e-5 to 1e-6. Scaling the matrix by a factor of 10 yields valid results.

Are there any requirements on the size of the elements for the eigenvalue decomposition?

4 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

Hello,

what is your IPP version? It seems we fixed similar issue in the past. Please note that the latest available package is IPP 5.3 update 4. (and IPP 6.0 beta 2)

Regards,
Vladimir

Hello Vladimir,

I am using IPP 5.3 update 4.

The behaviour can be reproduced for example by the eigenvalue decomposition of the following matrix (i.e. the identity matrix scaled by 1e-6).

A=[1e-6, 0, 0;

0, 1e-6, 0;

0, 0, 1e-6 ];

Setting this to 1e-5 produces a valid result.

Thanks,

Michael

Hi Michael,

then it seems to be a different issue. I've informed our engineers on that, they will work on investigation of this. Please submit your issue report to Intel Premier Support.

Regards,
Vladimir

Leave a Comment

Please sign in to add a comment. Not a member? Join today