https://github.com/cran/pracma
Tip revision: 7ba95550c2066826c77bc78b974b7cb8cd416301 authored by Hans W. Borchers on 11 January 2017, 00:43:13 UTC
version 1.9.9
version 1.9.9
Tip revision: 7ba9555
bsxfun.R
##
## b s x f u n . R
##
bsxfun <- function(func, x, y) {
stopifnot(is.numeric(x), is.numeric(y))
# fun <- match.fun(f)
# f <- function(x, y) fun(x, y, ...)
dx <- dim(x); dy <- dim(y)
if ( is.vector(x) && is.vector(y) ) {
z <- mapply(func, x, y)
} else if (is.array(x) && is.array(y) && all(dx == dy)) {
z <- mapply(func, x, y)
dim(z) <- dx
} else {
stop("Argument 'x', 'y' must be vectors or arrays of the same size.")
}
return(z)
}
arrayfun <- function(func, ...) {
# func <- match.fun(func)
dots <- list(...)
if (length(dots) < 1)
stop("Empty list of arrays: Rsult cannot be computed.")
d <- dim(dots[[1]]) # no test on array sizes to be fast
r <- mapply(func, ...) # no try ... catch, number of variables
dim(r) <- d
return(r)
}