Re: [Rcpp-devel] Parallel R::bessel_k function in Rcpp

2016-03-23 Thread Baptiste Auguie
Hi, Check also the Bessel package from Martin M., https://cran.r-project.org/web/packages/Bessel/index.html It wraps some C code; maybe a C++ interface wouldn't be too hard to write. There's also the GSL package apparently. Cheers, baptiste On 24 March 2016 at 13:20, Hoang Nguyen wrote: > Tha

Re: [Rcpp-devel] Parallel R::bessel_k function in Rcpp

2016-03-23 Thread Hoang Nguyen
Hi, I found in wiki that Modified Bessel Functions of the Second Kinds has different names and one of it is the third kind. Sorry for this trouble. Thank you. Best. On Thu, Mar 24, 2016 at 1:20 AM, Hoang Nguyen wrote: > Thank Dirk, > I check Boost and they only have the Modified Bessel Functions

Re: [Rcpp-devel] Parallel R::bessel_k function in Rcpp

2016-03-23 Thread Hoang Nguyen
Thank Dirk, I check Boost and they only have the Modified Bessel Functions of the First and Second Kinds. I have to read more to find how to derive the third kind. Best regards. On Thu, Mar 24, 2016 at 1:17 AM, Dirk Eddelbuettel wrote: > > PS Looks like Boost has Bessel functions: > > > http://w

Re: [Rcpp-devel] Parallel R::bessel_k function in Rcpp

2016-03-23 Thread Dirk Eddelbuettel
PS Looks like Boost has Bessel functions: http://www.boost.org/doc/libs/1_60_0/libs/math/doc/html/math_toolkit/bessel/bessel_over.html So maybe try that. We should have access to that in the BH package. Dirk -- http://dirk.eddelbuettel.com | @eddelbuettel | e...@debian.org __

Re: [Rcpp-devel] Parallel R::bessel_k function in Rcpp

2016-03-23 Thread Dirk Eddelbuettel
On 24 March 2016 at 01:02, Hoang Nguyen wrote: | Hi, | I am working with Rcpp and need to call function bessel_k from R::bessel_k in | parallel. I got usually crash in Rstudio, and still not find out which is the | mistake You _cannot call back to R_ from parallel code. | Here is my code in R:

[Rcpp-devel] Parallel R::bessel_k function in Rcpp

2016-03-23 Thread Hoang Nguyen
Hi, I am working with Rcpp and need to call function bessel_k from *R::bessel_k* in parallel. I got usually crash in Rstudio, and still not find out which is the mistake Here is my code in R: Sys.setenv("PKG_CXXFLAGS"="-fopenmp") Sys.setenv("PKG_LIBS"="-fopenmp") library(Rcpp) sourceCpp("Bessel.c