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

Question

我在 Eigen C++ 中有一个稀疏矩阵 A。现在我想将它对称化为另一个稀疏矩阵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())));