CUDA实现数组倒序

数组倒序，将在主机上初始化的数组传输到设备上，然后用CUDA并行倒序，此时在全局内存上操作，再将结果返回到主机并验证。

 #include <stdio.h>
 #include <assert.h>
 #include "cuda.h"
 #include "cuda_runtime.h"
 #include "device_launch_parameters.h"
 //检查CUDA运行时是否有错误
 void checkCUDAError(const char* msg);
 // Part3: implement the kernel
 /*
 blockDim块内的线程数
 blockIdx网格内的块索引
 gridDim网格内块个数
 threadIdx块内线程索引
 */
 __global__ void reverseArrayBlock(int *d_out, int *d_in)
 {
     int inOffset = blockDim.x * blockIdx.x;
     int outOffset = blockDim.x * (gridDim.x -  - blockIdx.x);
     int in = inOffset + threadIdx.x;
     int out = outOffset + (blockDim.x -  - threadIdx.x);
     d_out[out] = d_in[in];
 }
 /////////////////////////////////////////////////////////////////////
 //主函数
 /////////////////////////////////////////////////////////////////////
 int main(int argc, char** argv)
 {
     //指向主机的内存空间和大小
     int *h_a;
     int dimA =  * ; // 256K elements (1MB total)
     //指向设备的指针和大小
     int *d_b, *d_a;
     //定义网格和块大小，每个块的线程数量
     int numThreadsPerBlock = ;

     /*
     根据数组大小和预设的块大小来计算需要的块数
     */
     int numBlocks = dimA / numThreadsPerBlock;
     //申请主机及设备上的存储空间
     size_t memSize = numBlocks * numThreadsPerBlock * sizeof(int);
     //主机上的大小
     h_a = (int *)malloc(memSize);
     //设备上的大小
     cudaMalloc((void **)&d_a, memSize);
     cudaMalloc((void **)&d_b, memSize);
     //在主机上初始化输入数组
     for (int i = ; i < dimA; ++i)
     {
         h_a[i] = i;
     }
     //将主机数组拷贝到设备上，h_a-->d_a
     cudaMemcpy(d_a, h_a, memSize, cudaMemcpyHostToDevice);
     //启动内核
     dim3 dimGrid(numBlocks);
     dim3 dimBlock(numThreadsPerBlock);
     reverseArrayBlock <<< dimGrid,    dimBlock >>>(d_b, d_a);
     //阻塞，一直到设备完成计算
     cudaThreadSynchronize();
     //检查是否设备产生了错误
     //检查任何CUDA错误
     checkCUDAError("kernel invocation");
     //将结果从设备拷贝到主机，d_b-->h_a
     cudaMemcpy(h_a, d_b, memSize, cudaMemcpyDeviceToHost);
     //检查任何CUDA错误
     checkCUDAError("memcpy");
     //核对返回到主机上的结果是否正确
     for (int i = ; i < dimA; i++)
     {
         assert(h_a[i] == dimA -  - i);
     }
     //释放设备内存
     cudaFree(d_a);
     cudaFree(d_b);
     //释放主机内存
     free(h_a);
     printf("Correct!\n");
     return ;
 }
 void checkCUDAError(const char *msg)
 {
     cudaError_t err = cudaGetLastError();
     if (cudaSuccess != err)
     {
         fprintf(stderr, "Cuda error: %s: %s.\n", msg,cudaGetErrorString(err));
         exit(EXIT_FAILURE);
     }
 }