https://github.com/cran/pracma
Tip revision: 3fdb68cc842f2ab3b59608f8e04d495763584f78 authored by Hans W. Borchers on 27 November 2015, 12:07:10 UTC
version 1.8.8
version 1.8.8
Tip revision: 3fdb68c
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)
}