(CUDA 编程9).CUDA shared memory使用------GPU的革命

/********************************************************************

*shared_memory_test.cu

*ThisisaexampleoftheCUDAprogram.

*Author:zhao.kaiyong(at)gmail.com

*'target='_blank'>http://blog.csdn.net/openhero[/code]
*'target='_blank'>http://www.comp.hkbu.edu.hk/~kyzhao/[/code]
*********************************************************************/

#include<stdio.h>

#include<stdlib.h>

#include<cutil.h>

#include<cutil_inline.h>

//1,2,3,4,5,6,7.....16

#defineBANK_CONFLICT16


#defineTHREAD_SIZE16




/************************************************************************/

/*static*/

/************************************************************************/

__global__voidshared_memory_static(float*result,intnum,float*table_1)

{


__shared__floatsh_data[THREAD_SIZE];


intidx=threadIdx.x;

floatret=0.0f;


sh_data[idx]=table_1[idx];


for(inti=0;i<num;i++)

{

ret+=sh_data[idx%BANK_CONFLICT];

}


result[idx]=ret;


}


/************************************************************************/

/*dynamic*/

/************************************************************************/

extern__shared__chararray[];


__global__voidshared_memory_dynamic(float*result,intnum,float*table_1,intshared_size)

{


float*sh_data=(float*)array;//这里就让sh_data指向了sharedmemory的第一个地址，就可以动态分配空间

float*sh_data2=(float*)&sh_data[shared_size];//这里的shared_size的大小为sh_data的大小；

intidx=threadIdx.x;

floatret=0.0f;


sh_data[idx]=table_1[idx];


for(inti=0;i<num;i++)

{

ret+=sh_data[idx%BANK_CONFLICT];

}


result[idx]=ret;


}


/************************************************************************/

/*Bankconflict*/

/************************************************************************/

__global__voidshared_memory_bankconflict(float*result,intnum,float*table_1)

{


__shared__floatsh_data[THREAD_SIZE];


intidx=threadIdx.x;

floatret=0.0f;


sh_data[idx]=table_1[idx];


for(inti=0;i<num;i++)

{

ret+=sh_data[idx%BANK_CONFLICT];

}


result[idx]=ret;


}


/************************************************************************/

/*HelloCUDA*/

/************************************************************************/

intmain(intargc,char*argv[])

{


if(cutCheckCmdLineFlag(argc,(constchar**)argv,"device"))

{

cutilDeviceInit(argc,argv);

}else

{

intid=cutGetMaxGflopsDeviceId();

cudaSetDevice(id);

}


float*device_result=NULL;

floathost_result[THREAD_SIZE]={0};


CUDA_SAFE_CALL(cudaMalloc((void**)&device_result,sizeof(float)*THREAD_SIZE));



float*device_table_1=NULL;

floathost_table1[THREAD_SIZE]={0};


for(inti=0;i<THREAD_SIZE;i++)

{

host_table1[i]=rand()%RAND_MAX;

}

CUDA_SAFE_CALL(cudaMalloc((void**)&device_table_1,sizeof(float)*THREAD_SIZE));

CUDA_SAFE_CALL(cudaMemcpy(device_table_1,host_table1,sizeof(float)*THREAD_SIZE,cudaMemcpyHostToDevice));




unsignedinttimer=0;

CUT_SAFE_CALL(cutCreateTimer(&timer));

CUT_SAFE_CALL(cutStartTimer(timer));


shared_memory_static<<<1,THREAD_SIZE>>>(device_result,1000,device_table_1);

//shared_memory_dynamic<<<1,THREAD_SIZE>>>(device_result,1000,device_table_1,16);

//shared_memory_bankconflict<<<1,THREAD_SIZE>>>(device_result,1000,device_table_1);

CUT_CHECK_ERROR("Kernelexecutionfailed");


CUDA_SAFE_CALL(cudaMemcpy(host_result,device_result,sizeof(float)*THREAD_SIZE,cudaMemcpyDeviceToHost));

CUT_SAFE_CALL(cutStopTimer(timer));

printf("Processingtime:%f(ms)",cutGetTimerValue(timer));

CUT_SAFE_CALL(cutDeleteTimer(timer));


for(inti=0;i<THREAD_SIZE;i++)

{

printf("%f",host_result[i]);

}


CUDA_SAFE_CALL(cudaFree(device_result));

CUDA_SAFE_CALL(cudaFree(device_table_1));


cutilExit(argc,argv);


}

这里只是一个简单的demo，大家可以测试一下。下一章节会将一些sharedmemory的更多的特性，更深入的讲解sharedmemory的一些隐藏的性质；

再在接下来的章节会讲一些constant和texture的使用；

写的内容一直都是文字比较多，代码比较少，其实学习的过程更重要的思想，实践的代码，最好是自己写，唯一可以学习的是思想，学习更重要的也是思想的交流，知识的传播，最好的是思想的传播，代码，方法，都是只是一些工具而已。但是工具的熟练层度，就得靠自己下来多练习。