Skip to content

Hypergeometric function of a matrix argument (Koev & Edelman's algorithm)

Notifications You must be signed in to change notification settings

stla/HypergeoMat

Repository files navigation

HypergeoMat

R-CMD-check

library(HypergeoMat)
library(microbenchmark)
jhpq <- hypergeomPFQ_julia()
## Starting Julia ...
microbenchmark(
  Rcpp = hypergeomPFQ(m = 30, a = c(1, 2, 3), b = c(4, 5), x = c(0.1, 0.2, 0.3)),
  Julia = jhpq(m = 30, a = c(1, 2, 3), b = c(4, 5), x = c(0.1, 0.2, 0.3)),
  times = 10
)
## Unit: milliseconds
##   expr       min        lq      mean    median       uq      max neval
##   Rcpp 3171.9094 3266.6886 3386.5241 3317.9526 3531.426 3659.559    10
##  Julia  339.2297  355.8744  755.2628  377.9874  438.013 3866.369    10