特征3-对矩阵Xd进行切片以反转列顺序



我需要使用RcppEigen反转矩阵Xd的列的顺序。

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++中,使用本征I可以使用环路

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

尝试m1.rowwise().reverse()类似地使用colwise()进行行反转。

在上面的问题中,我忘记提到我只想反转最后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++中使用特征:

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)

最新更新