test_that("Contourmap, test ind", {
  data <- integration.testdata1()
  ind1 <- c(1, 2, 3, 4)
  ind2 <- c(4, 3, 2, 1)
  ind3 <- rep(FALSE, data$n)
  ind3[1:4] <- TRUE

  res1 <- contourmap(data$mu, data$Q,
    n.levels = 2, ind = ind1,
    seed = data$seed, alpha = 0.1, max.threads = 1
  )
  res2 <- contourmap(data$mu, data$Q,
    n.levels = 2, ind = ind2,
    seed = data$seed, alpha = 0.1, max.threads = 1
  )
  res3 <- contourmap(data$mu, data$Q,
    n.levels = 2, ind = ind3,
    seed = data$seed, alpha = 0.1, max.threads = 1
  )

  expect_equal(res1$F, res2$F, tolerance = 1e-7)
  expect_equal(res2$F, res3$F, tolerance = 1e-7)
})


test_that("Contourmap, P measures", {
  data <- integration.testdata1()

  res1 <- contourmap(data$mu, data$Q,
    n.levels = 4,
    seed = data$seed, alpha = 0.1, max.threads = 1,
    compute = list(F = FALSE, measures = c("P2", "P1"))
  )

  expect_equal(res1$P1, 0.9217417, tolerance = 1e-3)
  expect_equal(res1$P2, 0.405841, tolerance = 1e-3)
})