argv <- commandArgs(trailingOnly = TRUE)
+if (!is.null(options("matprod")[[1]])) options(matprod = "blas")
+
nfrom <- 128
nto <- 2048
nstep <- 128
loops <- as.numeric(argv[z])
}
}
-
}
p <- Sys.getenv("OPENBLAS_LOOPS")
loops <- as.numeric(p)
}
-
-cat(sprintf(
- "From %.0f To %.0f Step=%.0f Loops=%.0f\n",
- nfrom,
- nto,
- nstep,
- loops
-))
+cat(sprintf("From %.0f To %.0f Step=%.0f Loops=%.0f\n", nfrom, nto, nstep, loops))
cat(sprintf(" SIZE Flops Time\n"))
n <- nfrom
while (n <= nto) {
- A <- matrix(runif(n * n),
- ncol = n,
- nrow = n,
- byrow = TRUE)
- B <- matrix(runif(n * n),
- ncol = n,
- nrow = n,
- byrow = TRUE)
+ A <- matrix(runif(n * n), nrow = n)
+ B <- matrix(runif(n * n), nrow = n)
C <- 1
z <- system.time(for (l in 1:loops) {
l <- l + 1
})
- mflops <- (2.0 * n * n * n) * loops / (z[3] * 1.0e6)
+ mflops <- (2.0 * n * n * n) * loops / (z[3] * 1e+06)
st <- sprintf("%.0fx%.0f :", n, n)
cat(sprintf("%20s %10.2f MFlops %10.6f sec\n", st, mflops, z[3]))
n <- n + nstep
-
}