我需要使用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)
。