Question

我正在尝试使用赋值

设置3D循环

 C(i,j,k) = A(i,j,k) + B(i,j,k)

在我的GPU上使用Python。这是我的GPU：

http://www.geforce.com/hardware/desktop-gpus/geforce-gt-520/specifications

我正在寻找/比较的来源是：

http://nbviewer.ipython.org/gist/harrism/f5707335f40af9463c43

http://nbviewer.ipython.org/github/ContinuumIO/numbapro-examples/blob/master/webinars/2014_06_17/intro_to_gpu_python.ipynb

我可能导入的模块多于必要的模块。这是我的代码：

import numpy as np
import numbapro
import numba
import math
from timeit import default_timer as timer
from numbapro import cuda
from numba import *

@autojit
def myAdd(a, b):
  return a+b

myAdd_gpu = cuda.jit(restype=f8, argtypes=[f8, f8], device=True)(myAdd)

@cuda.jit(argtypes=[float32[:,:,:], float32[:,:,:], float32[:,:,:]])
def myAdd_kernel(a, b, c):
    tx = cuda.threadIdx.x
    ty = cuda.threadIdx.y
    tz = cuda.threadIdx.z
    bx = cuda.blockIdx.x
    by = cuda.blockIdx.y
    bz = cuda.blockIdx.z
    bw = cuda.blockDim.x
    bh = cuda.blockDim.y
    bd = cuda.blockDim.z
    i = tx + bx * bw
    j = ty + by * bh
    k = tz + bz * bd
    if i >= c.shape[0]:
      return
    if j >= c.shape[1]:
      return
    if k >= c.shape[2]:
      return
    for i in xrange(0,c.shape[0]):
      for j in xrange(0,c.shape[1]):
        for k in xrange(0,c.shape[2]):
          # c[i,j,k] = a[i,j,k] + b[i,j,k]
          c[i,j,k] = myAdd_gpu(a[i,j,k],b[i,j,k])

def main():
    my_gpu = numba.cuda.get_current_device()
    print "Running on GPU:", my_gpu.name
    cores_per_capability = {1: 8,2: 32,3: 192,}
    cc = my_gpu.compute_capability
    print "Compute capability: ", "%d.%d" % cc, "(Numba requires >= 2.0)"
    majorcc = cc[0]
    print "Number of streaming multiprocessor:", my_gpu.MULTIPROCESSOR_COUNT
    cores_per_multiprocessor = cores_per_capability[majorcc]
    print "Number of cores per mutliprocessor:", cores_per_multiprocessor
    total_cores = cores_per_multiprocessor * my_gpu.MULTIPROCESSOR_COUNT
    print "Number of cores on GPU:", total_cores

    N = 100
    thread_ct = my_gpu.WARP_SIZE
    block_ct = int(math.ceil(float(N) / thread_ct))

    print "Threads per block:", thread_ct
    print "Block per grid:", block_ct

    a = np.ones((N,N,N), dtype = np.float32)
    b = np.ones((N,N,N), dtype = np.float32)
    c = np.zeros((N,N,N), dtype = np.float32)

    start = timer()
    cg = cuda.to_device(c)
    myAdd_kernel[block_ct, thread_ct](a,b,cg)
    cg.to_host()
    dt = timer() - start
    print "Wall clock time with GPU in %f s" % dt
    print 'c[:3,:,:] = ' + str(c[:3,1,1])
    print 'c[-3:,:,:] = ' + str(c[-3:,1,1])


if __name__ == '__main__':
    main()

运行此操作的结果如下：

Running on GPU: GeForce GT 520
Compute capability:  2.1 (Numba requires >= 2.0)
Number of streaming multiprocessor: 1
Number of cores per mutliprocessor: 32
Number of cores on GPU: 32
Threads per block: 32
Block per grid: 4
Wall clock time with GPU in 1.104860 s
c[:3,:,:] = [ 2.  2.  2.]
c[-3:,:,:] = [ 2.  2.  2.]

当我在源代码中运行示例时，我看到了显着的加速。我不认为我的例子运行正常，因为挂钟时间比我预期的要长得多。我在第一个示例链接中主要通过“使用cuda python进行更大的加速”部分对此进行了建模。

我相信我已正确安全地编入索引。也许问题出在我的blockdim上？还是griddim？或许我的GPU使用了错误的类型。我想我读过它们必定是某种类型。我对此非常陌生，所以这个问题很可能是微不足道的！

非常感谢任何和所有帮助！

Answer 1

您正在正确创建索引，但是您忽略了它们。运行嵌套循环

for i in xrange(0,c.shape[0]):
    for j in xrange(0,c.shape[1]):
        for k in xrange(0,c.shape[2]):

强制所有线程循环遍历所有维度中的所有值，这不是您想要的。您希望每个线程计算一个块中的一个值，然后继续。

我认为这样的事情会更好......

i = tx + bx * bw
while i < c.shape[0]:
    j = ty+by*bh
    while j < c.shape[1]:
        k = tz + bz * bd
        while k < c.shape[2]:
            c[i,j,k] = myAdd_gpu(a[i,j,k],b[i,j,k])
            k+=cuda.blockDim.z*cuda.gridDim.z
        j+=cuda.blockDim.y*cuda.gridDim.y
    i+=cuda.blockDim.x*cuda.gridDim.x

尝试编译并运行它。还要确保验证它，因为我没有。

Answer 2

我没有看到你使用imshow或show，所以没有必要导入它们。

看起来好像你没有使用数学导入（我没有看到任何math.some_function的调用。

你从numba和numbapro的进口似乎是重复的。你的“来自numba import cuda”会覆盖你的“来自numbapro import cuda”，因为它是紧随其后的。你打电话给cuda使用numba中的cuda而不是numbapro。当你打电话给“来自numba import *”时，你从numba导入所有东西，而不仅仅是cuda，这似乎是你唯一使用的东西。另外，（我相信）导入numba.cuda相当于来自numba import cuda。为什么不用一个“来自numba import cuda”消除你从numba和numbapro的所有进口。

cuda python GPU numbapro 3d循环性能不佳

2 个答案: