eigen库selfadjointView问题

Question

每当我尝试使用 eigen 库使用任何矩阵或稀疏矩阵的 selfadjointView 属性时，我都会不断收到错误消息。下面是一个简单的代码来检查这一点。在我的程序中，我尝试使用自伴随矩阵：

#define EIGEN_YES_I_KNOW_SPARSE_MODULE_IS_NOT_STABLE_YET
#include <Eigen/Sparse>
#include <Eigen/Dense>
#include <Eigen/Core>
#include <iostream>

using namespace Eigen;
int main ()
{
SparseMatrix<float>  mat(3,3);
Matrix<float, 3, 1> vec;
std::cout<<mat.selfadjointView<>()*vec;
}

我得到的错误消息是： 错误：没有匹配的函数可用于调用“Eigen::SparseMatrix::selfadjointView()”

Answer 1

您必须指定模板参数，因此它应该读取 mat.selfadjointView() 或 mat.selfadjointView() 。第一个意味着它应该使用 mat 上三角部分的条目并填充下三角部分以使矩阵自伴。第二个则相反。