内存层次优化相关的几种方法,,一些内存层次优化相关的知识
admin
2024-03-14 14:20:32
0

引言

软件开发人员往往期望计算机硬件拥有无限容量、零访问延迟、无限带宽以及便宜的内存,但是现实却是内存容量越大,相应的访问时间越长;内存访问速度越快,价格也更贵;带宽越大,价格越贵。为了解决大容量、高速度、低成本之间的矛盾,基于程序访问的局部性原理,将更常用数据放在小容量的高速存储器中,多种速度不同的存储器分层级联,协调工作。

图 1 memory hierarchy for sever [1]

现代计算机的存储层次可以分几层。如图 1 所示,位于处理器内部的是寄存器;稍远一点的是一级 Cache,一级 Cache 一般能够保存 64k 字节,访问它大约需要 1ns,同时一级 Cache 通常划分为指令 Cache(处理器从指令 Cache 中取要执行的指令)和数据 Cache(处理器从数据 Cache 中存 / 取指令的操作数);然后是二级 Cache,通常既保存指令又保存数据,容量大约 256k,访问它大约需要 3-10ns;然后是三级 Cache,容量大约 16-64MB,访问它大约需要 10-20ns;再接着是主存、硬盘等。注意,CPU 和 Cache 是以 word 传输的,Cache 到主存以块(一般 64byte)传输的。

前文提到了程序的局部性原理,一般指的是时间局部性(在一定时间内,程序可能会多次访问同一内存空间)和空间局部性(在一定时间内,程序可能会访问附近的内存空间),高速缓存(Cache)的效率取决于程序的空间和时间的局部性性质。比如一个程序重复地执行一个循环,在理想情况下,循环的第一个迭代将代码取至高速缓存中,后续的迭代直接从高速缓存中取数据,而不需要重新从主存装载。因此,为了使程序获得更好的性能,应尽可能让数据访问发生在高速缓存中。但是如果数据访问在高速缓存时发生了冲突,也可能会导致性能下降。

篇幅原因,本文重点讨论编译器在 Cache 优化中可以做哪些工作,如果读者对其他内存层次优化感兴趣,欢迎留言。下面将介绍几种通过优化 Cache 使用提高程序性能的方法。

对齐和布局

现代编译器可以通过调整代码和数据的布局方式,提高 Cache 命中率,进而提升程序性能。本节主要讨论数据和指令的对齐、代码布局对程序性能的影响,大部分处理器中 Cache 到主存是以 Cache line(一般为 64Byte,也有地方称 Cache 块,本文统一使用 Cache line)传输的,CPU 从内存加载数据是一次一个 Cache line,CPU 往内存写数据也是一次一个 Cache line。假设处理器首次访问数据对象 A,其大小刚好为 64Byte,如果数据对象 A 首地址并没有进行对齐,即数据对象 A 占用两个不同 Cache line 的一部分,此时处理器访问该数据对象时需要两次内存访问,效率低。但是如果数据对象 A 进行了内存对齐,即刚好在一个 Cache line 中,那么处理器访问该数据时只需要一次内存访问,效率会高很多。编译器可以通过合理安排数据对象,避免不必要地将它们跨越在多个 Cache line 中,尽量使得同一对象集中在一个 Cache 中,进而有效地使用 Cache 来提高程序的性能。通过顺序分配对象,即如果下一个对象不能放入当前 Cache line 的剩余部分,则跳过这些剩余的部分,从下一个 Cache line 的开始处分配对象,或者将大小(size)相同的对象分配在同一个存储区,所有对象都对齐在 size 的倍数边界上等方式达到上述目的。

Cache line 对齐可能会导致存储资源的浪费,如图 2 所示,但是执行速度可能会因此得到改善。对齐不仅仅可以作用于全局静态数据,也可以作用于堆上分配的数据。对于全局数据,编译器可以通过汇编语言的对齐指令命令来通知链接器。对于堆上分配的数据,将对象放置在 Cache line 的边界或者最小化对象跨 Cache line 的次数的工作不是由编译器来完成的,而是由 runtime 中的存储分配器来完成的 [2]。

图 2 因块对齐可能会浪费存储空间

前文提到了数据对象对齐,可以提高程序性能。指令 Cache 的对齐,也可以提高程序性能。同时,代码布局也会影响程序的性能,将频繁执行的基本块的首地址对齐在 Cache line 的大小倍数边界上能增加在指令 Cache 中同时容纳的基本块数目,将不频繁执行的指令和频繁指令的指令放到不同的 Cache line 中,通过优化代码布局来提升程序性能。

利用硬件辅助

Cache 预取是将内存中的指令和数据提前存放至 Cache 中,达到加快处理器执行速度的目的。Cache 预取可以通过硬件或者软件实现,硬件预取是通过处理器中专门的硬件单元实现的,该单元通过跟踪内存访问指令数据地址的变化规律来预测将会被访问到的内存地址,并提前从主存中读取这些数据到 Cache;软件预取是在程序中显示地插入预取指令,以非阻塞的方式让处理器从内存中读取指定地址数据至 Cache。由于硬件预取器通常无法正常动态关闭,因此大部分情况下软件预取和硬件预取是并存的,软件预取必须尽力配合硬件预取以取得更优的效果。本文假设硬件预取器被关闭后,讨论如何利用软件预取达到性能提升的效果。

预取指令 prefech (x) 只是一种提示,告知硬件开始将地址 x 中的数据从主存中读取到 Cache 中。它并不会引起处理停顿,但若硬件发现会产生异常,则会忽略这个预取操作。如果 prefech (x) 成功,则意味着下一次取 x 将命中 Cache;不成功的预取操作可能会导致下次读取时发生 Cache miss,但不会影响程序的正确性 [2]。

数据预取是如何改成程序性能的呢?如下一段程序:

double a[n];
for (int i = 0; i < 100; i++)a[i] = 0;

假设一个 Cache line 可以存放两个 double 元素,当第一次访问 a [0] 时,由于 a [0] 不在 Cache 中,会发生一次 Cache miss,需要从主存中将其加载至 Cache 中,由于一个 Cache line 可以存放两个 double 元素,当访问 a [1] 时则不会发生 Cache miss。依次类推,访问 a [2] 时会发生 Cache miss,访问 a [3] 时不会发生 Cache miss,我们很容易得到程序总共发生了 50 次 Cache miss。

我们可以通过软件预取等相关优化,降低 Cache miss 次数,提高程序性能。首先介绍一个公式 [3]:

上述公式中 L 是 memory latency,S 是执行一次循环迭代最短的时间。iterationAhead 表示的是循环需要经过执行几次迭代,预取的数据才会到达 Cache。假设我们的硬件架构计算出来的 iterationAhead=6,那么原程序可以优化成如下程序:

double a[n];
for (int i = 0; i < 12; i+=2)   //prologue prefetch(&a[i]);
for (int i = 0; i < 88; i+=2) { // steady stateprefetch(&a[i+12]);a[i] = 0;a[i+1] = 0;
}
for (int i = 88; i < 100; i++) //epiloguea[i] = 0;

由于我们的硬件架构需要循环执行 6 次后,预取的数据才会到达 Cache。一个 Cache line 可以存放两个 double 元素,为了避免浪费 prefetch 指令,所以 prologue 和 steady state 循环都展开了,即执行 prefetch (&a [0]) 后会将 a [0]、a [1] 从主存加载至 Cache 中,下次执行预取时就无需再次将 a [1] 从主存加载至 Cache 了。prologue 循环先执行数组 a 的前 12 个元素的预取指令,等到执行 steady state 循环时,当 i = 0 时,a [0] 和 a [1] 已经被加载至 Cache 中,就不会发生 Cache miss 了。依次类推,经过上述优化后,在不改变语义的基础上,通过使用预取指令,程序的 Cache miss 次数从 50 下降至 0,程序的性能将会得到很大提升。

注意,预取并不能减少从主存储器取数据到高速缓存的延迟,只是通过预取与计算重叠而隐藏这种延迟。总之,当处理器有预取指令或者有能够用作预取的非阻塞的读取指令时,对于处理器不能动态重排指令或者动态重排缓冲区小于我们希望隐藏的具体 Cache 延迟,并且所考虑的数据大于 Cache 或者是不能够判断数据是否已在 Cache 中,预取是适用的。预取也不是万能,不当的预取可能会导致高速缓存冲突,程序性能降低。我们应该首先利用数据重用来减少延迟,然后才考虑预取。

除了软件预取外,ARMv8 还提供了 Non-temporal 的 Load/Store 指令,可以提高 Cache 的利用率。对于一些数据,如果只是访问一次,无需占用 Cache,可以使用这个指令进行访问,从而保护 Cache 中关键数据不被替换,比如 memcpy 大数据的场景下,使用该指令对于其关键业务而言,是有一定的收益的。

循环变换

重用 Cache 中的数据是最基本的高效使用 Cache 方法。对于多层嵌套循环,可以通过交换两个嵌套的循环(loop interchange)、逆转循环迭代执行的顺序(loop reversal)、将两个循环体合并成一个循环体(loop fusion)、循环拆分(loop distribution)、循环分块(loop tiling)、loop unroll and jam 等循环变换操作。选择适当的循环变换方式,既能保持程序的语义,又能改善程序性能。我们做这些循环变换的主要目的是为了实现寄存器、数据高速缓存以及其他存储层次使用方面的优化。

篇幅受限,本节仅讨论循环分块(loop tiling)如何改善程序性能,若对 loop interchange 感兴趣,请点击查阅。下面这个简单的循环:

for(int i = 0; i < m; i++) {for(int j = 0; j < n; j++) {x = x+a[i]+c*b[j];}
}

我们假设数组 a、b 都是超大数组,m、n 相等且都很大,程序不会出现数组越界访问情况发生。那么如果 b [j] 在 j 层循环中跨度太大时,那么被下次 i 层循环重用时数据已经被清出高速缓存。即程序访问 b [n-1] 时,b [0]、b [1] 已经被清出缓存,此时需要重新从主存中将数据加载至缓存中,程序性能会大幅下降。

我们如何通过降低 Cache miss 次数提升程序的性能呢?通过对循环做 loop tiling 可以符合我们的期望,即通过循环重排,使得数据分成一个一个 tile,让每一个 tile 的数据都可以在 Cache 中被 hint [4]。从内层循环开始 tiling,假设 tile 的大小为 t,t 远小于 m、n,t 的取值使得 b [t-1] 被访问时 b [0] 依然在 Cache 中,将会大幅地减少 Cache miss 次数。假设 n-1 恰好被 t 整除,此时 b 数组的访问顺序如下所示:

i=1; b[0]、b[1]、b[2]...b[t-1]
i=2; b[0]、b[1]、b[2]...b[t-1]
...
i=n; b[0]、b[1]、b[2]...b[t-1]
...
...
...
i=1; b[n-t]、b[n-t-1]、b[n-t-2]...b[n-1]
i=2; b[n-t]、b[n-t-1]、b[n-t-2]...b[n-1]
...
i=n; b[n-t]、b[n-t-1]、b[n-t-2]...b[n-1]

经过 loop tiling 后循环变换成:

for(int j = 0; j < n; j+=t) {for(int i = 0; i < m; i++) {for(int jj = j; jj < min(j+t, n); jj++) {x = x+a[i]+c*b[jj];}}
}

假设每个 Cache line 能够容纳 X 个数组元素,loop tiling 前 a 的 Cache miss 次数为 m/X,b 的 Cache miss 次数是 m*n/X,总的 Cache miss 次数为 m*(n+1)/x。loop tiling 后 a 的 Cache miss 次数为 (n/t)*(m/X),b 的 Cache miss 次数为 (t/X)*(n/t)=n/X,总的 Cache miss 次数为 n*(m+t)/xt。此时,由于 n 与 m 相等,那么 loop tiling 后 Cache miss 大约可以降低 t 倍 [4]。

前文讨论了 loop tiling 在小用例上如何提升程序性能,总之针对不同的循环场景,选择合适的循环交换方法,既能保证程序语义正确, 又能获得改善程序性能的机会。

小结

汝之蜜糖,彼之砒霜。针对不同的硬件,我们需要结合具体的硬件架构,利用性能分析工具,通过分析报告和程序,从系统层次和算法层次思考问题,往往会有意想不到的收获。本文简单地介绍了内存层次优化相关的几种方法,结合一些小例子深入浅出地讲解了一些内存层次优化相关的知识。纸上得来终觉浅,绝知此事要躬行,更多性能优化相关的知识需要我们从实践中慢慢摸索。

相关内容

热门资讯

linux入门---制作进度条 了解缓冲区 我们首先来看看下面的操作: 我们首先创建了一个文件并在这个文件里面添加了...
C++ 机房预约系统(六):学... 8、 学生模块 8.1 学生子菜单、登录和注销 实现步骤: 在Student.cpp的...
A.机器学习入门算法(三):基... 机器学习算法(三):K近邻(k-nearest neigh...
数字温湿度传感器DHT11模块... 模块实例https://blog.csdn.net/qq_38393591/article/deta...
有限元三角形单元的等效节点力 文章目录前言一、重新复习一下有限元三角形单元的理论1、三角形单元的形函数(Nÿ...
Redis 所有支持的数据结构... Redis 是一种开源的基于键值对存储的 NoSQL 数据库,支持多种数据结构。以下是...
win下pytorch安装—c... 安装目录一、cuda安装1.1、cuda版本选择1.2、下载安装二、cudnn安装三、pytorch...
MySQL基础-多表查询 文章目录MySQL基础-多表查询一、案例及引入1、基础概念2、笛卡尔积的理解二、多表查询的分类1、等...
keil调试专题篇 调试的前提是需要连接调试器比如STLINK。 然后点击菜单或者快捷图标均可进入调试模式。 如果前面...
MATLAB | 全网最详细网... 一篇超超超长,超超超全面网络图绘制教程,本篇基本能讲清楚所有绘制要点&#...
IHome主页 - 让你的浏览... 随着互联网的发展,人们越来越离不开浏览器了。每天上班、学习、娱乐,浏览器...
TCP 协议 一、TCP 协议概念 TCP即传输控制协议(Transmission Control ...
营业执照的经营范围有哪些 营业执照的经营范围有哪些 经营范围是指企业可以从事的生产经营与服务项目,是进行公司注册...
C++ 可变体(variant... 一、可变体(variant) 基础用法 Union的问题: 无法知道当前使用的类型是什...
血压计语音芯片,电子医疗设备声... 语音电子血压计是带有语音提示功能的电子血压计,测量前至测量结果全程语音播报࿰...
MySQL OCP888题解0... 文章目录1、原题1.1、英文原题1.2、答案2、题目解析2.1、题干解析2.2、选项解析3、知识点3...
【2023-Pytorch-检... (肆十二想说的一些话)Yolo这个系列我们已经更新了大概一年的时间,现在基本的流程也走走通了,包含数...
实战项目:保险行业用户分类 这里写目录标题1、项目介绍1.1 行业背景1.2 数据介绍2、代码实现导入数据探索数据处理列标签名异...
记录--我在前端干工地(thr... 这里给大家分享我在网上总结出来的一些知识,希望对大家有所帮助 前段时间接触了Th...
43 openEuler搭建A... 文章目录43 openEuler搭建Apache服务器-配置文件说明和管理模块43.1 配置文件说明...