test_that("gsDesign2:::h1() returns results as expected ", {
  # the design
  gstry <- gsDesign::gsDesign(
    k = 3,
    sfl = gsDesign::sfLDOF,
    delta = 0
  )
  # probabilities calculated based on function gsDesign2:::h1()
  upper.null <- sum(gsDesign2:::h1(
    theta = gstry$theta[1],
    info = gstry$n.I[1],
    a = gstry$upper$bound[1],
    b = Inf
  )$h)
  upper.alt <- sum(gsDesign2:::h1(
    theta = gstry$theta[2],
    info = gstry$n.I[1],
    a = gstry$upper$bound[1],
    b = Inf
  )$h)
  lower.null <- sum(gsDesign2:::h1(
    theta = gstry$theta[1],
    info = gstry$n.I[1],
    a = -Inf,
    b = gstry$lower$bound[1]
  )$h)
  lower.alt <- sum(gsDesign2:::h1(
    theta = gstry$theta[2],
    info = gstry$n.I[1],
    a = -Inf,
    b = gstry$lower$bound[1]
  )$h)
  # probabilities calculated based on function gsProbability
  x <- gsDesign::gsProbability(
    k = 3,
    a = gstry$lower$bound,
    b = gstry$upper$bound,
    n.I = gstry$n.I, theta = gstry$theta
  )
  expect_equal(object = as.numeric(c(upper.null, upper.alt)), expected = x$upper$prob[1, ], tolerance = 0.0001)
  expect_equal(object = as.numeric(c(lower.null, lower.alt)), expected = x$lower$prob[1, ], tolerance = 0.0001)
})

test_that("gsDesign2:::h1() returns probability almost zero for extreme case", {
  exmtest1 <- sum(gsDesign2:::h1(theta = 9, info = 0.5, a = -Inf, b = 0)$h)
  exmtest2 <- sum(gsDesign2:::h1(theta = 1, info = 0.5, a = 9, b = Inf)$h)
  expect_equal(object = as.numeric(c(exmtest1, exmtest2)), expected = c(0, 0), tolerance = 0.0001)
})