tests/testthat/test-errors.R

In stla/HypergeoMat3: Hypergeometric Function of a Matrix Argument

context("Errors")

test_that("errors on m", {
  expect_error(hypergeomPFQ(m = -1, a = c(1,2), b = c(3,4), x = c(5,6)))
  expect_error(hypergeomPFQ(m = c(1,2), a = c(1,2), b = c(3,4), x = c(5,6)))
})

test_that("errors on alpha", {
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = c(5,6),
                            alpha = "a"))
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = c(5,6),
                            alpha = -1))
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = c(5,6),
                            alpha = c(1,2)))
})

test_that("errors on x", {
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4),
                            x = cbind(c(1,1),c(2,3),c(4,5))))
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = "a"))
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = list(5,6)))
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = NULL))
  expect_error(hypergeomPFQ(m = 5, a = c(1,2), b = c(3,4), x = numeric(0L)))
})

test_that("errors lmvgamma", {
  expect_error(lmvgamma(-1+2i, 3))
  expect_error(lmvgamma(1+2i, c(2,3)))
  expect_error(lmvgamma(c(3,4), 3))
})

test_that("errors mvgamma", {
  expect_error(mvgamma(c(-1+2i,3), 3))
  expect_error(mvgamma(1+2i, c(2,3)))
  expect_error(mvgamma(-2, 3))
})

test_that("errors mvbeta", {
  expect_error(mvbeta(a = c(1,2), b = 3, p = 4))
  expect_error(mvbeta(a = 2, b = 3, p = c(3,4)))
})

test_that("errors IncBeta", {
  expect_error(IncBeta(m = 5, a = 0.5, b = 2, x = c(1,2,3)))
  expect_error(IncBeta(m = 5, a = c(1,5), b = 2, x = c(1,2,3)))
  expect_error(IncBeta(m = 5, a = 5, b = 2, x = diag(c(0.5,-0.5))))
  expect_error(IncBeta(m = 5, a = 5, b = 2, x = diag(c(1.1,0.5))))
})

test_that("errors IncGamma", {
  expect_error(IncGamma(m = 5, a = 1, x = toeplitz(3:1)))
})

test_that("errors BesselA", {
  expect_error(BesselA(m = 5, x = diag(c(1,2+2i)), nu = -1+1i))
})