test-pmax0.R
In hutilscpp: Miscellaneous Functions in C++

# test_that("init", {
expect_true(TRUE)
do_pmax0_radix_sorted_dbl <- hutilscpp:::do_pmax0_radix_sorted_dbl
do_pmax0_radix_sorted_int <- hutilscpp:::do_pmax0_radix_sorted_int
do_pmin0_radix_sorted_dbl <- hutilscpp:::do_pmin0_radix_sorted_dbl
do_pmin0_radix_sorted_int <- hutilscpp:::do_pmin0_radix_sorted_int

do_pmax0_bitwise <- hutilscpp:::do_pmax0_bitwise

hutilscpp_rev <- hutilscpp:::hutilscpp_rev
is_covr <- hutilscpp:::is_covr
firstNonNegativeRadix <- hutilscpp:::firstNonNegativeRadix

# test_that("pmax0 abs", {
expect_true(TRUE) # to get started
x <- c(-1, 0, 1, -1)
expect_identical(pmax.int(x, 0), pmax0(x))
# Need to test extrema
if (is.integer(-.Machine$integer.max)) {
  y <- c(-.Machine$integer.max, .Machine$integer.max, 0L)
  expect_identical(pmax0(y), c(0L, .Machine$integer.max, 0L))
}

# test_that("pmax0 radix", {
x <- as.double(seq(-1e6, 1e7, length.out = 3e3))
expect_identical(pmax0(x), do_pmax0_radix_sorted_dbl(x))
x <- hutilscpp_rev(x)
expect_identical(pmax0(x), do_pmax0_radix_sorted_dbl(x))
x <- hutilscpp_rev(x)
x <- as.integer(x)
expect_identical(pmax0(x), do_pmax0_radix_sorted_int(x))
x <- hutilscpp_rev(x)
expect_identical(pmax0(x), do_pmax0_radix_sorted_int(x))

x <- -5:6
x[1] <- -5L
pmax0(x, sorted = TRUE, in_place = TRUE)
expect_equal(x, pmax(-5:6, 0))
x <- (-5:6 + 0)
pmax0(x, sorted = TRUE, in_place = TRUE)
expect_equal(x, pmax(-5:6, 0))
x <- copy(5:-6)
pmax0(x, sorted = TRUE, in_place = TRUE)
expect_equal(x, pmax(5:-6, 0))
x <- 5:-6 + 0
pmax0(x, sorted = TRUE, in_place = TRUE)
expect_equal(x, pmax(5:-6, 0))


# test_that("pmax0 radix extrema", {
# Some systems cannot handle this much memory
can_do <- tryCatch(is.integer(every_int()), error = function(e) FALSE)

if (can_do && (at_home()) && hutilscpp:::is64bit()) {
  min_int <- -.Machine$integer.max
  max_int <- +.Machine$integer.max
  x <- min_int:max_int
  res <- pmax0(x)
  expect_equal(res[1], 0L)
  if (!hutilscpp:::is_covr()) {
    expect_equal(which_first(res > 0), max_int)
  }
  res <- NULL
  x <- NULL
  x <- max_int:min_int
  res <- pmax0(x)
  if (!hutilscpp:::is_covr()) {
    expect_equal(which_first(x == 0), max_int)
  }


  # test_that("firstNonnegativeRadix", {
  x <- c(-1, -1, 0, 0, 1, 1)
  expect_equal(firstNonNegativeRadix(x) + 1L, 3L)
  x <- hutilscpp_rev(x)
  expect_equal(firstNonNegativeRadix(x, desc = TRUE) + 1L, 3L)

  big <- seq(-99e3, 75e4, length.out = 1e4)
  expect_equal(which_first(big >= 0), firstNonNegativeRadix(big) + 1L)
}

# test_that("firstNonnegativeRadix desc", {
x <- 10:-1
expect_equal(firstNonNegativeRadix(x, desc = TRUE) + 1L, which_first(x <= 0))


# test_that("firstNonnegativeRadix corners", {
x <- 1:10
expect_true(firstNonNegativeRadix(x) <= 1)
expect_true(firstNonNegativeRadix(x, mini = 2) <= 2)
expect_true(firstNonNegativeRadix(-x, desc = TRUE) <= 1)
expect_true(firstNonNegativeRadix(-x, desc = TRUE, maxi = 5) <= 5)
x <- as.double(x)
expect_true(firstNonNegativeRadix(x) <= 1)
expect_true(firstNonNegativeRadix(x, mini = 2) <= 2)
expect_true(firstNonNegativeRadix(-x, desc = TRUE) <= 1)
expect_true(firstNonNegativeRadix(-x, desc = TRUE, maxi = 5) <= 5)

# Check bad arguments to mini
expect_equal(firstNonNegativeRadix(1:5, mini = -2L), 0)
expect_equal(firstNonNegativeRadix(1:5 + 0, mini = -2L), 0)



# test_that("Already nonnegative", {
x <- 1:100 + 0L
expect_equal(pmax0(x, sorted = TRUE), pmax(x, 0))
expect_equal(pmax0(x), 1:100)
expect_equal(do_pmax0_radix_sorted_int(x), 1:100)
x <- as.double(x)
expect_equal(pmax0(x, sorted = TRUE), 1:100)
expect_equal(pmax0(x), 1:100)
expect_equal(do_pmax0_radix_sorted_dbl(x), 1:100)
x <- x - 2L
expect_equal(pmax0(x), pmax(x, 0))
expect_equal(do_pmax0_radix_sorted_dbl(x), pmax(x, 0))
x <- as.integer(x)
expect_equal(pmax0(x, sorted = TRUE), pmax(x, 0))
expect_equal(pmax0(x), pmax(x, 0))
expect_equal(do_pmax0_radix_sorted_int(x), pmax(x, 0))
x <- hutilscpp_rev(x)
expect_equal(pmax0(x, sorted = TRUE), pmax(x, 0))
expect_equal(pmax0(x), pmax(x, 0))
expect_equal(do_pmax0_radix_sorted_int(x), pmax(x, 0))
x <- as.double(x)
expect_equal(pmax0(x, sorted = TRUE), pmax(x, 0))
expect_equal(pmax0(x), pmax(x, 0))
expect_equal(do_pmax0_radix_sorted_dbl(x), pmax(x, 0))


x <- c(x, -1L, x)
expect_equal(pmax0(x), pmax(x, 0L))
xd <- as.double(x)
expect_equal(pmax0(x), pmax(xd, 0))


# test_that("do_pmin0s", {
x <- c(0L, 6L, -4L, -2L, -1L, 7L, 9L, 4L, 8L, 3L, 10L, 5L, -3L, 1L, 2L)
expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_int(x), pmin(x, 0L))

x <- hutilscpp_rev(x)

expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_int(x), pmin(x, 0L))


x <- as.double(x)

expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_dbl(x), pmin(x, 0L))

x <- rev(x)

expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_dbl(x), pmin(x, 0L))



x <- pmin0(x)

# Now check already nonpositive
x <- as.double(x)
expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_dbl(x), pmin(x, 0L))

x <- as.integer(x)
expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_int(x), pmin(x, 0L))

# Now strictly positive
x <- x - 1L
expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_int(x), pmin(x, 0L))

x <- as.double(x)
expect_equal(pmin0(x), pmin(x, 0L))
x <- sort(x)
expect_equal(do_pmin0_radix_sorted_dbl(x), pmin(x, 0L))

x <- c(2, 1, 0, -1, -2)
expect_equal(do_pmin0_radix_sorted_dbl(x), pmin(x, 0))
x <- as.integer(x)
expect_equal(do_pmin0_radix_sorted_int(x), pmin(x, 0))



# test_that("in-place", {
abc <- -1:5 + 0L
def <- -1:5 + 0L
expect_equal(pmax0(abc, in_place = TRUE), pmax.int(def, 0L))
expect_equal(abc, pmax.int(def, 0L))



# test_that("pmax0 bitwise", {
expect_equal(do_pmax0_bitwise(-1:5), pmax.int(-1:5, 0L))
z <- c(1:10, 0L)
expect_equal(pmax0(z, in_place = TRUE), z)



# test_that("pmax0 sorted but all negative", {
expect_equal(pmax0(rep(-1L, 10), sorted = TRUE), integer(10))
expect_equal(pmax0(-10:-1, sorted = TRUE), integer(10))
expect_equal(pmax0(rep(-1, 10), sorted = TRUE), double(10))
expect_equal(pmax0(-10:-1 + 0, sorted = TRUE), double(10))
z <- c(-10:-1, 0L)
zd <- as.double(z)
expect_identical(pmax0(z, sorted = TRUE, in_place = TRUE), integer(11))
expect_equal(z, integer(11))
expect_identical(pmax0(zd, sorted = TRUE, in_place = TRUE), double(11))
expect_equal(zd, double(11))


# test_that("pmax0 altrep", {
expect_warning(pmax0(1:10, in_place = TRUE), "ALTREP")
expect_equal(pmax0(1:10), pmax(1:10, 0L))
expect_equal(pmax0(-1:-10), pmax(-1:-10, 0L))
expect_equal(pmax0(-1:10), pmax(-1:10, 0L))
expect_equal(pmax0(1:-10), pmax(1:-10, 0L))


# test_that("pmax0 sorted double", {
expect_equal(do_pmax0_radix_sorted_dbl(double(0)), double(0))
expect_equal(do_pmax0_radix_sorted_dbl(0.25), 0.25)
expect_equal(do_pmax0_radix_sorted_dbl(-0.25), 0)