2016-04-29 14:13:49 +02:00
|
|
|
library(rbenchmark)
|
|
|
|
|
|
|
|
f = function()
|
|
|
|
{
|
|
|
|
set.seed (1)
|
|
|
|
m <- 1000
|
|
|
|
n <- 500
|
|
|
|
A <- matrix (runif (m*n),m,n)
|
|
|
|
|
|
|
|
# Matrix multiply
|
|
|
|
B <- crossprod(A)
|
|
|
|
}
|
|
|
|
|
2016-05-29 17:45:39 +02:00
|
|
|
print(benchmark(f(), replications=100))
|
2016-04-29 14:13:49 +02:00
|
|
|
|