> library(gputools)
> 
> set.seed(5446)
> p <- 20
> X <- matrix(rnorm(2^p),ncol = 2^(p/2))
> 
> dtime <- system.time(d <- dist(X))
> gputime <- system.time(gpud <- gpuDist(X))
> 
> dtime
   user  system elapsed 
  4.995   0.002   5.002 
> gputime
   user  system elapsed 
  0.193   0.079   0.278 
> dtime/gputime
       user      system     elapsed 
25.88082902  0.02531646 17.99280576 
> max(abs(c(d) - c(gpud)))
[1] 6.05222e-06