library("gpuR")
ORDER = 10000

A = matrix(rnorm(ORDER^2), nrow=ORDER)
B = matrix(rnorm(ORDER^2), nrow=ORDER)
gpuA = gpuMatrix(A, type="double")
gpuB = gpuMatrix(B, type="double")

print("cpu time: ")
system.time(A %*% B)

print("gpu time: ")
system.time(gpuA %*% gpuB)