Оптимизация apply()

Question

Я хотел был бы оптимизировать следующий код. В настоящее время работает около 0,085 секунд на 2 ГГц двухъядерный машины с 2 Мб кэша L2, для M быть 2404 на 100 числовая матрица:

Rescale <- function(M = utility.mat){ 
    exp.M <- exp(M) 
    result <- apply(exp.M, 1, function(x) x/sum(x)) 
    result <- t(result) 
    return (result) 
} 

Я попытался заменить apply() с цикл, который дает примерно такую ​​же производительность. Любые другие идеи?

Answer 1

Это, кажется, бой в 6 раз быстрее, на моей машине:

Rescale1 <- function(M){ 
    M <- exp(M) 
    result <- M/rowSums(M) 
    return (result) 
} 

Я подумал, что, может быть, вы могли бы ускорить его дальше по телефону .Internal(rowSums()), но это не сработало. Хотя это может быть просто потому, что я не использую его правильно.