tests/testthat/test-mvgamma.R

In stla/HypergeoMat3: Hypergeometric Function of a Matrix Argument

context("Multivariate Gamma function")

test_that("Comparison with gamma", {
  expect_equal(mvgamma(-1.5,1), gamma(-1.5))
  expect_equal(mvgamma(-2.5,1), gamma(-2.5))
})

test_that("Complex value for p=1", {
  y <- 2
  o1 <- Mod(mvgamma(y*1i,1))
  o2 <- sqrt(pi/y/sinh(pi*y))
  expect_equal(o1, o2)
})