Eigen 3 - 切片 MatrixXd 以反转列顺序

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

在 C++ 中使用 Eigen 我可以用循环来完成

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 文档，但没有找到任何有用的信息。

谢谢，马可

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).

Eigen 3 - 切片 MatrixXd 以反转列顺序

Eigen 3 - Slicing a MatrixXd to reverse the columns order

c++

r

rcpp

eigen3