8.3.3MKL数学核心函数库实验.pdfVIP

8.3.3 MKL数学核心函数库实验 本实验学习如何使用英特尔数学核心函数库中的例程以及多线程能力。 （1）矩阵乘法 本实验通过分别用原始运算、DDOT、DGEMV、DGEMM 编译和运行一个矩阵乘法程序，演示了 原始计算、一级 BLAS、二级 BLAS、三级 BLAS的性能特征和它们之间的差异。请按照以下实 验步骤执行： 定位到文件夹 code\ MKL_Overview\DGEMM，打开文件 mkl_lab_solution.c，仔细查看 4段 不同方法实现同一矩阵乘法的代码； 检查文件 Makefile中的 include和 library路径，修改为与本机一致； 打开开始菜单-所有程序-Intel? Software Development Tools-Intel? C++ Compiler 9.0-Build Environment for IA-32 Applications，输入：nmake，编译完 成后目录下会生成一个 matrix.exe文件，输入：matrix.exe运行程序，会提示输入矩 阵的维数（输入的数字尽量大些，如 1000），记录运行结果并观察它们的差异： roll_your_own: ___________ DDOT: __________ DGEMV: ___________ DGEMM: ___________ MKL默认使用单线程，可以通过改变环境变量 OMP_NUM_THREADS来自定义 MKL使用的线程数 量，输入：set OMP_NUM_THREADS=2，重新运行程序观察结果； 改变不同线程数量，记录各自的运行结果。 问题： 当自定义的线程数量超过处理器核的数量时，会发生什么现象？ （2）蒙特卡罗法计算圆周率 本实验通过使用 MKL中矢量统计库（VSL）的随机数生成器来修改经典的圆周率计算法—— 蒙特卡罗法，在实验中会用到 VSL的多线程能力。请按照以下实验步骤执行： 定位到文件夹 code\ MKL_Overview\MonteCarloPi，打开文件 pimonte.c。理解函数 rand（） 实现的原理； 问题： 代码中的循环体能用多线程实现吗？ 打开文件 pimonte_VSL.c，观察用库函数调用替代函数 rand（）所产生的变化； 问题： 为什么不是 1：1的直接替换 rand（）函数？ 设置 blocksize的目的是什么？ 参数 BRNG 和 VSL_BRNG_MCG31的作用是什么？ 对这个程序来说，它们是最好的选择吗？ 能用多线程实现吗？ 检查文件 Makefile.VSL中的 include和 library路径，修改为与本机一致； 先编译运行原始算法程序，输入：nmake –f Makefile.rand，生成文件 pimonte_rand.exe， 运行它并记录结果；再编译运行 VSL改进算法程序，输入：nmake –f Makefile.VSL，生成 文件 pimonte_VSL.exe，运行它并记录结果； 改变环境变量 OMP_NUM_THREADS，设置不同的线程数量，观察结果有什么变化。