使用迭代器从双循环构建矢量

时间:2016-01-03 01:25:16

标签: c++ iterator iteration rcpp armadillo

我正在学习如何使用迭代器,但我无法概念化如何使用它们从双循环构造向量。

这是一个例子:假设我想编写一个输入两个向量的函数,并输出通过组合每个向量中的单个元素生成的所有可能和的向量(这是一个人为的例子)。这是一个工作版本:

#include <RcppArmadillo.h>
using namespace Rcpp;
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::export]]
arma::colvec test(arma::colvec x, arma::colvec y) {
  arma::colvec out(x.size()*y.size());
  arma::colvec::iterator ix,iy;
  int count=0;
  for (ix=x.begin();ix!=x.end();ix++) {
    for (iy=y.begin();iy!=y.end();iy++) {
      out(count)=*ix+*iy;
      count+=1;
    }
  }
  return out;
}

/*** R
test(c(1:15),c(15:1))
*/

这很好用,但我觉得使用count作为out向量的ad-hoc迭代器似乎不太优雅,我觉得必须有一个更好的解决方案。我希望这样的东西能起作用:

arma::colvec::iterator ix,iy,io;
for (ix=x.begin(),io=out.begin();ix!=x.end(),io!=out.end();ix++,io++) {
  for (iy=y.begin();iy!=y.end();iy++,io++) {
    *io=*ix+*iy;
  }
}

这会导致我的计算机崩溃,所以我猜我无法使每个循环进行io迭代。如果有人能指出在这种情况下的良好做法,我将非常感激。

此致

1 个答案:

答案 0 :(得分:5)

每次嵌套循环结束时,你连续两次做io++次。这会使你的迭代器更快地移动,最终超出范围。 ix也超出了范围。

从外循环的迭代表达式中删除io++。然后它将等同于:

*(io++) = *ix + *iy;

并且io没有其他增量。

现在,要注意循环条件中的逗号运算符(这是造成崩溃的原因):

ix != x.end(), io != out.end()

- 评估ix != x.end(),丢弃其值,使用io != out.end()的值。应该有逻辑AND(&&)。如果,你应该在嵌套循环的条件下检查io != out.end()io增加的地方)。

但是你自己构建了out

arma::colvec out(x.size() * y.size());

没有支票就安全了。

代码(C ++ 11):

auto io = out.begin();

for (auto x_elem : x)
  for (auto y_elem : y)
    *(io++) = x_elem + y_elem;