[Rcpp-devel] Rf_integrate?

[Silkworth,David J.]

I have a current project desire to move something like R's integrate function inside a loop in Rcpp code.  The performance hit to call back to R seems to kill the advantage of Rcpp in the first place.  Actually my integrand is t*pdf(t), very similar indeed to pweibull which integrates pdf(t).  It has been hard to find previous discussion on this topic since the title of articles and book include the word Integration in  another context.

I realize that eventually R's integrate function calls rdqags (over a definite interval), but there is a lot of memory management that is taken care of before the call.  This is over my head.

I could try to incorporate the GSL, but this too seems daunting (even with RcppGSL).  I think there may be integration support in the boost headers, but my head is too small for this yet.

Any ideas that could help me?

Dave Silkworth

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/rcpp-devel/attachments/20140730/68d627f5/attachment.html>