Rcpp 函数导致包崩溃

如何解决Rcpp 函数导致包崩溃

我有一个包 https://github.com/tfrostig/RSEE，其中包含几 (3) 个 RcppArmadillo 函数。该软件包在其他计算机上运行良好。当我构建包时没有出现错误，但是每当我调用任何 RCPP 函数时，它都会导致 R 崩溃。

当我尝试使用我的单元测试时，出现错误：“Exited with status -1073741819”。

如果我使用 Rcpp::sourceCpp() 然后调用函数，一切正常。其他带有 Rcpp 函数的包运行良好。

例如：

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


// [[Rcpp::export]]
arma::mat localRegression(arma::mat weightmat,arma::mat modelmat,arma::vec xtemp) {
  return inv(modelmat.t() * weightmat * modelmat) * modelmat.t() * weightmat * xtemp;
}

使用 RSEE:::localRegression 会导致它崩溃。如果我使用 sourceCpp 加载源代码，然后调用 localRegression 它工作正常。

什么会导致这种类型的问题？

The session info is: 
R version 4.0.3 (2020-10-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 19042)

Matrix products: default

locale:
[1] LC_COLLATE=English_Israel.1252  LC_CTYPE=English_Israel.1252    LC_MONETARY=English_Israel.1252
[4] LC_NUMERIC=C                    LC_TIME=English_Israel.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] RSEE_0.1.0

loaded via a namespace (and not attached):
[1] compiler_4.0.3 tools_4.0.3    Rcpp_1.0.6    

解决方法

查看您的包，我假设错误和崩溃来自 arma::mat iterLowess(...,double epsmed = 10^(-6)) { 中的 src/RCPP_LOWESS.cpp。请注意，^ 不是 Cpp 中的幂运算，而是 byte XOR operation。此外，10 是一个整数，而 10.0 是一个双精度数，因此当编译器“应该”进行自动转换时，它可能会失败。

试试看

library(Rcpp)
f <- cppFunction('double powww(double x = 10^(-6)){
                 double y = x^2;
                 return y;}')

您会注意到这会引发错误。

几乎不可能给你一个确切的答案，因为我们在这个问题中的信息非常有限，但是因为你提到它在调用函数时立即崩溃（我假设它在调试状态下也是如此）我们应该查看错误的函数定义。