Bugzilla – Attachment 796 Details for
Bug 1476
SelfAdjointEigenSolver gives wrong results for RowMajor matrices when using MKL
Home
|
New
|
Browse
|
Search
|
[?]
|
Reports
|
Requests
|
Help
|
Log In
[x]
|
Forgot Password
Login:
[x]
This bugzilla service is closed. All entries have been migrated to
https://gitlab.com/libeigen/eigen
Minimal failing example
test.cc (text/x-c++src), 953 bytes, created by
Björn Barz
on 2017-10-16 12:15:51 UTC
(
hide
)
Description:
Minimal failing example
Filename:
MIME Type:
Creator:
Björn Barz
Created:
2017-10-16 12:15:51 UTC
Size:
953 bytes
patch
obsolete
>#include <iostream> >#include <cstdlib> >#include <Eigen/Core> >#include <Eigen/Eigenvalues> >using namespace std; > >typedef Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> RMatrix; >typedef Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic, Eigen::ColMajor> CMatrix; > >int main() >{ > srand(0); > CMatrix cmat = CMatrix::Random(5, 5); > cmat *= cmat.transpose(); > RMatrix rmat = cmat; > > { > Eigen::SelfAdjointEigenSolver<CMatrix> eig(cmat.rows()); > eig.compute(cmat); > if (eig.info() == Eigen::Success) > cout << eig.eigenvectors() << endl << endl; > else > cerr << "Error" << endl; > } > > { > Eigen::SelfAdjointEigenSolver<RMatrix> eig(rmat.rows()); > eig.compute(rmat); > if (eig.info() == Eigen::Success) > cout << eig.eigenvectors() << endl << endl << endl; > else > cerr << "Error" << endl; > } > > return 0; >}
You cannot view the attachment while viewing its details because your browser does not support IFRAMEs.
View the attachment on a separate page
.
View Attachment As Raw
Actions:
View
Attachments on
bug 1476
: 796