我正在尝试将函数用作另一个函数的参数。
我正在构建一个包。
func1.cpp:
#include <Rcpp.h>
using namespace Rcpp;
// [[Rcpp::export]]
NumericVector timesTwo(NumericVector x) {
return x * 2;
}
main.cpp中:
#include <Rcpp.h>
#include <cmath>
#include <cstdio>
#include <cstdlib>
// [[Rcpp::plugins(cpp11)]]
using namespace Rcpp;
NumericVector timesTwo(NumericVector x);
typedef NumericVector (*timesTwo_call)(NumericVector x);
// [[Rcpp::export]]
XPtr<timesTwo_call> putFunPtrInXPtr(std::string fstr) {
if (fstr == "timesTwo")
return(XPtr<timesTwo_call>(new timesTwo_call(×Two)));
else
return XPtr<timesTwo_call>(R_NilValue); // runtime error as NULL no XPtr
}
// [[Rcpp::export]]
NumericVector testFunc(SEXP func, NumericVector x)
{
XPtr<timesTwo_call> xpfun(func);
timesTwo_call fun = *xpfun;
NumericVector tmp = fun(x);
const int N = tmp.size();
NumericVector result(N);
for (int i = 0; i < N; ++i)
{
result[i] = tmp[i] * 3;
}
return result;
}
/*** R
x <- c(1,2,3)
fun <- putFunPtrInXPtr("timesTwo")
result <- testFunc(fun, x)
*/
尝试构建,它在RcppExports.cpp中给出了一些错误,在行:
XPtr<timesTwo_call> putFunPtrInXPtr(std::string fstr);
错误:
timesTwo_call was not declared in this scope
template argument 1 is invalid
template argument 3 is invalid
invalid type in declaration before ' token
- 更新 -
我应该注意,即使我将timesTwo
的定义放在主文件中,我仍然会遇到相同的错误。
但是,如果(当我将定义放在main中时)作为sourceCpp运行,那么它可以工作!
- 更新2 ---
我在src文件夹下创建了文件test_types.h
(其中所有cpp文件都是)。
#ifndef TEST_TYPES_H
#define TEST_TYPES_H
#include <Rcpp.h>
using namespace Rcpp;
typedef NumericVector (*timesTwo_call)(NumericVector x);
NumericVector timesTwo(NumericVector x);
#endif
但我仍然遇到同样的错误。
答案 0 :(得分:3)
问题是您的typedef
未导出到RcppExports.cpp
。
请参阅 2.4部分。 Rcpp attribute vignette的生成代码类型。简而言之:将typedef放入名为:
src/<pkg>_types.h
src/<pkg>_types.hpp
inst/include/<pkg>_types.h
inst/include/<pkg>_types.hpp
inst/include/<pkg>.h
此文件将自动包含在RcppExports.cpp
中。您必须手动将其包含在main.cpp
。