rcpp 代码示例 cppFunction

Question

我正在尝试运行“无缝 R 和 C++ 与 Rcpp 集成”的代码（第 32 页，清单 2.10），但它给出了一个错误。有人可以向我解释为什么不工作吗？谢谢

Code <- ' 
#include <gsL/gsl_const_mksa.h>           // decl of constants 
std::vector<double> volumes() { 
std::vector<double> v(5); 
v[0] = GSL_CONST_MKSA_US_GALLON;       // 1 US gallon 
v[1] = GSL_CONST_MKSA_CANADIAN_GALLON; // 1 Canadian gallon 
v[2] = GSL_CONST_MKSA_UK_GALLON;       // 1 UK gallon 
v[3] = GSL_CONST_MKSA_QUART;           // 1 quart 
v[4] = GSL_CONST_MKSA_PINT;            // 1 pint 
return v; 
}' 

gslVolumes <- cppFunction(code, depends="RcppGSL") 

这是消息错误：

file16e2b6cb966.cpp: In function ‘SEXPREC* sourceCpp_52966_volumes()’: 
file16e2b6cb966.cpp:30: error: ‘__result’ was not declared in this scope 
make: *** [file16e2b6cb966.o] Error 1 
llvm-g++-4.2 -arch x86_64 -I/Library/Frameworks/R.framework/Resources/include -DNDEBUG -I/usr/local/include -I/Library/Frameworks/R.framework/Versions/3.0/Resources/library/Rcpp/include -I/usr/local/include  -I"/Library/Frameworks/R.framework/Versions/3.0/Resources/library/Rcpp/include" -I"/Library/Frameworks/R.framework/Versions/3.0/Resources/library/RcppGSL/include"    -fPIC  -mtune=core2 -g -O2  -c file16e2b6cb966.cpp -o file16e2b6cb966.o 
Erro em sourceCpp(code = code, env = env, rebuild = rebuild, showOutput = showOutput,  : 
  Error 1 occurred building shared library. 

Answer 1

你好像有错别字：

Code <- ' 
#include <gsL/gsl_const_mksa.h>           // decl of constants 

应该是code <- 和小写c，然后是#include <gsl/gsl_const_mksa.h> 和小写“ell”。

一般来说，我建议开启详细模式：

gslVolumes <- cppFunction(code, depends="RcppGSL", verbose=TRUE) 

这会告诉你的

1. object code not found 来自第一个错误，并且

2. file....cpp:10:63: fatal error: gsL/gsl_const_mksa.h: No such file or directory

关于缺少的标题。

但我现在确实看到，在当前版本中，我也得到了__result not declared。将要 调查。

编辑：这是一个错误/更改。我写章节的时候就成功了，现在你需要

1. 从code 分配中删除带有#include <gsl/gsl_const_mksa.h> 的行

2. 向cppFunction() 调用添加一个新的includes=... 参数，如下所示：

正确的调用：

 gslVolumes <- cppFunction(code, depends="RcppGSL",
                           includes="#include <gsl/gsl_const_mksa.h>")

Answer 2

除了 Dirk 所说的，我建议您将代码提升为 .cpp 文件。

// [[Rcpp::depends(RcppGSL)]]
#include <Rcpp.h>
#include <gsl/gsl_const_mksa.h>           // decl of constants 

// [[Rcpp::export]]
std::vector<double> volumes() { 
  std::vector<double> v(5); 
  v[0] = GSL_CONST_MKSA_US_GALLON;       // 1 US gallon 
  v[1] = GSL_CONST_MKSA_CANADIAN_GALLON; // 1 Canadian gallon 
  v[2] = GSL_CONST_MKSA_UK_GALLON;       // 1 UK gallon 
  v[3] = GSL_CONST_MKSA_QUART;           // 1 quart 
  v[4] = GSL_CONST_MKSA_PINT;            // 1 pint 
  return v; 
}   

然后，您可以sourceCpp 该文件。