New user self-registration is disabled due to spam. Please email eigen-core-team @ lists.tuxfamily.org if you need an account.
Before reporting a bug, please make sure that your Eigen version is up-to-date!
Bug 854 - Error in SelfAdjointEigenSolver::computeDirect
Summary: Error in SelfAdjointEigenSolver::computeDirect
Status: RESOLVED FIXED
Alias: None
Product: Eigen
Classification: Unclassified
Component: Eigenvalues (show other bugs)
Version: 3.2
Hardware: x86 - 64-bit Windows
: Normal Unknown
Assignee: Nobody
URL:
Whiteboard:
Keywords:
Depends on:
Blocks:
 
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.