@@ -2095,7 +2095,7 @@ rollfun = function(x, n, FUN, fill=NA_real_, na.rm=FALSE, nf.rm=FALSE, partial=F
20952095 ans
20962096}
20972097base_compare = function(x, n, funs=c("mean","sum","max","min","prod"), algos=c("fast","exact")) {
2098- num.step = 0.001
2098+ num.step = 0.0001
20992099 for (fun in funs) {
21002100 for (na.rm in c(FALSE, TRUE)) {
21012101 for (fill in c(NA_real_, 0)) {
@@ -2129,7 +2129,40 @@ base_compare = function(x, n, funs=c("mean","sum","max","min","prod"), algos=c("
21292129 }
21302130 }
21312131}
2132- num = 7001.0
2132+ num = 7000.0
2133+ x = rnorm(1e3); n = 50
2134+ base_compare(x, n)
2135+ x = rnorm(1e3+1); n = 50 ## uneven len
2136+ base_compare(x, n)
2137+ x = rnorm(1e3); n = 51 ## uneven window
2138+ base_compare(x, n)
2139+ x = rnorm(1e3+1); n = 51
2140+ base_compare(x, n)
2141+ x = sort(rnorm(1e3)); n = 50 ## inc
2142+ base_compare(x, n)
2143+ x = sort(rnorm(1e3+1)); n = 50
2144+ base_compare(x, n)
2145+ x = sort(rnorm(1e3)); n = 51
2146+ base_compare(x, n)
2147+ x = sort(rnorm(1e3+1)); n = 51
2148+ base_compare(x, n)
2149+ x = rev(sort(rnorm(1e3))); n = 50 ## desc
2150+ base_compare(x, n)
2151+ x = rev(sort(rnorm(1e3+1))); n = 50
2152+ base_compare(x, n)
2153+ x = rev(sort(rnorm(1e3))); n = 51
2154+ base_compare(x, n)
2155+ x = rev(sort(rnorm(1e3+1))); n = 51
2156+ base_compare(x, n)
2157+ x = rep(rnorm(1), 1e3); n = 50 ## const
2158+ base_compare(x, n)
2159+ x = rep(rnorm(1), 1e3+1); n = 50
2160+ base_compare(x, n)
2161+ x = rep(rnorm(1), 1e3); n = 51
2162+ base_compare(x, n)
2163+ x = rep(rnorm(1), 1e3+1); n = 51
2164+ base_compare(x, n)
2165+ num = 7100.0
21332166## random NA non-finite
21342167x = makeNA(rnorm(1e3), nf=TRUE); n = 50
21352168base_compare(x, n)
@@ -2185,7 +2218,7 @@ if (requireNamespace("zoo", quietly=TRUE)) {
21852218 }
21862219 }
21872220 }
2188- num = 7002 .0
2221+ num = 7200 .0
21892222 ## no NA
21902223 x = rnorm(1e3); n = 50 # x even, n even
21912224 zoo_compare(x, n)
@@ -2285,7 +2318,40 @@ afun_compare = function(x, n, funs=c("mean","sum","max","min","prod"), algos=c("
22852318 }
22862319 }
22872320}
2288- num = 7003.0
2321+ num = 7300.0
2322+ x = rnorm(1e3); n = sample(50, length(x), TRUE)
2323+ afun_compare(x, n)
2324+ x = rnorm(1e3+1); n = sample(50, length(x), TRUE) ## uneven len
2325+ afun_compare(x, n)
2326+ x = rnorm(1e3); n = sample(51, length(x), TRUE) ## uneven window
2327+ afun_compare(x, n)
2328+ x = rnorm(1e3+1); n = sample(51, length(x), TRUE)
2329+ afun_compare(x, n)
2330+ x = sort(rnorm(1e3)); n = sample(50, length(x), TRUE) ## inc
2331+ afun_compare(x, n)
2332+ x = sort(rnorm(1e3+1)); n = sample(50, length(x), TRUE)
2333+ afun_compare(x, n)
2334+ x = sort(rnorm(1e3)); n = sample(51, length(x), TRUE)
2335+ afun_compare(x, n)
2336+ x = sort(rnorm(1e3+1)); n = sample(51, length(x), TRUE)
2337+ afun_compare(x, n)
2338+ x = rev(sort(rnorm(1e3))); n = sample(50, length(x), TRUE) ## desc
2339+ afun_compare(x, n)
2340+ x = rev(sort(rnorm(1e3+1))); n = sample(50, length(x), TRUE)
2341+ afun_compare(x, n)
2342+ x = rev(sort(rnorm(1e3))); n = sample(51, length(x), TRUE)
2343+ afun_compare(x, n)
2344+ x = rev(sort(rnorm(1e3+1))); n = sample(51, length(x), TRUE)
2345+ afun_compare(x, n)
2346+ x = rep(rnorm(1), 1e3); n = sample(50, length(x), TRUE) ## const
2347+ afun_compare(x, n)
2348+ x = rep(rnorm(1), 1e3+1); n = sample(50, length(x), TRUE)
2349+ afun_compare(x, n)
2350+ x = rep(rnorm(1), 1e3); n = sample(51, length(x), TRUE)
2351+ afun_compare(x, n)
2352+ x = rep(rnorm(1), 1e3+1); n = sample(51, length(x), TRUE)
2353+ afun_compare(x, n)
2354+ num = 7400.0
22892355#### no NA
22902356x = rnorm(1e3); n = sample(50, length(x), TRUE) # x even, n even
22912357afun_compare(x, n)
0 commit comments