Question

我需要使用RcppEigen来反转MatrixXd列的顺序。

在R我只会做

> M = matrix(1:9, ncol = 3)
> M
     [,1] [,2] [,3]
[1,]    1    4    7
[2,]    2    5    8
[3,]    3    6    9
> M = M[, 3:1]
> M
     [,1] [,2] [,3]
[1,]    7    4    1
[2,]    8    5    2
[3,]    9    6    3

在使用Eigen的C++中，我可以使用循环

Eigen::MatrixXd m1(3, 3);
Eigen::MatrixXd m2(3, 3);
m1 << 1, 4, 7,
      2, 5, 8,
      3, 6, 9;

for (int i = 0; i < 3; i++){
    m2.col(i) =  m1.col(2 - i);
}

return m2;

有更好的方法吗？我还查看了Eigen::Array文档，但没有找到任何有用的内容。

谢谢，Marco

Answer 1

尝试m1.rowwise().reverse()同样使用colwise()进行行反转。

Answer 2

在上面的问题中，我忘了提到我想要反转最后 m 列的顺序。因此，3x3矩阵的最后2列的示例：
在R

    > M = matrix(1:9, ncol = 3)
    > M
         [,1] [,2] [,3]
    [1,]    1    4    7
    [2,]    2    5    8
    [3,]    3    6    9
    > M = M[, 3:2]
    > M
         [,1] [,2]
    [1,]    7    4
    [2,]    8    5
    [3,]    9    6

在C++中使用Eigen：

using Eigen::MatrixXd;
Eigen::MatrixXd m1(3, 3);
Eigen::MatrixXd m2(3, 2);
m1 << 1, 2, 3,
      4, 5, 6,
      7, 8, 9;
m2 = m1.block(0, 1, 3, 2).rowwise().reverse();
return m2;

诀窍。
m2 = m1.rowwise().reverse().block(0, 0, 3, 2);也有效，但可能效率较低对于VectorsX？可以使用tail(m)，head(m)或segment(i, m)。

特征3 - 切割MatrixXd以反转列顺序

2 个答案: