参考资料 |大数据意外崩溃

Question

我正在尝试运行以下逐行递归操作以获得norm_vec：

#include <Rcpp.h>
#ifdef _OPENMP
#include <omp.h> // OpenMP header
#endif

using namespace Rcpp;
using namespace std;

// [[Rcpp::export]]
NumericVector sim1(NumericVector X) {

  int T = X.length();
  NumericVector norm_vec = NumericVector(T);

  norm_vec[0] = 0;

  int i = 1;

  while (i <= T) {

    NumericVector sim_vec = NumericVector(i);
    NumericVector weight_vec = NumericVector(i);
    NumericVector norm = NumericVector(i);

    int j = 0;
    double sim_vec_sum = 0;

    while (j <= i) {

      sim_vec[j] = exp(-abs(X[i] - X[i-j-1]));
      sim_vec_sum += sim_vec[j];
      j++;

    }

    j = 0;
    double norm_sum = 0;

    while (j <= i) {

      weight_vec[j] = sim_vec[j]/sim_vec_sum;
      norm[j] = X[j]*weight_vec[j];
      norm_sum += norm[j];
      j++;

    }

    norm_vec[i] = norm_sum;

    i++;

  }

  return norm_vec;

}

用 sourceCpp("sim1.cpp") 调用上面的代码，我想得到 sim1(rnorm(n))。虽然代码在较小的 n 上工作得很好，但一旦 n 超过 19，R 就会完全关闭。

我对 cpp 不太满意 - 可能我犯了一个基本错误。非常感谢您的提前帮助！

Answer 1

R 中的循环从 1 到 T，但在 C/C++ 中它们从 0 到 T-1。

在您的循环中 i 上升到 T，并且 j 上升到 i。这样i可以达到T，比T-1.

的最大值还多

参考资料 |大数据意外崩溃

Rcpp | unexpected crashes with large data

c++

r

rcpp