cuda - Cuda C中的加法赋值运算符

Question

我在 Cuda C 中遇到了加法赋值运算符的问题。我收到以下错误：

kernel.cu(5): error: expression must have integral or enum type

我的代码是：

import pycuda.driver as drv
import pycuda.autoinit
from pycuda.compiler import SourceModule
import numpy as np

mod=SourceModule("""
__global__ void addition(float* a,float* b,float*c){
int i=threadIdx.x + blockIdx.x * blockDim.x;
c[a[i]]+=b[i];
}
""")

addition=mod.get_function("addition")
a=np.array([1,2,3,1,2,3,2,1]).astype(np.float32)
b=np.array([0.1,0.2,0.1,0.5,0.1,0.2,0.1,0.5]).astype(np.float32)
c=np.zeros_like(a)
addition(drv.Out(c),drv.In(a),drv.In(b),block=(32,1,1))
print c

我想要的输出是 c = [0,1.1,0.4,0.3,0,0,0,0]。任何人都可以提出解决方案吗？

Answer 1

问题出在您的内核中，您使用 A 在 C 中进行索引
。A 是浮点类型。

另请注意，您正在启动 32 个线程，但您只会索引 8 个位置，这意味着您将索引超出范围。

您将面临的最后一个问题是，由于 a 中的索引重复，多个线程试图更改 C 中的相同位置。解决它的一种方法是使用 AtomicAdd。

__global__ void addition(float* a,float* b,float*c, int n)
{
int i=threadIdx.x + blockIdx.x * blockDim.x;
if(i < n)
atomicAdd(&c[(int)a[i]],b[i]);
}

以相同的方式启动内核，但不要忘记传递 a 或 b 的大小的 n。
您还可以在启动内核时消除 n 并更改线程块尺寸。