Question

我正在尝试构建一个包含使用RCpp执行简单卷积的函数的包。代码看起来像

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector conv_filter(NumericVector x, NumericVector c){
  int T = x.size();
  int M = (c.size()-1)/2;
  int t, i;    
    NumericVector fx(T);
  for(t=0; t<T; t++){
        for(i=-M; i<M+1; i++){
            if(t+i>=0 && t+i<T){
        fx(t) += c(M+i)*x(t+i);
            }
        }       
    }
    return fx;
}

在源代码时工作正常，但是当内置到包中时，我不断收到错误，说“期望单个值”。我想我犯了非常基本的错误，但即使在阅读了相关主题后我也看不到它的来源。非常感谢你的帮助。

Answer 1

您的计划适合我：

R> Rcpp::sourceCpp("/tmp/hrcho.cpp")

R> conv_filter(1:4, 4:1)
[1]  7 16 25 24
R>

使用这个（只是缩进的）源代码：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector conv_filter(NumericVector x, NumericVector c){
  int T = x.size();
  int M = (c.size()-1)/2;
  int t, i;    
  NumericVector fx(T);
  for(t=0; t<T; t++){
    for(i=-M; i<M+1; i++){
      if(t+i>=0 && t+i<T){
        fx(t) += c(M+i)*x(t+i);
      }
    }       
  }
  return fx;
}

/*** R
conv_filter(1:4, 4:1)
*/

RCpp：期望单值错误

1 个答案: