tests/testthat/test-input-format.R

In geodist: Fast, Dependency-Free Geodesic Distance Calculations

# tibble not tested to avoid pacakge suggests, but just needs these lines
#x <- tibble::tibble (x = -180 + 360 * runif (n),
#                     y = -90 + 180 * runif (n))
#y <- tibble::tibble (x = -180 + 360 * runif (2 * n),
#                     y = -90 + 180 * runif (2 * n))

test_all <- identical (Sys.getenv ("MPADGE_LOCAL"), "true")

test_that("geodist with df", {
              n <- 1e2
              x <- data.frame (x = -180 + 360 * runif (n),
                               y = -90 + 180 * runif (n))
              y <- data.frame (x = -180 + 360 * runif (2 * n),
                               y = -90 + 180 * runif (2 * n))
              d1 <- geodist (x)
              expect_equal (dim (d1), c (n, n))
              diag (d1) <- Inf
              expect_true (all (d1 >= 0))

              d2 <- geodist (x, y)
              expect_equal (dim (d2), c (n, 2 * n))
              diag (d2) <- Inf
              expect_true (all (d2 >= 0))
})

test_that("geodist with matrix", {
              n <- 1e2
              x <- cbind (-180 + 360 * runif (n), -90 + 180 * runif (n))
              y <- cbind (-180 + 360 * runif (2 * n), -90 + 180 * runif (2 * n))
              expect_message (d1 <- geodist (x), "object has no named columns")
              colnames (x) <- c ("x", "y")
              if (test_all)
                  expect_message (d1 <- geodist (x),
                                  "Maximum distance is > 100km")
              else
                  d1 <- geodist (x)
              expect_equal (dim (d1), c (n, n))
              diag (d1) <- Inf
              expect_true (all (d1 >= 0))

              expect_message (d2 <- geodist (x, y),
                              "object has no named columns")
              expect_equal (dim (d2), c (n, 2 * n))
              diag (d2) <- Inf
              expect_true (all (d2 > 0))
              expect_true (all (d2 >= 0))
})

test_that ("other columns", {
               n <- 50
               x <- cbind (-10 + 20 * runif (n), -10 + 20 * runif (n))
               y <- cbind (-10 + 20 * runif (2 * n), -10 + 20 * runif (2 * n))
               colnames (x) <- colnames (y) <- c ("x", "y")
               d <- geodist (x, y)
               x1 <- cbind (newx = runif (n), newy = runif (n), x)
               y1 <- cbind (newx = runif (n), newy = runif (n), y)
               d1 <- geodist (x1, y1)
               expect_identical (d, d1)
})

test_that ("column names", {
               n <- 50
               x <- cbind (-10 + 20 * runif (n), -10 + 20 * runif (n))
               y <- cbind (-10 + 20 * runif (n), -10 + 20 * runif (n))
               colnames (x) <- colnames (y) <- c ("x", "x")
               expect_error (d <- geodist (x, y, paired = TRUE),
                             paste0 ("Unable to determine longitude and ",
                                     "latitude columns"))
               colnames (x) <- colnames (y) <- c ("_x_x_", "_x_y_")
               expect_error (d <- geodist (x, y, paired = TRUE),
                             paste0 ("Unable to determine longitude and ",
                                     "latitude columns"))
               colnames (x) <- colnames (y) <- c ("_x_x", "_x_y")
               if (test_all)
                   expect_message (d1 <- geodist (x, y, paired = TRUE),
                                  "Maximum distance is > 100km")
               else
                   d1 <- geodist (x, y, paired = TRUE)

               colnames (x) <- colnames (y) <- c ("_a_x_", "_a_y_")
               if (test_all)
                   expect_message (d2 <- geodist (x, y, paired = TRUE),
                                   "Maximum distance is > 100km")
               else
                   d2 <- geodist (x, y, paired = TRUE)
               expect_identical (d1, d2)
})