Rcpp 函数导致包崩溃

时间:2021-02-15 09:43:56

标签: r package rcpp rcpparmadillo

我有一个包 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    

1 个答案:

答案 0 :(得分:5)

查看您的包,我假设错误和崩溃来自 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;}')

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

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