对存储器山和矩阵分块算法的探讨.pdfVIP

对存储器山和矩阵分块算法的探讨 曾虹驰 郭悦 [摘要]本文对课堂上提及的存储器山的测试代码作了进一步的解释说明，并利用该代码绘 制出了笔者本机的存储器山的图示。此外，对于矩阵分块算法，本文修正了课堂上给出的 代码的错误，并进一步探讨了不同的分块方式对缓存不命中期望的影响。 [关键词]存储器山 矩阵分块 缓存不命中期望 一、存储器山（Memory Mountain） 1、测试代码 1 对于ICS 课本P656 的测试代码的完整版 ，包含如下文件： 而课本上的代码节选自其中的mountain.c 文件，对其进一步说明如下： 第1 行的MAXELEMS 定义为： 25 15 5 2 MAXBYTE 的值为2 B= 2 KB = 2 MB = 32MB，MAXELEMS 为data 数组的长度 。 1 代码下载地址：/public/mountain.tar 2 根据代码的注释，MAXBYTES 为64MB，但是将1 左移25 位的实际值只有32MB，疑为源代码标注 错误。 这一段以传入的参数stride 为步长依次读取data 数组中的元素，并且通过将sink 放在 寄存器中来尽可能避免编译器的一些优化指令。 24 行传入的参数Mhz 的定义为： 即Mhz 参数在运行时会被设置为当前机器的CPU 时钟频率 （主频）。 30 行fcyc2 函数在头文件中的定义为： 即30 行返回test(elems, stride)函数运行所用的CPU 时钟周期数。故31 行中“cycles / 3 6 6 Mhz”的计算结果即为函数运行所用的时间 （单位：10 s ）。注意到时间单位为 10 s，这是 在因为在该程序中，参数 Mhz 的值是直接以兆赫兹作为单位的（例如，对于 CPU 主频为 3 一个程序所花的CPU 时间为T = N / f，其中N 表示CPU 时钟周期数，f 为CPU 主频。 /view/91a36cd2b14e852458fb5772.html 2500MHz 的机器，该参数的值会被设置为2 500，而不是2 500 000 000 ），这使得cycles / Mhz = cycles / (hz * 10-6 ) = (cycles / hz) * 106 = t * 106 。又因前半部分size / stride 的单位是Byte ，将 6 6 4 “10 ”与之结合，即可得到(10 * size / stride) / (cycles / hz) 的单位是MB/s （兆字节/秒）。 2、运行结果 5 运行mountain.c ，将结果输出到文件，结果节选如下： 第 4 行表示在工作集大小为 32M 时，步长从 s1 递增到 s64 的读吞吐量测试结果，同 理，第5 行表示工作集为16MB 时的结果。 利用 Excel 对数据处理，可以得到笔者的两台均使用 Intel Core i5-3210M 处理器的机器 的存储器山，如图所示： 4 唐渊在上课时对于返回值的单位是 “MB/s”并未给出非常清晰的解释。 5 测试时没有在BIOS 中锁定主频，只是简单地在虚拟机中运行程序。此外，测试过程中采用Oracle VirtualBox 运行程序会出现不合理的错误结果，而使用VMware Workstation 运行程序才能得出合理结果。 因而下文中所有的数据均为VMware 下的测试结果。 大致估计，可知红--蓝区域为主存（16MB 附近），黄—红区域为三级缓存 （4MB 附近）， 紫--黄区域 （256KB 附近）为二级缓存，最右侧区域 （16KB 附近）为一级缓存。