我的第一个CUDA计划出了问题。在其中我尝试添加两个数组但是当我打印结果时完全错误。
这是我的代码:
#include <stdio.h>
const int N=10;
__global__ void add(int *a, int *b, int *c) {
int tid = threadIdx.x;
c[tid] = a[tid] + b[tid];
}
int main(){
int a[N], b[N], c[N];
int *dev_a, *dev_b, *dev_c;
cudaMalloc( (void**)&dev_a, N * sizeof(int) );
cudaMalloc( (void**)&dev_b, N * sizeof(int) );
cudaMalloc( (void**)&dev_c, N * sizeof(int) );
for (int i=0; i<N; i++) {
a[i] = -i; b[i] = i * i;
}
cudaMemcpy ( dev_a, a, N * sizeof(int), cudaMemcpyHostToDevice );
cudaMemcpy ( dev_b, b, N * sizeof(int), cudaMemcpyHostToDevice );
add<<<1,N>>>(dev_a, dev_b, dev_c);
cudaMemcpy(c, dev_c, N * sizeof(int),cudaMemcpyDeviceToHost );
for (int i=0; i<N; i++) {
printf("%d + %d = %d\n", a[i],b[i],c[i]);
}
cudaFree (dev_a); cudaFree (dev_b); cudaFree (dev_c);
return 0;
}
正如您将看到我尝试在C中添加A和B,但是,我不知道为什么,它不起作用。
这是输出:
0 + 0 = 0
-1 + 1 = 0
-2 + 4 = 1139031032
-3 + 9 = 32737
-4 + 16 = 0
-5 + 25 = 0
-6 + 36 = 1203890432
-7 + 49 = 32767
-8 + 64 = 6299112
-9 + 81 = 0
最后,我为我的英语道歉......在学习了cuda后,接下来我要做的就是改进它。
答案 0 :(得分:3)
您的机器出了问题。请error checking,你会知道它是什么。
当我编译并运行您的确切代码时,我得到以下内容:
$ ./t81
0 + 0 = 0
-1 + 1 = 0
-2 + 4 = 2
-3 + 9 = 6
-4 + 16 = 12
-5 + 25 = 20
-6 + 36 = 30
-7 + 49 = 42
-8 + 64 = 56
-9 + 81 = 72
$
以下是使用错误检查修改的代码。如果你编译并运行它,我相信你会发现你的机器有问题:
#include <stdio.h>
#define cudaCheckErrors(msg) \
do { \
cudaError_t __err = cudaGetLastError(); \
if (__err != cudaSuccess) { \
fprintf(stderr, "Fatal error: %s (%s at %s:%d)\n", \
msg, cudaGetErrorString(__err), \
__FILE__, __LINE__); \
fprintf(stderr, "*** FAILED - ABORTING\n"); \
exit(1); \
} \
} while (0)
const int N=10;
__global__ void add(int *a, int *b, int *c) {
int tid = threadIdx.x;
c[tid] = a[tid] + b[tid];
}
int main(){
int a[N], b[N], c[N];
int *dev_a, *dev_b, *dev_c;
cudaMalloc( (void**)&dev_a, N * sizeof(int) );
cudaMalloc( (void**)&dev_b, N * sizeof(int) );
cudaMalloc( (void**)&dev_c, N * sizeof(int) );
cudaCheckErrors("cudamalloc fail");
for (int i=0; i<N; i++) {
a[i] = -i; b[i] = i * i;
}
cudaMemcpy ( dev_a, a, N * sizeof(int), cudaMemcpyHostToDevice );
cudaMemcpy ( dev_b, b, N * sizeof(int), cudaMemcpyHostToDevice );
cudaCheckErrors("cuda memcpy fail");
add<<<1,N>>>(dev_a, dev_b, dev_c);
cudaMemcpy(c, dev_c, N * sizeof(int),cudaMemcpyDeviceToHost );
cudaCheckErrors("cudamemcpy or cuda kernel fail");
for (int i=0; i<N; i++) {
printf("%d + %d = %d\n", a[i],b[i],c[i]);
}
cudaFree (dev_a); cudaFree (dev_b); cudaFree (dev_c);
return 0;
}
This poster做了同样的事情,发现他们的机器出了问题。