访问eigen3中的特征值

Accessing eigenvalues in eigen3

本文关键字：特征值 eigen3 访问      更新时间：2023-10-16

我正在尝试设置特征值求解器，看起来我面临一些复杂情况。代码可以很好地编译，但在某些地方会出错。失败的代码段是

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;
}
相关文章：