https://github.com/cran/pracma
Tip revision: 1305bf51cc38adca02d0c8a834c61a4c7e038309 authored by Hans W. Borchers on 08 February 2015, 00:00:00 UTC
version 1.8.3
version 1.8.3
Tip revision: 1305bf5
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)
}