Rcpp 函数导致包崩溃

Question

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

当我尝试使用单元测试时，出现错误：“退出状态为 -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

Answer 1

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

举个例子

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

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

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

Rcpp 函数导致包崩溃

Rcpp functions cause crash in package

r

package

rcpp

rcpparmadillo