tests/testthat/test-tsearch-tsearchn-comparison.R

In geometry: Mesh Generation and Surface Tesselation

context("Comparison of tsearch and tsearchn")
test_that("tsearch and tsearchn give the same results", {
  set.seed(1)
  X <- runif(50)
  Y <- runif(50)
  T <- delaunayn(cbind(X, Y))
  XI <- runif(1000)
  YI <- runif(1000)

  out <- tsearch(X, Y, T, XI, YI)
  outn <- tsearchn(cbind(X, Y), T, cbind(XI, YI), fast=FALSE)
  
  expect_that(na.omit(out), equals(na.omit(outn$idx)))

  out <- tsearch(X, Y, T, XI, YI, TRUE)
  expect_that(na.omit(outn$p), equals(na.omit(out$p), tolerance=1e-12))
})