tests/testthat/test_mistakeInHistory.R

In gpuMagic: An openCL compiler with the capacity to compile R functions and run the code on GPU

context("A smart man will never fall at the same place twice.")
skip_if(gpuMagic:::getTotalDeviceNum()==0)
n=100
m=200
k=50

test_that("matrix multiplication",{
  testFunc<-function(ind,A,B){
    C=matrix(0,nrow(A),1)
    for(i in 1:nrow(A)){
      tmp=0
      for(k in 1:ncol(A)){
        tmp=tmp+A[i,k]*B[k,ind]
      }
      C[i]=tmp
    }
    return(C)
  }
  
  
  A=matrix(runif(n*m),n,m)
  B=matrix(runif(k*m),m,k)
  res_gpu=gpuSapply(1:k,testFunc,A,B)
  res_cpu=sapply(1:k,testFunc,A,B)
  
  expect_equal(res_gpu,res_cpu)
})


test_that("Mixture of two index sub and one index sub",{
  testFunc<-function(ind,A){
    A_col = subRef(A, , ind)
    i=1
    tmp=subRef(A_col,i)
    return(tmp)
  }
  A=matrix(runif(n*m),n,m)
  res_gpu=gpuSapply(1:k,testFunc,A)
  res_cpu=sapply(1:k,testFunc,A)
  
  expect_equal(res_gpu,res_cpu)
})

test_that("precision bump",{
  testFunc<-function(ind){
    A=matrix(0,1,1)
    A[1,1]=ind/7
    return(A)
  }
  res_gpu=gpuSapply(1:k,testFunc)
  res_cpu=sapply(1:k,testFunc)
  
  expect_equal(res_gpu,res_cpu)
})