CUDA编程问题3

问题描述

| 我有两个一维数组。现在,我必须对这两个数组执行“与”运算,如果数组中的元素相等,则必须将其放入2d数组中。 我已经编写了以下代码,但是我的代码无法正常工作。谁能告诉我我做错了什么?
#include<stdio.h>
#include<stdio.h>
#include<cuda.h>
#define height 3
#define width 2

__global__ void fun(unsigned char *c,unsigned short *s,unsigned char *u,size_t pitch)
{ 

    int tid =blockIdx.x * blockDim.x + threadIdx.x ;
    //int tidy=blockIdx.y * blockDim.y + threadIdx.y ;

    if(tid<6)
    {

        for (int r = 0; r < height; ++r)
        {        //float* row = (float*)((char*)u + r*pitch);

                for (int d = 0; d< width; ++d) 
                {
                    u[r*width+d] = c[tid] & s[tid];
                }
        }
    } 
}


int main()
{
    int i,j;
    unsigned char man[6]=\"manis\",*dev_c,jan[3][2],*dev_r;

    unsigned short a[6]={32,33,43,35,36,37},*dev_s;

    size_t pitch;

    cudamalloc((void**)&dev_c,sizeof( unsigned char)*6);       // memory allocation on device

    cudamalloc((void**)&dev_s,sizeof(unsigned short)*6);

    cudamallocPitch((void**)&dev_r,&pitch,width*sizeof(unsigned char),height);

    cudamemcpy(dev_c,man,6*sizeof(unsigned char),cudamemcpyHostToDevice);

    cudamemcpy(dev_s,a,6*sizeof(short int),cudamemcpyHostToDevice);

    fun<<<5,2>>>(dev_c,dev_s,dev_r,pitch);

    cudamemcpy2D(jan,3*sizeof(unsigned char),pitch,sizeof( unsigned char),2,cudamemcpyDevicetoHost);

    printf(\"\\nThe call is completd\");

    for ( i=0; i < width; i++)
    {
        for (j=0; j < height; j++)
        {
            printf(\"%d \",jan[j+i*width]);
        }
有人可以告诉我我做错了什么吗?我使用cudamallocpitch在设备上分配2d,然后使用cudamemcpy将数组复制回主机,但是每次运行它显示所有新值时,代码均无法正常工作。 任何人都可以解释我如何实现我的目标,即“将2维数组的元素相加并将结果放入2d”。 请忽略我的英语,并帮助我。     

解决方法

        该代码有太多错误,几乎不可能构成一个合理的答案。但是,没有特别的顺序,我可以看到以下内容: 设备内存中存在多个缓冲区溢出,主要是由于使用25个线程写入分配了4或5个字的数组 您的内核包含一个可怕的内存争用。线程将覆盖另一个人的输出,这会在CUDA中提供未​​定义的行为。 您内核中的while循环将 产生无限循环行为 而且内核可能永远不会终止 结果是。 您一直在询问2D阵列,但是 代码中没有任何地方是2D数组。 一切都是线性一维存储器 我所看到的。 那应该让你忙一阵子。 编辑:改天,代码的另一个完全不同的版本,所以一些注释: 您的内核现在完全是串行的。前6个线程将全部尝试覆盖具有不同结果的输出数组。这是一场内存竞赛。 您正在内核中混合
char
unsigned short
数据类型。这可能意味着您要在8到16位值之间执行and运算,然后将结果存储为8位类型。即使没有内存争用,我也可以肯定这不是您想要的。 大约第十次,我将重复自己的内容-ѭ3分配带有填充的线性内存,使其适合GPU内存控制器和纹理单元。它没有分配2D数组。因此,如果要访问该线性存储器,则需要使用
cudaMallocPitch
调用返回的音高。文档中提供了访问内存的公式,但是我将在这里重复:
T* pElement = (T*)((char*)BaseAddress + Row * pitch) + Column;
程序中的每个API调用都会返回一个错误代码。您应该检查它们中的每一个,以查看是否检测到任何运行时错误。