Re: [Rcpp-devel] Integrals with Rcpp

2015-01-09 Thread Federico Andreis
Dear Dirk, thanks for your answer, I'll keep studying and decide what way to follow. Dear Qiang, here's the code for the function (I already did some work to the math behind the density itself in order to quicken a bit) library(cubature) library(zipfR) dcbinom <- Vectorize(function(x,nu,pp) { #

Re: [Rcpp-devel] Integrals with Rcpp

2015-01-09 Thread Qiang Kou
Hi, The backend of cubature package is the cubature C library from MIT. The same author with nlopt. Can you give a little more detail about the double integral? Cubature should not be that slow from my experience. Best, KK On Fri, Jan 9, 2015 at 9:38 AM, Dirk Eddelbuettel wrote: > > On 9 J

Re: [Rcpp-devel] Integrals with Rcpp

2015-01-09 Thread Dirk Eddelbuettel
On 9 January 2015 at 15:10, Federico Andreis wrote: | Dear all, | | I was wondering what, in your opinion, would be the best way to handle the | computation of definite integrals using Rcpp. | | To be more precise, I need to evaluate a double integral that turns out to be | somewhat nasty (reall

[Rcpp-devel] Integrals with Rcpp

2015-01-09 Thread Federico Andreis
Dear all, I was wondering what, in your opinion, would be the best way to handle the computation of definite integrals using Rcpp. To be more precise, I need to evaluate a double integral that turns out to be somewhat nasty (really slow computation and presumably inaccurate results using the R cu

Re: [Rcpp-devel] Setting time zone on Rcpp::Datetime

2015-01-09 Thread Dirk Eddelbuettel
On 9 January 2015 at 11:50, janus Larsen wrote: | Hi, | How do I set the timezone on an Rcpp::Datetime? | Thanks in advance, | Sunaj | | This returns "1970-01-01 01:00:00 CET" (computer setting), but I want GMT... R does the formatting in its session based on your locale. That is generally the

[Rcpp-devel] Setting time zone on Rcpp::Datetime

2015-01-09 Thread janus Larsen
Hi, How do I set the timezone on an Rcpp::Datetime? Thanks in advance, Sunaj This returns "1970-01-01 01:00:00 CET" (computer setting), but I want GMT... // [[Rcpp::export]] Rcpp::Datetime test() { double d=0; return(d); } ___ Rcpp-devel mailing list