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
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
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
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
__
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:
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