C++ 从Rcpp函数接收结果后会话崩溃
下面的代码编译和执行正确，但每次我运行它时，我的R会话在完成后不久就会出现致命错误。我正在运行R版本3.3.2和Rtools 3.3 有什么我错过的吗？我怎样才能追踪导致坠机的原因C++ 从Rcpp函数接收结果后会话崩溃,c++,r,rcpp,C++,R,Rcpp,下面的代码编译和执行正确，但每次我运行它时，我的R会话在完成后不久就会出现致命错误。我正在运行R版本3.3.2和Rtools 3.3 有什么我错过的吗？我怎样才能追踪导致坠机的原因 #include<Rcpp.h> using namespace Rcpp; NumericMatrix dupCheckRcpp(NumericMatrix x) { int nrow, ncol; int i, j, k, m, n; bool flag; Numer
#include<Rcpp.h>
using namespace Rcpp;
NumericMatrix dupCheckRcpp(NumericMatrix x) {
int nrow, ncol;
int i, j, k, m, n;
bool flag;
NumericMatrix dupMat(300,ncol);
n = 0;
nrow = 0; ncol = 0;
nrow = x.nrow();
ncol = x.ncol();
for (i = 0; i < nrow - 1 ; ++i) {
for (j = i + 1; j < nrow; ++j) {
flag = TRUE;
for (k = 0; k < ncol; ++k) {
if (x(i,k) != x(j,k)) {
flag = FALSE;
break;
}
}
if (flag == TRUE) {
for (m = 0; m < ncol; ++m) {
dupMat(n,m) = x(i,m);
}
n = n + 1;
}
}
}
return dupMat;
}
#包括
使用名称空间Rcpp；
NumericMatrix dupCheckRcpp（NumericMatrix x）{
int nrow，ncol；
int i，j，k，m，n；
布尔旗；
数字矩阵dupMat（300，ncol）；
n=0；
nrow=0；ncol=0；
nrow=x.nrow（）；
ncol=x.ncol（）；
对于（i=0；ibool以下定义：
NumericMatrix dupMat(300, ncol);
ncolxnrowdupMatncolnrownrow = x.nrow();
ncol = x.ncol();
Rcpp::NumericMatrix dupMat(nrow, ncol);
<> >
flagtruetruefalsefalse有三种方法可以访问
numerimatrixi，j
：以这种方式访问元素将放弃边界检查和后续异常标志，如果该点不在范围内，则会发出警告。本质上，这种访问方法会导致UB，因为（i，j）
很容易超出行索引。可能是UB在稍后RStudio或R运行后台任务时造成了严重破坏，导致崩溃发生n=n+1
：这是首选方法，因为它提供边界检查并抛出一个漂亮的异常，例如.at（i，j）
if (flag == true) {
for (m = 0; m < ncol; ++m) {
Rcpp::Rcout << "dupMat (" << n << ","<< m << ")" << std::endl <<
"x (" << i << ","<< m << ")" << std::endl;
dupMat.at(n, m) = x.at(i, m);
}
n = n + 1; // able to exceed nrow.
}
if（标志==true）{
对于（m=0；mbool
，然后详细说明矩阵子集导致的未定义行为（UB）问题
以下定义：
NumericMatrix dupMat(300, ncol);
有两个问题：
它位于初始化ncol
之前
假设x
矩阵的nrow
固定为300
将dupMat
的实例化移动到初始化ncol
和nrow
之后。或者，将其移动到知道重复行的数量之后
nrow = x.nrow();
ncol = x.ncol();
Rcpp::NumericMatrix dupMat(nrow, ncol);
<> > BoOL C++中的值是用小写的。
也就是说，在设置flag
变量的值时，使用true
代替true
和false
代替false
有三种方法可以访问numerimatrix
中的单个元素，但是，我们将只关注其中两种使用i，j
索引的方法
（i，j）
：以这种方式访问元素会放弃边界检查和随后的异常标志，如果点不在范围内则会发出警告。本质上，这种访问方法会导致UB，因为n=n+1
很容易超出行索引。当RStudio或R运行backgr时，UB可能会在稍后点造成严重破坏导致崩溃发生的错误任务.at（i，j）
：这是首选方法，因为它提供边界检查并抛出一个漂亮的异常，例如
dupCheckRcpp（a）中出错：索引超出范围
由以下代码段触发：
if (flag == true) {
for (m = 0; m < ncol; ++m) {
Rcpp::Rcout << "dupMat (" << n << ","<< m << ")" << std::endl <<
"x (" << i << ","<< m << ")" << std::endl;
dupMat.at(n, m) = x.at(i, m);
}
n = n + 1; // able to exceed nrow.
}
if（标志==true）{
对于（m=0；m