cudaLab1矩阵相乘.pdfVIP

Lab1 矩阵相乘 1. 目标 熟悉基本的CUDA 程序架构以及如何调用相应的API 进行CUDA 编程 2. 过程 请结合PPT 中的内容，完成矩阵相乘的并行程序的实现（不用share memory 实现） 要求： 1）实现2 个矩阵（32*32 ）的相乘，M 矩阵的初始值全为2，N 矩阵的初始值全为5 。 同时用CPU 代码实现，比较两个代码的运行时间 2 ）实现2 个矩阵（1024*1024 ）的相乘，M 矩阵的初始值全为2，N 矩阵的初始值全为 5 。同时用CPU 代码实现，比较两个代码的运行时间 3 ）实现2 个矩阵（10000*10000 ）的相乘，M 矩阵的初始值全为2，N 矩阵的初始值全 为5 。同时用CPU 代码实现，比较两个代码的运行时间 vim matrixMulD.cu Nvcc -arch sm_20 –o matrixD matrixMulD.cu ./matrix 32 CPU 实现代码为 /** * C = A * B. */ #include stdio.h #include time.h #include malloc.h /** * Host code * */ void MatrixMulOnHost(float* M, float* N, float* P, int Width) { int i,j,k; double sum,a,b; for (i = 0; i Width; ++i) for (j = 0; j Width; ++j) { sum = 0; for (k = 0; k Width; ++k) { a = M[i * Width+ k]; b = N[k * Width+ j]; sum += a * b; } P[i * Width + j] = sum; } } /** * Host main routine */ int main(int argc,char * argv[]) { int Width=0; if(argc==1) { printf(Input the width\n); exit(1); } Width=atoi(argv[1]); int numElements,size,i; float *h_A=NULL,*h_B=NULL,*h_C=NULL,time1; clock_t start,end; start=clock(); numElements= Width*Width; size = numElements * sizeof(float); printf([matrix multiplication of two %d by %d matrix]\n,Width,Width); h_A = (float *)malloc(size); h_B = (float *)malloc(size); h_C = (float *)malloc(size); for (i = 0; i numElements; ++i) { h_A[i] = 2.0; h_B[i] = 5.0; } MatrixMulO