c++ - 如何在 Eigen C++ 中对称化稀疏矩阵？

Question

我A在 Eigen C++ 中有一个稀疏矩阵。现在我想将它对称到另一个稀疏矩阵Asym：

我希望它会像这样简单：

Eigen::SparseMatrix<FLOATDATA> A;
...
Eigen::SparseMatrix<FLOATDATA> Asym = 0.5*(A+A.transpose()); // error here

但是由于明显的原因，它给出了以下断言失败错误：

error: static assertion failed: THE_STORAGE_ORDER_OF_BOTH_SIDES_MUST_MATCH

我的问题是如何在 Eigen C++ 中巧妙地完成上述操作？

Answer 1

使代码编译的最简单方法是将转置矩阵评估为正确存储顺序的临时矩阵：

Eigen::SparseMatrix<FLOATDATA> Asym = 0.5*(A+Eigen::SparseMatrix<FLOATDATA>(A.transpose())));