高性能计算课程设计讲义.docxVIP

高性能课程设计 问题描述 用OpenMP完成查找MAX_VALU以内的所有“完全数”的程序，输出完数及其个数。并且将单线程串行计算与多线程并行计算进行比较，结合所选的硬件的平台对实验结果进行分析。 算法描述 概念：若一个数恰好等于其所有的真因子（即除了自身以外的 HYPERLINK /view/461750.htm \t _blank 约数）的和，则称该数为“完全数”。 举例：前5个完全数分别为6、28、496、8128 6 = 1+2+3、28 =1+2+3+4+5+6+7 算法设计：设计一个双重循环，外层循环次数为MAX_VALU，即对每一个小于MAX_VALU的数index进行验证。由于一个数index除了其本身外没有比index/2+1大的约数，所以内层循环次数设置为index/2+1，对于内层循环的每一个i值，若index能够被i整除，则说明i是index的因子。计算index的所有真因子的和s，若s与index相等，则index为“完全数”。该算法程序代码如下： for(index=1;index=MAX_VALUE;index++) { //外层循环MAX_VALUE次 int s=0; for(int i=1;i index/2+1;i++){ //内层循环 index/2+1次 if(index % i==0){ //判断i是否是index的约数 s = s + i; //计算index的因子的和 } } if(s == index) { //若和与该数相等，则index是“完全数” printf(“%d是完全数\n”,index); } } 算法流程图 程序流程图 OpenMP代码 #include stdafx.h #include stdio.h #include omp.h #include time.h #define MAX_VALUE 10000 #define THREAD 2 int main() { omp_set_num_threads(THREAD ); int index= 0; double startTime = 0.0; double endTime = 0.0; //串行部分 startTime = clock(); //开始计时 int count1 = 0; for(index=1;index=MAX_VALUE;index++) { //外层循环MAX_VALUE次 int s=0; for(int i=1;i index/2+1;i++){ //内层循环 index/2+1次 if(index % i==0){ //判断i是否是index的约数 s = s + i; //计算index的因子的和 } } if(s == index) { //若和与该数相等，则是“完全数”，个数加1并输出 count1++; printf(%d ,index); } } endTime = clock(); //串行计算结束 printf(\n串行：THREAD = %d , MAX_VALUE = %d , count1 = %d ,time1 = %.10f\n, 1 ,MAX_VALUE ,count1,(endTime - startTime)/CLOCKS_PER_SEC); //并行部分 startTime = clock(); //开始计时 int count2 = 0,id=0; #pragma omp parallel private(id,index)reduction(+:count2) //并行计算 { id = omp_get_thread_num(); //获取当前线程编号 for(index=id+1; index = MAX_VALUE; index=index+THREAD) { // 循环分配 int s=0; for(int i=1;i index/2+1;i++){ //内层循环 index/2+1次 if(index % i==0){ //判断i是否是index的约数 s = s + i; //计算index的因子的和