test_that("kerPIQ works", {
  expect_equal(kerPIQ(1,c(1,2),c(2,1)), 1/4)
})