Эффективное создание матрицы из значений функций

Что опытные разработчики R считают наиболее эффективным (но при этом читабельным) способом построения матрицы с заданным количеством строк и столбцов из заданной функции, такой что, например, A_ij = someFun(i,j) с 1 <= i <= rows, 1 <= j <= cols?

Поскольку я не смог найти что-то в документации, я придумал

initMatrix <- function(rows, cols, fn) {
    A <- matrix(nrow=rows, ncol=cols)
    for (i in 1:rows)
        for (j in 1:cols)
            A[i,j] <- fn(i,j)
    return(A)
}

что кажется мне глупым и медленным. Любые улучшения (особенно однострочные) приветствуются! :)