并行计算实验报告.pdfVIP

并行计算上机实验报告 课程名称 ：并行计算 实验名称 ：并行计算第一次实验 姓名 ：雷远 学号 ：P 上机日期 ：2015/5/7 一、实验题目 1.针对教材中求PI 的实例程序，请给出至少两种不同并行方式的OpenMP 实现。（划分数 = 1, 000, 000 ） 2.使用OPENMP 编写矩阵乘法程序。 3.用OpenMP 实现2 到1, 000, 000 素数的求解，并把最后的结果输出到一个文件中，同 时在屏幕中显示寻找到的素数的个数。（尽量尝试用利于并行的高效的确定性算法实现）。 二、实验环境（软、硬件平台） 处理器：Intel Core i5-3210 4 核 2.8GHz 安装内存(RAM) ：8.00GB 系统类型：Ubuntu 14.04.1 LTS 编译器：g++ 4.8.1 openmp 三、基本思想 1.给出计算PI 的近似公式，我们编写程序来模拟这一过程即可。串行程序不再说明思想 ， N 越大计算的误差越小。实验要求使用至少实现两种并行方式,我选择以下两种： (1)并行域并行化：如1-10 的循环，有两个线程，那么第一个线程计算迭代1 ，第二个线 程计算迭代2 ，第一个再处理迭代3…… ，在并行域外进行汇总即可。 (2)共享任务结构：如1-10 的循环，有两个线程，那么第一个线程处理1-5 ，第二线程处 理2-10 ，最后依然在并行域外进行汇总。 2.矩阵乘法采用了共享任务结构的并行方式，使用schedule 子句动态调度任务给不同的 线程，在两个线程下跑出了比较好的加速比。 3.这里采用的想法是这样:求1000000 以内的素数，通常使用的串行算法效率最高的就是 筛法了，但是普通的埃斯托尼斯筛法因为一个数会被他的所有质因子筛到,重复较多.时间复 杂度是O(n*loglogn).这里我们采用线性筛算法，我们可以在O(n)的时间内求出素数,保证每个 数只被他最小的质因子筛到。为了方便并行，我使用线性筛求出前sqrt(1000000) = 1000 的素 数，然后利用这些素数去排除剩余数中的合数，我们只要枚举这些质数的倍数即可.此时便可 以利用共享任务结构的办法去求解sqrt(1000000)=1000 之后的素数了，最后的并行时间应该 是两部分时间的总和. 四、运行结果与分析 (给出不同线程数的加速比图或加速比表格) 加速比=串行时间/并行时间（所有串行时间直接给出平均值） 1. 分析:由图表可以看出，规约并行还是得到了比较理想的加速比，在四核的环境下， 线程数为5 的时候加速比变低了，猜测和任务调度有关，线程数在5-8 区间内加速比成线 性增长 ，但是始终不会超过线程数为4 的加速比。 2. 分析:加速比随着线程数的增加而线性增长 ，超过4 之后基本不变，这可能就是共享 任务的特点 ，不会像规约并行那样线程数超过4 后加速比反倒是下降了. 3. 分析:这样的加速比, 是符合预测, 很好的. 附 ：（实验源码 ） 1 pi.cpp #include cstdio #include cstdlib #include cstring #include cctype #include cmath #include ctime #include cassert #include climits #include iostream #include iomanip #include string #include vector #include set #include map #include queue #include deque #include bitset #include algorithm #include omp.h #define MST(a, b) memset(a, b, sizeof(a)) #define REP(i, a) for (int i = 0; i int(a); i++) #define REPP(i, a, b) for (int i = int(a); i = int(b); i++) #define NUM_THREADS 4 using namespace std; const i