为什么'&'运算符在 RcppArmadillo 中不起作用

Question

我正在尝试从通过元素“&”连接的另外两个逻辑向量中获取逻辑向量：

//[[Rcpp::export]]
arma::uvec test1(arma::vec t1, double R1, double R2){
arma::uvec t = (t1 >= R1) & (t1 < R2);
return t;
}

当我尝试编译时它返回以下错误

error: no match for 'operator&' (operand types are 'arma::enable_if2<true, const arma::mtOp<unsigned int, arma::Col<double>, arma::op_rel_gteq_post> >::result {aka const arma::mtOp<unsigned int, arma::Col<double>, arma::op_rel_gteq_post>}' and 'arma::enable_if2<true, const arma::mtOp<unsigned int, arma::Col<double>, arma::op_rel_lt_post> >::result {aka const arma::mtOp<unsigned int, arma::Col<double>, arma::op_rel_lt_post>}')
arma::uvec t = (t1 >= R1) & (t1 < R2);
                          ^

我不知道发生了什么。我猜犰狳做事不同。但是我找不到任何资源来帮助我解决问题。任何帮助，将不胜感激！非常感谢！

Answer 1

我不知道发生了什么。我猜犰狳做事不同。但是我找不到任何资源来帮助我解决问题。

这里的最终来源是the Armadillo docs。如果您转到the section on operators，您会发现& 运算符不是“[o] Mat、Col、Row 和 的重载运算符之一Cube 类。”因此，如果您想要这样的操作员，您必须自己编写代码（或者查看其他人是否已经在互联网上发布了它）。 Rcpp::NumericVectors有这样一个算子：

#include <Rcpp.h>

// [[Rcpp::export]]
Rcpp::LogicalVector test1(const Rcpp::NumericVector& t1, double R1, double R2){
    return (t1 >= R1) & (t1 < R2);
}

test1(1:10, 3, 7)
# [1] FALSE FALSE  TRUE  TRUE  TRUE  TRUE FALSE FALSE
# [9] FALSE FALSE

当然，如果您的其余代码确实依赖于 Armadillo，那也无济于事。

更新：只需使用&&

正如 mtall 在 cmets 中指出的那样，&& 运算符实际上是可用的，即使它没有在 Armadillo 文档中讨论（也许它不像我那样是终极的来源想）。

所以，只需将您的代码更改为以下内容：

#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]

//[[Rcpp::export]]
arma::uvec test1(arma::vec t1, double R1, double R2){
    arma::uvec t = (t1 >= R1) && (t1 < R2);
    return t;
}

根据您对 cme​​ts 的问题和回复，我相信它会按照您希望的方式工作：

test1(1:10, 3, 7)
      [,1]
 [1,]    0
 [2,]    0
 [3,]    1
 [4,]    1
 [5,]    1
 [6,]    1
 [7,]    0
 [8,]    0
 [9,]    0
[10,]    0