Bugzilla – Attachment 947 Details for
Bug 1739
GeneralizedSelfAdjointEigenSolver not preallocating all internal data
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
Repro case
ReproEigenAllocation.cpp (text/plain), 769 bytes, created by
Flavio Protasio Ribeiro
on 2019-08-09 20:24:40 UTC
(
hide
)
Description:
Repro case
Filename:
MIME Type:
Creator:
Flavio Protasio Ribeiro
Created:
2019-08-09 20:24:40 UTC
Size:
769 bytes
patch
obsolete
>#define EIGEN_RUNTIME_NO_MALLOC 1 >#define EIGEN_NO_AUTOMATIC_RESIZING 1 > >#include <iostream> >#include <Eigen/Core> >#include <Eigen/Eigenvalues> > >using namespace Eigen; >using namespace std; > >int main(int argc, char **argv) >{ > int N = 5; > GeneralizedSelfAdjointEigenSolver< MatrixXcf > ei(N); > MatrixXcf A = MatrixXcf::Random(N, N); > MatrixXcf B = MatrixXcf::Random(N, N); > > A = (A + A.adjoint()).eval(); > A = (A * A.adjoint()).eval(); > B = (B + B.adjoint()).eval(); > B = (B * B.adjoint()).eval(); > A = A * A; > B = B * B; > > Eigen::internal::set_is_malloc_allowed(false); > ei.compute(A,B, ComputeEigenvectors | Ax_lBx); > Eigen::internal::set_is_malloc_allowed(true); > auto v = ei.eigenvalues(); > cout << v << endl << 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 1739
: 947 |
948