|
|
|
@@ -2,9 +2,7 @@ |
|
|
|
|
|
|
|
argv <- commandArgs(trailingOnly = TRUE) |
|
|
|
|
|
|
|
if (!is.null(options("matprod")[[1]])) { |
|
|
|
options(matprod = "blas") |
|
|
|
} |
|
|
|
if (!is.null(options("matprod")[[1]])) options(matprod = "blas") |
|
|
|
|
|
|
|
nfrom <- 128 |
|
|
|
nto <- 2048 |
|
|
|
@@ -42,11 +40,10 @@ while (n <= nto) { |
|
|
|
solve(A, B) |
|
|
|
}) |
|
|
|
|
|
|
|
mflops <- (2.0/3 * n * n * n + 2 * n * n * n) * loops/ (z[3] * 1e6) |
|
|
|
mflops <- (2.0/3 * n * n * n + 2 * 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 |
|
|
|
} |
|
|
|
|