我想计算数组的特征值/向量而不是矩阵。
我尝试了EigenSolver<ArrayXf>
,但这给出了编译错误。
我可以将数组复制到矩阵中。但那是浪费记忆。
以下代码给出了分段错误。
测试1:
#include <Eigen/Eigen>
using namespace Eigen;
int main() {
ArrayXf A = ArrayXf::Ones(3,3);
EigenSolver<MatrixXf> es(A);
}
结果:
<...>/Eigen/src/Core/util/XprHelper.h:130:
Eigen::internal::variable_if_dynamic<T, Value>::variable_if_dynamic(T)
[with T = long int; int Value = 1]: Assertion `v == T(Value)' failed.
Aborted (core dumped)
我也试过EigenSolver<MatrixXf> es(A.matrix())
。但这也行不通。
测试2:
#include <Eigen/Eigen>
using namespace Eigen;
int main() {
ArrayXf A = ArrayXf::Ones(3,3);
EigenSolver<MatrixXf> es(A.matrix());
}
结果:
<...>/XprHelper.h:130: Eigen::internal::variable_if_dynamic<T, Value>::variable_if_dynamic(T) [with T = long int; int Value = 1]: Assertion `v == T(Value)' failed.
Aborted (core dumped)
问题是ArrayXf
是一维数组,而你想要一个2D数据:ArrayXXf
。
一些历史:我们在引入Array之前提出了VectorXf / MatrixXf名称,没有自然名称来区分1D和2D,因此单个X与双XX ...