6个派生优化器的简单介绍及其实现
By 苏剑林 | 2019-11-25 | 49900位读者 |优化器可能是深度学习最“玄学”的一个模块之一了:有时候换一个优化器就能带来明显的提升,有时候别人说提升很多的优化器用到自己的任务上却一丁点用都没有,理论性质好的优化器不一定工作得很好,纯粹拍脑袋而来的优化器也未必就差了。但不管怎样,优化器终究也为热爱“深度炼丹”的同学提供了多一个选择。
近几年来,关于优化器的工作似乎也在慢慢增多,很多论文都提出了对常用优化器(尤其是Adam)的大大小小的改进。本文就汇总一些优化器工作或技巧,并统一给出了代码实现,供读者有需调用。
基本形式 #
所谓“派生”,就是指相关的技巧都是建立在已有的优化器上的,任意一个已有的优化器都可以用上这些技巧,从而变成一个新的优化器。
已有的优化器的基本形式为:
\begin{equation}\begin{aligned}\boldsymbol{g}_t =&\, \nabla_{\boldsymbol{\theta}} L\\
\boldsymbol{h}_t =&\, f(\boldsymbol{g}_{\leq t})\\
\boldsymbol{\theta}_{t+1} =&\, \boldsymbol{\theta}_t - \gamma \boldsymbol{h}_t
\end{aligned}\end{equation}
其中$\boldsymbol{g}_t$即梯度,而$\boldsymbol{g}_{\leq t}$指的是截止到当前步的所有梯度信息,它们经过某种运算$f$(比如累积动量、累积二阶矩校正学习率等)后得到$\boldsymbol{h}_t$,然后由$\boldsymbol{h}_t$来更新参数,这里的$\gamma$就是指学习率。
变式杂烩 #
下面介绍优化器的6个变式,也可以理解为用优化器时的一些技巧。这些技巧有时候会很有效,有时候也可能无效甚至反作用,不能一概而论,只能理解为多一种选择就多一种可能。
权重衰减 #
权重衰减指的是直接在优化器每一步更新后面加上一个衰减项:
\begin{equation}\begin{aligned}\boldsymbol{g}_t =&\, \nabla_{\boldsymbol{\theta}} L\\
\boldsymbol{h}_t =&\, f(\boldsymbol{g}_{\leq t})\\
\boldsymbol{\theta}_{t+1} =&\, \boldsymbol{\theta}_t - \gamma \boldsymbol{h}_t - \gamma \lambda \boldsymbol{\theta}_t
\end{aligned}\end{equation}
其中$\lambda$称为“衰减率”。在SGD中,权重衰减等价于往loss里边加入$l_2$正则项$\frac{1}{2}\lambda \Vert \boldsymbol{\theta}\Vert_2^2$,但在Adagrad、Adam等带有自适应学习率的优化器中,$f$变成了非线性,所以两者不在等价。《Decoupled Weight Decay Regularization》一文特别指出权重衰减的防过拟合能力优于对应的$l_2$正则,推荐大家使用权重衰减而不是$l_2$正则。
层自适应 #
在优化器中,最后的更新量由$\boldsymbol{h}_t$和学习率$\gamma$决定,有时候$\boldsymbol{h}_t$的模长会大于参数$\boldsymbol{\theta}_t$的模长,这可能会导致更新的不稳定。所以,一个直接的想法是:每一层的参数的更新幅度应该由$\boldsymbol{\theta}_t$的模长的模长来调控。这个直接的想法就导致了如下的优化器变体:
\begin{equation}\begin{aligned}\boldsymbol{g}_t =&\, \nabla_{\boldsymbol{\theta}} L\\
\boldsymbol{h}_t =&\, f(\boldsymbol{g}_{\leq t})\\
\boldsymbol{\theta}_{t+1} =&\, \boldsymbol{\theta}_t - \gamma \boldsymbol{h}_t\times \frac{\Vert\boldsymbol{\theta}_t\Vert_2}{\Vert\boldsymbol{h}_t\Vert_2}
\end{aligned}\end{equation}
如果基础优化器是Adam,那么上述优化器就是LAMB。论文《Large Batch Optimization for Deep Learning: Training BERT in 76 minutes》指出LAMB在batch size较大(成千上万)的时候比Adam效果要好。
分段线性学习率 #
学习率也是优化器中的一个迷之存在,通常来说精调学习率策略能够获得一定提升,而不恰当的学习率甚至可能导致模型不收敛。常见的学习率策略有warmup、指数衰减、断层式下降(比如某个epoch后直接降到原来的1/10)等,比较迷的还有cos型学习率策略、多项式型学习率策略等。
考虑到常见的函数都可以用分段线性函数逼近,所以笔者干脆引入了一个分段线性学习率的策略,供大家随便玩。形式如下:
\begin{equation}\begin{aligned}\boldsymbol{g}_t =&\, \nabla_{\boldsymbol{\theta}} L\\
\boldsymbol{h}_t =&\, f(\boldsymbol{g}_{\leq t})\\
\boldsymbol{\theta}_{t+1} =&\, \boldsymbol{\theta}_t - \gamma \rho_t\boldsymbol{h}_t
\end{aligned}\end{equation}
其中$\rho_t$是某个以步数$t$为自变量的分段线性函数。
梯度累积 #
梯度累积之前在《用时间换取效果:Keras梯度累积优化器》一文中也介绍过了,其实不算优化器的变式,但可以写到优化器中,通过小batch size达到大batch size的效果,实现时间换空间。更大的batch size有时候能提升效果,尤其是基准batch size过小的情况下(8以下?)。
重述《用时间换取效果:Keras梯度累积优化器》一文关于梯度下降的描述:
所谓梯度累积,其实很简单,我们梯度下降所用的梯度,实际上是多个样本算出来的梯度的平均值,以batch_size=128为例,你可以一次性算出128个样本的梯度然后平均,我也可以每次算16个样本的平均梯度,然后缓存累加起来,算够了8次之后,然后把总梯度除以8,然后才执行参数更新。当然,必须累积到了8次之后,用8次的平均梯度才去更新参数,不能每算16个就去更新一次,不然就是batch_size=16了。
Lookahead #
Lookahead优化器来自论文《Lookahead Optimizer: k steps forward, 1 step back》,在之前的文章《Keras实现两个优化器:Lookahead和LazyOptimizer》也介绍过。Lookahead的含义是用常用优化器预先往前摸索几步,然后根据摸索结果来更新,流程如下:
\begin{equation}\begin{aligned}&\boldsymbol{g}_t =\, \nabla_{\boldsymbol{\theta}} L\\
&\boldsymbol{h}_t =\, f(\boldsymbol{g}_{\leq t})\\
&\boldsymbol{\theta}_{t+1} =\, \boldsymbol{\theta}_t - \gamma\boldsymbol{h}_t\\
&\text{如果}t\,\text{mod}\,k = 0\text{:}\\
&\qquad\boldsymbol{\Theta}_{t+1} = \boldsymbol{\Theta}_t + \alpha (\boldsymbol{\theta}_{t+1}- \boldsymbol{\Theta}_t)\\
&\qquad\boldsymbol{\theta}_{t+1} = \boldsymbol{\Theta}_{t+1} \,(\text{即覆盖原来的}\boldsymbol{\theta}_{t+1})
\end{aligned}\end{equation}
其实这个优化器叫Lookaback也无妨,也就是每走几步就往后看看,跟几步前的权重做个插值。
Lazy优化器 #
Lazy优化器在刚才的文章《Keras实现两个优化器:Lookahead和LazyOptimizer》中也介绍过,其本意就是Embedding层的更新应当稀疏化一些,这有助于防止过拟合。(参考知乎讨论)
参考实现 #
前面的介绍比较简单,事实上这些变式本身确实不难理解,关键还是代码实现。从前面的介绍中大家可以发现,这6个变式并无矛盾之处,因此良好的实现应当能让我们能搭积木般组合其中一个或多个变式使用;另外,目前keras也有两个分支:纯keras和tf.keras,良好的实现应当能同时兼容它们(或者同时给出两种实现)。
终极移花接木 #
虽然部分变式在之前也实现过了,但这里还是用新的方式重新实现了它们。这里的实现方式源于笔者意外发现的一种移花接木技巧。
假设我们有这样一个类:
import numpy as np
class A(object):
def __init__(self):
self.a = np.ones(1)
self.b = np.ones(2)
self.c = np.ones(3)
然后假设我们要继承A类得到一个B类,B类是要将__init__
方法的所有np.ones
替换为np.zeros
,其余都不变。由于__init__
可能是一个很复杂的流程,如果将它完整复制过来然后改写显然太冗余了。
有没有直接写几行代码就能替换所有的呢?还真有!
class B(A):
def __init__(self):
_ = np.ones
np.ones = np.zeros
super(B, self).__init__()
np.ones = _
有了这个demo,我们就可以“魔改”已有优化器了。在keras中,参数的更新都是通过K.update
来实现的(参考keras的optimizers.py),我们只需要用上述方式重新定义一下K.update
就好。
如果是tf.keras呢?很遗憾,这种方式不可行,因为tf.keras中常用优化器的迭代流程都被写到C里边去了(参考tf.keras的adam.py),我们看不到代码,也就不能用这种方法改了。一个解决办法是我们自己重新实现一个优化器如Adam,将迭代流程暴露出来,这样我们就能用上述方式魔改了。
使用示例 #
根据上述思路统一实现的6个优化器变体,都被放在笔者的bert4keras项目中:bert4keras.optimizers。
所有函数会根据keras还是tf.keras来进行正确的导入,从而实现keras/tf.keras都能用同样的方式使用。里边自带了一个Adam实现,这个Adam是专门写给tf.keras的。对于tf.keras,如果想要实现上述变式,那只能用bert4keras自带的优化器(目前只有Adam),不能用tf.keras内置的优化器。
参考代码:
from bert4keras.optimizers import *
# 变成带权重衰减的Adam
AdamW = extend_with_weight_decay(Adam, 'AdamW')
optimizer = AdamW(learning_rate=0.001, weight_decay_rate=0.01)
# 变成带分段线性学习率的Adam
AdamLR = extend_with_piecewise_linear_lr(Adam, 'AdamLR')
# 实现warmup,前1000步学习率从0增加到0.001
optimizer = AdamLR(learning_rate=0.001, lr_schedule={1000: 1.})
# 变成带梯度累积的Adam
AdamGA = extend_with_gradient_accumulation(Adam, 'AdamGA')
optimizer = AdamGA(learning_rate=0.001, grad_accum_steps=10)
# 组合使用
AdamWLR = extend_with_piecewise_linear_lr(AdamW, 'AdamWLR')
# 带权重衰减和warmup的优化器
optimizer = AdamWLR(learning_rate=0.001,
weight_decay_rate=0.01,
lr_schedule={1000: 1.})
(注:一下实现这么多个优化器,而且同时要兼容keras和tf.keras,难免可能会有错漏之处,如有发现,万望不吝指正。)
写在最后 #
炼丹不易,且炼且珍惜。
转载到请包括本文地址:https://spaces.ac.cn/archives/7094
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Nov. 25, 2019). 《6个派生优化器的简单介绍及其实现 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/7094
@online{kexuefm-7094,
title={6个派生优化器的简单介绍及其实现},
author={苏剑林},
year={2019},
month={Nov},
url={\url{https://spaces.ac.cn/archives/7094}},
}
March 15th, 2020
执行优化器说 无learning_rate 参数
keras < 2.3.0 需要改learning_rate为lr,这些只需要去看看keras的源码就知道答案了。
July 4th, 2020
苏神,bert4keras里的这个extend_with_gradient_accumulation是不是有bug? 现在训一个unilm 优化器用AdamGA loss不会下降,但换成Adam或AdamW的时候都很正常
September 24th, 2020
感谢大佬分享!!
August 21st, 2021
苏神 请教一下 打比赛用Keras的时候,要求结果复现。
于是,我在主程序最前面加上了
random.seed(SEED)
os.environ['PYTHONHASHSEED'] = str(SEED)
np.random.seed(SEED)
tf.set_random_seed(SEED)
但是,结果仍然无法复现...需要怎么改吗?还是说keras的确存在无法复现的问题...而且经过测试优化器使用SGD可以完全复现,但是使用Adam或者RMSprop结果不能完全复现
具体环境如下:
TensorFlow backend (yes / no): yes
TensorFlow version: 1.14.0
Keras version:2.3.1
Python version: 3.6.13
CUDA/cuDNN version: Cuda compilation tools, release 10.0, V10.0.130
GPU model and memory: Tesla V100-DGXS-32GB
tf + gpu基本上都没什么完全的可复现性,不要做这个尝试了~
我也遇到这个问题,在pytorch多次实验发现,只有把Adam系的beta2改得很小(变回SGD),才能保证复现
May 4th, 2022
哈喽,请教一个问题, Adam优化器有权重衰减变种AdamW, 即将weight decay decouple开, AdamW优化器在很多任务中被采用.那么为啥LAMB优化器的权重衰减变种(假设叫做LAMBW,hhh)却没有看到大家使用呢? 这里头存在什么更深的原因吗?
你说的LAMBW,是指将weight decay放到LAMB外边?LAMB本身是想着通过参数的模长来控制更新量的模长,将weight decay放到LAMB外边不符合LAMB的设计初衷。
August 8th, 2022
[...]6个派生优化器的简单介绍及其实现 – 科学空间[...]