RcppArmadillo find() function 无法识别 Rcpp 中的两个相等值

Question

让我们假设我有一个值向量l和一个我称之为x的洗牌版本。 对于向量x的任何值i ，我可以应用arma::uvec index = find(l==x(i)) 。 但是，如果我的输入向量是x+mu并且我将其转换为 C++ function，那么find()并不总是能够识别相等性。 我给你一个代码的例子：

#include <RcppArmadilloExtensions/sample.h>
// [[Rcpp::depends(RcppArmadillo)]]
using namespace Rcpp;
using namespace arma;
using namespace std;


// [[Rcpp::export]]
arma::vec my_function(arma::vec x, double mu, arma::vec l){
    arma::vec results(x.n_rows);
    results.zeros();
    for(int i=0; i<x.n_rows; ++i){
      Rcout << i << "\n";
      double rescaled = x(i)-mu;
      arma::uvec index = find(l==rescaled);
      results(i) = index(0);
     }
    return results;
  }

以下 R 代码返回错误。

set.seed(1)
l <- seq(-4.975,5.025,length=201)
mu <- 3.1
x <- sample(l,length(l))+mu
my_function(x = x, mu = mu, l = l)