Rcpp内の他のパッケージのC関数を使用

Question

多次元統合を実行するためにC++関数のcubatureパッケージからCルーチンを呼び出そうとしています。

私が再現しようとしている基本的なRの例は、

library(cubature) 
integrand <- function(x) sin(x) 
adaptIntegrate(integrand, 0, pi) 

である私はちょうどthis recipe from the gallery次RcppからこのR関数を呼び出すことができますが、C/Cから前後に切り替えるにはいくつかのパフォーマンスのペナルティが存在することになる++ C++からC関数を直接呼び出す方が合理的です。

adapt_integrateルーチンCは、私はしかし、C++からそれを呼び出す方法を理解していない

// R_RegisterCCallable("cubature", "adapt_integrate", (DL_FUNC) adapt_integrate); 

でcubatureからエクスポートされます。ここで私の不自由な試みは、

sourceCpp(code = ' 
#include <Rcpp.h> 
using namespace Rcpp; 

// [[Rcpp::export]] 
double integrand(double x){ 
return(sin(x)); 
} 

// [[Rcpp::depends(cubature)]] 
// [[Rcpp::export]] 
Rcpp::List integratecpp(double llim, double ulim) 
{ 
    Rcpp::Function p_cubature = R_GetCCallable("cubature", "adapt_integrate"); 

    Rcpp::List result = p_cubature(integrand, llim, ulim); 
    return(result); 
} 
' 
) 

integratecpp(0, pi) 

これはコンパイルに失敗します。明らかに私は非常にばかげた何かをして、R_GetCCallableの出力をRcpp::Functionに変換する（またはそれを直接呼びますか？）いくつかの重要なステップがありません。関数ポインタを扱ういくつかの関連記事を読んだが、外部のC関数を使った例は見ていない。

Answer 1

残念ながらcubatureはinst/includeのヘッダーを出荷していないので、あなたは彼らからそれを借りて、あなたのコード内でこのような何かをする必要があります：

typedef void (*integrand) (unsigned ndim, const double *x, void *, 
      unsigned fdim, double *fval); 

int adapt_integrate(
    unsigned fdim, integrand f, void *fdata, 
    unsigned dim, const double *xmin, const double *xmax, 
    unsigned maxEval, double reqAbsError, double reqRelError, 
    double *val, double *err) 
{ 
    typedef int (*Fun)(unsigned,integrand,void*,unsigned, 
     const double*,const double*, unsigned, double, double, double*, double*) ; 
    Fun fun = (Fun) R_GetCCallable("cubature", "adapt_integrate") ;   
    return fun(fdim,f,fdata,dim,xmin,xmax,maxEval,reqAbsError, reqRelError,val,err); 
} 

のメンテナでnegociateするのは良い考えかもしれませんcubatureinst/includeに宣言が含まれているので、LinkingToを使用するだけです。

Answer 2

この質問は先に表示されていませんでした。@Romainが解決したようです。

完全性のために、すべての当事者が一緒に遊んでいるときにこれを行う方法の実例は、xtsおよびRcppXtsパッケージによって提供されています。 xtsでは、我々はこれを行う（ソース）中（約10機能のための）inst/include/xtsAPI.hファイル：

SEXP attribute_hidden xtsLag(SEXP x, SEXP k, SEXP pad) {  
    static SEXP(*fun)(SEXP,SEXP,SEXP) = NULL;   
    if (fun == NULL)         
     fun = (SEXP(*)(SEXP,SEXP,SEXP)) R_GetCCallable("xts","lagXts"); 
    return fun(x, k, pad);        
} 

をR_registerRoutinesとR_RegisterCCallableの通常のビジネスと一緒に。

はRcppXtsでは、これはかなりうまく機能

function("xtsLag", 
     &xtsLag,  
     List::create(Named("x"), Named("k"), Named("pad")), 
     "Extract the coredata from xts object"); 

として（Rcppモジュールで）ピックアップされます。誰かが私に返信して、xtsをもっとコンパクトに書いてください（if NULLは偽りです）。

Answer 3

この質問は今3歳ですが、私はRcppNumericalライブラリが利用可能であることを今Rcppと多次元統合が容易になることを指摘したい： https://github.com/yixuan/RcppNumerical

積分を計算するためのルーチンはトーマス・ハーンのキューバに基づいていますCRANのR2Cubaライブラリにもありますので、Cubatureの関数adaptIntegrate()でキューバのルーチンを使用することができれば、このパッケージに興味があります。