系统结构实验二_12121469.docx

计算机系统结构 实验报告 专业：计算机科学与技术 指导老师： 雷咏梅 姓名： 钱成浩 学号： 计算机工程与科学学院 实验二 多核环境下OpenMP并行编程 实验2-1. OpenMP程序的编译和运行 实验目的 1) 在Linux平台上编译和运行OpenMP程序； 2) 在Windows平台上编译和运行OpenMP程序。 3) 掌握OpenMP并行编程基础。 实验环境 1) 硬件环境：i5双核四线程pc机 2) 软件环境：centos6+gcc，以及win8.1+Visual Studio2013 运行内容 代码如下： #include omp.h #include stdio.h int main() { int nthreads,tid; omp_set_num_threads(8); #pragma omp parallel private(nthreads,tid) { tid=omp_get_thread_num(); printf(Hello World from OMP thread %d\n,tid); if(tid==0) { nthreads=omp_get_num_threads(); printf(Number of threads is %d\n,nthreads); } } } 运行结果 linux Windows 实验2-2 矩阵乘法的OpenMP实现及性能分析 实验目的 1) 用OpenMP实现最基本的数值算法“矩阵乘法” 2) 掌握for编译制导语句 3) 对并行程序进行简单的性能调优 实验要求 运行并测试OpenMP编写两个n阶的方阵a和b的相乘程序，结果存放在方阵c中，其中乘法用for编译制导语句实现并行化操作，并调节for编译制导中schedule的参数，使得执行时间最短。要求在window环境（不用虚拟机），在linux环境（用和不用虚拟机情况下）测试程序的性能，并写出详细的分析报告。 3. 实验内容 1. 用OpenMP编写两个n阶的方阵a和b的相乘程序，结果存放在方阵c中，其中乘法用for编译制导语句实现并行化操作，并调节for编译制导中schedule的参数，使得执行时间最短，写出代码。 方阵a和b的初始值如下： 输入： 方阵的阶n、并行域的线程数 输出： c中所有元素之和、程序的执行时间 提示： a,b,c的元素定义为int型，c中所有元素之和定义为long long型。 Windows计时: 用time.h中的clock_t clock( void )函数得到当前程序执行的时间 在windows下使用Microsofe Visual Studio编程，源代码如下： #includestdio.h #includeomp.h #includetime.h void comput(float* A, float* B, float* C)//两个矩阵相乘传统方法 { int x, y; for (y = 0; y4; y++) { for (x = 0; x4; x++) { C[4 * y + x] = A[4 * y + 0] * B[4 * 0 + x] + A[4 * y + 1] * B[4 * 1 + x] + A[4 * y + 2] * B[4 * 2 + x] + A[4 * y + 3] * B[4 * 3 + x]; } } } int main() { double duration; clock_t s, f; int x = 0; int y = 0; int n = 0; int k = 0; float A[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 }; float B[] = { 0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, 0.8f, 0.9f, 0.10f, 0.11f, 0.12f, 0.13f, 0.14f, 0.15f, 0.16f }; float C[16]; s = clock(); //#pragma omp parallel if(false) for (n = 0; n1000000; n++) { comput(A, B, C); } f = clock(); duration = (double)(f - s) / CLOCKS_PER_SEC; printf(s---1,000,000 :%f\n, duration); fo