c++ - 如何在 C++ 的 Eigen 中交换对称稀疏矩阵中的列和行？

Question

我正在尝试为稀疏对称矩阵实现以下密集矩阵运算:

    dm.col(j).swap(dm.col(i));
    dm.row(j).swap(dm.row(i));

为了交换稀疏对称矩阵中的两列(和两行)，我试图生成一个置换矩阵，以便我可以将它传递给 .twistedBy():

    Eigen::PermutationMatrix<Eigen::Dynamic, Eigen::Dynamic> P;
    P.setIdentity();
    P.col(j).swap(P.col(i));

    sm = sm.twistedBy(P);

不幸的是，置换矩阵不支持 .col().swap()。我该怎么办？

这个问题阐明了.twistedBy的使用，但没有解释如何构造置换矩阵: Permuting sparse matrices in Eigen

我很难理解 eigen 的文档: https://eigen.tuxfamily.org/dox/classEigen_1_1PermutationMatrix.html

也将不胜感激任何对此的一般帮助!

感谢您的宝贵时间!

Answer 1

您正在寻找applyTranspositionOnTheRight :

Eigen::PermutationMatrix<Eigen::Dynamic, Eigen::Dynamic >::PermutationMatrix P(n);
P.setIdentity();
P.applyTranspositionOnTheRight(j, i);