blas-benchmarks/benchmark-sample.R

runs = 10

f = function()
{
    set.seed (1)
    m <- 3500
    n <- 3500
    A <- matrix (runif (m*n),m,n)

    # Matrix multiply
    B <- crossprod(A)
}

cumulate = 0
for (i in 1:runs)
    cumulate = cumulate + as.numeric(system.time(f())[3])
    
cat(paste0(round(cumulate/runs, 3), "\n"))