tests/extremes.r

In float: 32-Bit Floats

suppressPackageStartupMessages(library(float))
set.seed(1234)
tol=1e-6

same = function(test, truth) stopifnot(all.equal(dbl(test), truth, tol=tol))

x1 = matrix(-4:5, 5)
s1 = fl(x1)

x2 = x1
x2[3, 1] = NA_integer_
x2[4, 2] = NA_integer_
s2 = fl(x2)

y = 1L
z = 1.0
s3 = flrunif(5)


### min
test = min(s1)
truth = min(x1)
same(test, truth)

test = min(s1, s3)
truth = min(x1, dbl(s3))
same(test, truth)

test = min(s1, y)
truth = min(x1, y)
same(test, truth)

# NA
test = sum(s2, na.rm=FALSE)
truth = sum(x2, na.rm=FALSE)
same(test, truth)

test = sum(s2, na.rm=TRUE)
truth = sum(x2, na.rm=TRUE)
same(test, truth)



### max
test = max(s1)
truth = max(x1)
same(test, truth)

test = max(s1, s3)
truth = max(x1, dbl(s3))
same(test, truth)

test = max(s1, y)
truth = max(x1, y)
same(test, truth)

# NA
test = sum(s2, na.rm=FALSE)
truth = sum(x2, na.rm=FALSE)
same(test, truth)

test = sum(s2, na.rm=TRUE)
truth = sum(x2, na.rm=TRUE)
same(test, truth)




# which
test = which.min(s1)
truth = which.min(x1)
same(test, truth)

test = which.min(s2)
truth = which.min(x2)
same(test, truth)

test = which.max(s1)
truth = which.max(x1)
same(test, truth)

test = which.max(s2)
truth = which.max(x2)
same(test, truth)