Rcpp

Question

diğer paketinden C işlevini kullanarak çok boyutlu birleştirme gerçekleştirmek için C++ işlevinde bir C yordamı çağırmak çalışıyorum.

ben yeniden oluşturmaya çalışıyorum temel Ar örneği sadece this recipe from the gallery aşağıdaki Rcpp bu R işlevi diyebiliriz ama ileri geri c/c geçiş bazı performans cezası olacağını

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

olduğunu ++ R'ye C++'dan direkt olarak C fonksiyonu çağırmak daha mantıklı görünüyor.

adapt_integrate rutin C ı ancak, C++ onu aramak için anlamıyorum

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

ile cubature ihraç edilmektedir. İşte bu benim topal girişimi,

Bu, derlenemiyor; Açıkçası çok aptalca bir şey yapıyorum ve R_GetCCallable çıktısını Rcpp::Function'a dönüştürmek için bazı önemli adımları atıyorum (veya doğrudan arayın?). İşlev işaretçileriyle ilgili birçok ilgili yazı okudum, ancak harici bir C işlevi kullanarak bir örnek görmedim.

Answer 1

Maalesef cubatureinst/include başlıkların birlikte gelmez, bu yüzden onlardan ödünç ve kodunuzda böyle bir şey yapmak zorunda:

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); 
} 

O sürdürücü ile negociate için iyi bir fikir olabilir cubature, inst/include numaralı belgede beyanları gönderir, böylece yalnızca LinkingTo'u kullanmanız gerekir.

Answer 2

Bu soruyu daha önce görmedik ve @Romain buna hitap ediyor gibi görünüyor.

Tüm taraflar birlikte çalışırken bunun nasıl yapılacağına dair bir çalışma örneği, tamlık açısından xts ve RcppXts paketleri tarafından sağlanır. xts, biz bunu (kaynak) (yaklaşık on fonksiyonlar için) inst/include/xtsAPI.h dosyası:

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 ve R_RegisterCCallable olağan iş ile birlikte.

RcppXts bu oldukça iyi çalışıyor

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

olarak (bir Rcpp Modülünde) alınır. Birisi bana xts tarafını daha sıkı bir şekilde yazmam için beni azarladı (if NULL sahte).

Answer 3

Bu soru şimdi üç yaşında ama RcppNumerical kütüphanesi kullanılabilir olduğunu şimdi Rcpp ile çok boyutlu entegrasyon daha kolay olabileceğini işaret etmek istiyorum: https://github.com/yixuan/RcppNumerical

integraller hesaplamak için rutinleri Thomas Hahn Küba dayanmaktadır paketi ve ayrıca CRAN'daki R2Cuba kütüphanesinde mevcuttur, bu yüzden Küba'dan adaptIntegrate() fonksiyonu üzerinden Küba rutinlerini kabul edebiliyorsanız, bu paket ilgi çekebilir.