tests/testthat/test-transformSeedPoint.R

In hitandrun: "Hit and Run" and "Shake and Bake" for Sampling Uniformly from Convex Shapes

context("seedPoint argument to har.init()")

test_that("it is transformed correctly", {
  constr <- list(
    constr=rbind(rep(1, 3), -diag(3)),
    rhs=c(1, rep(0, 3)),
    dir=c("=", rep("<=", 3)))
  x0 <- c(1/2,1/4,1/4)

  state <- har.init(constr, x0=x0)

  expect_equal(as.vector(state$transform %*% state$x0), x0)

  constr$dir <- rep("<=", 4)
  state <- har.init(constr, x0=x0)
  expect_equal(as.vector(state$x0), c(x0, 1))
  expect_equal(as.vector(state$transform %*% state$x0), x0)
})