线性Attention的探索:Attention必须有个Softmax吗?
By 苏剑林 | 2020-07-04 | 253551位读者 |众所周知,尽管基于Attention机制的Transformer类模型有着良好的并行性能,但它的空间和时间复杂度都是O(n2)级别的,n是序列长度,所以当n比较大时Transformer模型的计算量难以承受。近来,也有不少工作致力于降低Transformer模型的计算量,比如模型剪枝、量化、蒸馏等精简技术,又或者修改Attention结构,使得其复杂度能降低到O(nlogn)甚至O(n)。
前几天笔者读到了论文《Transformers are RNNs: Fast Autoregressive Transformers with Linear Attention》,了解到了线性化Attention(Linear Attention)这个探索点,继而阅读了一些相关文献,有一些不错的收获,最后将自己对线性化Attention的理解汇总在此文中。
Attention #
当前最流行的Attention机制当属Scaled-Dot Attention,形式为
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V}) = softmax\left(\boldsymbol{Q}\boldsymbol{K}^{\top}\right)\boldsymbol{V}\label{eq:std-att}\end{equation}
这里的\boldsymbol{Q}\in\mathbb{R}^{n\times d_k}, \boldsymbol{K}\in\mathbb{R}^{m\times d_k}, \boldsymbol{V}\in\mathbb{R}^{m\times d_v},简单起见我们就没显式地写出Attention的缩放因子了。本文我们主要关心Self Attention场景,所以为了介绍上的方便统一设\boldsymbol{Q}, \boldsymbol{K}, \boldsymbol{V}\in\mathbb{R}^{n\times d},一般场景下都有n > d甚至n\gg d(BERT base里边d=64)。相关解读可以参考笔者的《Attention is All You Need》浅读(简介+代码),以及它的一些改进工作也可以参考《突破瓶颈,打造更强大的Transformer》、《Google新作Synthesizer:我们还不够了解自注意力》,这里就不多深入介绍了。
摘掉Softmax #
读者也许想不到,制约Attention性能的关键因素,其实是定义里边的Softmax!事实上,简单地推导一下就可以得到这个结论。\boldsymbol{Q}\boldsymbol{K}^{\top}这一步我们得到一个n\times n的矩阵,就是这一步决定了Attention的复杂度是\mathcal{O}(n^2);如果没有Softmax,那么就是三个矩阵连乘\boldsymbol{Q}\boldsymbol{K}^{\top}\boldsymbol{V},而矩阵乘法是满足结合率的,所以我们可以先算\boldsymbol{K}^{\top}\boldsymbol{V},得到一个d\times d的矩阵,然后再用\boldsymbol{Q}左乘它,由于d \ll n,所以这样算大致的复杂度只是\mathcal{O}(n)(就是\boldsymbol{Q}左乘那一步占主导)。
也就是说,去掉Softmax的Attention的复杂度可以降到最理想的线性级别\mathcal{O}(n)!这显然就是我们的终极追求:Linear Attention,复杂度为线性级别的Attention。所以,本文的主题就是探究摘掉Softmax后的线形Attention。
一般的定义 #
问题是,直接去掉Softmax还能算是Attention吗?它还能有标准的Attention的效果吗?为了回答这个问题,我们先将Scaled-Dot Attention的定义\eqref{eq:std-att}等价地改写为(本文的向量都是列向量)
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V})_i = \frac{\sum\limits_{j=1}^n e^{\boldsymbol{q}_i^{\top}\boldsymbol{k}_j}\boldsymbol{v}_j}{\sum\limits_{j=1}^n e^{\boldsymbol{q}_i^{\top}\boldsymbol{k}_j}}\label{eq:std-att-2}\end{equation}
所以,Scaled-Dot Attention其实就是以e^{\boldsymbol{q}_i^{\top}\boldsymbol{k}_j}为权重对\boldsymbol{v}_j做加权平均。所以我们可以提出一个Attention的一般化定义
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V})_i = \frac{\sum\limits_{j=1}^n \text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)\boldsymbol{v}_j}{\sum\limits_{j=1}^n \text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)}\label{eq:gen-att}\end{equation}
也就是把e^{\boldsymbol{q}_i^{\top}\boldsymbol{k}_j}换成\boldsymbol{q}_i, \boldsymbol{k}_j的一般函数\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j),为了保留Attention相似的分布特性,我们要求\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)\geq 0恒成立。也就是说,我们如果要定义新式的Attention,那么要保留式\eqref{eq:gen-att}的形式,并且满足\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)\geq 0。
这种一般形式的Attention在CV中也被称为Non-Local网络,出自论文《Non-local Neural Networks》。
几个例子 #
如果直接去掉Softmax,那么就是\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j) = \boldsymbol{q}_i^{\top}\boldsymbol{k}_j,问题是内积无法保证非负性,所以这还不是一个合理的选择。下面我们简单介绍几种可取的方案。
值得指出的是,下面介绍的这几种Linear Attention,前两种来自CV领域,第三种是笔者自己构思的,所以都还没有在NLP任务上做过什么实验,各位做模型改进的NLPer们就有实验方向了(^_^)~~顺便说一下,CV领域有不少对Attention的改进工作(除了下面介绍的外,还有EMANet等),很多内容都值得做NLP的我们参考阅读。
核函数形式 #
一个自然的想法是:如果\boldsymbol{q}_i,\boldsymbol{k}_j的每个元素都是非负的,那么内积自然也就是非负的。为了完成这点,我们可以给\boldsymbol{q}_i,\boldsymbol{k}_j各自加个激活函数\phi,\varphi,即
\begin{equation}\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j) = \phi(\boldsymbol{q}_i)^{\top} \varphi(\boldsymbol{k}_j)\label{eq:gen-att-2}\end{equation}
其中\phi(\cdot),\varphi(\cdot)是值域非负的激活函数。本文开头提到的论文《Transformers are RNNs: Fast Autoregressive Transformers with Linear Attention》选择的是\phi(x)=\varphi(x)=\text{elu}(x)+1。
非要讲故事的话,式\eqref{eq:gen-att-2}可以联想到“核方法(kernal method)”,尤其是\phi=\varphi时\phi就相当于一个核函数,而\langle \phi(\boldsymbol{q}_i), \phi(\boldsymbol{k}_j)\rangle就是通过核函数所定义的内积。这方面的思考可以参考论文《Transformer dissection: An unified understanding for transformer’s attention via the lens of kernel》,此处不做过多延伸。
妙用Softmax #
另一篇更早的文章《Efficient Attention: Attention with Linear Complexities》则给出了一个更有意思的选择。它留意到在\boldsymbol{Q}\boldsymbol{K}^{\top}中,\boldsymbol{Q}, \boldsymbol{K}, \in\mathbb{R}^{n\times d},如果“\boldsymbol{Q}在d那一维是归一化的、并且\boldsymbol{K}在n那一维是归一化的”,那么\boldsymbol{Q}\boldsymbol{K}^{\top}就是自动满足归一化了,所以它给出的选择是:
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V}) = softmax_2\left(\boldsymbol{Q}\right)softmax_1(\boldsymbol{K})^{\top}\boldsymbol{V}\end{equation}
其中softmax_1、softmax_2分别指在第一个(n)、第二个维度(d)进行Softmax运算。也就是说,这时候我们是各自给\boldsymbol{Q},\boldsymbol{K}加Softmax,而不是\boldsymbol{Q}\boldsymbol{K}^{\top}算完之后才加Softmax。
如果直接取\phi(\boldsymbol{q}_i)=softmax(\boldsymbol{q}_i),\varphi(\boldsymbol{k}_j)=softmax(\boldsymbol{k}_j),那么很显然这个形式也是式\eqref{eq:gen-att-2}的一个特例。另外这个设计在CV中出现过不止一次,比如A2-Nets也包含了同样的做法。
自己的构思 #
在这里,笔者给出自己的一种构思。这个构思的出发点不再是式\eqref{eq:gen-att-2},而是源于我们对原始定义\eqref{eq:std-att-2}的近似。由泰勒展开我们有
\begin{equation}e^{\boldsymbol{q}_i^{\top}\boldsymbol{k}_j} \approx 1 + \boldsymbol{q}_i^{\top}\boldsymbol{k}_j\end{equation}
如果\boldsymbol{q}_i^{\top}\boldsymbol{k}_j\geq -1,那么就可以保证右端的非负性,而从可以让\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)=1 + \boldsymbol{q}_i^{\top}\boldsymbol{k}_j。到这里读者可能已经想到了,想要保证\boldsymbol{q}_i^{\top}\boldsymbol{k}_j\geq -1,只需要分别对\boldsymbol{q}_i,\boldsymbol{k}_j做l_2归一化。所以,笔者最终提出的方案就是:
\begin{equation}\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j) = 1 + \left( \frac{\boldsymbol{q}_i}{\Vert \boldsymbol{q}_i\Vert}\right)^{\top}\left(\frac{\boldsymbol{k}_j}{\Vert \boldsymbol{k}_j\Vert}\right)\end{equation}
这不同于形式\eqref{eq:gen-att-2},但理论上它更加接近原始的Scaled-Dot Attention。
相关工作 #
通过修改Attention的形式来降低它的计算复杂度,相关的工作有很多,这里简要列举一些。
稀疏Attention #
我们之前介绍过OpenAI的Sparse Attention,通过“只保留小区域内的数值、强制让大部分注意力为零”的方式,来减少Attention的计算量。经过特殊设计之后,Attention矩阵的大部分元素都是0,因此理论上它也能节省显存占用量和计算量。后续类似工作还有《Explicit Sparse Transformer: Concentrated Attention Through Explicit Selection》、《Longformer: The Long-Document Transformer》等。
但是很明显,这种思路有两个不足之处:
1、如何选择要保留的注意力区域,这是人工主观决定的,带有很大的不智能性;
2、它需要从编程上进行特定的设计优化,才能得到一个高效的实现,所以它不容易推广。
Reformer #
Reformer也是有代表性的改进工作,它将Attention的复杂度降到了\mathcal{O}(n\log n)。某种意义上来说,Reformer也是稀疏Attention的一种,只不过它的稀疏Pattern不是事先指定的,而是通过LSH(Locality Sensitive Hashing)技术(近似地)快速地找到最大的若干个Attention值,然后只去计算那若干个值。此外,Reformer通过构造可逆形式的FFN(Feedforward Network)替换掉原来的FFN,然后重新设计反向传播过程,从而降低了显存占用量。
所以,相比前述稀疏Attention,Reformer解决了它的第一个缺点,但是依然有第二个缺点:实现起来复杂度高。要实现LSH形式的Attention比标准的Attention复杂多了,对可逆网络重写反向传播过程对普通读者来说更是遥不可及~
Linformer #
跟本文所介绍的Linear Attention很相似的一个工作是Facebook最近放出来的Linformer,它依然保留原始的Scaled-Dot Attention形式,但在进行Attention之前,用两个m\times n的矩阵\boldsymbol{E},\boldsymbol{F}分别对\boldsymbol{K},\boldsymbol{V}进行投影,即变为
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V}) = softmax\left(\boldsymbol{Q}(\boldsymbol{E}\boldsymbol{K})^{\top}\right)\boldsymbol{F}\boldsymbol{V}\end{equation}
这样一来,\boldsymbol{Q}(\boldsymbol{E}\boldsymbol{K})^{\top}就只是一个n\times m的矩阵,而作者声称对于哪怕对于很大的序列长度n,m也可以保持为一个适中的常数,从而这种Attention也是线性的。跟Linformer类似的思路还出现在更早一些的CV论文《Asymmetric Non-local Neural Networks for Semantic Segmentation》中。
但是,笔者认为“对于超长序列m可以保持不变”这个结论是值得质疑的,对于长序列原论文只做了MLM任务,而很明显MLM并不那么需要长程依赖,所以这个实验没什么说服力。因此,Linformer是不是真的Linear,还有待商榷。
自回归生成 #
Linformer的另一个缺点是\boldsymbol{E}\boldsymbol{K},\boldsymbol{F}\boldsymbol{V}这两个运算直接把整个序列的信息给“糅合”起来了,所以它没法简单地把将来信息给Mask掉(Causal Masking),从而无法做语言模型、Seq2Seq等自回归生成任务,这也是刚才说的原作者只做了MLM任务的原因。相比之下,本文介绍的几种Linear Attention都能做到这一点。以式\eqref{eq:gen-att}和式\eqref{eq:gen-att-2}为例,如果要Mask掉未来信息,那么只需要把求和\sum\limits_{j=1}^n改为\sum\limits_{j=1}^i:
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V})_i = \frac{\sum\limits_{j=1}^i \left(\phi(\boldsymbol{q}_i)^{\top} \varphi(\boldsymbol{k}_j)\right)\boldsymbol{v}_j}{\sum\limits_{j=1}^i \phi(\boldsymbol{q}_i)^{\top} \varphi(\boldsymbol{k}_j)}=\frac{ \phi(\boldsymbol{q}_i)^{\top} \sum\limits_{j=1}^i\varphi(\boldsymbol{k}_j)\boldsymbol{v}_j^{\top}}{ \phi(\boldsymbol{q}_i)^{\top} \sum\limits_{j=1}^i\varphi(\boldsymbol{k}_j)}\end{equation}
实现上式有两种方式:第一方式是设\boldsymbol{S}_i=\sum\limits_{j=1}^i\varphi(\boldsymbol{k}_j)\boldsymbol{v}_j^{\top}以及\boldsymbol{z}_i=\sum\limits_{j=1}^i\varphi(\boldsymbol{k}_j),我们有
\begin{equation}Attention(\boldsymbol{Q},\boldsymbol{K},\boldsymbol{V})_i =\frac{ \phi(\boldsymbol{q}_i)^{\top} \boldsymbol{S}_i}{ \phi(\boldsymbol{q}_i)^{\top} \boldsymbol{z}_i},\quad \begin{aligned}&\boldsymbol{S}_i=\boldsymbol{S}_{i-1}+\varphi(\boldsymbol{k}_i)\boldsymbol{v}_i^{\top}\\
&\boldsymbol{z}_i=\boldsymbol{z}_{i-1}+\varphi(\boldsymbol{k}_i)
\end{aligned}\end{equation}
这说明这种Attention可以作为一个RNN模型用递归的方式实现,它的空间复杂度最低,但是要串性计算,适合预测解码时使用;第二种是直接将\varphi(\boldsymbol{K}),\boldsymbol{V}\in\mathbb{R}^{n\times d}做外积,得到一个n\times d\times d的矩阵,然后对n那一维执行\text{cumsum}运算,这样就一次性得到\boldsymbol{S}_1,\boldsymbol{S}_2,\dots,\boldsymbol{S}_n了,它的速度最快,但空间占用最大,适合训练时使用,不过很多时候都有d^2\gg n,一般情况下训练时都很难承受这个空间复杂度,因此多数还是用RNN形式。
下采样技术 #
从结果上来看,Linformer的\boldsymbol{E}\boldsymbol{K}, \boldsymbol{F}\boldsymbol{V}就是将序列变短(下采样)了,而将序列变短的一个最朴素的方法就是Pooling了,所以笔者之前也尝试过把Pooling技术引入到Transformer中去。近来也有类似的工作发出来,比如IBM的《PoWER-BERT: Accelerating BERT Inference via Progressive Word-vector Elimination》和Google的《Funnel-Transformer: Filtering out Sequential Redundancy for Efficient Language Processing》。除了Pooling之外,其实还有其他的下采样技术,比如可以通过stride > 1的一维卷积来实现,基于这个思路,或许我们可以把FFN里边的Position-Wise全连接换成stride > 1的一维卷积?总之这方面应该也能玩出很多花样来,不过跟Linformer一样,这样糅合之后做自回归生成就很难了。
文章小结 #
本文介绍了一些从结构上对Attention进行修改从而降低其计算复杂度的工作,其中最主要的idea是去掉标准Attention中的Softmax,就可以使得Attention的复杂度退化为理想的\mathcal{O}(n)级别(Linear Attention)。相比于其他类似的改进结构的工作,这种修改能在把复杂度降到\mathcal{O}(n)的同时,依然保留所有的“token-token“的注意力,同时还能保留用于做自回归生成的可能性。
转载到请包括本文地址:https://spaces.ac.cn/archives/7546
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Jul. 04, 2020). 《线性Attention的探索:Attention必须有个Softmax吗? 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/7546
@online{kexuefm-7546,
title={线性Attention的探索:Attention必须有个Softmax吗?},
author={苏剑林},
year={2020},
month={Jul},
url={\url{https://spaces.ac.cn/archives/7546}},
}
July 4th, 2020
可惜的是,式10在实际使用的时候如果使用并行策略,空间复杂度会达到n×d×d,比原先的n×n还大,显存爆炸。。。还是得依赖作者的cuda编程才行
我没看作者的cuda代码,它是用cuda实现了RNN版本吗?就是 CuDNNLSTM vs 普通的LSTM 一样,也能快上好多倍。
不是很清楚,我也没看cuda代码,但是n×d×d这个思路不可行。他cuda应该是不会占用这么大的显存
每个 step 只需要保留一个 dxd, 循环 n 次, 就可以模仿 RNN 写了.
每个step只需要一个dxd,循环n次,那这样在训练的时候会不会很慢啊
如果序列比较短,那么可以按照常规的attention思路,先算\boldsymbol{Q}\boldsymbol{K}^{\top},然后做好mask再乘以\boldsymbol{V}。训练完成后,在预测阶段再转为Linear Attention的RNN模式。
July 4th, 2020
前提假设是否有问题呢?有的BERT的d有768,但是n最大只有512。
(另外评论有bug。。输入的跟显示的不一样。。)
768要除以head数,也就是12或16
本文已经补充说了,对于BERT base来说,d=64而不是768,因为是Multi-Head Attention,768是12个head的总和,每个head是64。
(评论的bug是因为你刚好输入了符号 > ,如果输入符号 > 不能跟其他符号粘着,不然会被当作html解释。实际上符号\gg的latex代码是\gg.)
明白了,谢谢!@ewrfcas|comment-13733
@苏剑林|comment-13734
July 9th, 2020
謝謝博主,我是來自台灣的讀者,今天才發現您的網站,有種相見恨晚的感覺,目前是研一也是研究自然語言處理,但和您真是相差太多了,感謝博主的分享。
过奖了,欢迎常来。
July 12th, 2020
公式(2)应该有问题吧,分母应该还要对i遍历求和吧
没有问题。每个下标为i的分量独立做归一化。
还是不能理解公式(2)
先接受,后理解。
要的是每一个分量都满足归一化,不是整体满足。
August 2nd, 2020
苏神,我想到一个bug,BERT在mask的时候使用的是把softmax(x)里的x置为-∞,用了这个近似以后可能会有问题?
这算什么bug...你就不能换一种方式,直接对\text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)进行mask么...
transformer are rnn这篇文章中,mask确实好像没法用进去
transformer are rnn这篇文章就是直接使用rnn的方式来mask掉将来信息。
苏神,我理解对sim(\bf{q}_i, \bf{k}_j)添加mask应该是一个矩阵元素相乘操作,那怎么再应用矩阵结合律来避免O(n^2)复杂度呢?
如果你说的是mask掉padding部分,那么只需要将\boldsymbol{V}对应的padding部分置零就差不多了。
如果是一个自回归的语言模型呢?
认真看本文“自回归生成”一节
December 25th, 2020
摘掉 softmax 一节的 QK^T 的维度应该是 nxm ~
额,我知道了。苏神统一了维度。。抱歉~
February 3rd, 2021
苏神 bert4keras能集成一下, longformer不
看到有中文版放出,我了解一下~不过看上去高效实现并不容易。
看一些测试效果竟然比 roberta还好,期待
February 19th, 2021
请问下,应该是用kernel来approximate normalized e^{q^Tk} 而不是直接的e^{q^Tk}? 不然的话底下的normalization似乎被直接扔了?
你说的是Performer?没有把归一化因子扔掉。你认真看看(3)式,是有归一化的,而且计算量主要集中在\sum\limits_{j=1}^n \text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)\boldsymbol{v}_j,所以只需要把这一步的计算解决了,那么归一化因子也很容易计算(相当于\boldsymbol{v}_j=1),所以主要研究如何降低\sum\limits_{j=1}^n \text{sim}(\boldsymbol{q}_i, \boldsymbol{k}_j)\boldsymbol{v}_j的计算量。
March 2nd, 2021
请问(9)式分母的级数,应该是n还是i呢?谢谢~~
噢噢,是i,改过来了,谢谢指出。
March 12th, 2021
苏神你好,请问式7中,为什么说做l2归一化就可以保证qk>=-1呢?确实一时没想明白
两个模长为1的向量的内积等于它们的夹角余弦,这个知道吧?余弦值在[-1,1]之间,这个知道吧?