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 1476 - SelfAdjointEigenSolver gives wrong results for RowMajor matrices when using MKL
Summary: SelfAdjointEigenSolver gives wrong results for RowMajor matrices when using MKL
Status: NEW
Alias: None
Product: Eigen
Classification: Unclassified
Component: Eigenvalues (show other bugs)
Version: 3.3 (current stable)
Hardware: x86 - 64-bit Linux
: Normal Wrong Result
Assignee: Nobody
URL:
Whiteboard:
Keywords:
Depends on:
Blocks:
 
Reported: 2017-10-16 12:15 UTC by Björn Barz
Modified: 2017-10-16 12:15 UTC (History)
4 users (show)



Attachments
Minimal failing example (953 bytes, text/x-c++src)
2017-10-16 12:15 UTC, Björn Barz
no flags Details

Description Björn Barz 2017-10-16 12:15:51 UTC
Created attachment 796 [details]
Minimal failing example

When defining EIGEN_USE_MKL and EIGEN_USE_MKL_ALL and linking against Intel MKL 2018.0, SelfAdjointEigenSolver will give wrong results for RowMajor matrices.

Attached is a minimal test case. It computes eigenvectors two times for the same matrix, one time stored in column-major order, the second time in row-major order. When compiled without MKL support, both results are equal. With MKL enabled, however, the result for the row-major matrix is the transpose of the correct result.

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