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

Summary: SelfAdjointEigenSolver gives wrong results for RowMajor matrices when using MKL
Product: Eigen Reporter: Björn Barz <bjoern.barz>
Component: EigenvaluesAssignee: Nobody <eigen.nobody>
Status: NEW ---    
Severity: Wrong Result CC: chtz, gael.guennebaud, jacob.benoit.1, jitseniesen
Priority: Normal    
Version: 3.3 (current stable)   
Hardware: x86 - 64-bit   
OS: Linux   
Whiteboard:
Attachments:
Description Flags
Minimal failing example none

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.