Padding in Pitch memory

Accelerated Computing CUDA CUDA Programming and Performance
1

Hi,
I am having some trouble in understanding the way the padding works while using the pitch memory. I have implemented a small program using cudaMallocPitch() , but I am not able to access the entire matrix allocated in the pitch in the kernel.While printing out the values it looks like the threads are accessing the padding instead of the matrix. Could someone explain where I am going wrong?

The program:

global void mat_add(int a_d,size_t pitch_a,int width,int height)
{
int i = blockIdx.x10 + threadIdx.x;
int c;

int row=(int)((char*)a_d+i*pitch_a);
for(c=0;c<width;c++){
row[c]=row[c]+1;
printf(“i:%d val:%d\n”,i,row[c]);

}
}

int main(void)
{
int **a_h,*a_d;
int i,j;
size_t pitch_a;
int width=10;
int height=30;

a_h=(int **)calloc(height,sizeof(int *));
for(i=0;i<height;i++)
a_h[i]=(int *)calloc(width,sizeof(int));

// memory allocation for device
cudaMallocPitch((void**)&a_d,&pitch_a,width*sizeof(int),height);

printf(" pitch a %d \n",pitch_a);

// assigning values to host array
for(i=0;i<height;i++)
{
for(j=0;j<width;j++)
{
a_h[i][j] =10;
}
}
// copy data from host to device
cudaMemcpy2D(a_d,pitch_a,a_h,sizeof(int)*width,sizeof(int)*w
idth,height,cudaMemcpyHostToDevice);

// call device function
mat_add<<<3,10>>>(a_d,pitch_a,width,height);

free(a_h);

cudaFree(a_d);

}

Thanks,
Sudha

2

Any luck with this one yet?

3

Firstly, you can not copy a pointer of pointer to device memory. Memory is allocated with dynamic method will not allocate with a continuous space so when copy this data to device memory you may not get correct data. I recommend you should use the static allocate method to allocate 2D matrix on host memory. in case you had to use dynamic allocate method, you should allocate 1D matrix on host memory.

I saw many problems in your kernel function.

you had allocated a 2D array on device memory with 10X3, and you had created 10X3 threads ( 3 blocks and 10 threads per block), it means that you want 1 thread processes 1 element of you 2D array, is it right?

if yes, the indexing in the kernel function is not correct, and you don’t need to use the for loop inside kernel.

I had changed it.

[codebox]global void mat_add(int *a_d,size_t pitch_a,int width,int height)

{

int row=(int)((char*)a_d + blockIdx.x * pitch_a);

row[threadIdx.x] = row[threadIdx.x] + 1;

}

}[/codebox]

you can not use the printf() function to see the data allocated in device memory, you had to copy data back to host memory and then using printf() function to see your result.