Commit 00923a3e authored by Martin Maechler's avatar Martin Maechler

a bit more on platform-dependencies (BLAS/Lapack)

parent 1b5d7de9
Package: CLA
Version: 0.90-2
Date: 2018-02-08
Date: 2018-09-10
Title: Critical Line Algorithm in Pure R
Author: Yanhao Shi <syhelena@163.com>,
Martin Maechler <maechler@stat.math.ethz.ch>
......
......@@ -80,7 +80,9 @@ file.info("wtsn0.rds")$size # 27049
wtsn0.ref <- readRDS("wtsn0.rds")
## see on all platforms what we get
all.equal(target = wtsn0.ref, current = wts.non0, tol=0) # expect TRUE only on 64bit (Lnx)
all.equal(target = wtsn0.ref, current = wts.non0, tol=0)
# expect TRUE only on 64bit (Lnx, R's BLAS)
# 3.10416e-15 and 1.366427e-15 on other BLAS
stopifnot(all.equal(target = wtsn0.ref, current = wts.non0,
tol = 1e-13))
......@@ -134,17 +136,42 @@ if(all(non.0.assets == non.0.TARG)) { ## show differences:
cat("\n'non.0.assets' differing from non.0.TARG:\n")
cat("+1:\n"); dput(names(b64.n0)[non.0.assets - non.0.TARG == +1])
cat("-1:\n"); dput(names(b64.n0)[non.0.assets - non.0.TARG == -1])
## solaris :
## +1:
## c("AET", "BCR", "CI", "CL", "ED", "FE", "HAL", "MCD", "SII", "SYK")
## -1:
## c("AZO", "BAX", "CLX", "COST", "DGX", "DVN", "ESRX", "LMT", "MUR",
## "PEP", "RIG", "SYMC", "TYC", "UST")
## ATLAS :
## +1:
## c("AZO", "BRL", "CCE", "CLX", "INTU", "JNJ", "K", "LMT", "MUR",
## "PEP", "SSP", "TYC", "UST", "XTO")
## -1:
## c("AET", "BCR", "CI", "ED", "FE", "MCD", "NEM", "SII", "SYK", "WMT")
## MKL:
## +1:
## c("CLX", "INTU", "LH", "LLL", "LMT", "PBG", "SYMC", "TYC")
## -1:
## c("AMGN", "BCR", "BUD", "CL", "CTL", "ED", "HAL", "NEM", "XTO")
## OpenBLAS:
## +1:
## c("BAX", "COST", "HST", "JNJ", "MDT", "MUR", "NOC", "PDCO", "WAG")
## -1:
## c("BCR", "CI", "CTL", "ED", "EIX", "HAL", "MCD", "RAI", "SYK")
if(any(isB <- abs(non.0.assets - non.0.TARG) > 1)) {
cat("more different, showing differences:\n")
dput((non.0.assets - non.0.TARG)[isB])
}
## OpenBLAS (only!):
## c(XTO = -2L)
}
## They have the same names and only differ by +/- 1:
stopifnot(
identical(names(b64.n0), names(b32.n0))
## ______ ______
, if(b64) identical(non.0.assets, non.0.TARG)
, if(b64) identical(non.0.assets, non.0.TARG) # fails on ATLAS, MKL, OpenBLAS
else if(nonWindows) identical(non.0.assets, b32.n0)
else ## 32-bit Windows
TRUE ## for now
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment