swh:1:snp:a4c99a50dc49f82b591f268001b320f8c3ca0041
Tip revision: 698183c35adc80e5e08cb2ef09ed13dc3c445a97 authored by John M Chambers on 12 June 2014, 00:00:00 UTC
version 1.0-6
version 1.0-6
Tip revision: 698183c
binaryRep.R
binaryRep <-
function(data, m = .Machine$double.digits) {
x <- data
n <- length(x)
xSign <- sign(x)
x <- xSign * x
exponent <- ifelse(x > 0, floor(1+log(x, 2)), 0)
x <- x/2 ^ exponent
pwrs <- binaryRepPowers(n, m)
x <- matrix(x, n, m)
xpwrs <- x * pwrs
xrep <- trunc(xpwrs)/pwrs
bits <- (xrep %*% binaryRepA(m)) *pwrs
bits[] <- as.integer(bits[])
new("binaryRep", original = data,
sign = as.integer(xSign),
exponent = as.integer(exponent),
bits = binaryRepBits(bits))
}