tests/sum.r

In float: 32-Bit Floats

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

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

x1 = matrix(-4:5, 5)
s1 = fl(x1)
y = 1L
z = 1.0
s2 = flrunif(5)


test = dbl(sum(s1))
truth = sum(x1)
same(test, truth)

test = dbl(sum(s1, s2))
truth = sum(x1, dbl(s2))
same(test, truth)

test = dbl(sum(s1, y))
truth = sum(x1, y)
same(test, truth)



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

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

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