R 的革命

Question

自从最新的 Ubuntu 版本(karmic koala)以来，我注意到内部 R 包在启动时做广告 REvolution包裹。
它似乎是一个用于高性能矩阵计算的库集合。显然，它似乎真的有效。例如在带有 REvolution 的矩阵转置上:

> system.time(t(matrix(rnorm(10000000),ncol=1000)))
   user  system elapsed 
  1.280   0.150   1.556 

没有革命:

> system.time(t(matrix(rnorm(10000000),ncol=1000)))
   user  system elapsed 
  1.320   0.170   1.725 

有人使用它吗？它真的有效吗？它正在改进哪些特定类型的计算以及如何改进？有什么缺点吗？

谢谢 :-)

Answer 1

是的，在多核机器上，英特尔 MKL 库实现了 BLAS -- 并由包 revolution-mkl 提供并通过封装开启 r-revolution-revobase将并行处理线性代数问题，您应该会看到与仅使用 libblas* 包的基本情况有所不同。

但是，你上面的例子并不是那么重要，我经常做类似的事情

 mean(replicate(N, system.time( someStuffHere() )["elapsed"]), trim=0.05)

计算多次重复的修剪平均值。

更重要的是，请注意您的示例包括 RNG 绘制的时间
i) 昂贵，并且 ii) 对所使用的方法是不变的，因此您应该在 system.time() 之外生成它.

除了 MKL， revolution-r还引入了一些来自 CRAN 的 REvolution 编写的包，可用于并行执行。

(免责声明:我帮助 REvo 将其整合到 Ubuntu 9.10 中)