高性能计算导论课程设计讲述.docxVIP

高性能课程设计问题描述用OpenMP完成查找MAX_VALU以内的所有“完全数”的程序，输出完数及其个数。并且将单线程串行计算与多线程并行计算进行比较，结合所选的硬件的平台对实验结果进行分析。算法描述概念：若一个数恰好等于其所有的真因子（即除了自身以外的约数）的和，则称该数为“完全数”。举例：前5个完全数分别为6、28、496、81286 = 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代码#includestdafx.h#includestdio.h#includeomp.h#includetime.h#define MAX_VALUE 10000#define THREAD 2int 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的因子的和 }}if(s == index) {//若和与该数相等，则是“完全数”，个数加1并输出count2++; printf(%d ,index);}}} endTime = clock(); //并行计算结束 printf(\n并行：THREAD = %d , MAX_VALUE = %d , count2 = %d ,time1 = %.10f\n, THREAD ,MAX_VALUE ,count2,(endTime - startTime)/CLOCKS_PER_SEC);return 0;}并行化关键代码分析#pragma omp parallel private(id,index)redu