import pycuda.autoinit import pycuda.gpuarray import numpy as np from scikits.cuda import integrate import pycuda.gpuarray as gpuarray import pycuda.autoinit integrate.init() x = np.asarray(np.random.rand(10), np.float32) x_gpu = gpuarray.to_gpu(x) z = integrate.trapz(x_gpu, 0.1) np.allclose(np.trapz(x), z)
def setUp(self): np.random.seed(0) integrate.init()
def setUp(self): integrate.init()
import numpy as np import pycuda.autoinit import pycuda.gpuarray as gpuarray import scikits.cuda.integrate as integrate integrate.init() a = np.asarray(np.random.rand(10), np.float32) #x = np.asarray(np.random.rand(10, 10), np.float32) x = np.ones( (2000, 2000) ) * 10 x_gpu = gpuarray.to_gpu(x) a_gpu = gpuarray.to_gpu(a) b = integrate.trapz(a_gpu) z = integrate.trapz2d(x_gpu) print(z) np.allclose(np.trapz(np.trapz(x)), z)