Question

我正在尝试用RcppArmadillo做一些矩阵乘法。但是，我的代码显示使用RcppArmadillo 不会变得更快。

我正在使用Windows_10_Pro和R 3.2.4，以及RcppArmadillo 0.6.600.4.0

例如：

library(RcppArmadillo)
library(inline)

MCplus <- cxxfunction(signature(X_="numeric", Y_="numeric"),body ='
arma::mat X = Rcpp::as<arma::mat>(X_);
arma::mat Y = Rcpp::as<arma::mat>(Y_);
arma::mat ans = X * Y * X;
return(wrap(ans));
', plugin="RcppArmadillo")

A <- matrix(1:16000000,4000,4000)
C <- matrix(2:16000001,4000,4000)

R_M <- proc.time()
ans_R <- A%*%C%*%A # test with R
proc.time() - R_M

C_M <- proc.time()
ans_C <- MCplus(A,C) # test with RcppArmadillo
proc.time() - C_M

R输出：

user  system elapsed
106.75   0.24 106.98

RcppArmadillo输出：

user  system elapsed
108.28   0.23 108.56

有什么可以改进吗？

提前致谢！

Answer 1

R本身将其归于LAPACK / BLAS - 链接到R的代码通过LAPACK / BLAS进行调用。所以，是的，两种方法都将运行相同的代码，差异只是由于开销很小。

有许多教程可以告诉您如何更改LAPACK库。这当然取决于操作系统。可以使用R Installation and Administration手册及其附录开始。

使用RcppArmadillo进行矩阵乘法：为什么它不会更快？

1 个答案: