R Under development (unstable) (2023-12-09 r85665 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > 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) > > proc.time() user system elapsed 0.32 0.01 0.36