子集犰狳场

nop*_*eva 0 c++ armadillo rcpp

如果我理解正确,字段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);
Run Code Online (Sandbox Code Playgroud)

但得到错误

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

joh*_*ack 6

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

F(2,1) = C;
Run Code Online (Sandbox Code Playgroud)

我认为它应该是:

F(2,0) = C;
Run Code Online (Sandbox Code Playgroud)

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

arma::field<arma::mat> G=F.subfield(1,0,2,0);
Run Code Online (Sandbox Code Playgroud)

要么:

arma::field<arma::mat> G=F.rows(1,2);
Run Code Online (Sandbox Code Playgroud)

有关此页面上子字段视图的更多信息.