实验1 寻找完数的OpenMP并行程序.docxVIP

深圳大学实验报告课程名称：并行计算实验名称：寻找完数的OpenMP并行程序姓名：学号：班级：软件工程实验日期：2014.11.19一.实验目的1) 用OpenMP编写基本的并行程序2) 掌握for编译制导语句3) 对并行程序进行简单的性能分析二.实验环境1) 硬件环境：32核CPU、32G内存计算机；2) 软件环境：Microsoft Visual Studio 2008；3) 登录方式：通过远程桌面(mstsc)连接hpc.szu.edu.cn:804，用户名和初始密码都是自己的学号。三.实验内容1) 实验要求：用OpenMP语言编写程序，求小于等于n的所有完数。程序必须采用原始方法依次判断每个小于等于n的数是否完数，不可利用其他高阶原理。2) 评分标准：首先要保证程序的运行结果正确，否则不得分，然后主要根据执行时间、加速比和分析来评分。3)参考函数setint FinishedNumber(intn)，完成以下程序中的函数setint FinishedNumber(intn, intp)，其中n代表一个整数，p代表线程数，该函数返回小于等于n的所有完数的集合。#includeiostream#includemath.h#includevector#includeset#includeutility#includealgorithm#includetime.husingnamespace std;setint FinishedNumber(intn){setint s;for (inti = 2; i = n; i++){intmaxv = (int)(sqrt(double(i)) + 1);int sum = 1;for (int j = 2; jmaxv; j++)if (i%j == 0)sum += j + i / j;if (sum == i)s.insert(i);}return s;}setint FinishedNumber(intn, intp){setint s;return s;}int main(){vectorpairint, int vTest;for (inti = 1000000; i = 5000000; i += 1000000)for (int j = 1; j = 32; j*=2)vTest.push_back(pairint, int(i, j));for (vectorpairint, int ::iterator iter = vTest.begin(); iter != vTest.end(); ++iter){clock_t pt = 0;pt -= clock();setint ps = FinishedNumber(iter-first, iter-second);pt += clock();clock_t st = 0;st -= clock();setint ss = FinishedNumber(iter-first);st += clock();cout iter-first , iter-second :;/*for (setint::iteratoriter = ps.begin(); iter != ps.end(); ++iter){cout *iter;}*/cout , execution time: pt / 1000.0;cout , speedup: st*1.0/pt endl;if (ps != ss){cout error! endl;}}return 0;}4)测试并行程序并分析其执行时间和加速比表1 执行时间(秒)节点数n1248163210000002000000300000040000005000000表2加速比节点数n1248163210000002000000300000040000005000000