apply（）の最適化

Question

次のコードを最適化したいと思います。現在はMが2404 100の数値行列であるために、2メガバイトL2キャッシュと2GHzのデュアルコアマシン上の周りに0.085秒を実行します：

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())を呼び出すことによって、さらにそれをスピードアップすることができますが、それがうまくいかなかったことを考えました。私はそれを正しく使用していないので、単純にそれができます。