子集犰狳场

时间:2015-03-22 13:36:01

标签: c++ rcpp armadillo

如果我理解正确,Armadillo中的字段就像任意对象的列表一样。例如,一组不同大小的矩阵,或矩阵和向量。在文档中,我看到了可以与cube一起使用的类型slices,因此您可以使用它们进行子集化。但是,似乎没有特定的方法来对字段进行子集化。

我的代码的简化版本是:

arma::mat A = eye(2,2);
arma::mat B = eye(3,3)*3;
arma::mat C = eye(4,4)*4;
arma::field<arma::mat> F(3,1);
F(0,0) = A;
F(1,0) = B;
F(2,1) = C;

// to get matrices B and C
F.slices(1,2);

但得到错误

Error: field::slices(): indicies out of bounds or incorrectly used

1 个答案:

答案 0 :(得分:5)

首先,您提供的代码中存在一个小错误:

F(2,1) = C;

我认为它应该是:

F(2,0) = C;

其次,函数slices()仅对3D字段有效。但是,您的字段F只是一个2D字段,因为您只在构造函数中指定行和列。要访问矩阵B和C,您可以使用:

arma::field<arma::mat> G=F.subfield(1,0,2,0);

或:

arma::field<arma::mat> G=F.rows(1,2);

有关this page的子字段视图的详细信息。