New user self-registration is currently disabled. Please email eigen-core-team @ lists.tuxfamily.org if you need an account.
Bug 854 - Error in SelfAdjointEigenSolver::computeDirect
Error in SelfAdjointEigenSolver::computeDirect
Status: RESOLVED FIXED
Product: Eigen
Classification: Unclassified
Component: Eigenvalues
3.2
x86 - 64-bit Windows
: Normal Unknown
Assigned To: Nobody
:
Depends on:
Blocks:
  Show dependency treegraph
 
Reported: 2014-08-15 00:28 UTC by dario.sassi
Modified: 2014-08-21 10:50 UTC (History)
4 users (show)



Attachments

Description dario.sassi 2014-08-15 00:28:40 UTC
With a Matrix3d like

 850.961   51.966    0
  51.966  254.841    0
   0        0        0

using SelfAdjointEigenSolver::computeDirect the eigenvalues are good, the eigenvectors are good but in a wrong order.

With SelfAdjointEigenSolver::compute the order of eigenvectors is good, as explained in documentation.
Comment 1 Gael Guennebaud 2014-08-21 10:50:50 UTC
Thank you for the report.

https://bitbucket.org/eigen/eigen/commits/02b8a8c4089e/
Changeset:   02b8a8c4089e
User:        ggael
Date:        2014-08-21 10:49:09
Summary:     Bug 854: fix numerical issue in SelfAdjointEigenSolver::computeDirect for 3x3 matrices. The tolerance to detect stable cross products was too optimistic.
Add respective unit tests.


https://bitbucket.org/eigen/eigen/commits/cf7f702f6bbf/
Changeset:   cf7f702f6bbf
Branch:      3.2
User:        ggael
Date:        2014-08-21 10:49:09
Summary:     Bug 854: fix numerical issue in SelfAdjointEigenSolver::computeDirect for 3x3 matrices. The tolerance to detect stable cross products was too optimistic.
Add respective unit tests.
(grafted from 02b8a8c4089e62a36b7f72bbfc57c65e326848e9)

Note You need to log in before you can comment on or make changes to this bug.