多核行列式.docVIP

一、代码： 非并行： #includestdio.h #includewindows.h #includetime.h #includeiostream.h #define p1 100 void main() { clock_t start, stop; start = clock(); int a1[p1][p1],a2[p1][p1],a3[p1][p1]; for(int i=0;ip1;i++) { for(int j=0;jp1;j++) { a1[i][j]=(i+j); a2[i][j]=(i+j); a3[i][j]=0; } } for(int m=0;mp1;m++) { for(int n=0;np1;n++) { for(int z=0;zp1;z++) { a3[m][n]+=a1[m][z]*a2[z][n]; } } } for(int p=0;pp1;p++) { for(int q=0;qp1;q++) // printf(%d/n,a3[p][q]); couta3[p][q] ; coutendl; } stop = clock(); cout运行时间为：((double)(stop - start)/1000.0)秒endl; } 并行代码： #include stdio.h #include windows.h #include time.h #includeiostream.h #define p1 100 const int gNumThreads = 15; int a1[p1][p1],a2[p1][p1],a3[p1][p1]; CRITICAL_SECTION gCS; DWORD WINAPI threadFunction(LPVOID pArg) { int myNum = *((int *)pArg); for ( int i=myNum; ip1; i+=gNumThreads ) { for(int n=0;np1;n++) { for(int z=0;zp1;z++) { a3[i][n]+=a1[i][z]*a2[z][n]; } } } return 0; } void main() { HANDLE threadHandles[gNumThreads]; int tNum[gNumThreads]; clock_t start, stop; start = clock(); for(int i=0;ip1;i++) { for(int j=0;jp1;j++) { a1[i][j]=(i+j); a2[i][j]=(i+j); a3[i][j]=0; } } InitializeCriticalSection(gCS); for ( int ii=0; iigNumThreads; ++ii ) { tNum[ii] = ii; threadHandles[ii] = CreateThread( NULL, 0, threadFunction, (LPVOID)tNum[ii], 0, NULL); } WaitForMultipleObjects(gNumThreads, threadHandles, TRUE, INFINITE); DeleteCriticalSection(gCS); for(int p=0;pp1;p++) { for(int q=0;qp1;q++) couta3[p][q] ; coutendl; } stop = clock(); cout运行时间为：((double)(stop - start)/1000.0)秒endl; } 二、设计思想： 非多线程：根据行列式的特点，两个行列式相乘所形成的新行列式取第一个的行第二个的列，则我们可以先用两个循环控制行列式的行和列。最后再用一个循环控制进行将所有乘积相加即可得到新行列式的一个值。 多线程：由于行列式每行的计算互不影响，由此我们可以根据其这个特点将其以行来做成多线程，执行时由线程控制行，列不变，每行计算时不再控制行的增加。 三、运行结果： p1=100 gNumThreads = 15 非多线程： 多线程：