noC <- function(x) {
  getComment(x) <- NULL
  return(x)
}

test_that("ifelse method works correctly", {
  p <- maxample("pop")
  p[1,1,1] <- Inf
  p[2,1,1] <- -Inf
  p2 <- p[10:1,,]
  expect_identical(noC(ifelse(p2 > 1000, p, 0)), noC(ifelse(p2 > 1000, p2, 0))) 
  p[c(1,4,6)] <- NA
  p2 <- p[10:1,,]
  expect_identical(noC(ifelse(!is.na(p2) & !is.infinite(p2), p, 0)), noC(ifelse(!is.na(p2)  & !is.infinite(p2), p2, 0))) 
  expect_error(ifelse(p,1,9), "must only contain booleans")
})