eps.R
##
## e p s . R
##
eps <- function(x = 1.0) {
stopifnot(is.numeric(x))
x <- max(abs(x))
if (x < .Machine$double.xmin) {
e <- .Machine$double.xmin
} else {
e <- 2^floor(log2(x)) * .Machine$double.eps
}
e
}