tests/testthat/test-numerical-comparisons.R

In PredictiveEcology/fpCompare: Reliable Comparison of Floating Point Numbers

test_that("relational operators within tolerance", {
  x <- 0.5 - 0.3
  y <- 0.3 - 0.1
  expect_equal(x %==% y, TRUE)
  expect_equal(x %!=% y, FALSE)

  set.seed(123L)
  a <- jitter(1:10, 1e-3)
  b <- jitter(1:10, 1e-3)
  lessEq    <- a %<=% b
  less      <- a %<<% b
  greaterEq <- a %>=% b
  greater   <- a %>>% b
  equal     <- a %==% b
  notEqual  <- a %!=% b
  expect_equal(less,     !greater)
  expect_equal(notEqual, !equal)

  a <- jitter(1:10, 1e-7)
  b <- jitter(1:10, 1e-7)
  lessEq     <- a %<=% b
  less       <- a %<<% b
  greaterEq  <- a %>=% b
  greater    <- a %>>% b
  equal      <- a %==% b
  notEqual   <- a %!=% b
  ids <- c(5L, 8L)
  expect_equal(less[ids],  !greater[ids])
  expect_equal(less[-ids], greater[-ids])
  expect_equal(lessEq,     !greater)
  expect_equal(all(!equal[ids]), TRUE)
  expect_equal(all(equal[-ids]), TRUE)
})