我正在尝试设置 eigen3 特征值求解器,看起来我遇到了一些并发症。代码编译得很好,但在某些时候会失败。失败的那段代码是
EigenSolver<MatrixXd> eigensolver(A);
cout << "tell me something" << endl;
Matrix<complex<double>, -1, 1, 0, -1, 1> E = eigensolver.eigenvalues();
cout << "tell me something more" << endl;
cout << E;
以及带有错误消息的输出:
tell me something
tell me something more
(3.5,1.93649)
class_out: /usr/include/eigen3/Eigen/src/Core/DenseCoeffsBase.h:407: Eigen::internal::traits<Derived>::Scalar &Eigen::DenseCoeffsBase<Derived, 1>::operator()(Eigen::internal::traits<Derived>::Index) [with Derived = Eigen::Matrix<std::complex<double>, -1, 1, 0, -1, 1>]: Assertion `index >= 0 && index < size()' failed.
(3.5,-1.93649)Aborted
这看起来有点奇怪,但看起来不像是
<<
操作符的问题,像普通的一样cout << eigensolver.eigenvalues();
工作正常以及
cout
对于普通矩阵。另外,我用 Mathematica 检查了 A 的特征值,结果正好是 (3.5,1.93649) 和 (3.5,-1.93649)。有谁知道为什么会发生这种情况,或者有人可以建议其他方式来访问特征值?
最佳答案
解决方案似乎很简单:
EigenSolver<MatrixXd> eigensolver(A);
cout << "tell me something" << endl;
complex<double> E;
cout << "tell me something more" << endl;
for(int i = 0; i < A.rows(); ++i){
E = eigensolver.eigenvalues().col(0)[i];
cout << E << endl;
}
关于c++ - 访问 eigen3 中的特征值,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/15553766/