swh:1:snp:81eadaa089e8253d8469bcef66aa332632c6c669
Tip revision: 6bedfa382c05ccbdda30ab12438b32fd9ba13876 authored by HwB on 12 December 2012, 00:00:00 UTC
version 1.3.3
version 1.3.3
Tip revision: 6bedfa3
pow2.R
###
### POW2.R Raise 2 to some power
###
pow2 <- function(f, e) {
if (!is.numeric(f) && !is.complex(f))
stop("Argument 'f' must be numeric or complex.")
if (missing(e)) {
e <- f
f <- rep(1, length(e))
} else {
if (!is.numeric(e) && !is.complex(e))
stop("Argument 'e' must be numeric or complex.")
if (is.complex(f) || is.complex(e)) {
f <- Re(f)
e <- Re(e)
warning("Imaginary part of arguments ignored.")
}
}
if (length(f) != length(e))
stop("Arguments 'e' and 'f' must be of same length.")
return(f * 2^e)
}