आरसीपी

Question

में अन्य पैकेज से सी फ़ंक्शन का उपयोग करके मैं बहु-आयामी एकीकरण करने के लिए एक C++ फ़ंक्शन में क्यूबचर पैकेज से सी सी रूटीन कॉल करने का प्रयास कर रहा हूं।

बुनियादी अनुसंधान उदाहरण मैं पुन: पेश करने की कोशिश कर रहा हूँ

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

है मैं सिर्फ निम्नलिखित this recipe from the gallery Rcpp से आर समारोह कह सकते हैं, लेकिन आगे और पीछे से सी/सी स्विच करने में कुछ प्रदर्शन जुर्माना होगा ++ आर के लिए सी सी + से सीधे सी समारोह को कॉल करने के लिए यह अधिक समझदार लगता है।

सी दिनचर्या adapt_integrate

// 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 में परिवर्तित करने के लिए कुछ महत्वपूर्ण कदम उठा रहा हूं (या इसे सीधे कॉल करें?)। मैंने फ़ंक्शन पॉइंटर्स से निपटने वाली कई संबंधित पोस्ट पढ़ी हैं, लेकिन बाहरी सी फ़ंक्शन का उपयोग करके कोई उदाहरण नहीं देखा है।

Answer 1

दुर्भाग्य cubatureinst/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 हो सकता है cubature कि वह inst/include में घोषणाएं भेजता है ताकि आपको केवल LinkingTo का उपयोग करना पड़े।

Answer 2

पहले इस प्रश्न को नहीं देखा था, और ऐसा लगता है कि @Romain ने इसे संबोधित किया।

पूर्णता के लिए, यह सुनिश्चित करने के लिए एक कार्य उदाहरण है कि जब सभी पार्टियां खेलती हैं तो xts और RcppXts पैकेज द्वारा प्रदान की जाती है। R_registerRoutines और R_RegisterCCallable के सामान्य व्यापार के साथ साथ

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

: xts में, हम ऐसा करते हैं (स्रोत) में (लगभग दस कार्यों के लिए) फ़ाइल inst/include/xtsAPI.h।

RcppXts में इस रूप में

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

जो बहुत अच्छी तरह से काम करता है (एक Rcpp मॉड्यूल में) उठाया है। किसी ने मुझे xts पक्ष को अधिक कॉम्पैक्टली लिखने के लिए झगड़ा दिया (जैसा कि if NULL नकली है) जो मैं ... अंततः प्राप्त करूंगा।

Answer 3

यह सवाल तीन साल अब पुरानी है, लेकिन मैं बाहर बात करने के लिए अब Rcpp साथ बहुआयामी एकीकरण आसान हो सकता है कि कि RcppNumerical पुस्तकालय उपलब्ध है चाहता हूँ: https://github.com/yixuan/RcppNumerical

कंप्यूटिंग अभिन्न के लिए दिनचर्या थॉमस हैन के क्यूबा पर आधारित होते हैं पैकेज और सीआरएएन पर आर 2 सी क्यूबा लाइब्रेरी में भी उपलब्ध हैं, इसलिए यदि आप क्यूबचर से adaptIntegrate() फ़ंक्शन पर क्यूबा रूटीन का उपयोग कर स्वीकार कर सकते हैं, तो यह पैकेज ब्याज का हो सकता है।